/[PAMELA software]/DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp
ViewVC logotype

Diff of /DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.8 by mocchiut, Tue Jan 23 15:36:45 2007 UTC revision 1.9 by mocchiut, Wed Jan 24 10:08:58 2007 UTC
# Line 525  int OrbitalInfoCore(UInt_t run, TFile *f Line 525  int OrbitalInfoCore(UInt_t run, TFile *f
525        float jyear=0;        float jyear=0;
526    
527        if(atime >= gltle->GetToTime()) {        if(atime >= gltle->GetToTime()) {
528          gltle->Query(atime, dbc);          if ( !gltle->Query(atime, dbc) ){
529              
530          // Compute the magnetic dipole moment.            // Compute the magnetic dipole moment.
531          UInt_t year, month, day, hour, min, sec;            UInt_t year, month, day, hour, min, sec;
532              
533          TTimeStamp t = TTimeStamp(atime, kTRUE);            TTimeStamp t = TTimeStamp(atime, kTRUE);
534          t.GetDate(kTRUE, 0, &year, &month, &day);            t.GetDate(kTRUE, 0, &year, &month, &day);
535          t.GetTime(kTRUE, 0, &hour, &min, &sec);            t.GetTime(kTRUE, 0, &hour, &min, &sec);
536          jyear = (float) year            jyear = (float) year
537            + (month*31.+ (float) day)/365.              + (month*31.+ (float) day)/365.
538            + (hour*3600.+min*60.+(float)sec)/(24*3600*365.);              + (hour*3600.+min*60.+(float)sec)/(24*3600*365.);
539              
540          feldcof_(&jyear, &dimo); // get dipole moment for year            feldcof_(&jyear, &dimo); // get dipole moment for year
541            } else {
542              code = -56;
543              goto closeandexit;
544            };
545        }        }
546    
547        // Propagate the orbit from the tle time to atime, using SGP(D)4.        // Propagate the orbit from the tle time to atime, using SGP(D)4.
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23