[daip] AIPS antenna table
Eric Greisen
egreisen at nrao.edu
Mon Dec 6 16:15:19 EST 2010
Eric R. Keto wrote:
> Dear Eric,
>
> Thanks for your prompt response. I was preparing to
> send you some more information, but the case got more
> curious as I did some experiments. It took me a while to
> sort it out.
>
> The report of the error came from one of our staff.
> He showed me the output:
>
> LISTR1: Task LISTR (release of 31DEC10) begins
> LISTR1: Doing no flagging this time
> LISTR1: ANREFM: NO ANTENNA OR ARRAY RECOGNIZED
> LISTR1: CATDIR: CAN'T WRITE ON CATNO 2 NONZERO STATUS
> LISTR1: ANREFM: ERROR 7 CHANGING WRIT STATUS
> LISTR1: ANTINI: ERROR 7 RETURNED BY ANREFM
> LISTR1: GETANS: ERROR 7 READING AN FILE 1
> LISTR1: MATXUV: GETANS ERROR 7 SEARCHING ANTENNA TABLES
> LISTR1: Purports to die of UNNATURAL causes
> LISTR1: rtdc3 31DEC10 TST: Cpu= 0.0 Real= 0
> AIPS 1: Resumes
> AIPS 1: RETURN CODE 2 RECEIVED: STOPPING
>
> and I could see from running PRTAB that the antenna table
> format in the newer version of AIPS was different. In my
> first e-mail to you, I was wondering if there was a document with
> the list of changes to the antenna table that we could try.
>
> I was preparing an e-mail for you in response to your questions.
> When I ran LISTR to get the error messages to copy to you, I had
> no problem. The difference seems to be that I used optype='MATX',
> and LISTR was fine. But when he used optype='', he got the error
> above. The help file for LISTR says that the default optype is MATX.
> In my 2006 version of AIPS, we have no problems at all, and
> optype='' and optype='MATX' produce identical results.
>
> I'm not sure we have a problem with our antenna table or not.
> I suppose it would be better if AIPS recognized our array name,
> and we otherwise updated our antenna file, but we seem to be
> functional for all optypes in LISTR except ''.
>
> I presume CASA can read the older antenna table as well as other
> older VLA data?
>
> As for your questions, you can see that we used
> 'SAO SMA' (easily changed) as the array name, and AN01, etc
> as the antenna names.
>
> If you have time to look this over, I would be interested
> to know what you find out.
>
> Below I put in more information than you might want to know.
>
> Also attached is a small FITS file that I was using for these
> tests.
>
> Regards, Eric Keto
>
>
>
>
>
> Here is our Antenna Table function in the IDL language.
>
> The functions fits_wint, fits_wstr, are used in all the
> table writing functions and only the Antenna table is
> giving us trouble. Those function look like this:
>
> function fits_wint,name,value,comment
> ;
> ; write fits integer
> ;
> common fits,str_buf,unit
> res=fits_buff(string(name,value,comment,format='(a,T9,"=",T11,i20,T32,"/",T35,a,T81)'))
> return,1
> end
>
>
> function fits_buff,str,init=init,flush=flush
> ;
> ; add str to buffer str_buf and write out when 2880 bytes
> ;
> common fits,str_buf,unit
> if keyword_set(init) then begin
> str_buf=''
> return,1
> endif
> ;;print,str
> ;print,'new size of buffer ',strlen(str_buf+str)
> if keyword_set(flush) then begin
> while strlen(str_buf+str) lt 2880 do begin
> str=str+string('',format='(a,T81)')
> endwhile
> endif
> str_buf=str_buf+str
> if strlen(str_buf) eq 2880 then begin
> ;print,'Now writing out the 2880 byte str_buf'
> writeu,unit,str_buf
> str_buf=''
> endif
> return,1
> end
>
>
> function fits_want,ns,ver,stnum,stname,stx,sty,stz,freq0,datobs,site
> ;-----------------------------------------------------------------------
> ; Append AIPS antenna table as a FITS table extension to the FITS file.
> ;
> ; Arguments:
> ; ns (input, I*4) Number of antennas.
> ; ver (input, i*4) version number of table.
> ; stnum (input, I*4 array) Numbers of antennas
> ; stname (input, C*8 array) Names of antennas.
> ; stx (input, R*8 array) Antenna x-coordinates (m).
> ; sty (input, R*8 array) Antenna y-coordinates (m).
> ; stz (input, R*8 array) Antenna z-coordinates (m).
> ;-----------------------------------------------------------------------
> ;
> common fits,str_buf,unit
>
> obs_date=fix(strsplit(datobs,'-',/extract))
> mjd=uti_date2mjd(obs_date[0],obs_date[1],obs_date[2])
> Tu = (mjd - 51544.5d)/36525.
> gst=(24110.54841d + 8640184.812866d*Tu + .093104d*Tu^2 - 6.2e-6*Tu^3) mod 86400.d
> gst=gst >0 ? gst:gst+86400.
> gstia0=gst/240.d
>
> wbytes = 70
>
>
> res=fits_buff('',/init)
> result=fits_wstr('XTENSION','BINTABLE','Extension type')
> result=fits_wint('BITPIX ',8,'Binary data')
> result=fits_wint('NAXIS ',2,'Table is a matrix ')
> result=fits_wint ('NAXIS1', wbytes, 'WIDTH OF TABLE IN BYTES')
> result=fits_wint ('NAXIS2', ns, 'NUMBER OF ENTRIES IN TABLE')
> result=fits_wint ('PCOUNT', 0, 'NO RANDOM PARAMETERS')
> result=fits_wint ('GCOUNT', 1, 'GROUP COUNT')
> result=fits_wint ('TFIELDS', 12, 'NUMBER OF FIELDS IN EACH ROW')
> result=fits_wstr ('EXTNAME', 'AIPS AN', 'AIPS ANTENNA TABLE')
> ;result=fits_wstr ('EXTNAME', 'ARRAY_GEOMETRY', 'AIPS ANTENNA TABLE')
> result=fits_wint ('EXTVER',ver, 'VERSION NUMBER OF TABLE')
>
> result=fits_wstr('TFORM1','8A ', 'Fortran format of field 1')
> result=fits_wstr('TTYPE1','ANNAME ','Type (heading) of field 1')
> result=fits_wstr('TUNIT1',' ', 'Physical units of field 1')
>
> result=fits_wstr('TFORM2','3D ', 'Fortran format of field 2')
> result=fits_wstr('TTYPE2','STABXYZ ','Type (heading) of field 2')
> result=fits_wstr('TUNIT2','METERS ', 'Physical units of field 2')
>
> result=fits_wstr('TFORM3','0D ', 'Fortran format of field 3')
> result=fits_wstr('TTYPE3','ORBPARM ','Type (heading) of field 3')
> result=fits_wstr('TUNIT3',' ', 'Physical units of field 3')
>
> result=fits_wstr('TFORM4','1J ', 'Fortran format of field 4')
> result=fits_wstr('TTYPE4','NOSTA ','Type (heading) of field 4')
> result=fits_wstr('TUNIT4',' ', 'Physical units of field 4')
>
> result=fits_wstr('TFORM5','1J ', 'Fortran format of field 5')
> result=fits_wstr('TTYPE5','MNTSTA ','Type (heading) of field 5')
> result=fits_wstr('TUNIT5',' ', 'Physical units of field 5')
>
> result=fits_wstr('TFORM6','1E ', 'Fortran format of field 6')
> result=fits_wstr('TTYPE6','STAXOF ','Type (heading) of field 6')
> result=fits_wstr('TUNIT6','METERS ', 'Physical units of field 6')
>
> result=fits_wstr('TFORM7','1A ', 'Fortran format of field 7')
> result=fits_wstr('TTYPE7','POLTYA ','Type (heading) of field 7')
> result=fits_wstr('TUNIT7','DEGREES ', 'Physical units of field 7')
>
> result=fits_wstr('TFORM8','1E ', 'Fortran format of field 8')
> result=fits_wstr('TTYPE8','POLAA ','Type (heading) of field 8')
> result=fits_wstr('TUNIT8','DEGREES ', 'Physical units of field 8')
>
> result=fits_wstr('TFORM9','2E ', 'Fortran format of field 9')
> result=fits_wstr('TTYPE9','POLCALA ','Type (heading) of field 9')
> result=fits_wstr('TUNIT9',' ', 'Physical units of field 9')
>
> result=fits_wstr('TFORM10','1A ', 'Fortran format of field 10')
> result=fits_wstr('TTYPE10','POLTYB ','Type (heading) of field 10')
> result=fits_wstr('TUNIT10',' ', 'Physical units of field 10')
>
> result=fits_wstr('TFORM11','1E ', 'Fortran format of field 11')
> result=fits_wstr('TTYPE11','POLAB ','Type (heading) of field 11')
> result=fits_wstr('TUNIT11','DEGREES ', 'Physical units of field 11')
>
> result=fits_wstr('TFORM12','2E ', 'Fortran format of field 12')
> result=fits_wstr('TTYPE12','POLCALB ','Type (heading) of field 12')
> result=fits_wstr('TUNIT12',' ', 'Physical units of field 12')
>
>
>
> ; The array center coordinates have the X axis pointing along
> ; longitude 0.0 (Greenwich and East). This is different from the
> ; convention used in STABXYZ above where the X axis points along
> ; 0 hour angle and the Y axis points along -6 HA (East).
> ; The difference between the 2 systems is that X and Y are interchanged
> ; and the sign of Y is reversed.
>
> if (site eq 'saosma') then begin
> ; This is the SMA Pad 1 location.
> result=fits_wdble('ARRAYX',-5.46242840757d6,'')
> result=fits_wdble('ARRAYY',-2.49196045618d6,'')
> result=fits_wdble('ARRAYZ', 2.28652665343d6,'')
> endif
> if (site eq 'ovromma') then begin
> ; This is the OVRO VLBA site which is not quite the
> ; correct location for the reference point of the
> ; millimeter array.
> result=fits_wdble('ARRAYX',-2.40612582165d6,'')
> result=fits_wdble('ARRAYY',-4.47312856557d6,'')
> result=fits_wdble('ARRAYZ', 3.85876003303d6,'')
> endif
>
> result=fits_wdble('GSTIA0',gstia0,'')
> result=fits_wdble('DEGPDY',3.6098564497557D+02,'')
> result=fits_wdble('FREQ ',freq0,'')
> result=fits_wstr ('RDATE ',datobs,'')
> result=fits_wdble('POLARX',0.d0,'')
> result=fits_wdble('POLARY',0.d0,'')
> result=fits_wdble('UT1UTC',0.d0,'')
> result=fits_wdble('DATUTC',0.d0,'')
> result=fits_wstr ('TIMSYS','UTC ','')
> result=fits_wstr ('ARRNAM','SMA ','')
> result=fits_wint ('NUMORB',0,'')
> result=fits_wint ('NOPCAL',0,'')
> ;result=fits_wint ('NUMORB',0,'')
> result=fits_wint ('FREQID',-1,'')
> result=fits_wdble('IATUTC',0.d0,'')
>
>
> if (site eq 'saosma') then result=fits_wstr('ARRNAM ','SAO SMA ','')
> if (site eq 'ovromma') then result=fits_wstr('ARRNAM ','OVRO MMA','')
>
> result=fits_wend ()
>
> ;;print,'stnum ',stnum
> ;;print,'stname ',stname
>
> m = ' '
>
> for i=0,ns-1 do begin
>
> ;ANNAME
> nchars = strlen(stname[i])
> for k = 0,nchars-1 do begin
> writeu,unit,strupcase(strmid(stname[i],k,1))
> endfor
> if (nchars lt 8) then begin
> for k = nchars,7 do begin
> writeu,unit,m
> endfor
> endif
>
>
> ;STABXYZ
> d = double(stx[i]) & byteorder,d,/DTOXDR & writeu,unit,d
> d = double(sty[i]) & byteorder,d,/DTOXDR & writeu,unit,d
> d = double(stz[i]) & byteorder,d,/DTOXDR & writeu,unit,d
>
> ;NOSTA
> ; j = long(stnum[i]) & byteorder,j,/HTONL & writeu,unit,j
> j = i+1L & byteorder,j,/HTONL & writeu,unit,j
>
> ;MNTSTA
> j = 0L & byteorder,j,/HTONL & writeu,unit,j
>
> ;ORBPARM
> ; d = 0.d0 & byteorder,d,/DTOXDR & writeu,unit,d
>
> ;STXOF
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
>
> ;POLTYA
> writeu,unit,'R'
>
> ;POLAA
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
>
> ;POLCALA
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
>
> ;POLTYB
> writeu,unit,'L'
>
> ;POLAB
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
>
> ;POLCALB
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
> f = 0.0 & byteorder,f,/FTOXDR & writeu,unit,f
>
>
> endfor
>
>
>
> for i = 1,2880 - wbytes*ns do begin
> j = byte(0) & writeu,unit,j
> endfor
>
>
>
> end
>
>
>
>
>
> ----------------------------------------------------------------------------------------------------
> ----------------------------------------------------------------------------------------------------
> ----------------------------------------------------------------------------------------------------
>
> Here is some relevant output.
>
>
>
> ----------------------------------------------------------------------------------------------------
> ----------------------------------------------------------------------------------------------------
> ----------------------------------------------------------------------------------------------------
> 2006 Version of AIPS 2006 2006 2006
> ----------------------------------------------------------------------------------------------------
> Here is the PRTAB output of the antenna table
>
>
> W3OH_SMA .CONT . 1 Disk= 1 AN Table version 1
> Title: AIPS AN
> Created by FITLD on 05-NOV-2010 14:35:27
> Last written by FITLD on 05-NOV-2010 14:35:27
> Ncol 12 Nrow 8 Sort cols:
> Table has 18 keyword-value pairs:
> ARRAYX = -5.4624284075700D+06
> ARRAYY = -2.4919604561800D+06
> ARRAYZ = 2.2865266534300D+06
> GSTIA0 = 3.1785746980141D+02
> DEGPDY = 3.6098564497557D+02
> FREQ = 3.5449500489732D+11
> RDATE = 20080809
> POLARX = 0.0000000000000D+00
> POLARY = 0.0000000000000D+00
> UT1UTC = 0.0000000000000D+00
> DATUTC = 0.0000000000000D+00
> TIMSYS = UTC
> ARRNAM = SMA
> NUMORB = 0
> NOPCAL = 0
> FREQID = -1
> IATUTC = 0.0000000000000D+00
> ARRNAM = SAO SMA
> Table format incompatable with FITS ASCII tables
>
> COL. NO. 1 2 3 4 5 6 7 8 9
> ROW ANNAME STABXYZ ORBPARM NOSTA MNTSTA STAXOF POLTYA POLAA POLCALA
> NUMBER METERS METERS DEGREES DEGREES
> 1 AN01 0.000000D+00 1 0 0.000E+00 R 0.000E+00 0.000E+00
> 1 0.000000D+00 0.000E+00
> 1 0.000000D+00
> 2 AN02 1.175428D+02 2 0 0.000E+00 R 0.000E+00 0.000E+00
> 2 3.454445D+02 0.000E+00
> 2 -2.499002D+02
> 3 AN03 4.842663D+01 3 0 0.000E+00 R 0.000E+00 0.000E+00
> 3 4.074886D+02 0.000E+00
> 3 -5.952677D+01
> 4 AN04 -3.702497D+01 4 0 0.000E+00 R 0.000E+00 0.000E+00
> 4 4.055341D+02 0.000E+00
> 4 1.495775D+02
> 5 AN05 5.338816D+01 5 0 0.000E+00 R 0.000E+00 0.000E+00
> 5 1.293543D+02 0.000E+00
> 5 -1.140387D+02
> 6 AN06 9.102487D+01 6 0 0.000E+00 R 0.000E+00 0.000E+00
> 6 4.790408D+02 0.000E+00
> 6 -1.457231D+02
> 7 AN07 1.027643D+02 7 0 0.000E+00 R 0.000E+00 0.000E+00
> 7 4.253761D+02 0.000E+00
> 7 -1.864468D+02
> 8 AN08 -4.031359D+01 8 0 0.000E+00 R 0.000E+00 0.000E+00
> 8 2.274532D+02 0.000E+00
> 8 1.636865D+02
> Type Q to stop, just hit RETURN to continue
>
> laelaps PRTAB(31DEC06) 1346 02-DEC-2010 10:12:08 Page 2
> COL. NO. 10 11 12
> ROW POLTYB POLAB POLCALB
> NUMBER DEGREES
> 1 L 0.000E+00 0.000E+00
> 1 0.000E+00
> 2 L 0.000E+00 0.000E+00
> 2 0.000E+00
> 3 L 0.000E+00 0.000E+00
> 3 0.000E+00
> 4 L 0.000E+00 0.000E+00
> 4 0.000E+00
> 5 L 0.000E+00 0.000E+00
> 5 0.000E+00
> 6 L 0.000E+00 0.000E+00
> 6 0.000E+00
> 7 L 0.000E+00 0.000E+00
> 7 0.000E+00
> 8 L 0.000E+00 0.000E+00
> 8 0.000E+00
> PRTAB1: Appears to have ended successfully
> PRTAB1: laelaps 31DEC06 TST: Cpu= .0 Real= 11
> AIPS 1: Resumes
>
>
>
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
>
> 2010 version of AIPS 2010 2010 2010
>
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
>
>
>
> Now the PRTAB output of the same data from the 2010 TST version of AIPS
>
>
> PRTAB1: Task PRTAB (release of 31DEC10) begins
> rtdc3 PRTAB(31DEC10) 1346 02-DEC-2010 10:53:15 Page 1
> W3OH .UVDATA. 1 Disk= 1 AN Table version 1
> Title: AIPS AN
> Created by FITLD on 02-DEC-2010 10:52:36
> Last written by FITLD on 02-DEC-2010 10:52:36
> Ncol 12 Nrow 8 Sort cols:
> Table has 18 keyword-value pairs:
> ARRAYX = -5.4624284075700D+06
> ARRAYY = -2.4919604561800D+06
> ARRAYZ = 2.2865266534300D+06
> GSTIA0 = 3.1785746980141D+02
> DEGPDY = 3.6098564497557D+02
> FREQ = 3.5449500489732D+11
> RDATE = 20080809
> POLARX = 0.0000000000000D+00
> POLARY = 0.0000000000000D+00
> UT1UTC = 0.0000000000000D+00
> DATUTC = 0.0000000000000D+00
> TIMSYS = UTC
> ARRNAM = SMA
> NUMORB = 0
> NOPCAL = 0
> FREQID = -1
> IATUTC = 0.0000000000000D+00
> ARRNAM = SAO SMA
> Table format incompatable with FITS ASCII tables
>
> COL. NO. 1 2 3 4 5 6 7 8 9 10 11
> ROW ANNAME STABXYZ ORBPARM NOSTA MNTSTA STAXOF POLTYA POLAA POLCALA POLTYB POLAB
> NUMBER METERS METERS DEGREES DEGREES DEGREES
> 1 AN01 0.000000D+00 1 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 1 0.000000D+00 0.000E+00
> 1 0.000000D+00
> 2 AN02 1.175428D+02 2 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 2 3.454445D+02 0.000E+00
> 2 -2.499002D+02
> 3 AN03 4.842663D+01 3 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 3 4.074886D+02 0.000E+00
> 3 -5.952677D+01
> 4 AN04 -3.702497D+01 4 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 4 4.055341D+02 0.000E+00
> 4 1.495775D+02
> 5 AN05 5.338816D+01 5 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 5 1.293543D+02 0.000E+00
> 5 -1.140387D+02
> 6 AN06 9.102487D+01 6 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 6 4.790408D+02 0.000E+00
> 6 -1.457231D+02
> 7 AN07 1.027643D+02 7 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 7 4.253761D+02 0.000E+00
> 7 -1.864468D+02
> 8 AN08 -4.031359D+01 8 0 0.000E+00 R 0.000E+00 0.000E+00 L 0.000E+00
> 8 2.274532D+02 0.000E+00
> 8 1.636865D+02
> Type Q to stop, just hit RETURN to continue
>
> rtdc3 PRTAB(31DEC10) 1346 02-DEC-2010 10:53:18 Page 2
> COL. NO. 12
> ROW POLCALB
> NUMBER
> 1 0.000E+00
> 1 0.000E+00
> 2 0.000E+00
> 2 0.000E+00
> 3 0.000E+00
> 3 0.000E+00
> 4 0.000E+00
> 4 0.000E+00
> 5 0.000E+00
> 5 0.000E+00
> 6 0.000E+00
> 6 0.000E+00
> 7 0.000E+00
> 7 0.000E+00
> 8 0.000E+00
> 8 0.000E+00
> PRTAB1: Appears to have ended successfully
> PRTAB1: rtdc3 31DEC10 TST: Cpu= 0.0 Real= 3
> AIPS 1: Resumes
>>
>
>
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
>
>
> And now I run LISTR ....
>
> AIPS 1: LISTR: Task to print UV data and calibration tables.
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INNAME 'W3OH_SMA' UV data (name).
> AIPS 1: INCLASS 'CONT' UV data (class).
> AIPS 1: INSEQ 1 UV data (seq. #). 0 => high
> AIPS 1: INDISK 1 Disk unit #. 0 => any
> AIPS 1: OPTYPE 'MATX' List type:
> AIPS 1: 'MATX','LIST','GAIN','SCAN'
> AIPS 1: 'EFST','SEFD' special GAINs
> AIPS 1: INEXT 'AN' CL, SN or TY table for 'GAIN'
> AIPS 1: INVER 0 CL, Sn or TY table version
> AIPS 1: SOURCES *all ' ' Source list
> AIPS 1: CALCODE ' ' Calibrator code ' '=>all
> AIPS 1: TIMERANG *all 0 Time range to list
> AIPS 1: STOKES ' ' Stokes type to list.
> AIPS 1: SELBAND -1 Bandwidth to select (kHz)
> AIPS 1: SELFREQ -1 Frequency to select (MHz)
> AIPS 1: FREQID -1 Freq. ID to select.
> AIPS 1: None selected => 1.
> AIPS 1: BIF 0 Lowest IF number 0=1
> AIPS 1: EIF 0 Highest IF number
> AIPS 1: BCHAN 1 Low channel number 0=>1
> AIPS 1: ECHAN 0 High channel number
> AIPS 1: ANTENNAS *all 0 Antennas to list
> AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
> AIPS 1: UVRANGE 0 0 UV range in kilolambda
> AIPS 1: SUBARRAY 0 Subarray, 0=>1
> AIPS 1: Cal. info for input:
> AIPS 1: DOCALIB -1 > 0 calibrate data & weights
> AIPS 1: > 99 do NOT calibrate weights
> AIPS 1: GAINUSE 0 CAL (CL or SN) table to apply
> AIPS 1: DOPOL -1 If >0 correct polarization.
> AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
> AIPS 1: BLVER -1 BL table to apply.
> AIPS 1: FLAGVER 0 Flag table version
> AIPS 1: DOBAND -1 If >0 apply bandpass cal.
> AIPS 1: Method used depends on value
> AIPS 1: of DOBAND (see HELP file).
> AIPS 1: BPVER -1 Bandpass table version
> AIPS 1: SMOOTH *all 0 Smoothing function. See
> AIPS 1: HELP SMOOTH for details.
> AIPS 1: DPARM *all 0 Control info:
> AIPS 1: (1) 0=amp, 1=phase, 2=rms
> AIPS 1: 3=amp+rms, 4=phase+rms
> AIPS 1: 5=amp+phase, 6=delay,
> AIPS 1: 7=rate, 8=SNR, 9=par. ang
> AIPS 1: 10 = Tsys, 11=elev,
> AIPS 1: 12 = multiband delay
> AIPS 1: 13 = Tant 14 = weights
> AIPS 1: N.B. Options 3 and 4 only
> AIPS 1: valid for OPTYPE = 'MATX'
> AIPS 1: (2) 0=vec avg, 1=scalar
> AIPS 1: (3) No. col (4-10) def=4.
> AIPS 1: (4) avg. time min (MATX,LIST)
> AIPS 1: *** significant
> AIPS 1: (5) Type of gain listings:
> AIPS 1: 0 => all with same scale
> AIPS 1: like DEC-10 listings;
> AIPS 1: 1 => scan/source info at
> AIPS 1: source/scan boundaries.
> AIPS 1: (6) Matrix scaling control:
> AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
>
>
>
> LISTR1: Task LISTR (release of 31DEC10) begins
> LISTR1: Doing no flagging this time
> LISTR1: ANREFM: NO ANTENNA OR ARRAY RECOGNIZED
> rtdc3 LISTR(31DEC10) 1346 02-DEC-2010 10:57:05 Page 1
> File = W3OH .UVDATA. 1 Vol = 1 Userid = 1346 Chans= 1 - 1 IF = 1
> Freq=354.495004897 GHz Ncor= 1 No. vis= 8085
> Stokes = I Subarray = 1
>
>
> Time = 0/12:25:35 to 0/12:36:28 Source = W3OH : 0000
> Flux = 0.7448 Jy, Calcode = , Freq =354.495004897 GHz
> Amplitudes, 1000 = 10.000 Jy, averaging type = Vector
> IPOL in upper right
>
> Ant -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8
> 1| 16 8 2 122 10 10 39
> 2| 58 49 46 83 213 19
> 3| 35 48 137 223 36
> 4| 45 79 42 108
> 5| 15 10 74
> 6| 56 33
> 7| 7
> 8|
> IPO|
> Ampscalar average of matrix = 5.805E-01(1.082E-01) sigma = 5.624E-01
> Vector average of upper data= 1.116E-01 sigma = 8.980E-01
>
> Time = 0/12:43:53 to 0/12:55:16 Source = W3OH : 0000
> Flux = 0.7448 Jy, Calcode = , Freq =354.495004897 GHz
> Amplitudes, 1000 = 10.000 Jy, averaging type = Vector
> IPOL in upper right
>
> Ant -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8
> 1| 30 34 58 198 25 9 10
> 2| 86 33 78 71 119 24
> 3| 62 65 175 199 4
> 4| 22 30 10 125
> 5| 19 8 97
> 6| 79 26
> 7| 42
> 8|
> IPO|
> Ampscalar average of matrix = 6.207E-01(1.065E-01) sigma = 5.535E-01
> Vector average of upper data= 7.861E-02 sigma = 7.360E-01
>
> Time = 0/13:15:33 to 0/13:26:56 Source = W3OH : 0000
> Flux = 0.7448 Jy, Calcode = , Freq =354.495004897 GHz
> Amplitudes, 1000 = 10.000 Jy, averaging type = Vector
> IPOL in upper right
>
>
> Ant -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8
> 1| 22 18 17 18 11 7 48
> 2| 49 50 24 54 76 31
> 3| 98 48 124 95 46
> 4| 17 77 5 53
> 5| 21 34 75
> 6| 197 14
> 7| 45
> 8|
> IPO|
> Ampscalar average of matrix = 4.911E-01(7.906E-02) sigma = 4.108E-01
> Vector average of upper data= 8.072E-02 sigma = 6.370E-01
>
> Time = 0/13:34:21 to 0/13:45:43 Source = W3OH : 0000
> Flux = 0.7448 Jy, Calcode = , Freq =354.495004897 GHz
> Amplitudes, 1000 = 10.000 Jy, averaging type = Vector
> Type Q to stop, just hit RETURN to continue
>
> ---------------------------------------------------------------------------------------------------
> ---------------------------------------------------------------------------------------------------
>
> AIPS 1: LISTR: Task to print UV data and calibration tables.
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INNAME 'W3OH_SMA' UV data (name).
> AIPS 1: INCLASS 'CONT' UV data (class).
> AIPS 1: INSEQ 1 UV data (seq. #). 0 => high
> AIPS 1: INDISK 1 Disk unit #. 0 => any
> AIPS 1: OPTYPE ' ' List type:
> AIPS 1: 'MATX','LIST','GAIN','SCAN'
> AIPS 1: 'EFST','SEFD' special GAINs
> AIPS 1: INEXT ' ' CL, SN or TY table for 'GAIN'
> AIPS 1: INVER 0 CL, Sn or TY table version
> AIPS 1: SOURCES *all ' ' Source list
> AIPS 1: CALCODE ' ' Calibrator code ' '=>all
> AIPS 1: TIMERANG *all 0 Time range to list
> AIPS 1: STOKES ' ' Stokes type to list.
> AIPS 1: SELBAND -1 Bandwidth to select (kHz)
> AIPS 1: SELFREQ -1 Frequency to select (MHz)
> AIPS 1: FREQID -1 Freq. ID to select.
> AIPS 1: None selected => 1.
> AIPS 1: BIF 0 Lowest IF number 0=1
> AIPS 1: EIF 0 Highest IF number
> AIPS 1: BCHAN 1 Low channel number 0=>1
> AIPS 1: ECHAN 0 High channel number
> AIPS 1: ANTENNAS *all 0 Antennas to list
> AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
> AIPS 1: UVRANGE 0 0 UV range in kilolambda
> AIPS 1: SUBARRAY 0 Subarray, 0=>1
> AIPS 1: Cal. info for input:
> AIPS 1: DOCALIB -1 > 0 calibrate data & weights
> AIPS 1: > 99 do NOT calibrate weights
> AIPS 1: GAINUSE 0 CAL (CL or SN) table to apply
> AIPS 1: DOPOL -1 If >0 correct polarization.
> AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
> AIPS 1: BLVER -1 BL table to apply.
> AIPS 1: FLAGVER 0 Flag table version
> AIPS 1: DOBAND -1 If >0 apply bandpass cal.
> AIPS 1: Method used depends on value
> AIPS 1: of DOBAND (see HELP file).
> AIPS 1: BPVER -1 Bandpass table version
> AIPS 1: SMOOTH *all 0 Smoothing function. See
> AIPS 1: HELP SMOOTH for details.
> AIPS 1: DPARM *all 0 Control info:
> AIPS 1: (1) 0=amp, 1=phase, 2=rms
> AIPS 1: 3=amp+rms, 4=phase+rms
> AIPS 1: 5=amp+phase, 6=delay,
> AIPS 1: 7=rate, 8=SNR, 9=par. ang
> AIPS 1: 10 = Tsys, 11=elev,
> AIPS 1: 12 = multiband delay
> AIPS 1: 13 = Tant 14 = weights
> AIPS 1: N.B. Options 3 and 4 only
> AIPS 1: valid for OPTYPE = 'MATX'
> AIPS 1: (2) 0=vec avg, 1=scalar
> AIPS 1: (3) No. col (4-10) def=4.
> AIPS 1: (4) avg. time min (MATX,LIST)
> AIPS 1: *** significant
> AIPS 1: (5) Type of gain listings:
> AIPS 1: 0 => all with same scale
> AIPS 1: like DEC-10 listings;
> AIPS 1: 1 => scan/source info at
> AIPS 1: source/scan boundaries.
> AIPS 1: (6) Matrix scaling control:
>
>
>> go
> LISTR1: Task LISTR (release of 31DEC10) begins
> LISTR1: Doing no flagging this time
> LISTR1: ANREFM: NO ANTENNA OR ARRAY RECOGNIZED
> LISTR1: CATDIR: CAN'T WRITE ON CATNO 2 NONZERO STATUS
> LISTR1: ANREFM: ERROR 7 CHANGING WRIT STATUS
> LISTR1: ANTINI: ERROR 7 RETURNED BY ANREFM
> LISTR1: GETANS: ERROR 7 READING AN FILE 1
> LISTR1: MATXUV: GETANS ERROR 7 SEARCHING ANTENNA TABLES
> LISTR1: Purports to die of UNNATURAL causes
> LISTR1: rtdc3 31DEC10 TST: Cpu= 0.0 Real= 0
> AIPS 1: Resumes
> AIPS 1: RETURN CODE 2 RECEIVED: STOPPING
>
>
>
>
>
>
>
>
>
>
> On Mon, Nov 29, 2010 at 02:59:37PM -0700, Eric Greisen wrote:
>> Eric R. Keto wrote:
>>> Dear Eric:
>>>
>>> About 6 or 7 years ago I wrote a program to write
>>> SMA data into AIPS format. And you were critically
>>> helpful in ironing out a few bugs. We have been
>>> using that program quite successfully. Now we find
>>> that we have a problem with the Antenna table
>>> and the newest version of AIPS.
>>>
>>> >From some minor research into this, I see that there
>>> have been some changes to the Antenna table in
>>> the last year or so. Still, I am having a little trouble
>>> reverse-engineering the problem from an example
>>> of a new table.
>>>
>>> Is there any available documentation on the changes
>>> that would help us update the format of our Antenna
>>> table?
>>>
>>> Regards,
>>>
>>> Eric Keto
>> A number of things were added but a translation on the fly routine was
>> also added so the old format should be fine. The one fly in that
>> ointment is the routine that tries to recognize what antenna you are
>> using so that we can check the handedness of the coordinates. What
>> antenna names are you using (and array name)? The handedness of all
>> arrays is now forced to RIGHT.
>>
>> XYZHAND and FRAME were added as header keywords both 8 character.
>> Also NO_IF an integer giving the number of IFs.
>>
>> Columns DIAMETER and BEAMFWHM were added. The DIAMETER is a scalar for
>> each antenna in meters. The other has dimension NO_IF for each antenna
>> and is in degrees/meter wavelength and is the FWHM of the single dish.
>> Note that these are not filled in by anyone yet.
>>
>> Note that sending me the error messages would have been useful.
>>
>> Eric Greisen
LISTR had bad logic that caused it to put 2 READ flags on the file for
an unrecognized OPTYPE. ANREFM cleared one and hung up on the 2nd. I
have fixed LISTR's logic. I also added 'SAO SMA' as a recognized array
name. Note that your files have ARRNAM twice - the first with 'SAO' and
the second with 'SAO SMA'. Please drop the first to avoid confusion.
This fix will be in 31DEC10 and 31DEC11 by tomorrow's MNJ.
Thanks,
Eric Greisen
More information about the Daip
mailing list