Changeset 3938 in Sophya for trunk/AddOn/TAcq/brfitsrd.cc
- Timestamp:
- Jan 13, 2011, 11:42:21 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/AddOn/TAcq/brfitsrd.cc
r3909 r3938 57 57 // Extraction de qques parametres utiles depuis les fichiers FITS 58 58 string fkvs; 59 cpaqdeltatime_=0.; 59 60 if (fib==0) { 60 61 fkvs=mff_[fib].GetKeyValue("DATEOBS"); 61 if (fkvs.length()>0) { 62 cdateobs_.Set(fkvs); 63 cout << " BRMultiFitsReader/First file (for fiber 0) DATEOBS=" << fkvs << endl; 64 } 65 } 66 fkvs=mff_[fib].GetKeyValue("DATEOBS"); 62 if (fkvs.length()>0) cdateobs_.Set(fkvs); 63 fkvs=mff_[fib].GetKeyValue("TMSTART"); 64 if (fkvs.length()>0) { 65 ctmstart_.Set(fkvs); 66 cout << " BRMultiFitsReader/First file (for fiber 0) TMSTART=" << fkvs << endl; 67 fkvs=mff_[fib].GetKeyValue("TMEND"); 68 SOPHYA::TimeStamp tmend_=ctmstart_; 69 if (fkvs.length()>0) tmend_.Set(fkvs); 70 cpaqdeltatime_=((double)(tmend_.DaysPart()-cdateobs_.DaysPart())*86400.+ 71 (tmend_.SecondsPart()-cdateobs_.SecondsPart()))/(double)mff_[fib].NAxis2(); 72 } 73 } 74 fkvs=mff_[fib].GetKeyValue("FIBERNUM"); 67 75 memgr_.FiberId(fib) = atoi( fkvs.c_str() ); 68 76 … … 222 230 string fkvs=mff_[fib].GetKeyValue("DATEOBS"); 223 231 if (fkvs.length()>0) cdateobs_.Set(fkvs); 232 fkvs=mff_[fib].GetKeyValue("TMSTART"); 233 if (fkvs.length()>0) { 234 ctmstart_.Set(fkvs); 235 cout << " BRMultiFitsReader/First file (for fiber 0) TMSTART=" << fkvs << endl; 236 fkvs=mff_[fib].GetKeyValue("TMEND"); 237 SOPHYA::TimeStamp tmend_=ctmstart_; 238 if (fkvs.length()>0) tmend_.Set(fkvs); 239 cpaqdeltatime_=((double)(tmend_.DaysPart()-cdateobs_.DaysPart())*86400.+ 240 (tmend_.SecondsPart()-cdateobs_.SecondsPart()))/(double)mff_[fib].NAxis2(); 241 } 224 242 } 225 243 vfpos_[fib]=0; … … 254 272 for(size_t fib=0; fib<(size_t)memgr_.NbFibres(); fib++) 255 273 mmbufib_[fib]=memgr_.GetMemZone(mid_,fib); 274 // Definition temps pour la zone a remplir 275 memgr_.GetAuxData(mid_)->FillTime().Set(ctmstart_.ToDays()+cpaqdeltatime_*(double)vfpos_[0]/86400.); 256 276 return false; 257 277 }
Note:
See TracChangeset
for help on using the changeset viewer.