Changeset 17956 for trunk/Mars/fact


Ignore:
Timestamp:
08/13/14 12:34:21 (10 years ago)
Author:
tbretz
Message:
Replaced current prediction by the new central current prediction.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/fact/plots/quality.C

    r17664 r17956  
    9595    Double_t jd = time + 40587 + 2400000.5;
    9696
    97     // Sun properties
    98     Nova::EquPosn  sun  = Nova::GetSolarEquCoords(jd);
    99     Nova::ZdAzPosn hrzs = Nova::GetHrzFromEqu(sun, jd);
    100 
    10197    // Get source position
    10298    Nova::EquPosn pos = FindPointing(time);
    10399
    104     // Moon properties
    105     Nova::EquPosn moon = Nova::GetLunarEquCoords(jd, 0.01);
    106     Nova::HrzPosn hrzm = Nova::GetHrzFromEqu(moon, jd);
    107     double        disk = Nova::GetLunarDisk(jd);
    108 
    109     // Derived moon properties
    110     double angle = Nova::GetAngularSeparation(moon, pos);
    111     double edist = Nova::GetLunarEarthDist(jd)/384400;
    112 
    113     // Current prediction
    114     double sin_malt  = hrzm.alt<0 ? 0 : sin(hrzm.alt*TMath::DegToRad());
    115     double cos_mdist = cos(angle*TMath::DegToRad());
    116     double sin_szd   = sin(hrzs.zd*TMath::DegToRad());
    117 
    118     double c0 = pow(disk,      2.63);
    119     double c1 = pow(sin_malt,  0.60);
    120     double c2 = pow(edist,    -2.00);
    121     double c3 = exp(0.67*cos_mdist*cos_mdist*cos_mdist*cos_mdist);
    122     double c4 = exp(-97.8+105.8*sin_szd*sin_szd);
    123 
    124     double cur = 6.2 + 95.7*c0*c1*c2*c3 + c4;
    125 
    126     return cur;
     100    return FACT::PredictI(jd, pos);
    127101}
    128102
Note: See TracChangeset for help on using the changeset viewer.