Context Navigation

NROFITSDataset.cc@ 2777

Visit:

Last change on this file since 2777 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

Rev	Line
[1757]	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
	41	using namespace std ;
	42
	43	//header size (8*2880bytes)
	44	#define FITS_HEADER_SIZE 23040
	45
	46	// constructor
	47	NROFITSDataset::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_ ) ;
[2768]	65	sizes_.resize( numField_ ) ;
	66	offsets_.resize( numField_ ) ;
[1757]	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
	95	NROFITSDataset::~NROFITSDataset()
	96	{
	97	// close file
	98	close() ;
	99	}
	100
	101	// data initialization
	102	void 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
[1868]	182
	183	refFreq_.resize( ARYNM, 0.0 ) ;
[1757]	184	}
	185
	186	// fill data header
	187	int 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
	203	int 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
[2664]	623	//nro_debug_output( "RX", ARYNM, RX ) ;
[1757]	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
[2436]	630	// nro_debug_output( "HPBW", ARYNM, HPBW ) ;
[1757]	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
[2436]	637	// nro_debug_output( "EFFA", ARYNM, EFFA ) ;
[1757]	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
[2436]	644	// nro_debug_output( "EFFB", ARYNM, EFFB ) ;
[1757]	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
[2436]	651	// nro_debug_output( "EFFL", ARYNM, EFFL ) ;
[1757]	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
[2436]	658	// nro_debug_output( "EFSS", ARYNM, EFSS ) ;
[1757]	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
[2436]	665	// nro_debug_output( "GAIN", ARYNM, GAIN ) ;
[1757]	666	//
	667	if ( readColumn( HORN, "HORN" ) != 0 ) {
	668	os << LogIO::WARN << "Error while reading data HORN." << LogIO::POST ;
	669	return -1 ;
	670	}
	671	// DEBUG
[2436]	672	// nro_debug_output( "HORN", ARYNM, HORN ) ;
[1757]	673	//
	674	if ( readColumn( POLTP, "POLTP" ) != 0 ) {
	675	os << LogIO::WARN << "Error while reading data POLTP." << LogIO::POST ;
	676	return -1 ;
	677	}
	678	// DEBUG
[2436]	679	// nro_debug_output( "POLTP", ARYNM, POLTP ) ;
[1757]	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
[2436]	689	// nro_debug_output( "POLDR", ARYNM, POLDR ) ;
[1757]	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
[2436]	696	// nro_debug_output( "POLAN", ARYNM, POLAN ) ;
[1757]	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
[2436]	703	// nro_debug_output( "DFRQ", ARYNM, DFRQ ) ;
[1757]	704	//
	705	if ( readColumn( SIDBD, "SIDBD" ) != 0 ) {
	706	os << LogIO::WARN << "Error while reading data SIDBD." << LogIO::POST ;
	707	return -1 ;
	708	}
	709	// DEBUG
[2436]	710	// nro_debug_output( "SIDBD", ARYNM, SIDBD ) ;
[1757]	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
[2436]	717	// nro_debug_output( "REFN", ARYNM, REFN ) ;
[1757]	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
[2436]	724	// nro_debug_output( "IPINT", ARYNM, IPINT ) ;
[1757]	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
[2436]	731	// nro_debug_output( "MULTN", ARYNM, MULTN ) ;
[1757]	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
[2436]	738	// nro_debug_output( "MLTSCF", ARYNM, MLTSCF ) ;
[1757]	739	//
	740	if ( readColumn( LAGWIND, "LAGWIN" ) != 0 ) {
	741	os << LogIO::WARN << "Error while reading data LAGWIND." << LogIO::POST ;
	742	return -1 ;
	743	}
	744	// DEBUG
