Code covered by the BSD License
-
aberat (pos1, ve, tlight)
this function corrects position vector for aberration of light.
-
angles (pos)
this function converts a vector to angular quantities
-
applan1 (tjd, ujd, l, n, topo...
this function computes the apparent geocentric or topocentric place
-
atan3 (a, b)
four quadrant inverse tangent
-
brent (f, x1, x2, rtol)
solve for a single real root of a nonlinear equation
-
broot (f, x1in, x2in, factor,...
bracket a single root of a nonlinear equation
-
deg2dms (dd)
convert decimal degrees to degrees,
-
etilt1 (tjd)
this function computes quantities related to the orientation
-
findleap(jdate)
find number of leap seconds for utc julian date
-
funarg (t)
this subroutine computes fundamental arguments (mean elements)
-
gast2 (tjdh, tjdl, k)
this function computes the greenwich sidereal time
-
gdate (jdate)
convert Julian date to Gregorian (calendar) date
-
geocen (pos1, pe)
this function moves the origin of coordinates from the
-
getdate
interactive request and input of calendar date
-
getobs
interactive request of observer coordinates
-
hrs2hms (hrs)
convert decimal hours to hours,
-
jd2str(jdate)
convert Julian date to string equivalent
-
jplephem (et, ntarg, ncent)
reads the jpl planetary ephemeris and gives
-
julian (month, day, year)
Julian date
-
minima (f, a, b, tolm)
one-dimensional minimization
-
nutate1 (tjd, pos1)
this function nutates equatorial rectangular coordinates from
-
nutation (jdate)
nutation in longitude and obliquity
-
oevent1 (objfunc, prtfunc, ti...
predict minimization/root-finding orbital events
-
precess (tjd1, pos1, tjd2)
this function precesses equatorial rectangular coordinates from
-
readleap
read leap seconds data file
-
seclevnt (objfunc, prtfunc, t...
predict solar eclipse circumstances
-
seclfunc(x)
solar eclipse objective function
-
seclprt(iflag, jdutc)
print solar eclipse conditions
-
solsys (tjd, body, origin)
purpose
-
sunfld (pos1, pe)
this function corrects position vector for the deflection
-
tdtimes (tdbjd)
this function computes the terrestrial time (tt)
-
terra (glon, glat, ht, st)
this functions computes the position and velocity vectors of
-
utc2tdt (jdutc)
convert UTC julian date to TDT julian date
-
seclipse.m
-
View all files
A MATLAB Script for Predicting Solar Eclipses
by David Eagle
07 Dec 2012
(Updated 04 Feb 2013)
Predict local circumstances of solar eclipses.
|
Watch this File
|
| File Information |
| Description |
PDF document and a MATLAB script named seclipse.m that can be used to predict local circumstances of solar eclipses. This computer program provides the universal times and topocentric coordinates of the Sun and Moon at the beginning and end of the penumbra contacts, and the time and coordinates at maximum eclipse. The source ephemeris for this routine is a JPL binary ephemeris file. This application uses several functions ported to MATLAB from the Fortran version of the NOVAS (Naval Observatory Vector Astrometry Subroutines) source code developed at the United States Naval Observatory (www.usno.navy.mil/USNO/astronomical-applications/software-products/novas). JPL binary ephemeris files for Windows compatible computers can be downloaded at www.cdeagle.com. |
| MATLAB release |
MATLAB 7.12 (R2011a)
|
| Other requirements |
Requires platform specific version of MICE mex file and de421.bsp binary ephemeris file. Both are available at naif.jpl.nasa.gov/naif/toolkit_MATLAB.html |
|
Tags for This File
|
| Everyone's Tags |
|
| Tags I've Applied |
|
| Add New Tags |
Please login to tag files.
|
| Updates |
| 04 Feb 2013 |
Added MICE version of this script in the solar_eclipse_mice.zip archive. |
|
Contact us
