load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

begin

   ;## day when integration started
  day	= "07"
  
  
  ;## hour when integration started
 hour	= "12"


  ;## open work station
  timedate = systemfunc("date +$Y-%m-%d-%H-%M")
  wks = gsn_open_wks("x11","avo_"+timedate)         
  
    
  ;## path to wpp outfiles
  path	= "/import/wrkdir1/asemenov/WRF/WRFV3/10_25.03_res_30_10km/postprd2/"
  

  ;## file name
  files	= systemfunc("ls "+path+"WRFPRS_d02.*")

  do j 	= 0, dimsizes(files)-1
  fname	= files(j)
  

  ;## open file
  print("Opening file: "+fname+".nc")
;  f     = addfile ("/import/wrkdir1/asemenov/WRF_Alex/WRF/run-number-01/wrfout_d01_2009-10-10_12:00:00"+".nc", "r")
    ;## open file
  print("Opening file: "+fname+".grb")
  f     = addfile (fname+".grb", "r")
 
 
  gsn_define_colormap(wks,"BlueWhiteOrangeRed")
 
 
 
  res                       			= True 
  res2						= True
  
  ;TICKMARK RESOURCES (ON/OFF - FONT SIZE)
	res@pmTickMarkDisplayMode = "Always" 
	res@tmXTOn				= False
	res@tmXBOn				= True
	res@tmYLOn				= False
	res@tmYROn				= False
	res@tmXTLabelFontHeightF		= 0.025
	res@tmXBLabelFontHeightF		= 0.025
	res@tmYLLabelFontHeightF		= 0.025
	res@tmYRLabelFontHeightF		= 0.025
	
	res2@tmXTOn				= False
	res2@tmXBOn				= False
	res2@tmYLOn				= False
	res2@tmYROn				= False
	
  ; MAP RESOURCES.
	res@mpGridAndLimbOn			= True
	res@mpGridLineDashPattern		= 2
	res@mpGeophysicalLineColor		= "Black"
	res@mpGeophysicalLineThicknessF		= 1.5
	res@mpNationalLineColor			= "Black"	
	res@mpNationalLineThicknessF		= 1.5
	res@mpUSStateLineColor			= "Black"
	res@mpUSStateLineThicknessF		= 1.5
	res@mpPerimLineColor			= "Black"
            
  ;LABEL BAR RESOURCES
	res@lbLabelAutoStride			= True            
	res@lbOrientation			= "Vertical"
	res@lbLabelFontHeightF			= 0.025
	
  ;CONTOUR RESOURCES
	res@cnFillOn				= True             
	res@cnLinesOn				= False            
	res@cnLineLabelsOn			= False  
	res@cnLevelSelectionMode   		= "ManualLevels"
	res@cnMinLevelValF         		= 30.
	res@cnMaxLevelValF         		= -50.
	res@cnLevelSpacingF			= 5.
          
	res2@cnFillOn				= False             
	res2@cnLinesOn				= True            
	res2@cnLineLabelsOn			= True
	res2@cnLevelSelectionMode   		= "ManualLevels"
	res2@cnMinLevelValF         		= 940.
	res2@cnMaxLevelValF         		= 1044.
	res2@cnLevelSpacingF			= 4.
	res2@cnLineThicknessF			= 1.
	res2@cnInfoLabelOn			= True	
res2@cnInfoLabelOn = False
;res2@cnHighLabelsOn = True
 
   ; res2@cnLowLabelsOn = True
   ; res2@cnLowLabelBackgroundColor=0

   ;  res2@cnHighLabelFontHeightF = 0.02
   ;  res2@cnLowLabelFontHeightF =0.04

;res2@cnLowLabelBackgroundColor=-1

	
  ;GSN RESOURCES
	res@gsnSpreadColors			= True               
	res@gsnLeftString			= " "
	res@gsnRightString			= " "
	res@gsnDraw				= False
	res@gsnFrame				= False
;	res@gsnMaximize				= True
	
	res2@gsnDraw				= False
	res2@gsnFrame				= False
	res2@gsnMaximize			= True
	res2@gsnLeftString			= " "
	
  ;TITLE RESOURCES
	res@tiMainOn				= False 
	res@tiMainPosition			= "Center"
	res@tiMainFontHeightF			= 0.025
 
  ;PROJECTION RESOURCES
	res@tfDoNDCOverlay 			= True 
	res@gsnAddCyclic 			= False
 	res@mpProjection			= "Stereographic"
	res@mpLimitMode				= "Corners"
        res@mpLeftCornerLatF                  = 65.2
        res@mpLeftCornerLonF                  = 132.8014
        res@mpRightCornerLatF                 = 63.8
        res@mpRightCornerLonF                 = -45.0
  	res@mpCenterLatF 			= 90
  	res@mpCenterLonF 			= -180
	
	res2@tfDoNDCOverlay 			= True 
	res2@gsnAddCyclic 			= False
 	res2@mpProjection			= "Stereographic"
	res2@mpLimitMode				= "Corners"
        res2@mpLeftCornerLatF                  = 65.2
        res2@mpLeftCornerLonF                  = 132.8014
        res2@mpRightCornerLatF                 = 63.8
        res2@mpRightCornerLonF                 = -45.0
  	res2@mpCenterLatF 			= 90
  	res2@mpCenterLonF 			= -180


  var			= f->TMP_GDS5_HTGL_10
  var@long_name		= "temperature at 2m"
  printVarSummary	(var)
  var = (/ var-273 /) 


 var2                  = f->PRMSL_GDS5_MSL_10
 var2 = (/ var2/100 /) 

  res@gsnLeftString	= "Init time: "+var@initial_time+"~C~Fcst time: "+var@forecast_time+" hrs ~C~"+var@long_name+" ("+var@units+")"              

 plot1                         = gsn_csm_contour_map(wks,var,res)
  plot2                         = gsn_csm_contour(wks,var2,res2)
  overlay(plot1,plot2)
  draw(plot1)
  frame(wks)

end do

end