Proccal Operator Instructions



Steps in the PROCCAL calibration image processing:

cd to the IRAF directory, and start IRAF. Load the grasp and calplus packages.
cd to the control directory corresponding to the data tape containing the calibration episodes to be processed, then execute getcalonly.
If Proccal has not been run from this control directory before there should not be a log file (calibration.log) already there, but if there is the operator will be queried as to whether to quit or overwrite the existing log file (example).


Example (Operator inputs are in bold)
% cd ~/iraf

% cl

cl> grasp
           ******           GRASP: V2005.06           ******
           ****** GONG Reduction and Analysis Package ******
Bug reports and ideas for program changes (functions, algorithms, packaging,
user queries, help page explanations, task options, etc.) are more than welcome.

Any questions, comments, bug reports about GRASP should be E-mailed to:

                   grasp@noao.edu          Internet
               or  grasp                   Solarmail

     IRAF packages IMAGES, TV, PLOT and TABLES already been loaded.

      drift.      grhistory   gzip        model.      qasun.      velcal.
      ephem.      grtools.    inst.       movie       qcgong.     
      gongcor.    gunzip      merge.      pipeline.   suntrans.   
gr> calplus
 
      Load packages for calplus processing
 
      Load tasks for calplus processing
 
      amoeba_ML_mask          get_30_sec_cal          param_disp
      amoeba_mask             getavgdark              parps@
      cal_qa                  getavgdark_for_pin      plots
      check_guide             getcal_ML_plus          prep_dir
      chkcallog               getcalonly              rd_fl
      cleanup                 getcalonly_ML_noplot    solar_ML_plus
      conv_cal_fits_plus      getcalonly_noplot       solar_plus
      ftr_cleanup             getcalplus              solar_plus_30_sec
      ftr_cleanup20           getcalplus_30_sec       test_chkcal
      ftr_cont                getmask                 test_chkcal_2
      ftr_copy                getrms                  within_20
      ftr_copy_sandbox        getrms_ML_mask          
      ftr_tweak_copy          getrms_mask             
ca> cd /ftr/armet/ftr/field.tape.copies
ca> cd BB.060228
ca> getcalonly
 
    Current directory: /ftr/armet/ftr/field.tape.copies/BB.060228
    Proccal directory prefix is /ftr/armet/ftr/proccal
block_256 deleted
    GETCALPLUS, version 01.07.18 - Solaris
    run by armet on tarat on terminal 3 (05/10/06 10:44:31)
 
    paramfile = /home/armet/iraf/initplus.params/BB.par
    swpath    = /home/armet/iraf/tape.cal.plus/
 
.
.
.

calibration.log

    GETCALPLUS, version 01.07.18 - Solaris
    run by armet on tarat on terminal t (05/09/06 17:23:27)
 
    paramfile = /home/armet/iraf/initplus.params/CT.par
    swpath    = /home/armet/iraf/tape.cal.plus/
 
    Call getipermi with date = 060420, site = CT
    Back from getipermi with value of 3 , and PM date = 20060420
 
    The directory /ftr/armet/ftr/proccal/CT/cals.20060420 does not exist, create it.
 
              Site Dependent Parameters
              -------------------------
 
      Site                            Cerro Tololo
      First clampr column             15 
      Last clampr column              25 
      Phi                             -0.7
      Oscillation constant            -2061.3592
      X column start                  22 
      X column end                    225 
      Y column start                  4 
      Y column end                    242 
      Fndlmb minor axis estimate      370.
      Fndlmb major axis estimate      370.
      Fndlmb ellipse rotation angle   0.
      Fndlmb number of parameters     4 
      Fndlmb band width               8.
      Fndlmb threshold                0.
      Fndlmb percent of fit           1.
      Fndlmb x pixel length           1.
      Fndlmb y pixel length           1.
 
 
    Next dark id: CT060420125114
 
 
        CONV_CAL_FITS_PLUS (17:23:47)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060420125114	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060420125414	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060420125514	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060420125614	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060420125214	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060420125016	-> diff
 
        CONV_CAL_FITS_PLUS done(17:24:08)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:24:09)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.35 , yc = 522.646 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.214 , yc256 = 131.041 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:24:52)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.214
  y-coordinate of the center = 131.041
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.58  130.95   94.16   94.38    0.00   3.9E4  1.5E-3    4  747
blk_csun_6_band  111.58  130.95   94.16   94.38    0.00   3.9E4  1.5E-3    2  747
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:24:52)
 
 
        Call AMOEBA_MASK (17:24:58)
 
    AMOEBA_MASK done(17:25:59)
 
 
