Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/mooncolong.c
- Timestamp:
- Jun 15, 2004, 6:54:12 PM (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/XephemAstroLib/mooncolong.c
r1719 r2551 3 3 4 4 #include <stdio.h> 5 #include <stdlib.h> 5 6 #include <math.h> 6 #if defined(__STDC__) 7 #include <stdlib.h> 8 #endif 9 10 #include "P_.h" 7 11 8 #include "astro.h" 12 9 13 static void Librations P_((double RAD, double LAMH, double BH, double OM,14 double F, double L, double L1, double *L0, double *B0) );15 static void Moon P_((double RAD, double T, double T2, double LAM0, double R,10 static void Librations (double RAD, double LAMH, double BH, double OM, 11 double F, double L, double L1, double *L0, double *B0); 12 static void Moon (double RAD, double T, double T2, double LAM0, double R, 16 13 double M, double *F, double *L1, double *OM, double *LAM, double *B, 17 double *DR, double *LAMH, double *BH) );18 static void Sun P_((double RAD, double T, double T2, double *L, double *M,19 double *R, double *LAM0) );14 double *DR, double *LAMH, double *BH); 15 static void Sun (double RAD, double T, double T2, double *L, double *M, 16 double *R, double *LAM0); 20 17 21 18 /* given a Julian date and a lunar location, find selenographic colongitude of … … 26 23 */ 27 24 void 28 moon_colong (jd, lt, lg, cp, kp, ap, sp) 29 double jd; /* jd */ 30 double lt, lg; /* lat/long of location on moon, rads +N +E */ 31 double *cp; /* selenographic colongitude (-lng of rising sun), rads */ 32 double *kp; /* illuminated fraction of surface from Earth */ 33 double *ap; /* sun altitude at location, rads */ 34 double *sp; /* lunar latitude of subsolar point, rads */ 25 moon_colong ( 26 double jd, /* jd */ 27 double lt, /* lat of location on moon, rads +N +E */ 28 double lg, /* long of location on moon, rads +N +E */ 29 double *cp, /* selenographic colongitude (-lng of rising sun), rads */ 30 double *kp, /* illuminated fraction of surface from Earth */ 31 double *ap, /* sun altitude at location, rads */ 32 double *sp) /* lunar latitude of subsolar point, rads */ 35 33 { 36 34 double RAD = .0174533; … … 89 87 90 88 static void 91 Librations (RAD, LAMH, BH, OM, F, L, L1, L0, B0) 92 double RAD; 93 double LAMH; 94 double BH; 95 double OM; 96 double F; 97 double L; 98 double L1; 99 double *L0; 100 double *B0; 89 Librations (double RAD, double LAMH, double BH, double OM, double F, 90 double L, double L1, double *L0, double *B0) 101 91 { 102 92 double I, PSI, W, NUM, DEN, A, TEMP; … … 123 113 124 114 static void 125 Moon (RAD, T, T2, LAM0, R, M, F, L1, OM, LAM, B, DR, LAMH, BH) 126 double RAD; 127 double T; 128 double T2; 129 double LAM0; 130 double R; 131 double M; 132 double *F; 133 double *L1; 134 double *OM; 135 double *LAM; 136 double *B; 137 double *DR; 138 double *LAMH; 139 double *BH; 115 Moon (double RAD, double T, double T2, double LAM0, double R, double M, 116 double *F, double *L1, double *OM, double *LAM, double *B, double *DR, 117 double *LAMH, double *BH) 140 118 { 141 119 double T3, M1, D2, SUMR, SUML, DIST; … … 175 153 176 154 static void 177 Sun (RAD, T, T2, L, M, R, LAM0) 178 double RAD; 179 double T; 180 double T2; 181 double *L; 182 double *M; 183 double *R; 184 double *LAM0; 155 Sun (double RAD, double T, double T2, double *L, double *M, double *R, 156 double *LAM0) 185 157 { 186 158 double T3, C, V, E, THETA, OM; … … 238 210 if (ac != 2) { 239 211 fprintf (stderr, "%s: JD\n", av[0]); 240 exit (1);212 abort(); 241 213 } 242 214 … … 262 234 263 235 /* For RCS Only -- Do Not Edit */ 264 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 200 1-10-22 12:08:27 $ $Revision: 1.2$ $Name: not supported by cvs2svn $"};236 static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: mooncolong.c,v $ $Date: 2004-06-15 16:52:39 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};
Note:
See TracChangeset
for help on using the changeset viewer.