Changeset 2361 in Sophya for trunk/ArchTOIPipe/Kernel/fitsringwtr.cc
- Timestamp:
- Apr 17, 2003, 3:47:34 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/ArchTOIPipe/Kernel/fitsringwtr.cc
r2355 r2361 4 4 // Christophe Magneville 5 5 // Reza Ansari 6 // $Id: fitsringwtr.cc,v 1. 6 2003-03-28 16:06:40cecile Exp $6 // $Id: fitsringwtr.cc,v 1.7 2003-04-17 13:47:34 cecile Exp $ 7 7 8 8 #include "fitsringwtr.h" … … 40 40 41 41 void FITSRingWriter::init() { 42 declareInput("CircTheta");43 declareInput("CircPhi");44 42 declareInput("signal"); 45 43 } … … 59 57 fits_write_key(fptr, TLONG, "NCIRCLES", &nRings, "Number of rings", &fstatus); 60 58 fits_unlock(); 61 59 62 60 char* ttype = "DATA"; 63 61 char* tform = "1D"; … … 67 65 sprintf(xname, "CRING_%d", ring); 68 66 69 double theta = getData(0, ring*nSamples); 70 double phi = getData(1, ring*nSamples); 67 wontNeedBefore(ring*(nSamples+2)); 68 double theta = getData(0, ring*(nSamples+2)+0); 69 double phi = getData(0, ring*(nSamples+2)+1); 71 70 72 71 fits_lock(); … … 79 78 fits_write_key(fptr, TLONG, "NSAMPLES", &nSamples, "Number of samples on the circle", &fstatus); 80 79 fits_unlock(); 81 wontNeedBefore(ring*nSamples);82 80 81 cout << nSamples << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << endl; 82 83 83 for (int i=0; i<nSamples; i++) { 84 double data = getData(2, ring*nSamples+i); 84 double data = getData(0, ring*(nSamples+2)+2+i); 85 if (i > 1000) cout << i << " " << data << endl; 85 86 fits_lock(); 86 87 fits_write_col(fptr, TDOUBLE, 1,i+1,1,1, &data, &fstatus);
Note:
See TracChangeset
for help on using the changeset viewer.