;
;
;  This Script makes spatial plot of nearSurface Rain observed from radar 
;
;
    load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
    load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin

;  1. Read file and Print Meta Info

 FILES       = systemfunc("ls ./DATAF/2A25*.HDF")
                   print(FILES)
                   print(FILES(0))
 F           = addfiles(FILES,"r")
                   printVarSummary(F)
;               print(F)
;               print(F[1])
 nf          = dimsizes(FILES)
                   print(nf)

;;;;   The Dimesions of the above files is 9247X 49                           

;
; longitude       = F[0]->Longitude
; latitude        = F[0]->Latitude
; longitude@units = "degrees_east"
; latitude@units  = "degrees_north"
;                    printVarSummary(longitude)
;                    printVarSummary(latitude)


  Wks_Rain      =   gsn_open_wks("png","NearSurfRain")
                    gsn_define_colormap(Wks_Rain,"default")

  PLOT          = new(4,graphic)

     res                      =    True
     res@mpMinLatF            =   5.00
     res@mpMaxLatF            =   20.0
     res@mpMinLonF            =   72.0
     res@mpMaxLonF            =   85.0
     res@mpFillOn             =   False

     res@gsnDraw              =   False
     res@gsnFrame             =   False

     res@cnFillOn             =   True 
     res@cnLevelSelectionMode =   "ExplicitLevels"
     res@cnLevels             =   ispan(5,40,5)    
     NLevels                  =   dimsizes(res@cnLevels)
;;  
     res@tmXTOn               =   False
     res@tmYROn               =   False

;
     do   ii     =   0, nf-1
     For_Date    =   str_split(FILES(ii),".")
                     print(For_Date)
     longitude   =   F[ii]->Longitude
     latitude    =   F[ii]->Latitude
     Surf_Rain   =   F[ii]->nearSurfRain
     Surf_Rain@_FillValue   =   -99                     

          print(";;;; of Values< "+ res@cnLevels(0) + " = " +Surf_Rain.lt.res@cnLevels(0))
;                     asciiwrite("abc.txt",Surf_Rain)
;                       printVarSummary(Surf_Rain)
;                       printVarSummary(longitude)
;                       printVarSummary(latitude)

     Surf_Rain!0    =  "longitude"
     Surf_Rain!1    =  "latitude"
     Surf_Rain@lon2d =  longitude              
     Surf_Rain@lat2d =  latitude                       
     longitude@units = "degrees_east"
     latitude@units  = "degrees_north"
;                        printVarSummary(Surf_Rain)
     res@gsnRightString = For_Date(1)+"."+For_Date(2)
     PLOT(ii)        =    gsn_csm_contour_map(Wks_Rain,Surf_Rain,res) ; keep Long in X axis
     end do
     delete(Surf_Rain)
     delete(latitude)
     delete(longitude)
     delete(For_Date)
     
;;;;   Problem with this data file
;;;;   The Dimesions of this file (F20110410) is different 
;;;;   from the 2 above 9248x49
;;;;   The Dimesions of the above files is 9247X 49                           

     F_20110410  =  "2A25.20110410.76333.7.HDF"
     F20110410   =  addfile(F_20110410,"r")
     For_Date    =   str_split(F_20110410,".")
     longitude   =   F20110410->Longitude
     latitude    =   F20110410->Latitude
     Surf_Rain   =   F20110410->nearSurfRain
     Surf_Rain@_FillValue   =   -99                     
;                       printVarSummary(Surf_Rain)
;                       printVarSummary(longitude)
;                       printVarSummary(latitude)

     Surf_Rain!0    =  "longitude"
     Surf_Rain!1    =  "latitude"
     Surf_Rain@lon2d = longitude              
     Surf_Rain@lat2d = latitude                       
     longitude@units = "degrees_east"
     latitude@units  = "degrees_north"
;                        printVarSummary(Surf_Rain)
     res@gsnRightString = For_Date(1)+"."+For_Date(2)
     PLOT(3)        =    gsn_csm_contour_map(Wks_Rain,Surf_Rain,res) ; keep Long in X axis
                     
;   PANEL PLOTS
    Panel_Res    = True
                   gsn_panel(Wks_Rain,PLOT,(/1,4/),Panel_Res) ; keep Long in X axis
end