Services & ResourcesWolfram ForumsMathGroup Archive

[Date Index] [Thread Index] [Author Index]

First function debug help

• To: mathgroup at smc.vnet.net
• Subject: [mg106627] First function debug help
• From: Canopus56 <canopus56 at yahoo.com>
• Date: Tue, 19 Jan 2010 05:13:22 -0500 (EST)

I took a stab at writing my first function - converting a system formatted list date into a Modified Julian Day.  The function appears to be written properly, but does not return anything. 

Any help in debugging it would be appreciated. 

Ideally, I would like to send a date-time list in the form {y,m,d,h,m,s} to the function and have the Julian Day returned. 

Thanks for your help - Kurt 

(* 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. *)
(* fractionalize the day value *)

JulianDayModified[y_, m_, d_, h_, m_, s_] := (
 B = 2 - IntegerPart[y/100] + IntegerPart[IntegerPart[y/100]/4]; 
 d = d + ((m * 60 )/ ( 24 * 60 )) + ((s * 3600)/(24 * 3600));
 JD = IntegerPart[(365.25*(y + 4716))] + 
 IntegerPart[(30.6001*(m + 1 ))] + d + b - 1524.5;
 MJD = JD - 2400000.5;
 MJD
 )

JulianDayModified[2010, 1, 1, 12, 0, 0]

Returns the string " 
JulianDayModified[2010, 1, 1, 12, 0, 0]" 

and not the computed date

Also tried it this way with the Module statement with no change in the result: 

JulianDayModified[y_, m_, d_, h_, m_, s_] := Module[{B, JD, MJD},
 B = 2 - IntegerPart[y/100] + IntegerPart[IntegerPart[y/100]/4];
 d = d + ((m*60)/(24*60)) + ((s*3600)/(24*3600));
 JD = IntegerPart[(365.25*(y + 4716))] + 
 IntegerPart[(30.6001*(m + 1))] + d + b - 1524.5;
 MJD = JD - 2400000.5;
 MJD
 ]

• Follow-Ups:
  • Re: Re: First function debug help
    • From: Herman Kuun <oomkoos1@gmail.com>
  • Re: Re: First function debug help
    • From: DrMajorBob <btreat1@austin.rr.com>
  • Re: First function debug help
    • From: Herman Kuun <oomkoos1@gmail.com>