Changeset 17956 for trunk/Mars/fact
- Timestamp:
- 08/13/14 12:34:21 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Mars/fact/plots/quality.C
r17664 r17956 95 95 Double_t jd = time + 40587 + 2400000.5; 96 96 97 // Sun properties98 Nova::EquPosn sun = Nova::GetSolarEquCoords(jd);99 Nova::ZdAzPosn hrzs = Nova::GetHrzFromEqu(sun, jd);100 101 97 // Get source position 102 98 Nova::EquPosn pos = FindPointing(time); 103 99 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); 127 101 } 128 102
Note:
See TracChangeset
for help on using the changeset viewer.