IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 7, 2008, 3:31:01 PM (18 years ago)
Author:
eugene
Message:

big update from eam_branch_20081124 with updates to Opihi math

Location:
trunk/Ohana/src/skycalc/src
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/skycalc/src/astro.c

    r19681 r20936  
    106106
    107107  double incl;
    108   double xpr,ypr,zpr;
     108  double ypr,zpr;
    109109  double T;
    110110
    … …  
    141141double SC_etcorr (double jd) {
    142142
    143   double jd1900 = 2415019.5;
    144143  double dates[20] = {1900,1905,1910,1915,1920,1925,1930,1935,1940,1945,
    145144                      1950,1955,1960,1965,1970,1975,1980,1985,1990,1993};

  • trunk/Ohana/src/skycalc/src/dusktime.c

    r19681 r20936  
    1313}
    1414
    15 main (int argc, char **argv) {
     15int main (int argc, char **argv) {
    1616
    1717  struct SC_date_time date, tmpdate;
    … …  
    5959
    6060  if (VERBOSE) {
    61     jd_to_date (jdnow, &tmpdate);
     61    SC_jd_to_date (jdnow, &tmpdate);
    6262    fprintf (stderr, "%4d/%02d/%02d %02d:%02d:%02f\n", tmpdate.y, tmpdate.mo, tmpdate.d, tmpdate.h, tmpdate.mn, tmpdate.s);
    6363  }

  • trunk/Ohana/src/skycalc/src/geometry.c

    r19681 r20936  
    88   double mod;    /* modulus */
    99   double xy;     /* component in xy plane */
    10    short sign;    /* for determining quadrant */
    1110   double radian_ra, radian_dec;
    1211

  • trunk/Ohana/src/skycalc/src/moon.c

    r19681 r20936  
    368368  double geora, geodec, geodist;  /* geocent for moon, not used here.*/
    369369  double ramoon, decmoon, distmoon;
    370   double hamoonrise, hamoonset, tmoonrise, tmoonset, jdmoonrise, jdmoonset;
    371   short dow; /* day of week */
    372   struct SC_date_time date_midnight, tmpdate;
     370  double hamoonset, tmoonset, jdmoonset;
     371  struct SC_date_time date_midnight;
    373372  double dt, lst0, lst1, djd, horiz;
    374373
    … …  
    420419  double geora, geodec, geodist;  /* geocent for moon, not used here.*/
    421420  double ramoon, decmoon, distmoon;
    422   double hamoonrise, hamoonset, tmoonrise, tmoonset, jdmoonrise, jdmoonset;
    423   short dow; /* day of week */
    424   struct SC_date_time date_midnight, tmpdate;
     421  double hamoonset, tmoonrise, jdmoonrise;
     422  struct SC_date_time date_midnight;
    425423  double dt, lst0, lst1, djd, horiz;
    426424

  • trunk/Ohana/src/skycalc/src/moondata.c

    r19681 r20936  
    1212}
    1313
    14 main (int argc, char **argv) {
     14int main (int argc, char **argv) {
    1515
    16   struct SC_date_time date, tmpdate;
     16  struct SC_date_time date;
    1717  double longit, lat, elevsea, elev;
    1818  double jdnow, ra, dec, sid;

  • trunk/Ohana/src/skycalc/src/sun.c

    r19681 r20936  
    55void SC_lpsun (double jd, double *ra, double *dec) {
    66
    7   double n, L, g, lambda,epsilon,alpha,delta,x,y,z;
     7  double n, L, g, lambda,epsilon,x,y,z;
    88
    99  n = jd - J2000;
    … …  
    5454}
    5555
    56 /* Given site position, prints Sun info for the given night.
     56/* Given site position, prints Sun info for the given night. */
    5757/* dates are all in UT now */
    5858double SC_sunset_tonight (struct SC_date_time date, double lat, double longit, double elev) {
    5959
    6060  double jd, jdmid0, jdmid, stmid;
    61   double rasun, decsun, min_alt, max_alt, horiz;
    62   double hasunset, jdsunset, jdsunrise, sid;
     61  double rasun, decsun, horiz;
     62  double hasunset, jdsunset;
    6363  double dt, lst0, lst1, djd;
    64   struct SC_date_time date_midnight, tmpdate;
    65   short dow;
     64  struct SC_date_time date_midnight;
    6665
    6766  horiz = sqrt (2. * elev / 6378140.) * DEG_IN_RADIAN;
    … …  
    106105}
    107106
    108 /* Given site position, prints Sun info for the given night.
     107/* Given site position, prints Sun info for the given night. */
    109108/* dates are all in UT now */
    110109double SC_sunrise_tonight (struct SC_date_time date, double lat, double longit, double elev) {
    111110
    112111  double jd, jdmid, stmid;
    113   double rasun, decsun, min_alt, max_alt, horiz;
    114   double hasunset, jdsunset, jdsunrise, sid;
     112  double rasun, decsun, horiz;
     113  double hasunset, jdsunrise;
    115114  double dt, lst0, lst1, djd;
    116   struct SC_date_time date_midnight, tmpdate;
    117   short dow;
     115  struct SC_date_time date_midnight;
    118116
    119117  horiz = sqrt (2. * elev / 6378140.) * DEG_IN_RADIAN;

  • trunk/Ohana/src/skycalc/src/sundata.c

    r19731 r20936  
    1212}
    1313
    14 main (int argc, char **argv) {
     14int main (int argc, char **argv) {
    1515
    16   struct SC_date_time date, tmpdate;
     16  struct SC_date_time date;
    1717  double longit, lat, elevsea, elev;
    18   double jdnow, ra, dec, sid, alt, az, sind, sinh, cosh;
    19   double geora, geodec, geodist, dist;
     18  double jdnow, sid, alt, az, sind, sinh, cosh;
    2019  double RAo, DECo, abx, aby, abz, cs, theta;
    21   double Rsun, Dsun, days, Hsun;
     20  double Rsun, Dsun, Hsun;
    2221  time_t tzero;
    2322  struct tm *stm;

  • trunk/Ohana/src/skycalc/src/sunmoon.c

    r19730 r20936  
    99double angular_separation (double ra, double dec, double RA, double DEC);
    1010
    11 main (int argc, char **argv) {
     11int main (int argc, char **argv) {
    1212
    1313  int N;
    14   struct SC_date_time date, tmpdate;
     14  struct SC_date_time date;
    1515  double longitude, latitude, elevation;
    1616  double jdnow, sid;

  • trunk/Ohana/src/skycalc/src/time.c

    r19681 r20936  
    66  time_t t;
    77  struct tm *stm;
    8   double jd, jdb, jde;
    9   short dow;
    108 
    119  t = time (0);
    … …  
    109107double SC_lst (double jd, double longit) {
    110108 
    111   double t, ut, jdmid, jdint, jdfrac, sid_g, sid;
     109  double t, ut, jdmid, jdint, jdfrac, sid_g;
    112110  long sid_int;
    113111 
Note: See TracChangeset for help on using the changeset viewer.