[2436]	745	// nro_debug_output( "LAGWIND", ARYNM, LAGWIND ) ;
[1757]	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
[2436]	752	// nro_debug_output( "BEBW", ARYNM, BEBW ) ;
[1757]	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
[2436]	759	// nro_debug_output( "BERES", ARYNM, BERES ) ;
[1757]	760	//
	761	if ( readColumn( CHWID, "CHWID", sameEndian ) != 0 ) {
	762	os << LogIO::WARN << "Error while reading data CHWID." << LogIO::POST ;
	763	return -1 ;
	764	}
[2436]	765	// DEBUG
	766	// nro_debug_output( "CHWID", ARYNM, CHWID ) ;
[1757]	767	//
	768	if ( readARRY() != 0 ) {
	769	os << LogIO::WARN << "Error while reading data ARRY." << LogIO::POST ;
	770	return -1 ;
	771	}
	772	// DEBUG
[2436]	773	// nro_debug_output( "ARRY", NRO_FITS_ARYMAX, ARRY ) ;
[1757]	774	//
	775	if ( readColumn( NFCAL, "NFCAL", sameEndian ) != 0 ) {
	776	os << LogIO::WARN << "Error while reading data NFCAL." << LogIO::POST ;
	777	return -1 ;
	778	}
[2436]	779	// DEBUG
	780	// nro_debug_output( "NFCAL", ARYNM, NFCAL ) ;
[1757]	781	//
	782	if ( readColumn( F0CAL, "F0CAL", sameEndian ) != 0 ) {
	783	os << LogIO::WARN << "Error while reading data F0CAL." << LogIO::POST ;
	784	return -1 ;
	785	}
[2436]	786	// DEBUG
	787	// nro_debug_output( "F0CAL", ARYNM, F0CAL ) ;
[1757]	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
[2436]	800	// nro_debug_output( "FQCAL", ARYNM, 10, FQCAL ) ;
[1757]	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
[2436]	813	// nro_debug_output( "CHCAL", ARYNM, 10, CHCAL ) ;
[1757]	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
[2436]	826	// nro_debug_output( "CWCAL", ARYNM, 10, CWCAL ) ;
[1757]	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
[2436]	862	// nro_debug_output( "DSBFC", ARYNM, DSBFC ) ;
[1757]	863	//
	864
	865	show() ;
	866
	867	return 0 ;
	868	}
	869
	870	int 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 {
[1868]	903	sprintf( record_->LAVST, "%4d%02d%02d%02d%02d%02d.000", itmp[0], itmp[1], itmp[2], itmp[3], itmp[4], itmp[5] ) ;
[1757]	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 ) ;
[2664]	1025	for (int j = strlen(record_->ARRYT)-1 ; j >= 0 ; j--) {
	1026	if (record_->ARRYT[j] == ' ') record_->ARRYT[j] = '\0';
	1027	else break;
	1028	}
