Re: Re: First function debug help
- To: mathgroup at smc.vnet.net
- Subject: [mg106746] Re: [mg106629] Re: First function debug help
- From: DrMajorBob <btreat1 at austin.rr.com>
- Date: Fri, 22 Jan 2010 05:40:13 -0500 (EST)
- References: <201001191013.FAA29050@smc.vnet.net>
- Reply-to: drmajorbob at yahoo.com
JulianDayNumber[{2006, 1, 21, 15, 22, 36}]
AccountingForm[%, 10]
2.45376*10^6
2453755.141
Bobby
On Tue, 19 Jan 2010 04:13:47 -0600, Canopus56 <canopus56 at yahoo.com> wrote:
> 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
> ]
>
>
>
--
DrMajorBob at yahoo.com
- References:
- Re: First function debug help
- From: Canopus56 <canopus56@yahoo.com>
- Re: First function debug help