; 
; Flight track overviews
;
;*************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
;**************************************************
 
begin
;
; constants
;
  deg2rad = asin(1.)/180.
;
; specify data file:
;
  wks = gsn_open_wks("pdf","circle") 
  gsn_define_colormap (wks,"gui_default")
  stem = "/Volumes/Data/RICO/C130/LRT/rico_rf"

  flight_nmbr = "09"

  xyres = True
  xyres@gsnDraw  = False
  xyres@gsnFrame = False
  xyres@tmYLLabelStride  = 2 
  xyres@xyLineThicknessF = 1.
  xyres@vpHeightF = 0.15

;  xyres@tmXBBorderOn = False
  xyres@tmXBLabelsOn = False
  xyres@tmXTBorderOn = False
  xyres@tmXTOn = False
  xyres@tmYRBorderOn = False
  xyres@tmYROn = False
  xyres@tiXAxisFontHeightF = 0.010
  xyres@tiYAxisFontHeightF = 0.010
  xyres@tmXBLabelFontHeightF = 0.008
  xyres@tmYLLabelFontHeightF = 0.008

  flight_nmbrs = (/"01","02","03","04","05","06","07","08","09","10",\
  "11","12","13","14","15","16","17","18","19"/)
  rate = "LRT"
  plot = new((/6/),"graphic")

  do k=3,3;18
    flight_nmbr = flight_nmbrs(k)
    load "../Data/flights.mdat"
    ii = stringtointeger(flight_nmbr)-1
    if (.not.ismissing(iFT1(ii))) then
      l1 = iFT1(ii)-100
      l2 = min( (/dimsizes(flight->time_offset)-1,l1 + 2000/) )

      x1 = 0.01*(l1 + 100)
      xyres@gsnXRefLine = (/x1,x1+17.99/)
      xyres@gsnXRefLineColor = "red"
      xyres@gsnXRefLineDashPattern = 2
      xyres@gsnXRefLineThicknessF = 2.

      time = flight->time_offset(l1:l2)/100.
  
      galt = flight->GGALT(l1:l2)
      galt@long_name = "GPS Altitude [m]"
      galt@_FillValue = -32767
    
      roll = flight->ROLL(l1:l2)
      roll@long_name = "Roll Angle [deg]"
      roll@_FillValue = -32767

      ptch = flight->PITCH(l1:l2)
      ptch@long_name = "Pitch Angle [deg]"
      ptch@_FillValue = -32767

      thdg = flight->THDG(l1:l2)
      thdg@long_name = "Heading [deg]"
      thdg@_FillValue = -32767

      xyres@trYMaxF =  avg(galt) + 100.
      xyres@trYMinF = xyres@trYMaxF - 200.
      plot(0) = gsn_csm_xy(wks,time,galt,xyres)
      xyres@trYMaxF =  10.
      xyres@trYMinF = -10.
      plot(1) = gsn_csm_xy(wks,time,roll,xyres)
      plot(2) = gsn_csm_xy(wks,time,ptch,xyres)

      xyres@trYMaxF =   370.
      xyres@trYMinF =   -10.
      xyres@tmXBLabelsOn = True
      xyres@tiXAxisString = "Time since start of tape [hs]"
      plot(3) = gsn_csm_xy(wks,time,thdg,xyres)
	
      pres = True
      pres@gsnPanelBottom = 0.1
      pres@gsnPanelTop = 0.9
      pres@gsnFrame = False
      gsn_panel(wks,plot(0:3),(/4,1/),pres)
      frame(wks)
      if (.not. ismissing(iSF1(ii))) then
        print ("SF1  "+avg(flight->GGALT(iSF1(ii):iSF1(ii)+1799)))
      end if
      if (.not. ismissing(iSF2(ii))) then
        print (" "+avg(flight->GGALT(iSF2(ii):iSF2(ii)+1799)))
      end if
      if (.not. ismissing(iSC1(ii))) then
        print (" "+avg(flight->GGALT(iSC1(ii):iSC1(ii)+1799)))
      end if
      if (.not. ismissing(iSC2(ii))) then
        print (" "+avg(flight->GGALT(iSC2(ii):iSC2(ii)+1799)))
      end if
      if (.not. ismissing(iFT1(ii))) then
        print (" "+avg(flight->GGALT(iFT1(ii):iFT1(ii)+1799)))
        print ("SF1  "+avg(flight->GGSPD(iFT1(ii):iFT1(ii)+1799)))
      end if
      if (.not. ismissing(iFT2(ii))) then
        print (" "+avg(flight->GGALT(iFT2(ii):iFT2(ii)+1799)))
      end if
    end if
  end do
  
end