AIPS HELP file for CLCORPRM in 31DEC24
As of Sat Oct 5 15:38:29 2024
HELP SECTION
CLCORPRM
Type: Adverb (Real(20))
Use: The task specific parameters for CLCOR. The details depend on the
OPCODE selected.
Null value: Usually 0.
Tasks:
CLCOR....Applies user-selected corrections to the calibration CL
table. CLCORPRM usage depends on OPCODE - see below.
DFCOR....Applies differential corrections to the calibration CL table
RLCOR....Corrects a data set for R-L phase differences.
CLCORPRM provides right-left phase difference for IFs BIF
to EIF (as OPCODE 'POLR')
RLDIF....Determines Right minus Left phase difference, corrects cal
files. CLCORPRM is an output adverb giving right-left
phase differences if SPECTRAL<0 suitable for input to
CLCOR OPCODE 'POLR'.
XYDIF....Find/apply X minus Y linear polarization phase difference.
CLCORPRM is an output adverb giving X-Y phase differences if
SPECTRAL<0 suitable for input to CLCOR OPCODE 'POLR'.
Details for CLCOR: Contents are given as a function of the value of
OPCODE. Note that more details and examples are given in the CLCOR help
'POLR' => Modify Right-Left phase difference using phases in CLCORPRM
(deg); up to 20 IFs may be processed at a time.
'PHAS' => Rotate phase of residual gain by CLCORPRM (deg); up to 20
IFs may be processed at a time
'RATE' => Rotate phase of residual gain versus time. CLCORPRM(1)
degrees constant term. CLCORPRM(2) = rate of change of phase
(degrees/day). CLCORPRM(3) - (6) = day, hr, min, sec at
which the "zero" phase (CLCORPRM(1)) is specified.
'OPAC' => apply atmospheric opacity amplitude corrections using zenith
opacity of CLCORPRM(1) nepers.
'ADEL' => Correct phases, delays and rates for neutral atmospheric
delay.
CLCORPRM(1) = total pressure (mbars) at station, NOT at sea
level.
CLCORPRM(2) = partial pressure of water.
CLCORPRM(3) = Temperature (C)
CLCORPRM(4) = Tropospheric lapse rate (K/km) (should be
negative)
CLCORPRM(5) = Height of tropopause (km)
CLCORPRM(6) = Scale height of water vapor (km).
'GAIN' => Correct using polynominal gain curve for antenna
**VOLTAGE** gain as a function of the zenith angle (ZA) in
degrees.
correction = CLCORPRM(1) + ZA * CLCORPRM(2) +
ZA * ZA * CLCORPRM(3) etc.
'POGN' => Correct using polynominal gain curve for antenna
** POWER ** gain as a function of the zenith angle (ZA) in
degrees.
correction = CLCORPRM(1) + ZA * CLCORPRM(2) +
ZA * ZA * CLCORPRM(3) etc.
'CLOC' => Correct residual delay and model parms for the effects of a
linear clock drift at a particular antenna.
CLCORPRM(1) = rate of change of station clock (nanosec/day)
CLCORPRM(2) = clock value at the "zero" time specified by
CLCORPRM(3)-(6) (nanosec)
CLCORPRM(3) - (6) = day, hr, min, sec at which the "zero"
clock (CLCORPRM(2)) is specified.
CLCORPRM(7) : correction has three modes:
if = 0 then the clock drift is added as a small
correction and CLCORPRM(2) is ignored.
if = 1 then the total correction set by the CLCORPRMs is
added.
if = 2 then the values present in the CL table are
replaced by those defined by
CLCORPRM(1)-CLCORPRM(6).
'PANG' => Add or remove parallactic anglecorrections from CL table
entries. CLCORPRM(1) > 0 => Add corrections
CLCORPRM(1) =< 0 => Remove corrections
'PONT' => Correct for predictable pointing offset of an antenna.
CLCORPRM(1) is the linear rate of change of antenna gain as
the pointing drifts.
'IONS' => Make ionospheric Faraday rotation corrections using one of
several models. CLCORPRM(1) = Model type: 1 = Chiu model
CLCORPRM(3) = Sunspot no.
'ANTP' => Correct antenna and/or source position; antenna / source
corrections are values to be added to the old positions in
meters / sec of arc. The antenna (only one) must be
specified in the case of antenna correction. The source
(only one) must be specified in the case of source
correction.
!!!!!!!!!!!!!!!!!!! ATTENTION !!!!!!!!!!!!!!!!!!
Starting June 2001 CLCOR corrects the AN and/or the SU
table, if the relevant correction of the antenna or/and the
source is carried out. So the application of the corrected
CL table is required to match the data. Once CLCOR has been
run, then the CL table must not be deleted and the
correction can be undone only by doing it again with
opposite sign. If you might forget applying this, then use
the task SPLAT to apply it immediately.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1 = "X" correction in meters.
2 = "Y" correction in meters.
3 = "Z" correction in meters.
4 No longer used
5 = Correction in the picture plane towards the RA direction
in sec of arc. So it is RA correction, multiplied by
COS(DECL).
6 = Declination correction in sec of arc
7 No longer used
8 = drift of the source right ascension in mas/hr
9 = drift of the source declination in mas/hr
CLCORPRM(8,9) are used if INFILE is blank. The source
position is shifted by CLCORP(8,9) times the time where the
time is relative to zero hours IAT on the start date of the
observation.
'ANTC' => Correct antenna and/or source position; The same as ANTP
except that ANTP corrects the apparent source position while
ANTC corrects the position of epoch.
'PCAL' => MkIII manual phase cal; replace the gain correction in the
CL table with unit vectors with phases given in CLCORPRM.
CLCORPRM(1) corresponds to BIF etc. Phases are given in
degrees.
'PCFX' => Patch up missing phase cals. CLCORPRM gives the expected
relationship between phase cals and uses any non blanked
values.
'SBDL' => Add a delay to the IF residual delays. Values given in
CLCORPRM correspond to IFs BIF, BIF+1,... EIF in nanosec.
'MBDL' => Change the multi-band delay by introducing the corresponded
slope at phase vs IF frequency dependence. Values given in
CLCORPRM correspond to IFs BIF, BIF+1,... EIF in nanosec.
'SSLO' => Correct the phase only for an incorrect frequency used to
calculate the phase at the VLA. CLCORPRM(1) gives the
frequency error in MHz.
'ANAX' => Correct the delay, rate and phase for antenna axis offset.
CLCORPRM is an array of axis offsets in meters corresponding
to the antenna list given. Note that with 'ANAX' the antenna
table is also modified.
'ATMO' => apply atmospheric and clock delay corrections using the
information for the antennas and times given in the input
file. The file format:
The first line is number of the data rows. DELZN adds
the date and the JD in case the DELZN file has a
different reference date from the current input file.
The remaining lines specify the data with one line for
each antenna/time. Each line is formatted thus:
Column 1: antenna name; The name can be done either by
the antenna number (two digits) or by the
relevant antenna name (two symbols)
Column 2-5: day, hour, min, sec
Column 6: zenith atmosphere delay, in cm
Column 7: clock delay, in cm
Column 8: derivative of the zenith atmosphere delay, in
sec/sec*1.0E14
Column 9: derivative of the clock delay (clock drift), in
sec/sec*1.0E14
This format is in agreement with the DELZN output.
The values of the delays and its derivatives should
correspond to the desired correction. The corrections
for the interferometer delay and phase are added (with
sign) to the relevant columns of the CL table, which in
turn are added to the data when calibration is applied.
CLCORPRM(1): What correction to make in CL table?
0 => only atmosphere
1 => atmosphere + clocks
'TROP' => Similar to 'ATMO' except the zenith atmospheric (column 6 in
INFILE) is the TOTAL zenith delay, (typically between
150-400cm). This will be subtracted from the atmospheric
delay applied by the correlator, which is in the ATMOS
column in the CL table. Make sure there is something
sensible in this column (i.e. loosely between 5-15 nsec) or
you will get nonsense. NOTE THAT EXPERIMENTS CORRELATED
BETWEEN NOVEMBER 2002 AND FEBRUARY 2004 AT THE VLBA
CORRELATOR THE ATMOS COLUMN IN THE CL, MC, AND IM TABLES ARE
WRONG, SO YOU SHOULD NOT USE THIS OPTION. Other than this
'TROP' runs identically to 'ATMO'
'DISP' => Correction of dispersive delays as a function of elevation;
file contains:
The first line is number of the data rows.
The remaining lines specify the data with one line for
each antenna/time as:
Column 1: antenna number;
Column 2-5: day, hour, min, sec
Column 6: zenith dispersive delay, in 1/cm
Column 7: derivative of the zenith dispersive delay,
n sec/m^2/sec*1.0E14
'SUND' => Correction of the extra time delay caused by the bending of
the light ray passing through the gravitational field of the
solar system planets. The VLBA correlator (before September
2004) carries out this correction considering the source
locates at the infinity. This option removed the correlator
correction and put the correction taking into account the
actual position of the source inside of the solar system.
The correction is carried out only for the one source
(located at the solar system) given at the adverb SOURCE.
CLCORPRM(1): What correction to make in CL table?
0 => correct actual space craft position minus infinite
position. The infinite position correction carried
out by correlator is substituted by actual space
craft position.
1 => correct just actual space craft position
CLCORPRM(2): Print out the test data? 0 => yes, 1 => no
CLCORPRM(3): Print the additional test data?
0 => no: If CLCORPRM(2) equal 0 then delays for each
antenna/time are printed out
1 => yes: If CLCORPRM(2) equal 0 then positions of the
spacecraft, the first planet in the list,
baseline are printed out additionally
'EOPS' => Correction of UT1-UTC and the Earth's pole position. The
UT1-UTC and the Earth's pole position need to be corrected
if the wrong values used during correlation.
OPTYPE=EOPS can be used only for data sets with proper CT
tables, namely data from the old VLBA correlator, DiFX
correlators, and, since November 2015, the EVLA.
The actual EOP is interpolated for the time of observation
The default behavior, with CLCORPRM=0, is to use the same
sequence of days ([-2,2] or [-1,3]) as is used for the
majority of cases of the correlator for data sets that were
not concatenated in FITLD. For special cases (mainly for
more than one day observation), CLCORPRM can be set to
interpolate over a different range of dates.
CLCORPRM(1) Number of days preceding OBSDAT (given in the
data header) NEGATIVE NOT APPROPRIATE
<= 0.0 => 1 (Use 0.1 if you want 0 as would be
appropriate if the first data start at time ~ 1
day)
CLCORPRM(2) number of days starting from the first selected
day <= 0.0 => 5 (Required to have real data for
the full date range, but there is harm if you do
not get enough for concatenated data sets)
'IONO' => apply ionosphere and clock delay corrections using the
information for the antennas and times given in the input
file. Format of data at the input file is identical to
OPCODE = 'ATMO'. The differences with 'ATMO':
1.Phase-delay has the opposite sign
2.Different mapping functions
for 'ATMO' 1/SIN(EL)
for 'IONO' 1/COS(Z_i), where
SIN(Z_i) = (R/(R+h_i)) * COS(EL)
R is the Earth radius
h_i is the ionosphere height
See Thomson, Moran & Swenson 2nd edition page 560, eq
13.139 and 13.140
EXPLAIN SECTION