Changeset 19719
- Timestamp:
- Sep 23, 2008, 3:06:18 PM (18 years ago)
- Location:
- trunk/Ohana/src/skycalc/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/skycalc/src/sundata.c
r19681 r19719 1 1 # include <skycalc_internal.h> 2 2 # include <ohana.h> 3 # define VERBOSE 0 3 4 4 5 void SC_set_site (double *longit, double *lat, double *elevsea, double *elev) { … … 42 43 date.s = (float) (stm->tm_sec); 43 44 44 /* 45 fprintf (stderr, "%4d/%02d/%02d %02d:%02d:%02f\n", date.y, date.mo, date.d, date.h, date.mn, date.s); 46 */ 45 if (VERBOSE) fprintf (stderr, "%4d/%02d/%02d %02d:%02d:%02f\n", date.y, date.mo, date.d, date.h, date.mn, date.s); 47 46 48 47 SC_set_site (&longit, &lat, &elevsea, &elev); … … 52 51 sid = SC_lst(jdnow, longit); 53 52 SC_lpsun (jdnow, &Rsun, &Dsun); 53 54 if (VERBOSE) fprintf (stderr, "jdnow: %lf, sid: $lf\n", jdnow, sid); 54 55 55 56 /* dot product of unit vectors of (RAo,DECo) & (Rsun,Dsun) */ -
trunk/Ohana/src/skycalc/src/sunmoon.c
r19681 r19719 77 77 sid = SC_lst(jdnow, longitude); 78 78 79 if (VERBOSE) fprintf (stderr, "jdnow: %lf, sid: $lf\n", jdnow, sid); 80 79 81 // get the sun coordinates 80 82 SC_lpsun (jdnow, &sun_ra, &sun_dec); … … 84 86 sun_angle = angular_separation (sun_ra, sun_dec, RAo, DECo); 85 87 86 if (VERBOSE) fprintf (std out, "sun @ ra,dec = %f %f : alt, az = %f %f\n", sun_ra, sun_dec, sun_alt, sun_az);88 if (VERBOSE) fprintf (stderr, "sun @ ra,dec = %f %f : alt, az = %f %f\n", sun_ra, sun_dec, sun_alt, sun_az); 87 89 88 90 // get the moon coordintes … … 93 95 moon_angle = angular_separation (moon_ra, moon_dec, RAo, DECo); 94 96 95 if (VERBOSE) fprintf (std out, "moon @ ra,dec = %f %f : alt, az = %f %f\n", moon_ra, moon_dec, moon_alt, moon_az);97 if (VERBOSE) fprintf (stderr, "moon @ ra,dec = %f %f : alt, az = %f %f\n", moon_ra, moon_dec, moon_alt, moon_az); 96 98 97 99 phase = (sun_ra - moon_ra - 180)/360.0; … … 99 101 while (phase > 0.5) phase -= 1.0; 100 102 101 if (VERBOSE) fprintf (std out, "sun @ %f %f dist %f\n", sun_ra, sun_dec, sun_angle);102 if (VERBOSE) fprintf (std out, "moon @ %f %f angle %f phase %f\n", moon_ra, moon_dec, moon_angle, phase);103 if (VERBOSE) fprintf (stderr, "sun @ %f %f dist %f\n", sun_ra, sun_dec, sun_angle); 104 if (VERBOSE) fprintf (stderr, "moon @ %f %f angle %f phase %f\n", moon_ra, moon_dec, moon_angle, phase); 103 105 104 106 fprintf (stdout, "-sun_alt %f -sun_angle %f -moon_alt %f -moon_angle %f -moon_phase %f\n", sun_alt, sun_angle, moon_alt, moon_angle, phase);
Note:
See TracChangeset
for help on using the changeset viewer.
