* CALGAIN.001 v04jun2002 - S.T. Myers (NRAO) runfile to load procedures:
*    CALGAIN(I,TAU,FBAND) - see end of file. This procedure works on 
*      catalog entry I of INDISK, applies opacity correction for optical
*      depth TAU (if TAU>0.0), and then applies the U, K or Q band gain
*      curves for band FBAND (if FBAND = 'U', 'K', or 'Q').  This proc
*      uses the other procedures:
*    FIXKGAIN(I,J)
*    FIXQGAIN(I,J)
*    FIXUGAIN(I,J)
*
* Example usage: CALGAIN(1,0.05,'K')
*
Proc dummy
* Usage: do not call this procedure, it is only here to define the string
* fband for CALGAIN to use!
  string*1 fband
*
  Finish
* 
Proc FIXKGAIN(i,j)
* Usage: FIXKGAIN( CATNO, CLVER )
* I = catalog number of file; J = version number of CL table to correct
* Procedure to put in polynomial gain coefficients
* Compiled by B.Butler 2002-02-08 from Perley-Butler run 2001-11-10
* Same curves as in VLA.GAINS applied in FILLM
*
* Normalization close to 1 at maximum
*
* Back-propagate to June 2001, since that is the point at
* which the subreflector rotation v. elevation was turned on
* Use average gain curve for AN 4
*
  type 'Applying K band gain curves from Perley-Butler 2001-11-10'
  type 'Version BButler/SMyers 2002-02-08' 
  type 'Gain curves IF and Pol averaged'
  getn i
  gainver j
  opcode 'gain'
  source ''
*
* AVERAGE COEFF = 1.0429E+00, -2.8623E-03,  3.9044E-05, 0
*
 ANTEN 1 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0225E+00, -2.4650E-03,  4.2225E-05, 0; 
 GO CLCOR;

 ANTEN 2 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.5400E-01,  2.3688E-04,  2.2775E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 3 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0850E+00, -4.1050E-03,  4.1700E-05, 0; 
 WAIT CLCOR; GO CLCOR;

* AN4 was in barn during run, use average curve
 ANTEN 4 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0429E+00, -2.8623E-03,  3.9044E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 5 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0425E+00, -2.6525E-03,  3.5075E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 6 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0375E+00, -2.7100E-03,  3.7450E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 7 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0900E+00, -4.1675E-03,  3.8625E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 8 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0700E+00, -3.8175E-03,  4.3300E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 9 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0350E+00, -2.5700E-03,  3.7550E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 10 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0700E+00, -3.6375E-03,  3.9600E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 11 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.8800E-01, -1.0907E-03,  3.2150E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 12 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1325E+00, -6.1375E-03,  5.6625E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 13 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0550E+00, -2.9950E-03,  3.7325E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 14 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0975E+00, -4.6775E-03,  4.6750E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 15 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0400E+00, -4.0575E-03,  6.7950E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 16 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1125E+00, -5.2850E-03,  5.0975E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 17 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0600E+00, -3.1500E-03,  3.3675E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 18 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.6675E-01, -3.8510E-04,  2.8925E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 19 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0375E+00, -2.7425E-03,  3.9650E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 20 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0425E+00, -3.0425E-03,  4.2000E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 21 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.6275E-01, -1.2383E-04,  2.3875E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 22 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0800E+00, -4.1275E-03,  4.5425E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 23 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0125E+00, -1.7350E-03,  3.1000E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 24 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.8500E-01, -6.5701E-05,  1.2300E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 25 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0650E+00, -3.0625E-03,  3.1125E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 26 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.9300E-01, -1.0235E-03,  3.4800E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 27 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0875E+00, -4.6600E-03,  5.1800E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 28 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0350E+00, -3.0325E-03,  4.9525E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 WAIT CLCOR; 
*
 Finish
* 
Proc FIXQGAIN(i,j)
* Usage: FIXQGAIN( CATNO, CLVER )
* I = catalog number of file; J = version number of CL table to correct
* Procedure to put in polynomial gain coefficients
* Compiled by B.Butler 2002-02-08 from Perley-Butler run 2001-11-10
* Same curves as in VLA.GAINS applied in FILLM
*
* Normalization close to 1 at maximum
*
* Back-propagate to June 2001, since that is the point at
* which the subreflector rotation v. elevation was turned on
* Use average gain curve for AN 4, 20
* AN 7, 9, 15 still without Q-band Rx - no curves 
*
  type 'Applying Q band gain curves from Perley-Butler 2001-11-10'
  type 'Version BButler/SMyers 2002-02-08' 
  type 'Gain curves IF and Pol averaged'
  getn i
  gainver j
  opcode 'gain'
  source ''
*
* AVERAGE COEFF = 1.0787E+00, -5.5261E-03,  7.7206E-05, 0
*
 ANTEN 1 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1800E+00, -1.3900E-02,  2.1250E-04, 0; 
 GO CLCOR;

 ANTEN 2 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.8850E-01, -2.9875E-03,  7.3575E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 3 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0875E+00, -5.1525E-03,  6.3550E-05, 0; 
 WAIT CLCOR; GO CLCOR;