phi = -0.545 ,   sa = 0.027
        Amoeba output = avgphi=-0.54489456641798 avgsa=0.026720130133132
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:26:01)
 
        getcal("ceast","cwest","wlm","nbc",-0.545,0.027,3)
 
        GETCAL done (17:26:08)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 508.5077
 in radians [divided by -2061.3592 ] = -0.2466856334403
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.640724
 
 Should be close to pi, 2.640724  + -0.2466856334403  = 2.3940383665597
 
        VCOR done (17:26:20)
    Now that wlm, et al, have been calc'd, wfits them with 200604201251.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060420125114...Calibrated
 
    Next dark id: CT060421125214
 
 
        CONV_CAL_FITS_PLUS (17:26:31)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060421125214	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060421125514	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060421125614	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060421125714	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060421125314	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060421125116	-> diff
 
        CONV_CAL_FITS_PLUS done(17:26:51)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:26:52)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.338 , yc = 522.632 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.209 , yc256 = 131.071 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:27:32)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.209
  y-coordinate of the center = 131.071
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.62  130.89   94.14   94.37    0.00   4.2E4  1.6E-3    4  747
blk_csun_6_band  111.62  130.89   94.14   94.37    0.00   4.2E4  1.6E-3    2  747
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:27:33)
 
 
        Call AMOEBA_MASK (17:27:34)
 
    AMOEBA_MASK done(17:28:37)
 
 
phi = -0.549 ,   sa = 0.033
        Amoeba output = avgphi=-0.54896827906293 avgsa=0.033068968549681
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:28:40)
 
        getcal("ceast","cwest","wlm","nbc",-0.549,0.033,3)
 
        GETCAL done (17:28:47)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 511.5251
 in radians [divided by -2061.3592 ] = -0.24814942490372
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.64551
 
 Should be close to pi, 2.64551  + -0.24814942490372  = 2.3973605750963
 
        VCOR done (17:29:00)
    Now that wlm, et al, have been calc'd, wfits them with 200604211252.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060421125214...Calibrated
 
    Next dark id: CT060422125314
 
 
        CONV_CAL_FITS_PLUS (17:29:12)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060422125314	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060422125614	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060422125714	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060422125814	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060422125414	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060422125216	-> diff
 
        CONV_CAL_FITS_PLUS done(17:29:33)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:29:33)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.33 , yc = 522.627 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.246 , yc256 = 131.114 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:30:14)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.246
  y-coordinate of the center = 131.114
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.63  130.86   94.12   94.33    0.00   4.2E4  1.5E-3    4  749
blk_csun_6_band  111.63  130.86   94.12   94.33    0.00   4.2E4  1.5E-3    2  749
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:30:15)
 
 
        Call AMOEBA_MASK (17:30:17)
 
    AMOEBA_MASK done(17:31:22)
 
 
phi = -0.549 ,   sa = 0.039
        Amoeba output = avgphi=-0.54895585363071 avgsa=0.039399782168524
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:31:25)
 
        getcal("ceast","cwest","wlm","nbc",-0.549,0.039,3)
 
        GETCAL done (17:31:31)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 508.5103
 in radians [divided by -2061.3592 ] = -0.24668689474401
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.645622
 
 Should be close to pi, 2.645622  + -0.24668689474401  = 2.398935105256
 
        VCOR done (17:31:45)
    Now that wlm, et al, have been calc'd, wfits them with 200604221253.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060422125314...Calibrated
 
    Next dark id: CT060423125314
 
 
        CONV_CAL_FITS_PLUS (17:31:58)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060423125314	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060423125614	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060423125714	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060423125814	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060423125414	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060423125216	-> diff
 
        CONV_CAL_FITS_PLUS done(17:32:18)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:32:19)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.329 , yc = 522.623 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.249 , yc256 = 131.109 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:32:59)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.249
  y-coordinate of the center = 131.109
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.61  130.92   94.09   94.30    0.00   4.4E4  1.6E-3    4  746
blk_csun_6_band  111.61  130.92   94.09   94.30    0.00   4.4E4  1.6E-3    2  746
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:32:59)
 
 
        Call AMOEBA_MASK (17:33:01)
 
    AMOEBA_MASK done(17:34:16)
 
 
