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

Last change on this file since 2774 was 2768, checked in by Takeshi Nakazato, 12 years ago

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: test_sdsave

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

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