* AN4 was in barn during run, use average curve
 ANTEN 4 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0787E+00, -5.5261E-03,  7.7206E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 5 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.5075E-01,  8.0000E-04,  1.6003E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 6 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0875E+00, -5.2225E-03,  6.4925E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 8 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1425E+00, -9.1000E-03,  1.1847E-04, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 10 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1025E+00, -5.5650E-03,  6.0525E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 11 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0225E+00, -3.7450E-03,  7.1025E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 12 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1425E+00, -7.7475E-03,  8.9300E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 13 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0850E+00, -5.6100E-03,  7.3550E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 14 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0825E+00, -5.6175E-03,  7.7150E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 16 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0325E+00, -1.2130E-03,  8.1949E-06, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 17 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1100E+00, -7.3550E-03,  9.9325E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 18 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0475E+00, -5.5825E-03,  9.7275E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 19 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1150E+00, -6.5225E-03,  7.7475E-05, 0; 
 WAIT CLCOR; GO CLCOR;

* AN20 had severe pointing problems during run, use average curve
 ANTEN 20 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0787E+00, -5.5261E-03,  7.7206E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 21 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.7425E-01, -1.2132E-03,  4.2825E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 22 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0825E+00, -5.3000E-03,  6.8150E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 23 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0275E+00, -2.4933E-03,  4.0450E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 24 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0825E+00, -7.0600E-03,  1.1338E-04, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 25 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1125E+00, -6.4275E-03,  7.5775E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 26 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1300E+00, -7.6500E-03,  8.9100E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 27 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.1275E+00, -6.6525E-03,  7.1850E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 28 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0975E+00, -5.7825E-03,  7.1375E-05, 0;
 WAIT CLCOR; GO CLCOR;

 WAIT CLCOR; 
*
 Finish
* 
Proc FIXUGAIN(i,j)
* Usage: FIXUGAIN( CATNO, CLVER )
* I = catalog number of file; J = version number of CL table to correct
* Procedure to put in polynomial gain coefficients
* Compiled by B.Butler 2002-02-08 from Perley-Butler run 2001-11-10
* Same curves as in VLA.GAINS applied in FILLM
*
* Normalization close to 1 at maximum
*
  type 'Applying U band gain curves from Perley-Butler 2001-11-10'
  type 'Version BButler/SMyers 2002-02-08' 
  type 'Gain curves IF and Pol averaged'
  getn i
  gainver j
  opcode 'gain'
  source ''
*
* AVERAGE COEFF = 1.0324E+00, -1.4352E-03,  1.5845E-05, 0
*
 ANTEN 1 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.6250E-01,  8.5425E-04,  5.7450E-06, 0; 
 GO CLCOR;

 ANTEN 2 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.9250E-01, -4.5000E-04,  1.2125E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 3 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0400E+00, -1.6825E-03,  1.3150E-05, 0; 
 WAIT CLCOR; GO CLCOR;

* AN4 was in barn during run, use average curve
 ANTEN 4 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0324E+00, -1.4352E-03,  1.5845E-05, 0;
 WAIT CLCOR; GO CLCOR;

 ANTEN 5 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0350E+00, -1.5950E-03,  1.3700E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 6 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0200E+00, -1.1850E-03,  1.4475E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 7 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0425E+00, -2.0425E-03,  1.6450E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 8 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0375E+00, -1.5900E-03,  1.4575E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 9 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0150E+00, -9.2050E-04,  1.2700E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 10 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0500E+00, -1.9775E-03,  1.6775E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 11 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0175E+00, -1.3487E-03,  1.6382E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 12 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0600E+00, -2.3675E-03,  1.7750E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 13 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0175E+00, -9.8350E-04,  1.2347E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 14 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0425E+00, -1.8375E-03,  1.6550E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 15 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0150E+00, -6.3875E-04,  9.7800E-06, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 16 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0550E+00, -2.3100E-03,  1.8600E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 17 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0450E+00, -1.8975E-03,  1.8900E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 18 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.9575E-01, -8.5100E-04,  2.0100E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 19 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0525E+00, -2.4175E-03,  2.2200E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 20 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.2050E+00, -3.1300E-03,  3.2375E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 21 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.9600E-01, -4.0625E-04,  1.0275E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 22 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0450E+00, -1.8825E-03,  1.6675E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 23 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0225E+00, -8.4950E-04,  7.6725E-06, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 24 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0200E+00, -1.0663E-03,  1.3138E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 25 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0500E+00, -2.1575E-03,  1.7025E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 26 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 9.9300E-01, -1.5057E-03,  3.1075E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 27 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0350E+00, -1.6850E-03,  1.6950E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 ANTEN 28 0; STOKES = 'RL'; BIF=1; EIF=2; 
 CLCOR 1.0125E+00, -8.2600E-04,  1.0335E-05, 0; 
 WAIT CLCOR; GO CLCOR;

 WAIT CLCOR; 
*
 Finish
*
Proc calgain(i,x,fband)
* 
* Procedure to calibrate antennas for gain curve
* Copies CL1 to CL2 and applies opacity correction using optical
* depth given in scalar x, if x>0
* Then applies gain curve for U,K,Q band (input band in fband)
* using procedures FIXUGAIN,FIXKGAIN,FIXQGAIN provided above.
*
  task 'tacop'
  getn i
  inext 'cl'
  outdi indi
  geto i
  inver 1 
  outver 2; ncount 1
  keyword ''; keyval 0
  go tacop; wait tacop
*
  task 'clcor'
  if x>0.0 then opcode 'opac';
     type ' Correcting for opacity: tau =  ', x;
     clcor x,0; gainver 2; stokes ''; source '';
     go clcor; wait clcor;
  end
*
  if fband='K' then fixkgain(i,2); type ' Applied K gains'; end
  if fband='Q' then fixqgain(i,2); type ' Applied Q gains'; end
  if fband='U' then fixugain(i,2); type ' Applied U gains'; end
*
  type ' Completed CALGAIN '
*
  Finish
*