phi = -0.557 ,   sa = 0.036
        Amoeba output = avgphi=-0.55699296044574 avgsa=0.036171200460133
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:34:19)
 
        getcal("ceast","cwest","wlm","nbc",-0.557,0.036,3)
 
        GETCAL done (17:34:25)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 514.9913
 in radians [divided by -2061.3592 ] = -0.24983093679161
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.655197
 
 Should be close to pi, 2.655197  + -0.24983093679161  = 2.4053660632084
 
        VCOR done (17:34:38)
    Now that wlm, et al, have been calc'd, wfits them with 200604231253.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060423125314...Calibrated
 
    Next dark id: CT060424125414
 
 
        CONV_CAL_FITS_PLUS (17:34:52)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060424125414	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060424125714	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060424125814	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060424125914	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060424125514	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060424125316	-> diff
 
        CONV_CAL_FITS_PLUS done(17:35:14)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:35:14)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.337 , yc = 522.619 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.159 , yc256 = 131.106 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:35:56)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.159
  y-coordinate of the center = 131.106
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.63  130.90   94.07   94.27    0.00   4.2E4  1.5E-3    4  745
blk_csun_6_band  111.63  130.90   94.07   94.27    0.00   4.2E4  1.5E-3    2  745
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:35:57)
 
 
        Call AMOEBA_MASK (17:35:58)
 
    AMOEBA_MASK done(17:37:13)
 
 
phi = -0.557 ,   sa = 0.025
        Amoeba output = avgphi=-0.55664126974249 avgsa=0.025022085226859
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:37:16)
 
        getcal("ceast","cwest","wlm","nbc",-0.557,0.025,3)
 
        GETCAL done (17:37:22)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 516.7422
 in radians [divided by -2061.3592 ] = -0.25068032781477
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.657665
 
 Should be close to pi, 2.657665  + -0.25068032781477  = 2.4069846721852
 
        VCOR done (17:37:35)
    Now that wlm, et al, have been calc'd, wfits them with 200604241254.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060424125414...Calibrated
 
    Next dark id: CT060425125514
 
 
        CONV_CAL_FITS_PLUS (17:37:48)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060425125514	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060425125814	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060425125914	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060425130014	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060425125614	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060425125416	-> diff
 
        CONV_CAL_FITS_PLUS done(17:38:09)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:38:09)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.348 , yc = 522.615 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.212 , yc256 = 131.109 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:38:50)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.212
  y-coordinate of the center = 131.109
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.60  130.96   94.04   94.25    0.00   4.2E4  1.5E-3    4  744
blk_csun_6_band  111.60  130.96   94.04   94.25    0.00   4.2E4  1.5E-3    2  744
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:38:51)
 
 
        Call AMOEBA_MASK (17:38:52)
 
    AMOEBA_MASK done(17:39:56)
 
 
phi = -0.553 ,   sa = 0.02
        Amoeba output = avgphi=-0.5530129269417 avgsa=0.019567327803058
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:39:59)
 
        getcal("ceast","cwest","wlm","nbc",-0.553,0.02,3)
 
        GETCAL done (17:40:05)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 506.6861
 in radians [divided by -2061.3592 ] = -0.24580194465865
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.650739
 
 Should be close to pi, 2.650739  + -0.24580194465865  = 2.4049370553413
 
        VCOR done (17:40:18)
    Now that wlm, et al, have been calc'd, wfits them with 200604251255.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060425125514...Calibrated
 
    Next dark id: CT060426125614
 
 
        CONV_CAL_FITS_PLUS (17:40:31)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060426125614	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060426125914	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060426130014	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060426130114	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060426125714	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060426125516	-> diff
 
        CONV_CAL_FITS_PLUS done(17:40:52)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:40:52)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.348 , yc = 522.604 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.166 , yc256 = 131.099 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:41:36)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.166
  y-coordinate of the center = 131.099
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.58  130.99   94.04   94.22    0.00   4.3E4  1.6E-3    4  747
blk_csun_6_band  111.58  130.99   94.04   94.22    0.00   4.3E4  1.6E-3    2  747
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:41:37)
 
 
        Call AMOEBA_MASK (17:41:38)
 
    AMOEBA_MASK done(17:42:54)
 
 
