[daip] VLBA + Arecibo cal transfer
Shami Chatterjee
shami at astro.cornell.edu
Wed Sep 18 17:09:48 EDT 2002
Hi Joan,
This is an aips question, but since it involves ANTAB, I thought you might
be able to comment on it as well. Briefly: I'm reducing data for a VLBA +
Arecibo experiment. (And it's good data, according to POSSM etc!) VLBA
came with cal transfer, which I have mergecal'ed to TY 2, GC 2. For
Arecibo, Tapasi sent me a file with SEFD values for my run, with the
following note:
# These are actually SEFD values. With Gain=1, these could be used
# as Tsys for VLBI calibration. As Arecibo Telescope's gain is both
# Az and Za dependent, and AIPS can not handle such function, we
# have decided to provide the users with SEFD values.
# Shami, for your Expt. these numbers have been derived from
# the continuum record at 1405 MHz, two pols. have been
# averaged. I think, you can use it in ANTAB to creat
# a Tsys table.
Gain = 1
Tsys /
109 7:10:20.00 3.66561
109 7:10:30.00 3.66012
<etc>
Looking at the VLBA TY values, I have:
----
ROW TIME TIME INTER SOURCE ID ANTENNA SUBARRAY FREQ ID
1 07:05:11.0 00:01:39.0 5 2 1 1
2 07:07:34.0 00:00:53.0 1 2 1 1
...
84 09:27:49.5 00:01:00.0 3 2 1 1
----
for each antenna.
Following notes from Emmanuel Momjian, Tapasi and mostly from
"EXPLAIN ANTAB", I created a fake header and Tsys input table:
----
CONTROL INDEX = 'R1', 'L1', 'R2', 'L2', 'R3', 'L3', 'R4', 'L4' /
!
GAIN AR ALTAZ DPFU = 1.0 POLY = 1.0 /
! (Flat gain curve for AR, all IF's and polarizations,
! zenith point source sensitivity of 1 K/Jy)
!
TSYS AR FT=1.0 INDEX = 'R1', 'L1', 'R2', 'L2', 'R3', 'L3', 'R4', 'L4' /
!
109 7:10:20.00 3.665610 3.665610 3.665610 3.665610 3.665610 3.665610
3.665610 3.665610
109 7:10:30.00 3.660120 3.660120 3.660120 3.660120 3.660120 3.660120
3.660120 3.660120
<etc> - one entry every 10 seconds, till
109 9:20:50.00 3.407700 3.407700 3.407700 3.407700 3.407700 3.407700
3.407700 3.407700
----
where I just repeated the SEFD value Tapasi sent me 8 times, once for each
IF and poln, as I specify in the CONTROL card.
Now I load it:
mingus> ANTAB1: Task ANTAB (release of 31DEC02) begins
mingus> ANTAB1: ZTXOP2: using translated file name = [...]
mingus> ANTAB1: Gain curve read for station: AR
mingus> ANTAB1: 0 unidentified records skipped on input
mingus> ANTAB1: Appears to have ended successfully
Setting prtlev 1 echoes the full input table, so no format problems.
But when I use PRTAB, there is no new information in the table. If I send
it to a new TY table, it gets created, but is empty:
mingus> SNPLT1: Task SNPLT (release of 31DEC02) begins
mingus> SNPLT1: ERROR: SELECTED TABLE IS EMPTY
I can see several potential problems, in order of likelihood:
1) I made a mistake in the header.
(I've tried it with and without the CONTROL card: no difference)
2) The timestamps (every 10 seconds) do not match the exact scan times
aips wants in the table. (How can that be?)
3) Some confusion regarding the day number.
4) ...?
Is there an obvious problem, and how do I fix it?
My apologies for the bother, and many thanks for your help!
Shami
--
Shami Chatterjee
http://www.astro.cornell.edu/~shami/
Preudhomme's Law of Window Cleaning: It's on the other side.
More information about the Daip
mailing list