source: trunk/external-alma/atnf/PKSIO/NROFITSDataset.cc @ 2783

Last change on this file since 2783 was 2783, checked in by Takeshi Nakazato, 11 years ago

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No?

Module(s): Module Names change impacts.

Description: Describe your changes here...

more refactoring on NRO filler.
File size: 63.5 KB
Line 
1//#---------------------------------------------------------------------------
2//# NROFITSDataset.cc: Class for NRO 45m FITS dataset.
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2000-2006
5//# Associated Universities, Inc. Washington DC, USA.
6//#
7//# This library is free software; you can redistribute it and/or modify it
8//# under the terms of the GNU Library General Public License as published by
9//# the Free Software Foundation; either version 2 of the License, or (at your
10//# option) any later version.
11//#
12//# This library is distributed in the hope that it will be useful, but WITHOUT
13//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14//# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
15//# License for more details.
16//#
17//# You should have received a copy of the GNU Library General Public License
18//# along with this library; if not, write to the Free Software Foundation,
19//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
20//#
21//# Correspondence concerning AIPS++ should be addressed as follows:
22//#        Internet email: aips2-request@nrao.edu.
23//#        Postal address: AIPS++ Project Office
24//#                        National Radio Astronomy Observatory
25//#                        520 Edgemont Road
26//#                        Charlottesville, VA 22903-2475 USA
27//#
28//# $Id$
29//#---------------------------------------------------------------------------
30//# Original: 2009/02/27, Takeshi Nakazato, NAOJ
31//#---------------------------------------------------------------------------
32
33#include <atnf/PKSIO/NROFITSDataset.h>
34#include <scimath/Mathematics/InterpolateArray1D.h>
35
36#include <iostream>
37#include <fstream>
38#include <casa/math.h>
39#include <casa/iomanip.h>
40
41using namespace std ;
42
43//header size (8*2880bytes)
44#define FITS_HEADER_SIZE 23040
45
46// constructor
47NROFITSDataset::NROFITSDataset( string name )
48  : NRODataset( name )
49{}
50
51// destructor
52NROFITSDataset::~NROFITSDataset()
53{
54  // close file
55  close() ;
56}
57
58// data initialization
59void NROFITSDataset::initialize()
60{
61  LogIO os( LogOrigin( "NROFITSDataset", "initialize()", WHERE ) ) ;
62
63  // open file
64  if ( open() )
65    os << LogIO::SEVERE << "error while opening file " << filename_ << LogIO::EXCEPTION ;
66
67  // field names, units, and sizes
68  readHeader( numField_, "TFIELDS", same_ ) ;
69  names_.resize( numField_ ) ;
70  units_.resize( numField_ ) ;
71  sizes_.resize( numField_ ) ;
72  offsets_.resize( numField_ ) ;
73  getField() ;
74
75  // check endian
76  // FITS file is always BIG_ENDIAN, but it is not true for NRO data
77  vector<int> itmp( 6 ) ;
78  if ( readTable( itmp, "LAVST", true, 0 ) != 0 ) {
79    os << LogIO::WARN << "Error while checking endian." << LogIO::POST ;
80    return ;
81  }
82//   os << "itmp = " << itmp[0] << " "
83//      << itmp[1] << " " << itmp[2] << " "
84//      << itmp[3] << " " << itmp[4] << " "
85//      << itmp[5] << LogIO::POST ;
86  // check endian by looking month value in LAVST (start time)
87  if ( itmp[1] > 0 && itmp[1] < 13 ) {
88    same_ = 1 ;
89    os << LogIO::NORMAL << "same endian " << LogIO::POST ;
90  }
91  else {
92    same_ = 0 ;
93    os << LogIO::NORMAL << "different endian " << LogIO::POST ;
94  }
95
96  // initialization
97  int status = 0 ;
98  status = readHeader( ARYNM, "ARYNM", same_ ) ;
99  if ( status != 0 )
100    ARYNM = 1 ;
101  readHeader( rowNum_, "NAXIS2", same_ ) ;
102  readHeader( scanLen_, "NAXIS1", same_ ) ;
103  status = 0 ;
104  scanNum_ = rowNum_ / ARYNM ;
105  int nchan = 0 ;
106  if ( readTable( nchan, "NCH", same_ ) != 0 ) {
107    os << LogIO::WARN << "Error while checking maximum channel number." << LogIO::POST ;
108    return ;
109  }
110  chmax_ = nchan ;
111  datasize_ = sizeof( int ) * chmax_ ;
112  //record_->JDATA.resize( chmax_ ) ;
113  JDATA.resize( chmax_ ) ;
114  // zero clear
115  for ( uInt i = 0 ; i < JDATA.size() ; i++ )
116    JDATA[i] = 0 ;
117
118  RX.resize( ARYNM ) ;
119  HPBW.resize( ARYNM ) ;
120  EFFA.resize( ARYNM ) ;
121  EFFB.resize( ARYNM ) ;
122  EFFL.resize( ARYNM ) ;
123  EFSS.resize( ARYNM ) ;
124  GAIN.resize( ARYNM ) ;
125  HORN.resize( ARYNM ) ;
126  POLTP.resize( ARYNM ) ;
127  POLDR.resize( ARYNM ) ;
128  POLAN.resize( ARYNM ) ;
129  DFRQ.resize( ARYNM ) ;
130  SIDBD.resize( ARYNM ) ;
131  REFN.resize( ARYNM ) ;
132  IPINT.resize( ARYNM ) ;
133  MULTN.resize( ARYNM ) ;
134  MLTSCF.resize( ARYNM ) ;
135  LAGWIND.resize( ARYNM ) ;
136  BEBW.resize( ARYNM ) ;
137  BERES.resize( ARYNM ) ;
138  CHWID.resize( ARYNM ) ;
139  ARRY.resize( NRO_FITS_ARYMAX ) ; 
140  NFCAL.resize( ARYNM ) ;
141  F0CAL.resize( ARYNM ) ;
142  FQCAL.resize( ARYNM ) ;
143  CHCAL.resize( ARYNM ) ;
144  CWCAL.resize( ARYNM ) ;
145  DSBFC.resize( ARYNM ) ;
146
147  ARYTP.resize( ARYNM ) ;
148  arrayid_.resize( ARYNM ) ;
149  for ( int i = 0 ; i < ARYNM ; i++ )
150    arrayid_[i] = -1 ;
151
152  for ( int i = 0 ; i < ARYNM ; i++ ) {
153    FQCAL[i].resize( 10 ) ;
154    CHCAL[i].resize( 10 ) ;
155    CWCAL[i].resize( 10 ) ;
156  }
157
158  datasize_ += sizeof( char ) * ARYNM * 16 // RX
159    + sizeof( double ) * ARYNM * 6         // HPBW, EFFA, EFFB, EFFL, EFSS GAIN
160    + sizeof( char ) * ARYNM * 4           // HORN
161    + sizeof( char ) * ARYNM * 4           // POLTP
162    + sizeof( double ) * ARYNM * 3         // POLDR, POLAN, DFRQ
163    + sizeof( char ) * ARYNM * 4           // SIDBID
164    + sizeof( int ) * ARYNM * 3            // REFN, IPINT, MULTN
165    + sizeof( double ) * ARYNM             // MLTSCF
166    + sizeof( char ) * ARYNM * 8           // LAGWIND
167    + sizeof( double ) * ARYNM * 3         // BEBW, BERES, CHWID
168    + sizeof( int ) * NRO_FITS_ARYMAX      // ARRY
169    + sizeof( int ) * ARYNM                // NFCAL
170    + sizeof( double ) * ARYNM             // F0CAL
171    + sizeof( double ) * ARYNM * 10 * 3    // FQCAL, CHCAL, CWCAL
172    + sizeof( char ) * 180 ;               // CDMY1
173
174  refFreq_.resize( ARYNM, 0.0 ) ;
175
176  // NRODataRecord
177  record_ = new NRODataRecord() ;
178  record_->LDATA = NULL ;
179}
180
181int NROFITSDataset::fillHeader( int sameEndian )
182{
183  LogIO os( LogOrigin( "NROFITSDataset", "fillHeader()", WHERE ) ) ;
184
185  // fill array type
186  fillARYTP() ;
187
188  // read data header
189  float ftmp = 0.0 ;
190  if ( readHeader( LOFIL, "LOFIL" ) != 0 ) {
191    os << LogIO::NORMAL << "LOFIL set to 'FITS'." << LogIO::POST ;
192    LOFIL = "FITS" ;
193  }
194  // DEBUG
195  //cout << "LOFIL = \'" << LOFIL << "\'" << endl ;
196  //
197  if ( readHeader( VER, "VER" ) != 0 ) {
198    if ( readHeader( VER, "HISTORY NEWSTAR VER" ) != 0 ) {
199      os << LogIO::NORMAL << "VER set to 'V000'." << LogIO::POST ;
200      VER = "V000" ;
201    }
202  }
203  // DEBUG
204  //cout << "VER = \'" << VER << "\'" << endl ;
205  //
206  if ( readHeader( GROUP, "GROUP" ) != 0 ) {
207    if ( readHeader( GROUP, "HISTORY NEWSTAR GROUP" ) != 0 ) {
208      os << LogIO::NORMAL << "GROUP set to 'GRP0'." << LogIO::POST ;
209      GROUP = "GROUP0" ;
210    }
211  }
212  // DEBUG
213  //cout << "GROUP = \'" << GROUP << "\'" << endl ;
214  //
215  if ( readHeader( PROJ, "PROJECT" ) != 0 ) {
216    if ( readHeader( PROJ, "HISTORY NEWSTAR PROJECT" ) != 0 ) {
217      os << LogIO::NORMAL << "PROJ set to 'PROJ0'." << LogIO::POST ;
218      PROJ = "PROJECT0" ;
219    }
220  }
221  // DEBUG
222  //cout << "PROJ = \'" << PROJ << "\'" << endl ;
223  //
224  if ( readHeader( SCHED, "SCHED" ) != 0 ) {
225    if ( readHeader( SCHED, "HISTORY NEWSTAR SCHED" ) != 0 ) {
226      os << LogIO::NORMAL << "SCHED set to 'SCHED0'." << LogIO::POST ;
227      SCHED = "SCHED0" ;
228    }
229  }
230  // DEBUG
231  //cout << "SCHED = \'" << SCHED << "\'" << endl ;
232  //
233  if ( readHeader( OBSVR, "OBSERVER" ) != 0 ) {
234    os << LogIO::NORMAL << "OBSVR set to 'SOMEONE'" << LogIO::POST ;
235    OBSVR = "SOMEONE" ;
236  }
237  // DEBUG
238  //cout << "OBSVR = \'" << OBSVR << "\'" << endl ;
239  //
240  if ( readHeader( LOSTM, "STRSC" ) != 0 ) {
241    os << LogIO::WARN << "Error while reading data LOSTM." << LogIO::POST ;
242    return -1 ;
243  }
244  if ( LOSTM[0] == '9' ) {
245    LOSTM = "19" + LOSTM ;
246  }
247  else if ( LOSTM[0] == '0') {
248    LOSTM = "20" + LOSTM ;
249  }
250  // DEBUG
251  //cout << "LOSTM = \'" << LOSTM << "\'" << endl ;
252  //
253  if ( readHeader( LOETM, "STPSC" ) != 0 ) {
254    os << LogIO::WARN << "Error while reading data LOETM." << LogIO::POST ;
255    return -1 ;
256  }
257  if ( LOETM[0] == '9' ) {
258    LOETM = "19" + LOETM ;
259  }
260  else if ( LOETM[0] == '0') {
261    LOETM = "20" + LOETM ;
262  }   
263  // DEBUG
264  //cout << "LOETM = \'" << LOETM << "\'" << endl ;
265  //
266  if ( readHeader( NSCAN, "NAXIS2", sameEndian ) != 0 ) {
267    os << LogIO::WARN << "Error while reading data NSCAN." << LogIO::POST ;
268    return -1 ;
269  }
270  // DEBUG
271  //cout << "NSCAN = " << NSCAN << endl ;
272  //
273  string subt ;
274  if ( readHeader( TITLE, "TITLE" ) != 0 ) {
275    int stat1 = readHeader( TITLE, "HISTORY NEWSTAR TITLE1" ) ;
276    int stat2 = readHeader( subt, "HISTORY NEWSTAR TITLE2" ) ;
277    if ( stat1 != 0 && stat2 != 0 ) {
278      os << LogIO::NORMAL << "TITLE set to 'NOTITLE'." << LogIO::POST ;
279      TITLE = "NOTITLE" ;
280    }
281    else {
282      //cout << "TITLE = \'" << TITLE << "\'" << endl ;
283      //cout << "subt  = \'" << subt << "\'" << endl ;
284      TITLE = TITLE + subt ;
285    }
286  }
287  // DEBUG
288  //cout << "TITLE = \'" << TITLE << "\'" << endl ;
289  //
290  if ( readHeader( OBJ, "OBJECT" ) != 0 ) {
291    os << LogIO::NORMAL << "OBJ set to 'SOMEWHERE'." << LogIO::POST ;
292    OBJ = "NOOBJ" ;
293  }
294  // DEBUG
295  //cout << "OBJ = \'" << OBJ << "\'" << endl ;
296  //
297  if ( readHeader( ftmp, "EPOCH", sameEndian ) != 0 ) {
298    os << LogIO::WARN << "Error while reading data EPOCH." << LogIO::POST ;
299    return -1 ;
300  }
301  if ( ftmp == 1950.0 )
302    EPOCH = "B1950" ;
303  else if ( ftmp == 2000.0 )
304    EPOCH = "J2000" ;
305  else {
306    os << LogIO::WARN << "Unidentified epoch. set to 'XXXXX'" << LogIO::POST ;
307    EPOCH = "XXXXX" ;
308  }
309  // DEBUG
310  //cout << "EPOCH = \'" << EPOCH << "\'" << endl ;
311  //
312  string stmp ;
313  if ( readHeader( stmp, "RA" ) != 0 ) {
314    os << LogIO::WARN << "Error while reading data RA0." << LogIO::POST ;
315    return -1 ;
316  }
317  RA0 = radRA( stmp.c_str() ) ;
318  // DEBUG
319  //cout << "RA0 = " << RA0 << endl ;
320  //
321  if ( readHeader( stmp, "DEC" ) != 0 ) {
322    os << LogIO::WARN << "Error while reading data DEC0." << LogIO::POST ;
323    return -1 ;
324  }
325  DEC0 = radDEC( stmp.c_str() ) ;
326  // DEBUG
327  //cout << "DEC0 = " << DEC0 << endl ;
328  //
329  if ( readHeader( GLNG0, "GL0", sameEndian ) != 0 ) {
330    os << LogIO::WARN << "Error while reading data GLNG0." << LogIO::POST ;
331    return -1 ;
332  }
333  // DEBUG
334  //cout << "GLNG0 = " << GLNG0 << endl ;
335  //
336  if ( readHeader( GLAT0, "GB0", sameEndian ) != 0 ) {
337    os << LogIO::WARN << "Error while reading data GLAT0." << LogIO::POST ;
338    return -1 ;
339  }
340  // DEBUG
341  //cout << "GLAT0 = " << GLAT0 << endl ;
342  //
343  if ( readHeader( NCALB, "NCALB", sameEndian ) != 0 ) {
344    os << LogIO::WARN << "Error while reading data NCALB." << LogIO::POST ;
345    return -1 ;
346  }
347  // DEBUG
348  //cout << "NCALB = " << NCALB << endl ;
349  //
350  if ( readHeader( SCNCD, "SCNCD", sameEndian ) != 0 ) {
351    os << LogIO::NORMAL << "SCNCD set to 0 (RADEC)." << LogIO::POST ;
352    SCNCD = 0 ;
353  }
354  // DEBUG
355  //cout << "SCNCD = " << SCNCD << endl ;
356  //
357  if ( readHeader( SCMOD, "SCMOD1" ) != 0 ) {
358    os << LogIO::WARN << "Error while reading data SCMOD." << LogIO::POST ;
359    return -1 ;
360  }
361  string::size_type pos = SCMOD.find( ' ' ) ;
362  if ( pos != string::npos ) {
363    SCMOD = SCMOD.substr( 0, pos ) ;
364    SCMOD = SCMOD + " " ;
365  }
366  //cout << "SDMOD1 = \'" << SCMOD << "\'" << endl ;
367  if ( readHeader( stmp, "SCMOD2" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
368    if ( ( pos = stmp.find( ' ' ) ) != string::npos )
369      stmp = stmp.substr( 0, pos ) ;
370    SCMOD = SCMOD + stmp + " " ;
371    //cout << "SCMOD2 = \'" << SCMOD << "\'" << endl ;
372  }
373  if ( readHeader( stmp, "SCMOD3" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
374    if ( ( pos = stmp.find( ' ' ) ) != string::npos )
375      stmp = stmp.substr( 0, pos ) ;
376    SCMOD = SCMOD + stmp + " " ;
377    //cout << "SCMOD3 = \'" << SCMOD << "\'" << endl ;
378  }
379  if ( readHeader( stmp, "SCMOD4" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
380    if ( ( pos = stmp.find( ' ' ) ) != string::npos )
381      stmp = stmp.substr( 0, pos ) ;
382    SCMOD = SCMOD + stmp + " " ;
383    //cout << "SCMOD4 = \'" << SCMOD << "\'" << endl ;
384  }
385  if ( readHeader( stmp, "SCMOD5" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
386    if ( ( pos = stmp.find( ' ' ) ) != string::npos )
387      stmp = stmp.substr( 0, pos ) ;
388    SCMOD = SCMOD + stmp + " " ;
389    //cout << "SCMOD5 = \'" << SCMOD << "\'" << endl ;
390  }
391  if ( readHeader( stmp, "SCMOD6" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
392    if ( ( pos = stmp.find( ' ' ) ) != string::npos )
393      stmp = stmp.substr( 0, pos ) ;
394    SCMOD = SCMOD + stmp + " " ;
395    //cout << "SCMOD6 = \'" << SCMOD << "\'" << endl ;
396  }
397  // DEBUG
398  //cout << "SCMOD = \'" << SCMOD << "\'" << endl ;
399  //
400  if ( readHeader( URVEL, "VEL", sameEndian ) != 0 ) {
401    os << LogIO::WARN << "Error while reading data URVEL." << LogIO::POST ;
402    return -1 ;
403  }
404  // DEBUG
405  //cout << "URVEL = " << URVEL << endl ;
406  //
407  if ( readHeader( VREF, "VREF" ) != 0 ) {
408    os << LogIO::WARN << "Error while reading data VREF." << LogIO::POST ;
409    return -1 ;
410  }
411  // DEBUG
412  //cout << "VREF = \'" << VREF << "\'" << endl ;
413  //
414  if ( readHeader( VDEF, "VDEF" ) != 0 ) {
415    os << LogIO::WARN << "Error while reading data VDEF." << LogIO::POST ;
416    return -1 ;
417  }
418  // DEBUG
419  //cout << "VDEF = \'" << VDEF << "\'" << endl ;
420  //
421  if ( readHeader( SWMOD, "SWMOD" ) != 0 ) {
422    os << LogIO::WARN << "Error while reading data SWMOD." << LogIO::POST ;
423    return -1 ;
424  }
425  // DEBUG
426  //cout << "SWMOD = \'" << SWMOD << "\'" << endl ;
427  //
428  if ( readHeader( FRQSW, "FRQSW", sameEndian ) != 0 ) {
429    os << LogIO::NORMAL << "FRQSW set to 0." << LogIO::POST ;
430    FRQSW = 0.0 ;
431  }
432  // DEBUG
433  //cout << "FRQSW = " << FRQSW << endl ;
434  //
435  if ( readHeader( DBEAM, "DBEAM", sameEndian ) != 0 ) {
436    os << LogIO::WARN << "Error while reading data DBEAM." << LogIO::POST ;
437    return -1 ;
438  }
439  // DEBUG
440  //cout << "DBEAM = " << DBEAM << endl ;
441  //
442  if ( readHeader( MLTOF, "MLTOF", sameEndian ) != 0 ) {
443    os << LogIO::NORMAL << "MLTOF set to 0." << LogIO::POST ;
444    MLTOF = 0.0 ;
445  }
446  // DEBUG
447  //cout << "MLTOF = " << MLTOF << endl ;
448  //
449  if ( readHeader( CMTQ, "CMTQ", sameEndian ) != 0 ) {
450    os << LogIO::WARN << "Error while reading data CMTQ." << LogIO::POST ;
451    return -1 ;
452  }
453  // DEBUG
454  //cout << "CMTQ = " << CMTQ << endl ;
455  //
456  if ( readHeader( CMTE, "CMTE", sameEndian ) != 0 ) {
457    os << LogIO::WARN << "Error while reading data CMTE." << LogIO::POST ;
458    return -1 ;
459  }
460  // DEBUG
461  //cout << "CMTE = " << CMTE << endl ;
462  //
463  if ( readHeader( CMTSOM, "CMTSOM", sameEndian ) != 0 ) {
464    os << LogIO::WARN << "Error while reading data CMTSOM." << LogIO::POST ;
465    return -1 ;
466  }
467  // DEBUG
468  //cout << "CMTSOM = " << CMTSOM << endl ;
469  //
470  if ( readHeader( CMTNODE, "CMTNODE", sameEndian ) != 0 ) {
471    os << LogIO::WARN << "Error while reading data CMTNODE." << LogIO::POST ;
472    return -1 ;
473  }
474  // DEBUG
475  //cout << "CMTNODE = " << CMTNODE << endl ;
476  //
477  if ( readHeader( CMTI, "CMTI", sameEndian ) != 0 ) {
478    os << LogIO::WARN << "Error while reading data CMTI." << LogIO::POST ;
479    return -1 ;
480  }
481  // DEBUG
482  //cout << "CMTI = " << CMTI << endl ;
483  //
484  if ( readHeader( CMTTM, "CMTTM" ) != 0 ) {
485    os << LogIO::WARN << "Error while reading data CMTTM." << LogIO::POST ;
486    return -1 ;
487  }
488  // DEBUG
489  //cout << "CMTTM = \'" << CMTTM << "\'" << endl ;
490  //
491  if ( readHeader( SBDX, "SDBX", sameEndian ) != 0 ) {
492    os << LogIO::WARN << "Error while reading data SBDX." << LogIO::POST ;
493    return -1 ;
494  }
495  // DEBUG
496  //cout << "SBDX = " << SBDX << endl ;
497  //
498  if ( readHeader( SBDY, "SDBY", sameEndian ) != 0 ) {
499    os << LogIO::WARN << "Error while reading data SBDY." << LogIO::POST ;
500    return -1 ;
501  }
502  // DEBUG
503  //cout << "SBDY = " << SBDY << endl ;
504  //
505  if ( readHeader( SBDZ1, "SDBZ1", sameEndian ) != 0 ) {
506    os << LogIO::WARN << "Error while reading data SBDZ1." << LogIO::POST ;
507    return -1 ;
508  }
509  // DEBUG
510  //cout << "SBDZ1 = " << SBDZ1 << endl ;
511  //
512  if ( readHeader( SBDZ2, "SDBZ2", sameEndian ) != 0 ) {
513    os << LogIO::WARN << "Error while reading data SBDZ2." << LogIO::POST ;
514    return -1 ;
515  }
516  // DEBUG
517  //cout << "SBDZ2 = " << SBDZ2 << endl ;
518  //
519  if ( readHeader( DAZP, "DAZP", sameEndian ) != 0 ) {
520    os << LogIO::NORMAL << "DAZP set to 0." << LogIO::POST ;
521    DAZP = 0.0 ;
522  }
523  // DEBUG
524  //cout << "DAZP = " << DAZP << endl ;
525  //
526  if ( readHeader( DELP, "DELP", sameEndian ) != 0 ) {
527    os << LogIO::NORMAL << "DELP set to 0." << LogIO::POST ;
528    DELP = 0.0 ;
529  }
530  // DEBUG
531  //cout << "DELP = " << DELP << endl ;
532  //
533  if ( readHeader( CHBIND, "CHBIND", sameEndian ) != 0 ) {
534    os << LogIO::NORMAL << "CHBIND set to 1." << LogIO::POST ;
535    CHBIND = 1 ;
536  }
537  // DEBUG
538  //cout << "CHBIND = " << CHBIND << endl ;
539  //
540  if ( readHeader( NUMCH, "NCH", sameEndian ) != 0 ) {
541    if ( readTable( NUMCH, "NCH", sameEndian ) != 0 ) {
542      os << LogIO::NORMAL << "NUMCH set to " << chmax_ << "." << LogIO::POST ;
543      NUMCH = chmax_ ;
544    }
545  }
546  // DEBUG
547  //cout << "NUMCH = " << NUMCH << endl ;
548  //
549  if ( readHeader( CHMIN, "CHMIN", sameEndian ) != 0 ) {
550    os << LogIO::NORMAL << "CHMIN set to 1." << LogIO::POST ;
551    CHMIN = 1 ;
552  }
553  // DEBUG
554  //cout << "CHMIN = " << CHMIN << endl ;
555  //
556  if ( readHeader( CHMAX, "CHMAX", sameEndian ) != 0 ) {
557    os << LogIO::NORMAL << "CHMAX set to " << chmax_ << "." << LogIO::POST ;
558    CHMAX = chmax_ ;
559  }
560  // DEBUG
561  //cout << "CHMAX = " << CHMAX << endl ;
562  //
563  if ( readHeader( ALCTM, "ALCTM", sameEndian ) != 0 ) {
564    if ( readTable( ALCTM, "ALCTM", sameEndian ) != 0 ) {
565      os << LogIO::WARN << "Error while reading data ALCTM." << LogIO::POST ;
566      return -1 ;
567    }
568  }
569  // DEBUG
570  //cout << "ALCTM = " << ALCTM << endl ;
571  //
572  int itmp ;
573  if ( readHeader( itmp, "INTEG", sameEndian ) != 0 ) {
574    if ( readTable( itmp, "INTEG", sameEndian ) != 0 ) {
575      os << LogIO::WARN << "Error while reading data IPTIM." << LogIO::POST ;
576      return -1 ;
577    }
578  }
579  IPTIM = (double)itmp ;
580  // DEBUG
581  //cout << "IPTIM = " << IPTIM << endl ;
582  //
583  if ( readHeader( PA, "PA", sameEndian ) != 0 ) {
584    if ( readTable( PA, "PA", sameEndian ) != 0 ) {
585      os << LogIO::WARN << "Error while reading data PA." << LogIO::POST ;
586      return -1 ;
587    }
588  }
589  // DEBUG
590  //cout << "PA = " << PA << endl ;
591  //
592
593  // find data index for each ARYTP
594  findData() ;
595  vector<char *> v( ARYNM ) ;
596  if ( readColumn( RX, "RX" ) != 0 ) {
597    os << LogIO::WARN << "Error while reading data RX." << LogIO::POST ;
598    return -1 ;
599  }
600  // DEBUG
601  //nro_debug_output( "RX", ARYNM, RX ) ;
602  //
603  if ( readColumn( HPBW, "HPBW", sameEndian ) != 0 ) {
604    os << LogIO::WARN << "Error while reading data HPBW." << LogIO::POST ;
605    return -1 ;
606  }
607  // DEBUG
608//   nro_debug_output( "HPBW", ARYNM, HPBW ) ;
609  //
610  if ( readColumn( EFFA, "EFFA", sameEndian ) != 0 ) {
611    os << LogIO::WARN << "Error while reading data EFFA." << LogIO::POST ;
612    return -1 ;
613  }
614  // DEBUG
615//   nro_debug_output( "EFFA", ARYNM, EFFA ) ;
616  //
617  if ( readColumn( EFFB, "EFFB", sameEndian ) != 0 ) {
618    os << LogIO::WARN << "Error while reading data EFFB." << LogIO::POST ;
619    return -1 ;
620  }
621  // DEBUG
622//   nro_debug_output( "EFFB", ARYNM, EFFB ) ;
623  //
624  if ( readColumn( EFFL, "EFFL", sameEndian ) != 0 ) {
625    os << LogIO::WARN << "Error while reading data EFFL." << LogIO::POST ;
626    return -1 ;
627  }
628  // DEBUG
629//   nro_debug_output( "EFFL", ARYNM, EFFL ) ;
630  //
631  if ( readColumn( EFSS, "EFSS", sameEndian ) != 0 ) {
632    os << LogIO::WARN << "Error while reading data EFSS." << LogIO::POST ;
633    return -1 ;
634  }
635  // DEBUG
636//   nro_debug_output( "EFSS", ARYNM, EFSS ) ;
637  //
638  if ( readColumn( GAIN, "GAIN", sameEndian ) != 0 ) {
639    os << LogIO::WARN << "Error while reading data GAIN." << LogIO::POST ;
640    return -1 ;
641  }
642  // DEBUG
643//   nro_debug_output( "GAIN", ARYNM, GAIN ) ;
644  //
645  if ( readColumn( HORN, "HORN" ) != 0 ) {
646    os << LogIO::WARN << "Error while reading data HORN." << LogIO::POST ;
647    return -1 ;
648  }
649  // DEBUG
650//   nro_debug_output( "HORN", ARYNM, HORN ) ;
651  //
652  if ( readColumn( POLTP, "POLTP" ) != 0 ) {
653    os << LogIO::WARN << "Error while reading data POLTP." << LogIO::POST ;
654    return -1 ;
655  }
656  // DEBUG
657//   nro_debug_output( "POLTP", ARYNM, POLTP ) ;
658  //
659  vector<int> ipoldr( ARYNM, 0 ) ;
660  if ( readColumn( ipoldr, "POLDR", sameEndian ) != 0 ) {
661    os << LogIO::WARN << "Error while reading data POLDR." << LogIO::POST ;
662    return -1 ;
663  }
664  for ( int i = 0 ; i < ARYNM ; i++ )
665    POLDR[i] = (double)ipoldr[i] ;
666  // DEBUG
667//   nro_debug_output( "POLDR", ARYNM, POLDR ) ;
668  //
669  if ( readColumn( POLAN, "POLAN", sameEndian ) != 0 ) {
670    os << LogIO::WARN << "Error while reading data POLAN." << LogIO::POST ;
671    return -1 ;
672  }
673  // DEBUG
674//   nro_debug_output( "POLAN", ARYNM, POLAN ) ;
675  //
676  if ( readColumn( DFRQ, "DFRQ", sameEndian ) != 0 ) {
677    os << LogIO::WARN << "Error while reading data DFRQ." << LogIO::POST ;
678    return -1 ;
679  }
680  // DEBUG
681//   nro_debug_output( "DFRQ", ARYNM, DFRQ ) ;
682  //
683  if ( readColumn( SIDBD, "SIDBD" ) != 0 ) {
684    os << LogIO::WARN << "Error while reading data SIDBD." << LogIO::POST ;
685    return -1 ;
686  }
687  // DEBUG
688//   nro_debug_output( "SIDBD", ARYNM, SIDBD ) ;
689  //
690  if ( readColumn( REFN, "REFN", sameEndian ) != 0 ) {
691    os << LogIO::WARN << "Error while reading data REFN." << LogIO::POST ;
692    return -1 ;
693  }
694  // DEBUG
695//   nro_debug_output( "REFN", ARYNM, REFN ) ;
696  //
697  if ( readColumn( IPINT, "IPINT", sameEndian ) != 0 ) {
698    os << LogIO::WARN << "Error while reading data IPINT." << LogIO::POST ;
699    return -1 ;
700  }
701  // DEBUG
702//   nro_debug_output( "IPINT", ARYNM, IPINT ) ;
703  //
704  if ( readColumn( MULTN, "MULTN", sameEndian ) != 0 ) {
705    os << LogIO::WARN << "Error while reading data MULTN." << LogIO::POST ;
706    return -1 ;
707  }
708  // DEBUG
709//   nro_debug_output( "MULTN", ARYNM, MULTN ) ;
710  //
711  if ( readColumn( MLTSCF, "MLTSCF", sameEndian ) != 0 ) {
712    os << LogIO::WARN << "Error while reading data MLTSCF." << LogIO::POST ;
713    return -1 ;
714  }
715  // DEBUG
716//   nro_debug_output( "MLTSCF", ARYNM, MLTSCF ) ;
717  //
718  if ( readColumn( LAGWIND, "LAGWIN" ) != 0 ) {
719    os << LogIO::WARN << "Error while reading data LAGWIND." << LogIO::POST ;
720    return -1 ;
721  }
722  // DEBUG
723//   nro_debug_output( "LAGWIND", ARYNM, LAGWIND ) ;
724  //
725  if ( readColumn( BEBW, "BEBW", sameEndian ) != 0 ) {
726    os << LogIO::WARN << "Error while reading data BEBW." << LogIO::POST ;
727    return -1 ;
728  }
729  // DEBUG
730//   nro_debug_output( "BEBW", ARYNM, BEBW ) ;
731  //
732  if ( readColumn( BERES, "BERES", sameEndian ) != 0 ) {
733    os << LogIO::WARN << "Error while reading data BERES." << LogIO::POST ;
734    return -1 ;
735  }
736  // DEBUG
737//   nro_debug_output( "BERES", ARYNM, BERES ) ;
738  //
739  if ( readColumn( CHWID, "CHWID", sameEndian ) != 0 ) {
740    os << LogIO::WARN << "Error while reading data CHWID." << LogIO::POST ;
741    return -1 ;
742  }
743  // DEBUG
744//   nro_debug_output( "CHWID", ARYNM, CHWID ) ;
745  //
746  if ( readARRY() != 0 ) {
747    os << LogIO::WARN << "Error while reading data ARRY." << LogIO::POST ;
748    return -1 ;
749  }
750  // DEBUG
751//   nro_debug_output( "ARRY", NRO_FITS_ARYMAX, ARRY ) ;
752  //
753  if ( readColumn( NFCAL, "NFCAL", sameEndian ) != 0 ) {
754    os << LogIO::WARN << "Error while reading data NFCAL." << LogIO::POST ;
755    return -1 ;
756  }
757  // DEBUG
758//   nro_debug_output( "NFCAL", ARYNM, NFCAL ) ;
759  //
760  if ( readColumn( F0CAL, "F0CAL", sameEndian ) != 0 ) {
761    os << LogIO::WARN << "Error while reading data F0CAL." << LogIO::POST ;
762    return -1 ;
763  }
764  // DEBUG
765//   nro_debug_output( "F0CAL", ARYNM, F0CAL ) ;
766  //
767  for ( int i= 0 ; i < 10 ; i++) {
768    vector<double> vv( ARYNM, 0 ) ;
769    if ( readColumn( vv, "FQCAL", sameEndian, i ) != 0 ) {
770      os << LogIO::WARN << "Error while reading data FQCAL." << LogIO::POST ;
771      return -1 ;
772    }
773    for ( int j = 0 ; j < ARYNM ; j++ ) {
774      FQCAL[j][i] = vv[j] ;
775    }
776  }
777  // DEBUG
778//   nro_debug_output( "FQCAL", ARYNM, 10, FQCAL ) ;
779  //
780  for ( int i= 0 ; i < 10 ; i++) {
781    vector<double> vv( ARYNM, 0 ) ;
782    if ( readColumn( vv, "CHCAL", sameEndian, i ) != 0 ) {
783      os << LogIO::WARN << "Error while reading data CHCAL." << LogIO::POST ;
784      return -1 ;
785    }
786    for ( int j = 0 ; j < ARYNM ; j++ ) {
787      CHCAL[j][i] = vv[j] ;
788    }
789  }
790  // DEBUG
791//   nro_debug_output( "CHCAL", ARYNM, 10, CHCAL ) ;
792  //
793  for ( int i= 0 ; i < 10 ; i++) {
794    vector<double> vv( ARYNM, 0 ) ;
795    if ( readColumn( vv, "CWCAL", sameEndian, i ) != 0 ) {
796      os << LogIO::WARN << "Error while reading data CWCAL." << LogIO::POST ;
797      return -1 ;
798    }
799    for ( int j = 0 ; j < ARYNM ; j++ ) {
800      CWCAL[j][i] = vv[j] ;
801    }
802  }
803  // DEBUG
804//   nro_debug_output( "CWCAL", ARYNM, 10, CWCAL ) ;
805  //
806  if ( readHeader( SCNLEN, "NAXIS1", sameEndian ) != 0 ) {
807    os << LogIO::WARN << "Error while reading data SCNLEN." << LogIO::POST ;
808    return -1 ;
809  }
810  // DEBUG
811  //cout << "SCNLEN = " << SCNLEN << endl ;
812  //
813  if ( readHeader( SBIND, "SBIND", sameEndian ) != 0 ) {
814    os << LogIO::NORMAL << "SBIND set to 0." << LogIO::POST ;
815    SBIND = 0 ;
816  }
817  // DEBUG
818  //cout << "SBIND = " << SBIND << endl ;
819  //
820  if ( readHeader( IBIT, "IBIT", sameEndian ) != 0 ) {
821    os << LogIO::NORMAL << "IBIT set to 8." << LogIO::POST ;
822    IBIT = 8 ;  // 8 bit? 12 bit?
823  }
824  // DEBUG
825  //cout << "IBIT = " << IBIT << endl ;
826  //
827  if ( readHeader( SITE, "TELESCOP" ) != 0 ) {
828    os << LogIO::WARN << "Error while reading data SITE." << LogIO::POST ;
829    return -1 ;
830  }
831  // DEBUG
832  //cout << "SITE = \'" << SITE << "\'" << endl ;
833  //
834  if ( readColumn( DSBFC, "DSBFC", sameEndian ) != 0 ) {
835    os << LogIO::NORMAL << "All DSBFC elements set to 1." << LogIO::POST ;
836    for ( int i = 0 ; i < ARYNM ; i++ )
837      DSBFC[i] = 1.0 ;
838  }
839  // DEBUG
840//   nro_debug_output( "DSBFC", ARYNM, DSBFC ) ;
841  //
842
843  return 0 ;
844}
845
846int NROFITSDataset::fillRecord( int i )
847{
848  LogIO os( LogOrigin( "NROFITSDataset", "fillRecord()", WHERE ) ) ;
849
850  int status = 0 ;
851  string str4( 4, ' ' ) ;
852  string str8( 8, ' ' ) ;
853  string str24( 24, ' ' ) ;
854
855  strcpy( record_->LSFIL, str4.c_str() ) ;
856  status = readTable( record_->LSFIL, "LSFIL", 4, i ) ;
857  if ( status ) {
858    os << LogIO::WARN << "Error while reading LSFIL." << LogIO::POST ;
859    return status ;
860  }
861  // DEBUG
862  //cout << "LSFIL(" << i << ") = " << record_->LSFIL << endl ;
863  //
864  status = readTable( record_->ISCAN, "ISCN", same_, i ) ;
865  if ( status ) {
866    os << LogIO::WARN << "Error while reading ISCAN." << LogIO::POST ;
867    return status ;
868  }
869  // DEBUG
870  //cout << "ISCAN(" << i << ") = " << record_->ISCAN << endl ;
871  //
872  vector<int> itmp( 6, 0 ) ;
873  status = readTable( itmp, "LAVST", same_, i ) ;
874  if ( status ) {
875    os << LogIO::WARN << "Error while reading LAVST." << LogIO::POST ;
876    return status ;
877  }
878  else {
879    sprintf( record_->LAVST, "%4d%02d%02d%02d%02d%02d.000", itmp[0], itmp[1], itmp[2], itmp[3], itmp[4], itmp[5] ) ;
880  }
881  // DEBUG
882  //cout << "LAVST(" << i << ") = " << record_->LAVST << endl ;
883  //
884  strcpy( record_->SCANTP, str8.c_str() ) ;
885    status = readTable( record_->SCANTP, "SCNTP", strlen(record_->SCANTP), i ) ;
886  if ( status ) {
887    os << LogIO::WARN << "Error while reading SCANTP." << LogIO::POST ;
888    return status ;
889  }
890  // DEBUG
891  //cout << "SCANTP(" << i << ") = " << record_->SCANTP << endl ;
892  //
893  char *name1 = "" ;
894  char *name2 = "" ;
895  if ( SCNCD == 0 ) {
896    name1 = "DRA" ;
897    name2 = "DDEC" ;
898  }
899  else if ( SCNCD == 1 ) {
900    name1 = "DGL" ;
901    name2 = "DGB" ;
902  }
903  else {
904    name1 = "DAZ" ;
905    name2 = "DEL" ;
906  }
907  status = readTable( record_->DSCX, name1, same_, i ) ;
908  if ( status ) {
909    os << LogIO::WARN << "Error while reading DSCX." << LogIO::POST ;
910    return status ;
911  }
912  // DEBUG
913  //cout << "DSCX(" << i << ") = " << record_->DSCX << endl ;
914  //
915  status = readTable( record_->DSCY, name2, same_, i ) ;
916  if ( status ) {
917    os << LogIO::WARN << "Error while reading DSCY." << LogIO::POST ;
918    return status ;
919  }
920  // DEBUG
921  //cout << "DSCY(" << i << ") = " << record_->DSCY << endl ;
922  //
923  if ( SCNCD == 0 ) {
924    name1 = "RA" ;
925    name2 = "DEC" ;
926  }
927  else if ( SCNCD == 1 ) {
928    name1 = "GL" ;
929    name2 = "GB" ;
930  }
931  else {
932    name1 = "AZ" ;
933    name2 = "EL" ;
934  }
935  status = readTable( record_->SCX, name1, same_, i ) ;
936  if ( status ) {
937    os << LogIO::WARN << "Error while reading SCX." << LogIO::POST ;
938    return status ;
939  }
940  // DEBUG
941  //cout << "SCX(" << i << ") = " << record_->SCX << endl ;
942  //
943  status = readTable( record_->SCY, name2, same_, i ) ;
944  if ( status ) {
945    os << LogIO::WARN << "Error while reading SCY." << LogIO::POST ;
946    return status ;
947  }
948  // DEBUG
949  //cout << "SCY(" << i << ") = " << record_->SCY << endl ;
950  //
951  status = readTable( record_->PAZ, "PAZ", same_, i ) ;
952  if ( status ) {
953    os << LogIO::WARN << "Error while reading PAZ." << LogIO::POST ;
954    return status ;
955  }
956  // DEBUG
957  //cout << "PAZ(" << i << ") = " << record_->PAZ << endl ;
958  //
959  status = readTable( record_->PEL, "PEL", same_, i ) ;
960  if ( status ) {
961    os << LogIO::WARN << "Error while reading PEL." << LogIO::POST ;
962    return status ;
963  }
964  // DEBUG
965  //cout << "PEL(" << i << ") = " << record_->PEL << endl ;
966  //
967  status = readTable( record_->RAZ, "RAZ", same_, i ) ;
968  if ( status ) {
969    os << LogIO::WARN << "Error while reading RAZ." << LogIO::POST ;
970    return status ;
971  }
972  // DEBUG
973  //cout << "RAZ(" << i << ") = " << record_->RAZ << endl ;
974  //
975  status = readTable( record_->REL, "REL", same_, i ) ;
976  if ( status ) {
977    os << LogIO::WARN << "Error while reading REL." << LogIO::POST ;
978    return status ;
979  }
980  // DEBUG
981  //cout << "REL(" << i << ") = " << record_->REL << endl ;
982  //
983  status = readTable( record_->XX, "XX", same_, i ) ;
984  if ( status ) {
985    os << LogIO::WARN << "Error while reading XX." << LogIO::POST ;
986    return status ;
987  }
988  // DEBUG
989  //cout << "XX(" << i << ") = " << record_->XX << endl ;
990  //
991  status = readTable( record_->YY, "YY", same_, i ) ;
992  if ( status ) {
993    os << LogIO::WARN << "Error while reading YY." << LogIO::POST ;
994    return status ;
995  }
996  // DEBUG
997  //cout << "YY(" << i << ") = " << record_->YY << endl ;
998  //
999  strcpy( record_->ARRYT, str4.c_str() ) ;
1000  status = readTable( record_->ARRYT, "ARRYT", strlen(record_->ARRYT), i ) ;
1001  for (int j = strlen(record_->ARRYT)-1 ; j >= 0 ; j--) {
1002    if (record_->ARRYT[j] == ' ') record_->ARRYT[j] = '\0';
1003    else break;
1004  }
1005  if ( status ) {
1006    os << LogIO::WARN << "Error while reading ARRYT." << LogIO::POST ;
1007    return status ;
1008  }
1009  // DEBUG
1010  //cout << "ARRYT(" << i << ") = " << record_->ARRYT << endl ;
1011  //
1012  double dtmp ;
1013  status = readTable( dtmp, "TEMP", same_, i ) ;
1014  if ( status ) {
1015    os << LogIO::WARN << "Error while reading TEMP." << LogIO::POST ;
1016    return status ;
1017  }
1018  else {
1019    record_->TEMP = dtmp ;
1020  }
1021  // DEBUG
1022  //cout << "TEMP(" << i << ") = " << record_->TEMP << endl ;
1023  //
1024  status = readTable( dtmp, "PATM", same_, i ) ;
1025  if ( status ) {
1026    os << LogIO::WARN << "Error while reading PATM." << LogIO::POST ;
1027    return status ;
1028  }
1029  else {
1030    record_->PATM = dtmp ;
1031  }
1032  // DEBUG
1033  //cout << "PATM(" << i << ") = " << record_->PATM << endl ;
1034  //
1035  status = readTable( dtmp, "PH2O", same_, i ) ;
1036  if ( status ) {
1037    os << LogIO::WARN << "Error while reading PH2O." << LogIO::POST ;
1038    return status ;
1039  }
1040  else {
1041    record_->PH2O = dtmp ;
1042  }
1043  // DEBUG
1044  //cout << "PH2O(" << i << ") = " << record_->PH2O << endl ;
1045  //
1046  status = readTable( dtmp, "VWIND", same_, i ) ;
1047  if ( status ) {
1048    os << LogIO::WARN << "Error while reading VWIND." << LogIO::POST ;
1049    return status ;
1050  }
1051  else {
1052    record_->VWIND = dtmp ;
1053  }
1054  // DEBUG
1055  //cout << "VWIND(" << i << ") = " << record_->VWIND << endl ;
1056  //
1057  status = readTable( dtmp, "DWIND", same_, i ) ;
1058  if ( status ) {
1059    os << LogIO::WARN << "Error while reading DWIND." << LogIO::POST ;
1060    return status ;
1061  }
1062  else {
1063    record_->DWIND = dtmp ;
1064  }
1065  // DEBUG
1066  //cout << "DWIND(" << i << ") = " << record_->DWIND << endl ;
1067  //
1068  status = readTable( dtmp, "TAU", same_, i ) ;
1069  if ( status ) {
1070    os << LogIO::WARN << "Error while reading TAU." << LogIO::POST ;
1071    return status ;
1072  }
1073  else {
1074    record_->TAU = dtmp ;
1075  }
1076  // DEBUG
1077  //cout << "TAU(" << i << ") = " << record_->TAU << endl ;
1078  //
1079  status = readTable( dtmp, "TSYS", same_, i ) ;
1080  if ( status ) {
1081    os << LogIO::WARN << "Error while reading TSYS." << LogIO::POST ;
1082    return status ;
1083  }
1084  else {
1085    record_->TSYS = dtmp ;
1086  }
1087  // DEBUG
1088  //cout << "TSYS(" << i << ") = " << record_->TSYS << endl ;
1089  //
1090  status = readTable( dtmp, "BATM", same_, i ) ;
1091  if ( status ) {
1092    os << LogIO::WARN << "Error while reading BATM." << LogIO::POST ;
1093    return status ;
1094  }
1095  else {
1096    record_->BATM = dtmp ;
1097  }
1098  // DEBUG
1099  //cout << "BATM(" << i << ") = " << record_->BATM << endl ;
1100  //
1101  status = readTable( record_->VRAD, "VRAD", same_, i ) ;
1102  if ( status ) {
1103    os << LogIO::WARN << "Error while reading TEMP." << LogIO::POST ;
1104    return status ;
1105  }
1106  // DEBUG
1107  //cout << "VRAD(" << i << ") = " << record_->VRAD << endl ;
1108  //
1109  status = readTable( record_->FREQ0, "FRQ0", same_, i ) ;
1110  if ( status ) {
1111    os << LogIO::WARN << "Error while reading FREQ0." << LogIO::POST ;
1112    return status ;
1113  }
1114  // DEBUG
1115  //cout << "FREQ0(" << i << ") = " << record_->FREQ0 << endl ;
1116  //
1117  status = readTable( record_->FQTRK, "FQTRK", same_, i ) ;
1118  if ( status ) {
1119    os << LogIO::WARN << "Error while reading FQTRK." << LogIO::POST ;
1120    return status ;
1121  }
1122  // DEBUG
1123  //cout << "FQTRK(" << i << ") = " << record_->FQTRK << endl ;
1124  //
1125  status = readTable( record_->FQIF1, "FQIF1", same_, i ) ;
1126  if ( status ) {
1127    os << LogIO::WARN << "Error while reading FQIF1." << LogIO::POST ;
1128    return status ;
1129  }
1130  // DEBUG
1131  //cout << "FQIF1(" << i << ") = " << record_->FQIF1 << endl ;
1132  //
1133  status = readTable( record_->ALCV, "ALCV", same_, i ) ;
1134  if ( status ) {
1135    os << LogIO::WARN << "Error while reading ALCV." << LogIO::POST ;
1136    return status ;
1137  }
1138  // DEBUG
1139  //cout << "ALCV(" << i << ") = " << record_->ALCV << endl ;
1140  //
1141  record_->IDMY0 = 0 ;
1142  status = readTable( record_->DPFRQ, "DPFRQ", same_, i ) ;
1143  if ( status ) {
1144    //os << LogIO::WARN << "Error  DPFRQ set to 0." << LogIO::POST ;
1145    record_->DPFRQ = 0.0 ;
1146  }
1147  // DEBUG
1148  //cout << "DPFRQ(" << i << ") = " << record_->DPFRQ << endl ;
1149  //
1150  status = readTable( record_->SFCTR, "SFCTR", same_, i ) ;
1151  if ( status ) {
1152    os << LogIO::WARN << "Error while reading SFCTR." << LogIO::POST ;
1153    return status ;
1154  }
1155  // DEBUG
1156  //cout << "SFCTR(" << i << ") = " << record_->SFCTR << endl ;
1157  //
1158  status = readTable( record_->ADOFF, "ADOFF", same_, i ) ;
1159  if ( status ) {
1160    os << LogIO::WARN << "Error while reading ADOFF." << LogIO::POST ;
1161    return status ;
1162  }
1163  // DEBUG
1164  //cout << "ADOFF(" << i << ") = " << record_->ADOFF << endl ;
1165  //
1166  //status = readTable( record_->JDATA, "LDATA", same_, i ) ;
1167  status = readTable( JDATA, "LDATA", same_, i ) ;
1168  if ( status ) {
1169    os << LogIO::WARN << "Error while reading JDATA." << LogIO::POST ;
1170    return status ;
1171  }
1172  // DEBUG
1173//   for ( int i = 0 ; i < chmax_ ; i++ )
1174//     //cout << "JDATA[" << i << "] = " << JDATA[i] << " " ;
1175//   //cout << endl ;
1176  //
1177
1178
1179  // Update IPTIM since it depends on the row for NROFITS
1180  int integ ;
1181  status = readTable( integ, "INTEG", same_, i ) ;
1182  if ( !status ) {
1183    IPTIM = (double)integ ;
1184  }
1185
1186  return status ;
1187}
1188
1189vector< vector<double> > NROFITSDataset::getSpectrum()
1190{
1191  vector< vector<double> > spec;
1192
1193  for ( int i = 0 ; i < rowNum_ ; i++ ) {
1194    spec.push_back( getSpectrum( i ) ) ;
1195  }
1196
1197  return spec ;
1198}
1199
1200vector<double> NROFITSDataset::getSpectrum( int i )
1201{
1202  vector<double> spec( chmax_, 0.0 ) ;
1203  vector<double> specout( chmax_, 0.0 ) ;
1204  const NRODataRecord *record = getRecord( i ) ;
1205  double scale = record->SFCTR ;
1206  double offset = record->ADOFF ;
1207  double dscale = MLTSCF[getIndex( i )] ;
1208  //vector<int> ispec = record->JDATA ;
1209  vector<int> ispec = JDATA ;
1210  for ( int ii = 0 ; ii < chmax_ ; ii++ ) {
1211    spec[ii] = (double)( ispec[ii] * scale + offset ) * dscale ;
1212  }
1213
1214  // for AOS, re-gridding is needed
1215  if ( strncmp( record->ARRYT, "H", 1 ) == 0
1216       || strncmp( record->ARRYT, "W", 1 ) == 0
1217       || strncmp( record->ARRYT, "U", 1 ) == 0 ) {
1218
1219    string arryt = string( record->ARRYT ) ;
1220    uInt ib = getArrayId( arryt ) ;
1221    vector<double> fqcal = getFQCAL()[ib] ;
1222    vector<double> chcal = getCHCAL()[ib] ;
1223    int ncal = getNFCAL()[ib] ;
1224
1225//     //cout << "NRODataset::getFrequencies()  ncal = " << ncal << endl ;
1226    while ( ncal < (int)fqcal.size() ) {
1227      fqcal.pop_back() ;
1228      chcal.pop_back() ;
1229    }
1230    Vector<Double> xin( chcal ) ;
1231    Vector<Double> yin( fqcal ) ;
1232    int nchan = getNUMCH() ;
1233    Vector<Double> xout( nchan ) ;
1234    indgen( xout ) ;
1235    Vector<Double> yout ;
1236    InterpolateArray1D<Double, Double>::interpolate( yout, xout, xin, yin, InterpolateArray1D<Double,Double>::cubic ) ;
1237    // debug
1238    //cout << "i=" << i << endl ;
1239//     if ( i == 16 ) {
1240//       ofstream ofs0( "spgrid0.dat" ) ;
1241//       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
1242//         ofs0 << xout[ii] << "," ;
1243//       ofs0 << endl ;
1244//       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
1245//         ofs0 << setprecision(16) << record->FREQ0+yout[ii] << "," ;
1246//       ofs0 << endl ;
1247//       ofs0.close() ;
1248//     }
1249    //
1250    Vector<Double> z( nchan ) ;
1251    Double bw = abs( yout[nchan-1] - yout[0] ) ;
1252    bw += 0.5 * abs( yout[nchan-1] - yout[nchan-2] + yout[1] - yout[0] ) ;
1253    Double dz = bw / (Double)nchan ;
1254    if ( yout[0] > yout[nchan-1] )
1255      dz = - dz ;
1256    z[0] = yout[0] - 0.5 * ( yout[1] - yout[0] - dz ) ;
1257    for ( int ii = 1 ; ii < nchan ; ii++ )
1258      z[ii] = z[ii-1] + dz ;
1259    Vector<Double> zi( nchan+1 ) ;
1260    Vector<Double> yi( nchan+1 ) ;
1261    zi[0] = z[0] - 0.5 * dz ;
1262    zi[1] = z[0] + 0.5 * dz ;
1263    yi[0] = yout[0] - 0.5 * ( yout[1] - yout[0] ) ;
1264    yi[1] = yout[0] + 0.5 * ( yout[1] - yout[0] ) ;
1265    for ( int ii = 2 ; ii < nchan ; ii++ ) {
1266      zi[ii] = zi[ii-1] + dz ;
1267      yi[ii] = yi[ii-1] + 0.5 * ( yout[ii] - yout[ii-2] ) ;
1268    }
1269    zi[nchan] = z[nchan-1] + 0.5 * dz ;
1270    yi[nchan] = yout[nchan-1] + 0.5 * ( yout[nchan-1] - yout[nchan-2] ) ;
1271//     // debug
1272//     //cout << "nchan=" << nchan << ", bw=" << bw << ", dz=" << dz
1273//          << ", y[1]-y[0]=" << yout[1]-yout[0] << endl ;
1274//     //cout << "z: " << z[0] << " - " << z[nchan-1]
1275//          << ", zi: " << zi[0] << " - " << zi[nchan] << endl ;
1276//     //cout << "y: " << yout[0] << " - " << yout[nchan-1]
1277//          << ", yi: " << yi[0] << " - " << yi[nchan] << endl ;
1278//     ofstream ofs1( "spgrid1.dat", ios::out | ios::app ) ;
1279//     ofs1 << "spid=" << i << ", ARRYT=" << record->ARRYT << endl ;
1280//     ofs1 << "z[0]=" << z[0] << ", yout[0]=" << yout[0] << endl ;
1281//     for ( int ii = 1; ii < nchan ; ii++ ) {
1282//       ofs1 << "               dz=" << z[ii]-z[ii-1] << ", dy=" << yout[ii]-yout[ii-1] << endl ;
1283//       ofs1 << "z[" << ii << "]=" << z[ii] << ", yout[" << ii << "]=" << yout[ii] << endl ;
1284//     }
1285//     ofs1.close() ;
1286//     ofstream ofs2( "spgrid2.dat", ios::out | ios::app ) ;
1287//     ofs2 << "spid=" << i << ", ARRYT=" << record->ARRYT << endl ;
1288//     for ( int ii = 0 ; ii < nchan+1 ; ii++ )
1289//       ofs2 << "zi[" << ii << "]=" << zi[ii] << ", yi[" << ii << "]=" << yi[ii] << endl ;
1290//     ofs2.close() ;
1291//     //
1292    int ichan = 0 ;
1293    double wsum = 0.0 ;
1294    // debug
1295    //ofstream ofs3( "spgrid3.dat", ios::out | ios::app ) ;
1296    if ( dz > 0.0 ) {
1297      for ( int ii = 0 ; ii < nchan ; ii++ ) {
1298        double zl = zi[ii] ;
1299        double zr = zi[ii+1] ;
1300        for ( int j = ichan ; j < nchan ; j++ ) {
1301          double yl = yi[j] ;
1302          double yr = yi[j+1] ;
1303          if ( yl <= zl ) {
1304            if ( yr <= zl ) {
1305              continue ;
1306            }
1307            else if ( yr <= zr ) {
1308              specout[ii] += spec[j] * ( yr - zl ) ;
1309              wsum += ( yr - zl ) ;
1310            }
1311            else {
1312              specout[ii] += spec[j] * dz ;
1313              wsum += dz ;
1314              ichan = j ;
1315              break ;
1316            }
1317          }
1318          else if ( yl < zr ) {
1319            if ( yr <= zr ) {
1320              specout[ii] += spec[j] * ( yr - yl ) ;
1321              wsum += ( yr - yl ) ;
1322            }
1323            else {
1324              specout[ii] += spec[j] * ( zr - yl ) ;
1325              wsum += ( zr - yl ) ;
1326              ichan = j ;
1327              break ;
1328            }
1329          }
1330          else {
1331            ichan = j - 1 ;
1332            break ;
1333          }
1334        }
1335        specout[ii] /= wsum ;
1336        wsum = 0.0 ;
1337      }
1338    }
1339    else if ( dz < 0.0 ) {
1340      for ( int ii = 0 ; ii < nchan ; ii++ ) {
1341        double zl = zi[ii] ;
1342        double zr = zi[ii+1] ;
1343        for ( int j = ichan ; j < nchan ; j++ ) {
1344          double yl = yi[j] ;
1345          double yr = yi[j+1] ;
1346          if ( yl >= zl ) {
1347            if ( yr >= zl ) {
1348              continue ;
1349            }
1350            else if ( yr >= zr ) {
1351              specout[ii] += spec[j] * abs( yr - zl ) ;
1352              wsum += abs( yr - zl ) ;
1353            }
1354            else {
1355              specout[ii] += spec[j] * abs( dz ) ;
1356              wsum += abs( dz ) ;
1357              ichan = j ;
1358              break ;
1359            }
1360          }
1361          else if ( yl > zr ) {
1362            if ( yr >= zr ) {
1363              specout[ii] += spec[j] * abs( yr - yl ) ;
1364              wsum += abs( yr - yl ) ;
1365            }
1366            else {
1367              specout[ii] += spec[j] * abs( zr - yl ) ;
1368              wsum += abs( zr - yl ) ;
1369              ichan = j ;
1370              break ;
1371            }
1372          }
1373          else {
1374            ichan = j - 1 ;
1375            break ;
1376          }
1377        }
1378        specout[ii] /= wsum ;
1379        wsum = 0.0 ;
1380      }
1381    }
1382    //specout = spec ;
1383    //ofs3.close() ;
1384  }
1385  else {
1386    specout = spec ;
1387  }
1388
1389  return specout ;
1390}
1391
1392int NROFITSDataset::getIndex( int irow )
1393{
1394  const NRODataRecord *record = getRecord( irow ) ;
1395  string str = record->ARRYT ;
1396  string::size_type pos = str.find( " " ) ;
1397  if ( pos != string::npos )
1398    str = str.substr( 0, pos ) ;
1399  int index = -1 ;
1400  for ( int i = 0 ; i < ARYNM ; i++ ) {
1401    if ( str.compare( 0, 3, ARYTP[i] ) == 0 ) {
1402      index = i ;
1403      break ;
1404    }
1405  }
1406  return index ;
1407}
1408
1409double NROFITSDataset::radRA( string ra )
1410{
1411  int pos1 = ra.find( ':' ) ;
1412  int pos2 ;
1413  string ch = ra.substr( 0, pos1 ) ;
1414  //cout << "ch = \'" << ch << "\'" << endl ;
1415  pos2 = pos1 + 1 ;
1416  pos1 = ra.find( ':', pos2 ) ;
1417  string cm = ra.substr( pos2, pos1 - pos2 ) ;
1418  //cout << "cm = \'" << cm << "\'" << endl ;
1419  pos2 = pos1 + 1 ;
1420  pos1 = ra.size() ;
1421  string cs = ra.substr( pos2, pos1 - pos2 ) ;
1422  //cout << "cs = \'" << cs << "\'" << endl ;
1423  double h ;
1424  if ( ra[0] != '-' )
1425    h = atof( ch.c_str() ) + atof( cm.c_str() ) / 60.0 + atof( cs.c_str() ) / 3600.0 ;
1426  else
1427    h = atof( ch.c_str() ) - atof( cm.c_str() ) / 60.0 - atof( cs.c_str() ) / 3600.0 ;
1428  double rra = h * M_PI / 12.0 ;
1429  return rra ;
1430}
1431
1432double NROFITSDataset::radDEC( string dec )
1433{
1434  int pos1 = dec.find( ':' ) ;
1435  int pos2 ;
1436  string cd = dec.substr( 0, pos1 ) ;
1437  //cout << "cd = \'" << cd << "\'" << endl ;
1438  pos2 = pos1 + 1 ;
1439  pos1 = dec.find( ':', pos2 ) ;
1440  string cm = dec.substr( pos2, pos1 - pos2 ) ;
1441  //cout << "cm = \'" << cm << "\'" << endl ;
1442  pos2 = pos1 + 1 ;
1443  pos1 = dec.size() ;
1444  string cs = dec.substr( pos2, pos1 - pos2 ) ;
1445  //cout << "cs = \'" << cs << "\'" << endl ;
1446  double h ;
1447  if ( dec[0] != '-' )
1448    h = atof( cd.c_str() ) + atof( cm.c_str() ) / 60.0 + atof( cs.c_str() ) / 3600.0 ;
1449  else
1450    h = atof( cd.c_str() ) - atof( cm.c_str() ) / 60.0 - atof( cs.c_str() ) / 3600.0 ;
1451  double rdec = h * M_PI / 180.0 ;
1452  return rdec ;
1453}
1454
1455void NROFITSDataset::getField()
1456{
1457  long offset = 0;
1458  for ( int i = 0 ; i < numField_ ; i++ ) {
1459    char key1[9] ;
1460    char key2[9] ;
1461    char key3[9] ;
1462    if ( i < 9 ) {
1463      sprintf( key1, "TFORM%d  ", i+1 ) ;
1464      sprintf( key2, "TTYPE%d  ", i+1 ) ;
1465      sprintf( key3, "TUNIT%d  ", i+1 ) ;
1466      //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1467    }
1468    else if ( i < 99 ) {
1469      sprintf( key1, "TFORM%2d ", i+1 ) ;
1470      sprintf( key2, "TTYPE%2d ", i+1 ) ;
1471      sprintf( key3, "TUNIT%2d ", i+1 ) ;
1472      //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1473    }
1474    else {
1475      sprintf( key1, "TFORM%3d", i+1 ) ;
1476      sprintf( key2, "TTYPE%3d", i+1 ) ;
1477      sprintf( key3, "TUNIT%3d", i+1 ) ;
1478      //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1479    }
1480    //char tmp[9] ;
1481    string tmp ;
1482    //strcpy( tmp, "         " ) ;
1483    if ( readHeader( tmp, key1 ) != 0 ) {
1484      cerr << "Error while reading field keyword for scan header." << endl ;
1485      return ;
1486    }
1487    string form = tmp ;
1488    string::size_type spos = form.find( " " ) ;
1489    if ( spos != string::npos )
1490      form = form.substr( 0, spos ) ;
1491    //strcpy( tmp, "         " ) ;
1492    if ( readHeader( tmp, key2 ) != 0 ) {
1493      cerr << "Error while reading field type for scan header." << endl ;
1494      return ;
1495    }
1496    //names_[i] = string( tmp ) ;
1497    names_[i] = tmp ;
1498    spos = names_[i].find( " " ) ;
1499    if ( spos != string::npos )
1500      names_[i] = names_[i].substr( 0, spos ) ;
1501    //strcpy( tmp, "         " ) ;
1502    if ( form.find( "A" ) != string::npos ) {
1503      //cout << "skip to get unit: name = " << form << endl ;
1504      //strcpy( tmp, "none    " ) ;
1505      tmp = "none" ;
1506    }
1507    else {
1508      //cout << "get unit: name = " << form << endl ;
1509      if ( readHeader( tmp, key3 ) != 0 ) {
1510        //strcpy( tmp, "none    " ) ;
1511        tmp = "none" ;
1512      }
1513    }
1514    //units_[i] = string( tmp ) ;
1515    units_[i] = tmp ;
1516    spos = units_[i].find( " " ) ;
1517    if ( spos != string::npos )
1518      units_[i] = units_[i].substr( 0, spos ) ;
1519    //cout << "i = " << i << ": name=" << form << " type=" << names_[i] << " unit=" << units_[i] << endl ;
1520
1521    offsets_[i] = offset ;
1522    string substr1 = form.substr( 0, form.size()-1 ) ;
1523    string substr2 = form.substr( form.size()-1, 1 ) ;
1524    //cout << "substr1 = " << substr1 << ", substr2 = " << substr2 << endl ;
1525    int o1 = atoi( substr1.c_str() ) ;
1526    int o2 = 0 ;
1527    if ( substr2 == "A" )
1528      o2 = sizeof(char) ;
1529    else if ( substr2 == "J" )
1530      o2 = sizeof(int) ;
1531    else if ( substr2 == "F" )
1532      o2 = sizeof(float) ;
1533    else if ( substr2 == "D" )
1534      o2 = sizeof(double) ;
1535    sizes_[i] = o1 * o2 ;
1536    offset += sizes_[i] ;
1537  } 
1538}
1539
1540void NROFITSDataset::fillARYTP()
1541{
1542  string arry ;
1543  int count = 0 ;
1544  string arry1 ;
1545  string arry2 ;
1546  string arry3 ;
1547  string arry4 ;
1548  char arytp[4];
1549  if ( readHeader( arry, "ARRY1" ) == 0 )
1550    arry1 = arry ;
1551  else
1552    arry1 = "00000000000000000000" ;
1553  for ( int i = 0 ; i < 20 ; i++ ) {
1554    if ( arry1[i] == '1' ) {
1555      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1556      sprintf( arytp, "H%d", i+1 ) ;
1557      ARYTP[count++] = string( arytp ) ;
1558    }
1559  }
1560  if ( readHeader( arry, "ARRY2" ) == 0 )
1561    arry2 = arry ;
1562  else
1563    arry2 = "00000000000000000000" ;
1564  for ( int i = 0 ; i < 10 ; i++ ) {
1565    if ( arry2[i] == '1' ) {
1566      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1567      sprintf( arytp, "W%d", i+1 ) ;
1568      ARYTP[count++] = string( arytp ) ;
1569    }
1570  }
1571  for ( int i = 10 ; i < 15 ; i++ ) {
1572    if ( arry2[i] == '1' ) {
1573      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1574      sprintf( arytp, "U%d", i-9 ) ;
1575      ARYTP[count++] = string( arytp ) ;
1576    }
1577  }
1578  for ( int i = 15 ; i < 20 ; i++ ) {
1579    if ( arry2[i] == '1' ) {
1580      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1581      sprintf( arytp, "X%d", i-14 ) ;
1582      ARYTP[count++] = string( arytp ) ;
1583    }
1584  }
1585  if ( readHeader( arry, "ARRY3" ) == 0 )
1586    arry3 = arry ;
1587  else
1588    arry3 = "00000000000000000000" ;
1589  for ( int i = 0 ; i < 20 ; i++ ) {
1590    if ( arry3[i] == '1' ) {
1591      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1592      sprintf( arytp, "A%d", i+1 ) ;
1593      ARYTP[count++] = string( arytp ) ;
1594    }
1595  }
1596  if ( readHeader( arry, "ARRY4" ) == 0 )
1597    arry4 = arry ;
1598  else
1599    arry4 = "00000000000000000000" ;
1600  for ( int i = 0 ; i < 20 ; i++ ) {
1601    if ( arry4[i] == '1' ) {
1602      for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
1603      sprintf( arytp, "A%d", i+21 ) ;
1604      ARYTP[count++] = string( arytp ) ;
1605    }
1606  }     
1607  //nro_debug_output("ARYTP", ARYTP.size(), ARYTP);
1608}
1609
1610int NROFITSDataset::readARRY()
1611{
1612  LogIO os( LogOrigin( "NROFITSDataset", "readARRY()", WHERE ) ) ;
1613
1614  string arry1 ;
1615  string arry2 ;
1616  string arry3 ;
1617  string arry4 ;
1618  int status = readHeader( arry1, "ARRY1" ) ;
1619  if ( status ) {
1620    os << LogIO::SEVERE << "Error while reading ARRY1" << LogIO::POST ;
1621    return status ;
1622  }
1623  status = readHeader( arry2, "ARRY2" ) ;
1624  if ( status ) {
1625    os << LogIO::SEVERE << "Error while reading ARRY2" << LogIO::POST ;
1626    return status ;
1627  }
1628  status = readHeader( arry3, "ARRY3" ) ;
1629  if ( status ) {
1630    os << LogIO::SEVERE << "Error while reading ARRY3" << LogIO::POST ;
1631    return status ;
1632  }
1633  status = readHeader( arry4, "ARRY4" ) ;
1634  if ( status ) {
1635    os << LogIO::SEVERE << "Error while reading ARRY4" << LogIO::POST ;
1636    return status ;
1637  }
1638  int index = 0 ;
1639  for ( int i = 0 ; i < 20 ; i++ ) {
1640    // AOSH
1641    if ( arry1[i] == '1' )
1642      ARRY[index] = 1 ;
1643    else
1644      ARRY[index] = 0 ;
1645    // AOSW, AOSU, FX
1646    if ( arry2[i] == '1' )
1647      ARRY[index+20] = 1 ;
1648    else
1649      ARRY[index+20] = 0 ;
1650    // AC45 (1-20)
1651    if ( arry3[i] == '1' )
1652      ARRY[index+40] = 1 ;
1653    else
1654      ARRY[index+40] = 0 ;
1655    // AC45 (21-35)
1656    if ( i < 15 ) {
1657      if ( arry4[i] == '1' )
1658        ARRY[index+60] = 1 ;
1659      else
1660        ARRY[index+60] = 0 ;
1661    }
1662    index++ ;
1663  }
1664  return status ;
1665}
1666
1667void NROFITSDataset::findData()
1668{
1669  LogIO os( LogOrigin( "NROFITSDataset", "findData()", WHERE ) ) ;
1670
1671  // skip header
1672  fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
1673
1674  // get offset
1675  long offset = getOffset( "ARRYT" ) ;
1676  if ( offset == -1 ) {
1677    //cerr << "Error, ARRYT is not found in the name list." << endl ;
1678    return ;
1679  }
1680  //cout << "offset for ARRYT is " << offset << " bytes." << endl ;
1681  fseek( fp_, offset, SEEK_CUR ) ;
1682  int count = 0 ;
1683  int index = 0 ;
1684  while ( count < ARYNM && index < rowNum_ ) {
1685    char ctmp[5] ;
1686    fread( ctmp, 1, 4, fp_ ) ;
1687    ctmp[4] = '\0' ;
1688    //cout << "ctmp = " << ctmp << endl ;
1689    for ( int i = 0 ; i < ARYNM ; i++ ) {
1690      if ( arrayid_[i] != -1 )
1691        continue ;
1692      else if ( strncmp( ctmp, ARYTP[i].c_str(), ARYTP[i].size() ) == 0 ) {
1693        //cout << "matched: i = " << i << ", ARYTP = " << ARYTP[i] << ", ctmp = " << ctmp << endl ;
1694        arrayid_[i] = index ;
1695        count++ ;
1696      }
1697    }
1698    fseek( fp_, scanLen_-4, SEEK_CUR ) ;
1699    index++ ;
1700  }
1701
1702  if ( count != ARYNM ) {
1703    os << LogIO::WARN << "NROFITSDataset::findData()  failed to find rows for " ;
1704    for ( int i = 0 ; i < ARYNM ; i++ ) {
1705      if ( arrayid_[i] == -1 ) {
1706        os << LogIO::WARN << ARYTP[i] << " " ;
1707      }
1708    }
1709    os.post() ;
1710  }
1711
1712//   for ( int i = 0 ; i < ARYNM ; i++ )
1713//     //cout << "arrayid_[" << i << "] = " << arrayid_[i] << endl ;
1714  }
1715
1716long NROFITSDataset::getOffset( char *name )
1717{
1718  long offset = 0 ;
1719  string sname( name ) ;
1720  long j = -1 ;
1721  for ( long i = 0 ; i < numField_ ; i++ ) {
1722    // escape if name is found
1723    //cout << "names_[" << i << "] = " << names_[i] << "  sname = " << sname << endl ;
1724    if ( names_[i] == sname ) {
1725      j = i ;
1726      break ;
1727    }
1728  }
1729
1730  offset = (j >= 0) ? offsets_[j] : j ;
1731
1732  return offset ;
1733}
1734
1735int NROFITSDataset::getPolarizationNum()
1736{
1737  int npol = 0 ;
1738
1739  vector<char> type( 2 ) ;
1740  type[0] = 'C' ;
1741  type[1] = 'L' ;
1742  vector<double> crot ;
1743  vector<double> lagl ;
1744
1745  for ( int i = 0 ; i < ARYNM ; i++ ) {
1746    if ( POLTP[i][0] == type[0] ) {
1747      // circular polarization
1748      if( count( crot.begin(), crot.end(), POLDR[i] ) != 0 ) {
1749        crot.push_back( POLDR[i] ) ;
1750        npol++ ;
1751      }
1752    }
1753    else if ( POLTP[i][0] == type[1] ) {
1754      // linear polarization
1755      if ( count( lagl.begin(), lagl.end(), POLAN[i] ) != 0 ) {
1756        lagl.push_back( POLAN[i] ) ;
1757        npol++ ;
1758      }
1759    }
1760  }
1761
1762  if ( npol == 0 )
1763    npol = 1 ;
1764
1765
1766  return npol ;
1767}
1768
1769int NROFITSDataset::readHeader( string &v, char *name )
1770{
1771  //
1772  // Read 'name' attribute defined as char from the FITS Header
1773  // 
1774  int status = 0 ;
1775 
1776  char buf[81] ;
1777  strcpy( buf, "     " ) ;
1778  fseek( fp_, 0, SEEK_SET ) ;
1779  int count = 0 ;
1780  while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
1781    fread( buf, 1, 80, fp_ ) ;
1782    buf[80] = '\0' ;
1783    count++ ;
1784  }
1785  if ( strncmp( buf, "END", 3 ) == 0 ) {
1786    //cerr << "NROFITSDataset::readHeader()  keyword " << name << " not found." << endl ;
1787    //cerr << "count = " << count << endl ;
1788    status = -1 ;
1789    return status ;
1790  }
1791  string str( buf ) ;
1792  int pos1 = str.find( '\'' ) + 1 ;
1793  int pos2 = str.find( '\'', pos1 ) ;
1794  unsigned int clen = pos2 - pos1 ;
1795  //cout << "string: " << str << endl ;
1796  //cout << "value:  " << str.substr( pos1, clen ).c_str() << endl ;
1797  //cout << "clen = " << clen << endl ;
1798  v = str.substr( pos1, clen ) ;
1799  //cout << "v = \'" << v << "\'" << endl ;
1800
1801  return status ;
1802}
1803
1804int NROFITSDataset::readHeader( int &v, char *name, int /*b*/ )
1805{
1806  //
1807  // Read 'name' attribute defined as int from the FITS Header
1808  //
1809  int status = 0 ;
1810
1811  char buf[81] ;
1812  strcpy( buf, "     " ) ;
1813  fseek( fp_, 0, SEEK_SET ) ;
1814  while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
1815    fread( buf, 1, 80, fp_ ) ;
1816    buf[80] = '\0' ;
1817    //char bufo[9] ;
1818    //strncpy( bufo, buf, 8 ) ;
1819    //bufo[8] = '\0' ;
1820    //cout << "header: " << bufo << endl ;
1821  }
1822  if ( strncmp( buf, "END", 3 ) == 0 ) {
1823    //cerr << "NROFITSDataset::readHeader()  keyword " << name << " not found." << endl ;
1824    status = -1 ;
1825    return status ;
1826  }
1827  string str( buf ) ;
1828  int pos1 = str.find( '=' ) + 1 ;
1829  int pos2 = str.find( '/' ) ;
1830  //cout << "string: " << str << endl ;
1831  //cout << "value:  " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
1832  v = atoi( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
1833  //cout << "v = " << v << endl ;
1834
1835  //cout << "NROFITSDataset::readHeader() end to read" << endl ;
1836  return status ;
1837}
1838
1839
1840int NROFITSDataset::readHeader( float &v, char *name, int /*b*/ )
1841{
1842  //
1843  // Read 'name' attribute defined as float from the FITS Header
1844  //
1845  int status = 0 ;
1846
1847  char buf[81] ;
1848  strcpy( buf, "     " ) ;
1849  fseek( fp_, 0, SEEK_SET ) ;
1850  while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
1851    fread( buf, 1, 80, fp_ ) ;
1852    buf[80] = '\0' ;
1853    //char bufo[9] ;
1854    //strncpy( bufo, buf, 8 ) ;
1855    //bufo[8] = '\0' ;
1856    //cout << "header: " << bufo << endl ;
1857  }
1858  if ( strncmp( buf, "END", 3 ) == 0 ) {
1859    //cerr << "NROFITSDataset::readHeader()  keyword " << name << " not found." << endl ;
1860    status = -1 ;
1861    return status ;
1862  }
1863  string str( buf ) ;
1864  int pos1 = str.find( '=' ) + 1 ;
1865  int pos2 = str.find( '/' ) ;
1866  //cout << "string: " << str << endl ;
1867  //cout << "value:  " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
1868  v = atof( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
1869  //cout << "v = " << v << endl ;
1870
1871  return status ;
1872}
1873
1874int NROFITSDataset::readHeader( double &v, char *name, int /*b*/ )
1875{
1876  //
1877  // Read 'name' attribute defined as double from the FITS Header
1878  //
1879  int status = 0 ;
1880
1881  char buf[81] ;
1882  strcpy( buf, "     " ) ;
1883  fseek( fp_, 0, SEEK_SET ) ;
1884  while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
1885    fread( buf, 1, 80, fp_ ) ;
1886    buf[80] = '\0' ;
1887    char bufo[9] ;
1888    strncpy( bufo, buf, 8 ) ;
1889    bufo[8] = '\0' ;
1890    //cout << "header: \'" << bufo << "\' bufo = \'" << bufo << "\' ";
1891    //cout << strncmp( buf, name, strlen(name) ) << endl ;
1892  }
1893  if ( strncmp( buf, "END", 3 ) == 0 ) {
1894    //cerr << "NROFITSDataset::readHeader()  keyword " << name << " not found." << endl ;
1895    status = -1 ;
1896    return status ;
1897  }
1898  string str( buf ) ;
1899  int pos1 = str.find( '=' ) + 1 ;
1900  int pos2 = str.find( '/' ) ;
1901  //cout << "string: " << str << endl ;
1902  //cout << "value:  " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
1903  v = atof( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
1904  //cout << "v = " << v << endl ;
1905
1906  return status ;
1907}
1908
1909int NROFITSDataset::readTable( char *v, char *name, int clen, int idx )
1910{
1911  //
1912  // Read 'name' attribute defined as char from the idx-th row
1913  // of the FITS Scan Record
1914  //
1915  int status = movePointer( name, idx ) ;
1916  if ( status < 0 )
1917    return status ;
1918
1919  // get length of char
1920  int index = -1 ;
1921  for ( int i = 0 ; i < numField_ ; i++ ) {
1922    if ( names_[i] == name ) {
1923      index = i ;
1924      break ;
1925    }
1926  }
1927
1928  int xsize = sizes_[index] ;
1929
1930  // read data
1931  if ( xsize < clen ) {
1932    fread( v, 1, xsize, fp_ ) ;
1933    //v[xsize] = '\0' ;
1934  }
1935  else {
1936    fread( v, 1, clen-1, fp_ ) ;
1937    //v[clen-1] = '\0' ;
1938  }
1939
1940  return status ;
1941}
1942
1943int NROFITSDataset::readTable( int &v, char *name, int b, int idx )
1944{
1945  //
1946  // Read 'name' attribute defined as int from the idx-th row
1947  // of the FITS Scan Record
1948  //
1949  int status = movePointer( name, idx ) ;
1950  if ( status < 0 )
1951    return status ;
1952
1953  // read data
1954  fread( &v, sizeof(int), 1, fp_ ) ;
1955  if ( b == 0 )
1956    convertEndian( v ) ;
1957 
1958  return status ;
1959}
1960
1961int NROFITSDataset::readTable( float &v, char *name, int b, int idx )
1962{
1963  //
1964  // Read 'name' attribute defined as float from the idx-th row
1965  // of the FITS Scan Record
1966  //
1967  int status = movePointer( name, idx ) ;
1968  if ( status < 0 )
1969    return status ;
1970
1971  // read data
1972  fread( &v, sizeof(float), 1, fp_ ) ;
1973  if ( b == 0 )
1974    convertEndian( v ) ;
1975
1976  return status ;
1977}
1978
1979int NROFITSDataset::readTable( double &v, char *name, int b, int idx )
1980{
1981  //
1982  // Read 'name' attribute defined as double from the idx-th row
1983  // of the FITS Scan Record
1984  //
1985  int status = movePointer( name, idx ) ;
1986  if ( status < 0 )
1987    return status ;
1988
1989  // read data
1990  fread( &v, sizeof(double), 1, fp_ ) ;
1991  if ( b == 0 )
1992    convertEndian( v ) ; 
1993
1994  return status ;
1995}
1996
1997int NROFITSDataset::readTable( vector<char *> &v, char *name, int idx )
1998{
1999  //
2000  // Read 'name' attribute defined as char array from the FITS Scan Record
2001  //
2002  int status = movePointer( name, idx ) ;
2003  if ( status < 0 )
2004    return status ;
2005
2006  // get length of char
2007  int index = -1 ;
2008  for ( int i = 0 ; i < numField_ ; i++ ) {
2009    if ( names_[i] == name ) {
2010      index = i ;
2011      break ;
2012    }
2013  }
2014
2015  int xsize = sizes_[index] ;
2016
2017  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2018    int clen = strlen( v[i] ) ;
2019    if ( clen > xsize ) {
2020      fread( v[i], 1, xsize, fp_ ) ;
2021      //v[i][xsize] = '\0' ;
2022    }
2023    else {
2024      fread( v[i], 1, clen, fp_ ) ;
2025      //v[i][clen-1] = '\0' ;
2026    }
2027    //cout << "v[" << i << "] = " << v[i] << endl ;
2028  }
2029
2030  return status ;
2031}
2032
2033int NROFITSDataset::readTable( vector<int> &v, char *name, int b, int idx )
2034{
2035  //
2036  // Read 'name' attribute defined as int array from the FITS Scan Record
2037  //
2038  int status = movePointer( name, idx ) ;
2039  if ( status < 0 )
2040    return status ;
2041
2042  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2043    fread( &v[i], 1, sizeof(int), fp_ ) ;
2044    if ( b == 0 )
2045      convertEndian( v[i] ) ;
2046    //cout << "v[" << i << "] = " << v[i] << endl ;
2047  }
2048
2049  return status ;
2050}
2051
2052int NROFITSDataset::readTable( vector<float> &v, char *name, int b, int idx )
2053{
2054  //
2055  // Read 'name' attribute defined as float array from the FITS Scan Record
2056  //
2057  int status = movePointer( name, idx ) ;
2058  if ( status < 0 )
2059    return status ;
2060
2061  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2062    fread( &v[i], 1, sizeof(float), fp_ ) ;
2063    if ( b == 0 )
2064      convertEndian( v[i] ) ;
2065    //cout << "v[" << i << "] = " << v[i] << endl ;
2066  }
2067
2068  return status ;
2069}
2070
2071int NROFITSDataset::readTable( vector<double> &v, char *name, int b, int idx )
2072{
2073  //
2074  // Read 'name' attribute defined as double array from the FITS Scan Record
2075  //
2076  int status = movePointer( name, idx ) ;
2077  if ( status < 0 )
2078    return status ;
2079
2080  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2081    fread( &v[i], 1, sizeof(double), fp_ ) ;
2082    if ( b == 0 )
2083      convertEndian( v[i] ) ;
2084    //cout << "v[" << i << "] = " << v[i] << endl ;
2085  }
2086
2087  return status ;
2088}
2089
2090int NROFITSDataset::readColumn( vector<string> &v, char *name, int idx )
2091{
2092  //
2093  // Read idx-th column of ARRYTP-dependent 'name' attributes
2094  // defined as char array from the FITS Scan Record
2095  //
2096  int status = movePointer( name ) ;
2097  if ( status < 0 )
2098    return status ;
2099
2100  // get length of char
2101  int index = -1 ;
2102  for ( int i = 0 ; i < numField_ ; i++ ) {
2103    if ( names_[i] == name ) {
2104      index = i ;
2105      break ;
2106    }
2107  }
2108
2109  int xsize = sizes_[index] ;
2110
2111  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2112    int offset = scanLen_ * arrayid_[i] + xsize * idx ;
2113    fseek( fp_, offset, SEEK_CUR ) ;
2114//     int clen = (int)strlen( v[i] ) ;
2115//     if ( clen > xsize ) {
2116//       fread( v[i], 1, xsize, fp_ ) ;
2117//       //v[i][xsize] = '\0' ;
2118//     }
2119//     else {
2120//       fread( v[i], 1, clen-1, fp_ ) ;
2121//       //v[i][clen-1] = '\0' ;
2122//     }
2123    char c[xsize+1] ;
2124    fread( c, 1, xsize, fp_ ) ;
2125    c[xsize] = '\0' ;
2126    v[i] = string( c ) ;
2127    //cout << "v[" << i << "] = \'" << v[i] << "\'" << endl ;
2128    fseek( fp_, -xsize-offset, SEEK_CUR ) ;
2129  }
2130
2131  return status ;
2132}
2133
2134int NROFITSDataset::readColumn( vector<int> &v, char *name, int b, int idx )
2135{
2136  //
2137  // Read idx-th column of ARRYTP-dependent 'name' attributes
2138  // defined as int array from the FITS Scan Record
2139  //
2140  int status = movePointer( name ) ;
2141  if ( status < 0 )
2142    return status ;
2143
2144  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2145    int offset = scanLen_ * arrayid_[i] + sizeof(int) * idx ;
2146    fseek( fp_, offset, SEEK_CUR ) ;
2147    fread( &v[i], 1, sizeof(int), fp_ ) ;
2148    if ( b == 0 )
2149      convertEndian( v[i] ) ;
2150    //cout << "v[" << i << "] = " << v[i] << endl ;
2151    fseek( fp_, -sizeof(int)-offset, SEEK_CUR ) ;
2152  }
2153
2154  return status ;
2155}
2156
2157int NROFITSDataset::readColumn( vector<float> &v, char *name, int b, int idx )
2158{
2159  //
2160  // Read idx-th column of ARRYTP-dependent 'name' attributes
2161  // defined as float array from the FITS Scan Record
2162  //
2163  int status = movePointer( name ) ;
2164  if ( status < 0 )
2165    return status ;
2166
2167  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2168    int offset = scanLen_ * arrayid_[i] + sizeof(float) * idx ;
2169    fseek( fp_, offset, SEEK_CUR ) ;
2170    fread( &v[i], 1, sizeof(float), fp_ ) ;
2171    if ( b == 0 )
2172      convertEndian( v[i] ) ;
2173    //cout << "v[" << i << "] = " << v[i] << endl ;
2174    fseek( fp_, -sizeof(float)-offset, SEEK_CUR ) ;
2175  }
2176
2177  return status ;
2178}
2179
2180int NROFITSDataset::readColumn( vector<double> &v, char *name, int b, int idx )
2181{
2182  //
2183  // Read idx-th column of ARRYTP-dependent 'name' attributes
2184  // defined as double array from the FITS Scan Record
2185  //
2186  int status = movePointer( name ) ;
2187  if ( status < 0 )
2188    return status ;
2189
2190  for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2191    int offset = scanLen_ * arrayid_[i] + sizeof(double) * idx ;
2192    fseek( fp_, offset, SEEK_CUR ) ;
2193    fread( &v[i], 1, sizeof(double), fp_ ) ;
2194    if ( b == 0 )
2195      convertEndian( v[i] ) ;
2196    //cout << "offset = " << offset << ", v[" << i << "] = " << v[i] << endl ;
2197    fseek( fp_, -sizeof(double)-offset, SEEK_CUR ) ;
2198  }
2199
2200//   //cout << "v: " << endl ;
2201//   for ( vector<double>::iterator i = v.begin() ; i != v.end() ; i++ )
2202//     //cout << *i << " " ;
2203//   //cout << endl ;
2204 
2205  return status ;
2206}
2207
2208uInt NROFITSDataset::getArrayId( string type )
2209{
2210  uInt ib = 99;
2211  uInt len0 = type.size();
2212  for (uInt i = 0 ; i < arrayid_.size() ; i++) {
2213    uInt len = ARYTP[i].size();
2214    if ( len0 == len && type.compare( 0, len, ARYTP[i], 0, len ) == 0 ) {
2215      ib = i ;
2216      break ;
2217    }
2218  }
2219  return ib ;
2220}
2221
2222double NROFITSDataset::getStartIntTime( int i )
2223{
2224  double v ;
2225  readTable( v, "MJDST", same_, i ) ;
2226  return v/86400.0 ;
2227}
2228
2229double NROFITSDataset::getScanTime( int i )
2230{
2231  double startTime = getStartIntTime( i ) ;
2232  double interval = getIPTIM() ;
2233  interval /= 86400.0 ;
2234  return startTime+0.5*interval ;
2235}
2236
2237uInt NROFITSDataset::getPolNo( int irow )
2238{
2239  char rx[9] ;
2240  readTable( rx, "RX", 8, irow ) ;
2241  rx[8] = '\0' ;
2242  //cout << rx << endl ;
2243  return polNoFromRX( rx ) ;
2244}
2245
2246int NROFITSDataset::movePointer( char *name, int idx )
2247{
2248  // find offset
2249  long offset = getOffset( name ) ;
2250  if ( offset == -1 ) {
2251    //cerr << "Error, " << name << " is not found in the name list." << endl ;
2252    return -1 ;
2253  }
2254
2255  offset += (long)(idx * scanLen_) ;
2256
2257  //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2258  fseek( fp_, FITS_HEADER_SIZE+offset, SEEK_SET ) ;
2259
2260  return 0 ;
2261}
2262
2263// double NROFITSDataset::toLSR( double v, double t, double x, double y )
2264// {
2265//   return v ;
2266// }
Note: See TracBrowser for help on using the repository browser.