[1757]	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	//
[1868]	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
[1757]	1210	return status ;
	1211	}
	1212
	1213	vector< 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
	1224	vector<double> NROFITSDataset::getSpectrum( int i )
	1225	{
	1226	vector<double> spec( chmax_, 0.0 ) ;
	1227	vector<double> specout( chmax_, 0.0 ) ;
[2766]	1228	const NRODataRecord *record = getRecord( i ) ;
[1757]	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 ;
[1869]	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	// }
[1757]	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
	1416	int NROFITSDataset::getIndex( int irow )
	1417	{
[2766]	1418	const NRODataRecord *record = getRecord( irow ) ;
[1757]	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
	1433	double 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
	1456	double 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
	1479	void NROFITSDataset::getField()
	1480	{
[2768]	1481	long offset = 0;
[1757]	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	}
[2768]	1511	string form = tmp ;
	1512	string::size_type spos = form.find( " " ) ;
[1757]	1513	if ( spos != string::npos )
[2768]	1514	form = form.substr( 0, spos ) ;
[1757]	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, " " ) ;
[2768]	1526	if ( form.find( "A" ) != string::npos ) {
	1527	//cout << "skip to get unit: name = " << form << endl ;
[1757]	1528	//strcpy( tmp, "none " ) ;
	1529	tmp = "none" ;
	1530	}
	1531	else {
[2768]	1532	//cout << "get unit: name = " << form << endl ;
[1757]	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 ) ;
[2768]	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] ;
[1757]	1561	}
	1562	}
	1563
	1564	void NROFITSDataset::fillARYTP()
	1565	{
	1566	string arry ;
	1567	int count = 0 ;
	1568	string arry1 ;
	1569	string arry2 ;
	1570	string arry3 ;
	1571	string arry4 ;
[2664]	1572	char arytp[4];
[1757]	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' ) {
[2664]	1579	for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
[1757]	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" ;
[2664]	1588	for ( int i = 0 ; i < 10 ; i++ ) {
[1757]	1589	if ( arry2[i] == '1' ) {
[2664]	1590	for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
	1591	sprintf( arytp, "W%d", i+1 ) ;
	1592	ARYTP[count++] = string( arytp ) ;
[1757]	1593	}
	1594	}
[2664]	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	}
[1757]	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' ) {
[2664]	1615	for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
[1757]	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' ) {
[2664]	1626	for (int j = 0 ; j < 4 ; j++) arytp[j] = '\0';
[1757]	1627	sprintf( arytp, "A%d", i+21 ) ;
	1628	ARYTP[count++] = string( arytp ) ;
	1629	}
	1630	}
[2664]	1631	//nro_debug_output("ARYTP", ARYTP.size(), ARYTP);
[1757]	1632	}
	1633
	1634	int 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
	1691	void 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
[2768]	1699	long offset = getOffset( "ARRYT" ) ;
[1757]	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
[2768]	1740	long NROFITSDataset::getOffset( char *name )
[1757]	1741	{
[2768]	1742	long offset = 0 ;
[1757]	1743	string sname( name ) ;
[2768]	1744	long j = -1 ;
	1745	for ( long i = 0 ; i < numField_ ; i++ ) {
[1757]	1746	// escape if name is found
	1747	//cout << "names_[" << i << "] = " << names_[i] << " sname = " << sname << endl ;
	1748	if ( names_[i] == sname ) {
[2768]	1749	j = i ;
[1757]	1750	break ;
	1751	}
	1752	}
	1753
[2768]	1754	offset = (j >= 0) ? offsets_[j] : j ;
[1757]	1755
	1756	return offset ;
	1757	}
	1758
	1759	int 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
	1793	int 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
[2763]	1828	int NROFITSDataset::readHeader( int &v, char name, int /b*/ )
[1757]	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
[2763]	1864	int NROFITSDataset::readHeader( float &v, char name, int /b*/ )
[1757]	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
[2763]	1898	int NROFITSDataset::readHeader( double &v, char name, int /b*/ )
[1757]	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
	1933	int 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	//
[2440]	1939	int status = movePointer( name, idx ) ;
	1940	if ( status < 0 )
	1941	return status ;
[1757]	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
[2768]	1952	int xsize = sizes_[index] ;
	1953
[1757]	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
	1967	int 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	//
[2440]	1973	int status = movePointer( name, idx ) ;
	1974	if ( status < 0 )
	1975	return status ;
[1757]	1976
	1977	// read data
	1978	fread( &v, sizeof(int), 1, fp_ ) ;
	1979	if ( b == 0 )
	1980	convertEndian( v ) ;
	1981
	1982	return status ;
	1983	}
	1984
	1985	int 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	//
[2440]	1991	int status = movePointer( name, idx ) ;
	1992	if ( status < 0 )
	1993	return status ;
[1757]	1994
	1995	// read data
	1996	fread( &v, sizeof(float), 1, fp_ ) ;
	1997	if ( b == 0 )
	1998	convertEndian( v ) ;
	1999
	2000	return status ;
	2001	}
	2002
	2003	int 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	//
[2440]	2009	int status = movePointer( name, idx ) ;
	2010	if ( status < 0 )
	2011	return status ;
