Re: First function debug help
- To: mathgroup at smc.vnet.net
- Subject: [mg106629] Re: First function debug help
- From: Canopus56 <canopus56 at yahoo.com>
- Date: Tue, 19 Jan 2010 05:13:47 -0500 (EST)
Please disregard the prior message. I was able to debug the function and get it running by applying a little more brain power. The final functions are appended. I could use some help in coaxing the JulianDayNumber output from an exponential form, e.g. - 2.43612*10^6 into the form: 2436116.31 Thanks again for all your help. - Kurt JulianDayNumber[{y_, m_, d_, h_, mins_, s_}] := Module[ {b, f, JD, MJD}, (* This function computes the Modified Julian Day from a system \ formatted date. Domain is restricted to Greogorian dates. Source: Meeus. 1998. Chap. 7. Astronomical Alogrithms. p. 61. Test data is for launch of Sputnik on {1957,10,4,19,26, 24} which should yield JD2436116 .31 *) b = 2 - IntegerPart[y/100] + IntegerPart[IntegerPart[y/100]/4]; (* fractionalize the day value *) f = d + (h/24) + (mins / ( 24 * 60 )) + (s /(24 * 3600)); JD = IntegerPart[(365.25 *(y + 4716))] + IntegerPart[(30.6001*(m + 1))] + f + b - 1524.5; MJD = JD - 2400000.5; (* NumberForm[JD,{15,10}] *) JD ] JulianDayModified[{y_, m_, d_, h_, mins_, s_}] := Module[ {b, f, JD, MJD}, (* This function computes the Modified Julian Day from a system \ formatted date. Domain is restricted to Greogorian dates. Source: Meeus. 1998. Chap. 7. Astronomical Alogrithms. p. 61. Test data is for launch of Sputnik on {1957,10,4,19,26, 24} which should yield MJD36115 .81 *) b = 2 - IntegerPart[y/100] + IntegerPart[IntegerPart[y/100]/4]; (* fractionalize the day value *) f = d + (h/24) + (mins / ( 24 * 60 )) + (s /(24 * 3600)); JD = IntegerPart[(365.25 *(y + 4716))] + IntegerPart[(30.6001*(m + 1))] + f + b - 1524.5; MJD = JD - 2400000.5; (* NumberForm[MJD,{15,5}] *) MJD ]
- Follow-Ups:
- Re: Re: First function debug help
- From: DrMajorBob <btreat1@austin.rr.com>
- Re: Re: First function debug help