Changeset 2320 in Sophya
- Timestamp:
- Feb 11, 2003, 3:26:40 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/ArchTOIPipe/Kernel/fitsringwtr.cc
r2281 r2320 4 4 // Christophe Magneville 5 5 // Reza Ansari 6 // $Id: fitsringwtr.cc,v 1. 1 2002-11-28 14:41:12 aubourgExp $6 // $Id: fitsringwtr.cc,v 1.2 2003-02-11 14:26:40 cecile Exp $ 7 7 8 8 #include "fitsringwtr.h" … … 21 21 fits_unlock(); 22 22 name = "ringwriter"; 23 } 23 cout << "FITSRINGWriter:: FileName=" << name << endl; 24 } 24 25 25 26 FITSRingWriter::~FITSRingWriter() { … … 64 65 char xname[100]; 65 66 sprintf(xname, "CRING_%d", ring); 67 66 68 fits_lock(); 67 69 fits_create_tbl(fptr,BINARY_TBL,nSamples,1,&ttype,&tform,NULL,xname,&fstatus); 68 70 strcpy(buf, "RINGDATA"); 69 71 fits_write_key(fptr, TSTRING, "PDMTYPE", &buf, "Planck Data Model Type", &fstatus); 70 double theta = getData(0, ring *nSamples);72 double theta = getData(0, ring);//*nSamples); 71 73 fits_write_key(fptr, TDOUBLE, "CIRTHETA", &theta, "Theta angle in radians", &fstatus); 72 double phi = getData(1, ring*nSamples); 74 double phi = getData(1, ring);//*nSamples); 75 wontNeedBefore(ring);//*nSamples); 73 76 fits_write_key(fptr, TDOUBLE, "CIRPHI", &phi, "Phi angle in radians", &fstatus); 74 77 fits_write_key(fptr, TDOUBLE, "CIRAPER", &aperture, "Aperture angle in radians", &fstatus); 75 78 fits_write_key(fptr, TLONG, "NSAMPLES", &nSamples, "Number of samples on the circle", &fstatus); 76 79 fits_unlock(); 80 77 81 for (int i=0; i<nSamples; i++) { 78 82 double data = getData(2, ring*nSamples+i); 83 wontNeedBefore(ring*nSamples+i); 79 84 fits_lock(); 80 85 fits_write_col(fptr, TDOUBLE, 1,i,1,1, &data, &fstatus); 81 86 fits_unlock(); 82 87 } 88 83 89 fits_lock(); 84 90 fits_close_file(fptr, &fstatus);
Note:
See TracChangeset
for help on using the changeset viewer.