********************************CIFEX03***************************************
****************************** April 2004 ************************************
******************************************************************************
Note: Slight adjustment in the mean values can be done to get norm=1.000000

;       BEAM-ANGLE (unit vectors in AC coord. system) CALIBRATION
; (statistics from surface returns during AC/WCR calibration maneuvers)

                25 Mar 04                HiCu03
Side beam     mean       std         mean       std 
 AC_x     -0.012633034           -0.0122113   0.00158659
 AC_y      0.99991505             0.999898    0.00005636
-AC_z      0.0032026578           0.00317781  0.00624242

Nadir
 AC_x     -0.054603764           -0.0554711   0.00101259  
 AC_y      0.0018075706           0.0019240   0.0068228 
-AC_z      0.99850614             0.998436    0.0000750 

******************************************************************************
******************************************************************************
; AC data file

acncf='/net/owl/data/kingair_data/cifex04/processed/20040325_2.c25.nc'

;############################################################################
; Side/up antenna beam calibration
;############################################################################

openwcrf,'/sova4/cifex04/mar25/W2004-03-25-19-35-55'

.run commonwcrset
readwcrdata

  ;--------------------------------------------------------------------------
  ;side beam
  ;--------------------------------------------------------------------------

  t=9.00
  report='calsidebeam'+num2str(wcrdate)+'.txt'
  beam=0 & wcrbeamfind,beam,5,status=status,acncf=acncf,ground=-30.,timeoff=t,$
                       /nowcr,/disp,report=report,epstop=1.e-4

  print,beam,t

  ;    forward      right wing    down       timeoffset

sidebeam=[$

  [-0.012633034, 0.99991505, 0.0032026578, 0.000],$ ;<- 3/25/2004
  [-0.012320710, 0.99992380,-0.00077003439      ],$ ;<- 3/25/2004,toff=-0.12
]

  wcrtimeoff,toff,tcov,beam,5,ppsel='HH',acncf=acncf,/nowcr,ground=-30.,xplot=0

  print,toff,tcov  

  ; timeoffset, mean Doppler, std Doppler from ground in m/s

sidetimedv=[$

  [-0.120000,    0.0158514,     0.167368  ],$ ; <- 3/25/2005
  [-0.120000,    0.0242464      0.168008  ],$ ; <- 3/25/2005, toff=-0.12
]

;############################################################################
; Down antenna beam calibration
;############################################################################

.run commonwcrset

openwcrf,'/sova4/cifex04/mar25/W2004-03-25-19-35-55'
readwcrdata

acncf='/net/owl/data/kingair_data/cifex04/processed/20040413.c25.nc'
openwcrf,'/sova4/cifex04/apr13/'
readwcrdata,193100,193300,/timein

  ;--------------------------------------------------------------------------
  ;down beam
  ;--------------------------------------------------------------------------

  t=9.00
  report='caldownbeam'+num2str(wcrdate)+'.txt'
  beam=3 & wcrbeamfind,beam,status=status,acncf=acncf,ground=-30.,timeoff=t,$
                       /nowcr,/disp,report=report

  print,beam,t
 
  ;  forward        right wing    down      timeoffset

nadirbeam=[$

  [-0.054603764, 0.0018075706,  0.99850614, 0.00],$ ;<- 3/25/2004
  [-0.054563419, 0.0000000,     0.99851064      ],$ ;<- 3/25/2004, toff=-0.12
  [-0.055213629, 0.013803381,   0.99837958      ],$ ;<- 4/13/2004
]

  wcrtimeoff,toff,tcov,beam,ppsel='VV',acncf=acncf,/nowcr,ground=-30.,xplot=0

  print,toff,tcov  

  ; timeoffset, mean Doppler, std Doppler from ground in m/s

nadirtimedv=[$

  [-0.120000,  0.000906878,     0.192339],$ ;<- 3/25/2004
  [-0.120000,  0.000208201,     0.192730],$ ;<- 3/25/2004, toff=-0.12
  [ 0.0400000 -0.00668451       0.422000],$ ;<- 4/13/2004
]


;===========================================================================
;Side beam -------------------------------------------------

 print,[mean(sidebeam[0,*]),std(sidebeam[0,*])],$
       [mean(sidebeam[1,*]),std(sidebeam[1,*])],$
       [mean(sidebeam[2,*]),std(sidebeam[2,*])]

                 All flights     
 AC_x   -0.0122113    0.00158659 
 AC_y    0.999898     0.00005636 
-AC_z    0.00317781   0.00624242 
   
print,[mean(sidetimedv[0,*]),std(sidetimedv[0,*])],$
       [mean(sidetimedv[1,*]),std(sidetimedv[1,*])],$
       [mean(sidetimedv[2,*]),std(sidetimedv[2,*])]

timeoff [s]    
DV_grnd [m/s]  
DV_stdev[m/s]  

; Nadir beam -------------------------------------------------

 print,[mean(nadirbeam[0,*]),std(nadirbeam[0,*])],$
       [mean(nadirbeam[1,*]),std(nadirbeam[1,*])],$
       [mean(nadirbeam[2,*]),std(nadirbeam[2,*])]

                   All flights 
 AC_x   -0.0554711  0.00101259 
 AC_y    0.0019240  0.0068228  
-AC_z    0.998436   0.0000750  

 print,[mean(nadirtimedv[0,*]),std(nadirtimedv[0,*])],$
       [mean(nadirtimedv[1,*]),std(nadirtimedv[1,*])],$
       [mean(nadirtimedv[2,*]),std(nadirtimedv[2,*])]

timeoff [s]   
DV_grnd [m/s] 
DV_stdev[m/s]