[1757]	2012
	2013	// read data
	2014	fread( &v, sizeof(double), 1, fp_ ) ;
	2015	if ( b == 0 )
	2016	convertEndian( v ) ;
	2017
	2018	return status ;
	2019	}
	2020
	2021	int 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	//
[2440]	2026	int status = movePointer( name, idx ) ;
	2027	if ( status < 0 )
	2028	return status ;
[1757]	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
[2768]	2039	int xsize = sizes_[index] ;
	2040
[1757]	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
	2057	int 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	//
[2440]	2062	int status = movePointer( name, idx ) ;
	2063	if ( status < 0 )
	2064	return status ;
[1757]	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
	2076	int 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	//
[2440]	2081	int status = movePointer( name, idx ) ;
	2082	if ( status < 0 )
	2083	return status ;
[1757]	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
	2095	int 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	//
[2440]	2100	int status = movePointer( name, idx ) ;
	2101	if ( status < 0 )
	2102	return status ;
[1757]	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
	2114	int 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	//
[2440]	2120	int status = movePointer( name ) ;
	2121	if ( status < 0 )
	2122	return status ;
[1757]	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
[2768]	2133	int xsize = sizes_[index] ;
	2134
[1757]	2135	for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
[2440]	2136	int offset = scanLen_ * arrayid_[i] + xsize * idx ;
[1757]	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
	2158	int 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	//
[2440]	2164	int status = movePointer( name ) ;
	2165	if ( status < 0 )
	2166	return status ;
[1757]	2167
	2168	for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
[2440]	2169	int offset = scanLen_ * arrayid_[i] + sizeof(int) * idx ;
[1757]	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
	2181	int 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	//
[2440]	2187	int status = movePointer( name ) ;
	2188	if ( status < 0 )
	2189	return status ;
[1757]	2190
	2191	for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
[2440]	2192	int offset = scanLen_ * arrayid_[i] + sizeof(float) * idx ;
[1757]	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
	2204	int 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	//
[2440]	2210	int status = movePointer( name ) ;
	2211	if ( status < 0 )
	2212	return status ;
[1757]	2213
	2214	for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
[2440]	2215	int offset = scanLen_ * arrayid_[i] + sizeof(double) * idx ;
[1757]	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
	2232	uInt NROFITSDataset::getArrayId( string type )
	2233	{
[2664]	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 ) {
[1757]	2239	ib = i ;
	2240	break ;
	2241	}
	2242	}
	2243	return ib ;
	2244	}
[1868]	2245
	2246	double NROFITSDataset::getStartIntTime( int i )
	2247	{
	2248	double v ;
	2249	readTable( v, "MJDST", same_, i ) ;
	2250	return v/86400.0 ;
	2251	}
	2252
[2156]	2253	double NROFITSDataset::getScanTime( int i )
	2254	{
	2255	double startTime = getStartIntTime( i ) ;
[2203]	2256	double interval = getIPTIM() ;
[2156]	2257	interval /= 86400.0 ;
	2258	return startTime+0.5*interval ;
	2259	}
	2260
[2434]	2261	uInt 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
[2440]	2270	int NROFITSDataset::movePointer( char *name, int idx )
	2271	{
	2272	// find offset
[2768]	2273	long offset = getOffset( name ) ;
[2440]	2274	if ( offset == -1 ) {
	2275	//cerr << "Error, " << name << " is not found in the name list." << endl ;
	2276	return -1 ;
	2277	}
	2278
[2768]	2279	offset += (long)(idx * scanLen_) ;
[2440]	2280
	2281	//cout << "offset for " << name << " is " << offset << " bytes." << endl ;
[2441]	2282	fseek( fp_, FITS_HEADER_SIZE+offset, SEEK_SET ) ;
[2440]	2283
	2284	return 0 ;
	2285	}
	2286
[1868]	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.

Context Navigation

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

Download in other formats: