Re: Astronomical question
- To: mathgroup at smc.vnet.net
- Subject: [mg115359] Re: Astronomical question
- From: Robert Zimmerman <rlz at uoregon.edu>
- Date: Sun, 9 Jan 2011 02:16:32 -0500 (EST)
Here is a possible general template for lunar phases that will get
you started:
AD=AstronomicalData;
(* Enter the data in the form {Time, Magnitude}*)
mag={ DateList[{2010,1,#}],AD["Moon",{"ApparentMagnitude",DateList
[{2010,1,#}]}]}&/@Range[13,13+ 30 ] ;
(* date and magnitude for the full and new moon*)
fullnew= Take[SortBy[mag,Last] ,#]& /@{1,-1} //Flatten[#,1]&
(* constellations *)
AD["Moon",{"Constellation",#}]&/@{{2010,1,15},{2010,1,23},{2010,1,30},
{2010,2,5}}
(* options for graphics*)
optionsphase={Joined->True,GridLines->
{Automatic,Automatic} ,GridLinesStyle->Black,PlotStyle->
{Thick,Black},Filling-> Bottom,FillingStyle-> Opacity[.
4,Black] ,PlotLabel-> {"Full and New",fullnew } ,Epilog->{Text["New
\nSagittarius",{{2010,1,15},- 6}], Text["First\nAries",
{{2010,1,23},- 8}], Text["Full\nCancer",{{2010,1,30},- 12}], Text
["Last\nVirgo",{{2010,2,5},- 8}] } } ;
(* Graph*)
DateListPlot[mag ,optionsphase ]
> Mokambo wrote
> Dear all,
> -------
> I want to write a program in Mathematica using AstronomicalData to
> be able to trace the moon from a given point (given coordinates)
> and a certain period of time. I want to learn by doing so any
> links to useful formulas or a general template would be greatly
> appreciated.
> Thanks.