Martha starts 060615

** There were two drifts for this night (planned)
**    22p2  +052112   1034 to 1316
**    16p2  +035336   1319 to 1631
** So the calibration will have to be split.
** This is a spring drift....
** The pos file names are therefore:    pos_060516_+052112s.sav
**                                      pos_060516_+035336s.sav
** Therefore, the calibration will be split.

-- Working in:  /home/adjuntas2/haynes/idl/work/06.05.16
   cp /home/dorado13/galaxy/idlraw/06.05.16/*list* .
   gedit cal060516_1.list   [make triplets]
   gedit drift060516_1.list [add directory to name]
   and likewise for the second set

   >idl
   >@wasinit2
   >@alfinit

-- Perform the calibration, in 2 steps:
   > callist='cal060516_1.list'               [assuming you already editted it...]
   > dir='/home/dorado13/galaxy/idlraw/06.05.16'     [where the raw IDL files are]
   > calib1,callist,dcalON,dcalOFF,dir=dir
   > help,dcalon,/st
   > calib2,dcalON,dcalOFF,ncalib
   > Any bad boards in system?
   n
** There are 16 cal values for this run ***
** The calibration for this run looks ok **
   > save, ncalib,file='ncalib_060516_1.sav'
   > help,ncalib,/st

-- Repeat for the second set:
   > callist='cal060516_2.list'               [assuming you already editted it...]
   > dir='/home/dorado13/galaxy/idlraw/06.05.16'     [where the raw IDL files are]
   > calib1,callist,dcalON,dcalOFF,dir=dir
   > help,dcalon,/st
   > calib2,dcalON,dcalOFF,ncalib
   > Any bad boards in system?
   n
** There are 19 cal values for this run ***
** The calibration for this run looks ok **
   > save, ncalib,file='ncalib_060516_2.sav'
   > help,ncalib,/st

-- Process the first set
   > bpdgui
   Enter the new calsession: calsession_060516_1.sav
   Enter the new runpos: pos_060516_+052112s.sav                  
   CRAT set to "fit"
**The interp regions here are [3476 to 3516] **

-- restore,'pos_060516_+052112s.sav'
   help,pos
** There are 16 drifts for this run **

-- Next, process the second set
   > bpdgui
   Enter the new calsession: calsession_060516_2.sav
   Enter the new runpos: pos_060516_+035336s.sav                  
   CRAT set to "fit"
**The interp regions here are [3480 to 3520] **

-- restore,'pos_060516_+035336s.sav'
   help,pos
** There are 19 drifts for this run **

-- restore,'pos_060516_+052112s.sav'
-- restore,'d214400+294831.526180333.sav'
-- flagbb,dred,cont_pt,mask,pos,GAUSAV=11,HAN=3,AGC=1,n1=0,n2=6
   save,pos,file='pos_060516_+052112s.sav'

--------------------------------------------------------------------------------
restore,'pos_060516_+052112s.sav'
  d103431+052112.613669822.sav  q1  1 AGC
  d104437+052112.613670427.sav  q1  1 AGC
  d105444+052112.613671032.sav  q1  2 AGC   GPS burst ~50+15 secs rec 38-219; fnce 1.0
  d110450+052112.613671637.sav  q1  3 AGC   fnce 1.0
  d111457+052112.613672242.sav  q1  1 AGC   ditto
  d112504+052112.613672847.sav  q1  1 AGC   ditto   
  d113510+052112.613673452.sav  q1  2 AGC   ditto  
  d114517+052112.613674057.sav  q1  2 AGC   ditto 
  d115524+052112.613674662.sav  q1  2 AGC   ditto
  d120530+052112.613675267.sav  q1  3 AGC   ditto
  d121537+052112.613675872.sav  q1  5 AGC   ditto
  d122544+052112.613676477.sav  q1  2 AGC   ditto
  d123550+052112.613677082.sav  q1  1 AGC + posvel HVC    cont 3,5; 1.0 looks ok
  d124557+052112.613677687.sav  q1  2 AGC                 GPS recs 360-490
  d125604+052112.613678292.sav  q1  2 AGC
  d130610+052113.613678897.sav  q1  2 AGC
save,pos,file='pos_060516_+052112s.sav'
Above; Beam 1.0 shows fence mod; setting status to 2 or 3, middle drifts (see notes)

flagbb,dred,cont_pt,mask,pos,GAUSAV=11,HAN=3,AGC=1,n1=0,n2=6
restore,'pos_060516_+035336s.sav'

  d131901+035336.613679666.sav  q1  1 AGC
  d132908+035336.613680271.sav  q1  1 AGC  GPS recs 476-534
  d133915+035336.613680876.sav  q1  3 AGC  sw 4.1
  d134921+035336.613681481.sav  q1  1 AGC  sw 4.1; 5.1
  d135928+035336.613682086.sav  q1  1 AGC  rfi5 really strong 2,1       
  d140934+035336.613682691.sav  q1  1 AGC  
  d141941+035336.613683296.sav  q1  3 AGC  
  d142948+035337.613683901.sav  q1  3 AGC 
  d143954+035337.613684506.sav  q1  1 AGC 
  d145001+035337.613685111.sav  q2  nada   strong GPS burst rec 338-582
  d150008+035337.613685716.sav  q1  nada 
  d151014+035337.613686321.sav  q1  1 AGC
  d152021+035336.613700526.sav  q1  1 AGC
  d153028+035337.613701131.sav  q1  1 AGC
  d154034+035337.613701736.sav  q1  2 AGC
  d155041+035336.613702341.sav  q1  1 AGC
  d160048+035337.613702946.sav  q1  1 AGC
  d161054+035337.613703551.sav  q1  1 AGC
  d162101+035337.613704156.sav  q1  1 AGC

save,pos,file='pos_060516_+035336s.sav'
 
** Be sure to save your pos file regularly!
   
 strip_pv,dred,cont_pt,msmooth,n1=0,n2=6,showpol=3,gausav=11,han=5,rfimod=0,units=1
 inspect,msmooth,showpol=3,nstrip=4,agc=1

 reviewbb,dred,cont_pt,mask,pos,GAUSAV=11,HAN=3,AGC=1,n1=0,n2=6

> print,where(pos.badbox[0,0,0,2] eq 0)

** Overall this is a good dataset, with good baselines and not much rfi **