phi = -0.55 ,   sa = 0.032
        Amoeba output = avgphi=-0.54999436701757 avgsa=0.032140057969112
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:42:57)
 
        getcal("ceast","cwest","wlm","nbc",-0.55,0.032,3)
 
        GETCAL done (17:43:03)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 505.7997
 in radians [divided by -2061.3592 ] = -0.24537193711799
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.659368
 
 Should be close to pi, 2.659368  + -0.24537193711799  = 2.413996062882
 
        VCOR done (17:43:16)
    Now that wlm, et al, have been calc'd, wfits them with 200604261256.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060426125614...Calibrated
 
    Next dark id: CT060427125714
 
 
        CONV_CAL_FITS_PLUS (17:43:28)
 
        Convert the following fits image...
 
        /ftr/armet/ftr/proccal/CT/dark/CT060427125714	-> dark
        /ftr/armet/ftr/proccal/CT/east.cal/CT060427130014	-> east
        /ftr/armet/ftr/proccal/CT/west.cal/CT060427130114	-> west
        /ftr/armet/ftr/proccal/CT/pinhole/CT060427130214	-> pin
        /ftr/armet/ftr/proccal/CT/unrotated/CT060427125814	-> sun
        /ftr/armet/ftr/proccal/CT/diff/CT060427125616	-> diff
 
        CONV_CAL_FITS_PLUS done(17:43:50)
    Now call solar_plus with:
       parps = /home/armet/iraf/initplus.params/CT.par
       ip    = 3
       docf  = /ftr/armet/ftr/field.tape.copies/CT.060420/calibration.log
       swp   = /home/armet/iraf/tape.cal.plus/
  
        SOLAR_PLUS (17:43:50)
 
        Call getavgdark
 
        Edit east and west headers
 
        threshold = 1700. , xc = 443.33 , yc = 522.603 , number of points = 64 
 
        threshold = 750. , 0xc256 = 111.149 , yc256 = 131.097 , number of points = 64 
 
 
 Call clampc for east, west, sun and diff to create
 ceast_6_band, cwest_6_band, csun_6_band, cdiff_6_band
 
        Call FNDLMB for block averaged 6 band csun (17:44:31)
 
 Standard output for fndlmb should follow:

######################
MODULE NAME = FNDLMB
Update the limb parameters in the header
Initial estimate for the limb ellipse parameters
  x-coordinate of the center = 111.149
  y-coordinate of the center = 131.097
  semi-minor axis length     = 92.5
  semi-major axis length     = 92.5
  ellipse rotation angle     = 0.
No. of parameters to fit = 4
Width of the search zone = 8.
Threshold for second derivative search = 0.
Percent, used to exclude outlying points = 1.
X-direction camera pixel length = 1.
Y-direction camera pixel length = 1.
No. of times the algorithm is invoked = 2

image           xcenter ycenter minorax majorax   angle  thresh     rms  no.  no.
name              (pix)   (pix)   (pix)   (pix)   (deg) (2nd d)   error iter zero
blk_csun_6_band  111.59  131.04   94.00   94.19    0.00   3.8E4  1.5E-3    4  747
blk_csun_6_band  111.59  131.04   94.00   94.19    0.00   3.8E4  1.5E-3    2  747
 
 majorax = 370.
 minorax = 370.
 (minorax+majorax/2) = 555.
 0.75 * 2. * 3.1415926 * (minorax+majorax/2)/4 = 653.843959875
 
 
    Minimum number of zero crossings, min_zero, = 653.843959875
 
        FNDLMB, block averaged csun,  done(17:44:31)
 
 
        Call AMOEBA_MASK (17:44:32)
 
    AMOEBA_MASK done(17:45:59)
 
 
phi = -0.555 ,   sa = 0.003
        Amoeba output = avgphi=-0.55529039414538 avgsa=0.0030423979482597
        Create flat image
 
  
        Call GETCAL in 1024 directory(17:46:01)
 
        getcal("ceast","cwest","wlm","nbc",-0.555,0.003,3)
 
        GETCAL done (17:46:08)
 
 
        Call VCOR to create vmi of csun, A band, ipermi: 3 , cosc: -2061.3592
 
 
        Call VCOR to create vmi of csun, B band, ipermi: 3 , cosc: -2061.3592
 
 
 
 Spatially averaged value, vmi[250:650,312:712,1], = 514.4436
 in radians [divided by -2061.3592 ] = -0.2495652383146
 
 Spatially averaged value, nbc[250:650,312:712,1], = 2.654145
 
 Should be close to pi, 2.654145  + -0.2495652383146  = 2.4045797616854
 
        VCOR done (17:46:20)
    Now that wlm, et al, have been calc'd, wfits them with 200604271257.fits extension.
    Now delete all the local cal images
 
    Current id:   CT060427125714...Calibrated
 
    master_gap_list is in /ftr/armet/ftr/proccal/CT/cals.20060420
    200604201251 added to master_gap_list
    200604211252 added to master_gap_list
    200604221253 added to master_gap_list
    200604231253 added to master_gap_list
    200604241254 added to master_gap_list
    200604251255 added to master_gap_list
    200604261256 added to master_gap_list
    200604271257 added to master_gap_list
 
    Call cleanup to remove all the files created
 
    Files and images deleted
 
    GETCALPLUS completed (05/09/06 17:46:30)