Ignore:
Timestamp:
Feb 11, 2003, 3:26:40 PM (23 years ago)
Author:
cecile
Message:

nsamples ? wontneedbefore ?

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ArchTOIPipe/Kernel/fitsringwtr.cc

    r2281 r2320  
    44//                               Christophe Magneville
    55//                               Reza Ansari
    6 // $Id: fitsringwtr.cc,v 1.1 2002-11-28 14:41:12 aubourg Exp $
     6// $Id: fitsringwtr.cc,v 1.2 2003-02-11 14:26:40 cecile Exp $
    77
    88#include "fitsringwtr.h"
     
    2121  fits_unlock();
    2222  name = "ringwriter";
    23 }
     23  cout << "FITSRINGWriter:: FileName=" << name << endl;
     24 }
    2425
    2526FITSRingWriter::~FITSRingWriter() {
     
    6465    char xname[100];
    6566    sprintf(xname, "CRING_%d", ring);
     67
    6668    fits_lock();
    6769    fits_create_tbl(fptr,BINARY_TBL,nSamples,1,&ttype,&tform,NULL,xname,&fstatus);
    6870    strcpy(buf, "RINGDATA");
    6971    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);
    7173    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);
    7376    fits_write_key(fptr, TDOUBLE, "CIRPHI", &phi, "Phi angle in radians", &fstatus);
    7477    fits_write_key(fptr, TDOUBLE, "CIRAPER", &aperture, "Aperture angle in radians", &fstatus);
    7578    fits_write_key(fptr, TLONG, "NSAMPLES", &nSamples, "Number of samples on the circle", &fstatus);
    7679    fits_unlock();
     80
    7781    for (int i=0; i<nSamples; i++) {
    7882      double data = getData(2, ring*nSamples+i);
     83      wontNeedBefore(ring*nSamples+i);
    7984      fits_lock();
    8085      fits_write_col(fptr, TDOUBLE, 1,i,1,1, &data, &fstatus);
    8186      fits_unlock();
    8287    }
     88
    8389    fits_lock();
    8490    fits_close_file(fptr, &fstatus);
Note: See TracChangeset for help on using the changeset viewer.