Changeset 1678 in Sophya for trunk/SophyaExt/XAstroPack/xastropack.h

Timestamp:

Oct 9, 2001, 11:47:57 PM (24 years ago)

Author:

cmv

Message:

EqHToHor et reverse cmv 9/10/01

File:

• trunk/SophyaExt/XAstroPack/xastropack.h (modified) (3 diffs)

trunk/SophyaExt/XAstroPack/xastropack.h

@@ -46 +46 @@
 double MJDfrYear(double yr);
 void YDfrMJD(double mjd,double *dy,int *yr);
+int IsLeapYear(int y);
+int DayOrder(double mjd,int *dow);
+int DaysInMonth(double mjd);
+double MJDat0hFrMJD(double mjd);
+double HfrMJD(double mjd);
 double GSTfrUTC(double mjd0,double utc);
 double UTCfrGST(double mjd0,double gst);
 double GST0(double mjd0);
+double LSTfrGST(double gst,double geolng);
+double LSTfrMJD(double mjd,double geolng);
 void Precess(double mjd1,double mjd2,double ra1,double dec1,double *ra2,double *dec2);
 double AirmassfrAlt(double alt);
-double HafrRaTS(double gst,double ra);
+double HafrRaTS(double lst,double ra);
+double RafrHaTS(double lst,double ha);
+double HelioCorr(double jd,double ra,double dec);
 void HMSfrHdec(double hd,int *h,int *mn,double *s);
 double HdecfrHMS(int h,int mn,double s);
@@ -58 +67 @@
 void EqtoGal(double mjd,double ra,double dec,double *glng,double *glat);
 void GaltoEq(double mjd,double glng,double glat,double *ra,double *dec);
-void EqtoHor(double geolat,double ha,double dec,double *az,double *alt);
-void HortoEq(double geolat,double az,double alt,double *ha,double *dec);
+void EqHtoHor(double geolat,double ha,double dec,double *az,double *alt);
+void HortoEqH(double geolat,double az,double alt,double *ha,double *dec);
+void EqtoHor(double geolat,double lst,double ra,double dec,double *az,double *alt);
+void HortoEq(double geolat,double lst,double az,double alt,double *ra,double *dec);
 void EqtoEcl(double mjd,double ra,double dec,double *eclng,double *eclat);
 void EcltoEq(double mjd,double eclng,double eclat,double *ra,double *dec);
@@ -69 +80 @@
 int  CoordConvertToStd(TypAstroCoord typ,double& coord1,double& coord2);
 
+/*!
+\brief Pour remettre la valeur "val" dans la dynamique [0.,range[.
+       Si "vmax" different de "range", c'est la borne superieure
+       qui peut etre atteinte
+       (si elle est depassee, on soustrait "range").
+\verbatim
+  r>0 vmax>0
+  r=24. vmax=24.   -> mets dans [  0,+24[ borne sup exclue
+  r=24. vmax=12.   -> mets dans ]-12,+12] borne inf exclue
+\endverbatim
+*/
+inline void InRange(double *val,double range)
+             {*val-=range*floor(*val/range);}
+inline void InRange(double *val,double range,double vmax)
+         {InRange(val,range); if(*val>vmax) *val-=range;}
+
 #endif