00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041 #include "qsastime_extra.h"
00042
00043 static double MJDtoJD = 2400000.5;
00044 static double SecInDay = 86400;
00045
00046 int setFromISOstring( const char* ISOstring, MJDtime *MJD, int forceJulian )
00047 {
00048 double seconds;
00049 int y, m, d, h, min;
00050 int startAt = 0;
00051 int len = strlen( ISOstring );
00052
00053
00054
00055
00056
00057 y = strtol( &( ISOstring[startAt] ), NULL, 10 );
00058 if ( ISOstring[startAt] == '-' || ISOstring[startAt] == '+' )
00059 startAt++;
00060 startAt += 5;
00061 if ( startAt > len )
00062 return 1;
00063
00064 m = strtol( &( ISOstring[startAt] ), NULL, 10 );
00065 startAt += 3;
00066 if ( startAt > len )
00067 return 1;
00068
00069 d = strtol( &( ISOstring[startAt] ), NULL, 10 );
00070 startAt += 3;
00071 if ( startAt > len )
00072 return 1;
00073
00074 h = strtol( &( ISOstring[startAt] ), NULL, 10 );
00075 startAt += 3;
00076 if ( startAt > len )
00077 return 1;
00078
00079 min = strtol( &( ISOstring[startAt] ), NULL, 10 );
00080 startAt += 3;
00081 if ( startAt > len )
00082 return 1;
00083
00084 seconds = strtod( &( ISOstring[startAt] ), NULL );
00085 setFromUT( y, m - 1, d, h, min, seconds, MJD, forceJulian );
00086
00087 return 0;
00088 }
00089
00090
00091 void setFromDOY( int year, int doy, int hour, int min, double sec, MJDtime *MJD, int forceJulian )
00092 {
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102 int leaps, lastyear, extraDays;
00103
00104
00105
00106
00107
00108
00109 if ( forceJulian && year <= 0 )
00110 {
00111
00112
00113
00114
00115 leaps = ( year - 4 ) / 4;
00116 MJD->base_day = year * 365 + leaps + doy - 678943;
00117 }
00118 else if ( forceJulian )
00119 {
00120
00121
00122
00123 leaps = ( year - 1 ) / 4;
00124 MJD->base_day = year * 365 + leaps + doy - 678943;
00125 }
00126 else
00127 {
00128
00129
00130
00131
00132 lastyear = year - 1;
00133 leaps = lastyear / 4 - lastyear / 100 + lastyear / 400;
00134 MJD->base_day = year * 365 + leaps + doy - 678941;
00135 }
00136
00137 MJD->time_sec = sec + ( (double) min + (double) hour * 60. ) * 60.;
00138
00139 if ( MJD->time_sec >= SecInDay )
00140 {
00141 extraDays = (int) ( MJD->time_sec / SecInDay );
00142 MJD->base_day += extraDays;
00143 MJD->time_sec -= extraDays * SecInDay;
00144 }
00145
00146 return;
00147 }
00148
00149
00150 void setFromBCE( int yearBCE, int month, int day, int hour, int min, double sec, MJDtime *MJD, int forceJulian )
00151 {
00152
00153
00154 int year = 1 - yearBCE;
00155 setFromUT( year, month, day, hour, min, sec, MJD, forceJulian );
00156 }
00157
00158 void setFromMJD( double ModifiedJulianDate, MJDtime *MJD )
00159 {
00160
00161 MJD->base_day = (int) ModifiedJulianDate;
00162 MJD->time_sec = ( ModifiedJulianDate - MJD->base_day ) * SecInDay;
00163 }
00164
00165 void setFromJD( double JulianDate, MJDtime *MJD )
00166 {
00167
00168
00169
00170 MJD->base_day = (int) ( JulianDate - MJDtoJD );
00171 MJD->time_sec = ( JulianDate - MJDtoJD - (double) MJD->base_day ) * SecInDay;
00172 }
00173
00174 void setFromCDFepoch( double cdfepoch, MJDtime *MJD )
00175 {
00176
00177
00178
00179 double seconds = cdfepoch * 0.001;
00180
00181 MJD->base_day = (int) ( seconds / 86400.0 );
00182 MJD->time_sec = seconds - MJD->base_day * SecInDay;
00183 MJD->base_day -= 678941;
00184 }
00185
00186 double getCDFepoch( MJDtime *MJD )
00187 {
00188
00189
00190
00191 int days = MJD->base_day + 678941;
00192 double seconds = days * SecInDay + MJD->time_sec;
00193 return seconds * 1000.;
00194 }
00195
00196 double getMJD( MJDtime *MJD )
00197 {
00198
00199 return (double) MJD->base_day + MJD->time_sec / SecInDay;
00200 }
00201
00202 double getJD( MJDtime *MJD )
00203 {
00204
00205 double JD = getMJD( MJD ) + MJDtoJD;
00206 return JD;
00207 }
00208
00209 double getDiffDays( MJDtime *MJD1, MJDtime *MJD2 )
00210 {
00211
00212 double diff = (double) ( MJD1->base_day - MJD2->base_day ) + ( MJD1->time_sec - MJD2->time_sec ) / SecInDay;
00213 return diff;
00214 }
00215
00216 double getDiffSecs( MJDtime *MJD1, MJDtime *MJD2 )
00217 {
00218
00219 double diff = (double) ( MJD1->base_day - MJD2->base_day ) * SecInDay + ( MJD1->time_sec - MJD2->time_sec );
00220 return diff;
00221 }
00222
00223 const char * getISOString( MJDtime* MJD, int delim, int forceJulian )
00224 {
00225
00226
00227
00228
00229
00230
00231
00232 static char DateTime[50];
00233 int y, m, d, hour, min;
00234 int sec1, ysign;
00235 double sec;
00236 int slen;
00237 char * ptr;
00238
00239 breakDownMJD( &y, &m, &d, &hour, &min, &sec, MJD, forceJulian );
00240
00241 if ( y < 0 )
00242 {
00243 ysign = 1;
00244 y = -y;
00245 }
00246 else
00247 ysign = 0;
00248
00249 sec1 = (int) sec / 10;
00250 sec -= (double) sec1 * 10;
00251
00252 if ( delim == 1 )
00253 {
00254 if ( ysign == 0 )
00255 sprintf( DateTime, "%04d-%02d-%02dT%02d:%02d:%01d%-11.10f", y, m + 1, d, hour, min, sec1, sec );
00256 else
00257 sprintf( DateTime, "-%04d-%02d-%02dT%02d:%02d:%01d%-11.10f", y, m + 1, d, hour, min, sec1, sec );
00258
00259
00260 while ( ( ptr = strrchr( &( DateTime[0] ), ' ' ) ) != NULL )
00261 ptr[0] = '\0';
00262 strcat( &( DateTime[0] ), "Z" );
00263 }
00264 else
00265 {
00266 if ( ysign == 0 )
00267 sprintf( DateTime, "%04d-%02d-%02d %02d:%02d:%01d%-11.10f", y, m + 1, d, hour, min, sec1, sec );
00268 else
00269 sprintf( DateTime, "-%04d-%02d-%02d %02d:%02d:%01d%-11.10f", y, m + 1, d, hour, min, sec1, sec );
00270
00271
00272 slen = strlen( DateTime ) - 1;
00273 while ( DateTime[slen] == ' ' )
00274 {
00275 DateTime[slen] = '\0';
00276 slen--;
00277 }
00278 }
00279 return &( DateTime[0] );
00280 }