;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "WRFOptions.ncl"     ; set basic plot options here
;load "gsn_code.ncl"
;load "WRFPlot.ncl"
;load "WRFUserARW.ncl"
;load "SkewTFunc.ncl"
external MY_SORT "./my_sort.so"
 begin
  DATADir = "/data1/wrf/01_23_2016/GEFS"
   cmd = "find " + DATADir + "/GEFS* -type d -print"
   dirWRF = systemfunc (cmd)
   print(dirWRF)
  n_dirWRF = dimsizes(dirWRF)
  print("n_dirWRa" + n_dirWRF)



   filename = "/data1/wrf/GEFS/GEFS_01/wrfout_d02_2015-01-26_00:00:00"
   a = addfile(filename,"r")
;   mpdatabaseversion= "Earth..4"

; Set some basic resources
  res = True
  res@MainTitle                   = "REAL-TIME WRF"
;

; find ntimes
  times  = wrf_user_list_times(a)  ; get times in the file
  ntimes = dimsizes(times)         ; number of times in the file
  print ("ntimes = " + ntimes)
  ter = wrf_user_getvar(a,"HGT",0)
  ter = ter - 10
   printMinMax(ter,0)
  dims=dimsizes(ter)
  print("dims(0) " + dims(0))
  print("dims(1) " + dims(1))
  n_loops = 9
; wtype = "png"
  wtype = "pdf"
; for png --> gives "error" for pdf plot
  wtype@wkWidth  = 2500    ; produces a better looking PNG
  wtype@wkHeight = 2500
; wks = gsn_open_wks(wtype,"wrfrefl_01_26_2015_1.3km")
  wks = gsn_open_wks(wtype,"wrfrefl_01_23_2016_12.0km_pm")
   its=24
   ite=24
  do it = its,ite,1
   i_cnt = 0
   do n=0,n_dirWRF- 1,20      ; loop over directories
   i_loop = n*20/20
  print("i_loop" + i_loop)
  print("dirWRF(n) = " + dirWRF(n))
   filename = addfiles(dirWRF(n)+"/"+"wrfout_d02_2016-01-23_00:00:00", "r")
  print("filename = " + filename)

;  filename = "/data1/wrf/GEFS/GEFS_00/wrfout_d02_2015-01-26_00:00:00"
   print("filename = " + filename)
   a = addfile(filename,"r")
    time=it
    print("Working on time: " + times(it) )
   end do  ;       end of the time loop
  end do  ;       end of the model loop
end