Changeset 14885 for trunk/FACT++
- Timestamp:
- 02/09/13 16:35:05 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FACT++/src/makedata.cc
r14784 r14885 86 86 if (!conf.DoParse(argc, argv)) 87 87 return 127; 88 88 /* 89 89 if (!conf.Has("source-name")) 90 90 { … … 92 92 return 1; 93 93 } 94 94 */ 95 95 // ------------------ Eval config --------------------- 96 96 … … 106 106 time.SetFromStr(conf.Get<string>("date-time")); 107 107 108 const string source_name = conf.Get<string>("source-name");109 108 const double max_current = conf.Get<double>("max-current"); 110 109 const double max_zd = conf.Get<double>("max-zd"); … … 128 127 cout << jd1 << "\n"; 129 128 130 const string fDatabase = conf.Get<string>("source-database"); 129 if (!conf.Has("source-name")) 130 return 1; 131 132 const string source_name = conf.Get<string>("source-name"); 133 const string fDatabase = conf.Get<string>("source-database"); 131 134 132 135 // ------------------ Precalc moon coord --------------------- … … 181 184 182 185 // Current prediction 183 const double lc = angle*hrzm.alt*pow(disk, 6)/360/360; 184 const double cur = lc>0 ? 7.7+4942*lc : 7.7; 186 const double cang = sin(angle *M_PI/180); 187 const double calt = sin(hrzm.alt*M_PI/180); 188 189 const double lc = cang==0 ? -1 : calt*pow(disk, 3)/sqrt(cang); 190 const double cur = lc>0 ? 8.1+94.6*lc : 8.1; 185 191 186 192 // Relative energy threshold prediction … … 205 211 206 212 if (no_limits || (cur<max_current && 90-hrz.alt<max_zd)) 207 cout << ratio*cur/ 7.7;213 cout << ratio*cur/8.1; 208 214 cout << ", "; 209 215
Note:
See TracChangeset
for help on using the changeset viewer.