  filename = "/wrfout_d03.nc"
  f1       = addfile(filename,"r")

  q1 = f1->QCLOUD(73:77,0:8,:,:)
  time = f1->XTIME(73:77)
  dp1 = wrf_user_getvar(f1,"p",-1)
  dp  = dp1(73:77,0:8,:,:)
 
  lat2d1  = f1->XLAT(73:77,:,:)
  lon2d1  = f1->XLONG(73:77,:,:)
  z    = f1->ZNU(73:77,0:8)

  q1!0 = "time"
  q1!1 = "lev"
  q1!2 = "lat"
  q1!3 = "lon"
   
  printVarSummary(q1)
 
  dp!0 = "time"
  dp!1 = "lev"
  dp!2 = "lat"
  dp!3 = "lon"

  printVarSummary(dp)

  lwp1 = prcwater_dp (q1(time|:,lat|:,lon|:,lev|:), dp(time|:,lat|:,lon|:,lev|:))

  lwp = lwp1*1000
  ;lwp@long_name = "total column precipitable water"
  ;lwp@units     = "g/m2"

  copy_VarCoords(q1(:,0,:,:),lwp) ; add coordinate information

  printVarSummary(lwp)

  lat2d  = dim_avg_n_Wrap(lat2d1,0)
  lon2d  = dim_avg_n_Wrap(lon2d1,0)
 
  wks = gsn_open_wks("pdf","test01")

  dim  = dimsizes(lwp)
  nlat = dim(1)
  nlon = dim(2)

  res = True ; resource varb
  res@cnFillOn = True ; turn on color
  res@cnLinesOn = False ; no cn lines
  res@cnLineLabelsOn = False
  ;res@cnInfoLabelOn = False
  ;res@cnConstFEnableFill = True

res@cnFillPalette="GMT_seis"
;res@cnLevelSelectionMode = "Explicitlevels"

;res@cnLevels = (/0,10,20,30,40,50,60,70,80,90,100/)
;res@cnLevels = (/0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1/)


;res@mpGridAndLimbOn = True
res@mpLimitMode = "Corners"

res@mpLeftCornerLatF = lat2d(0,0)
res@mpLeftCornerLonF = lon2d(0,0)
res@mpRightCornerLatF = lat2d(nlat-1, nlon-1)
res@mpRightCornerLonF = lon2d(nlat-1, nlon-1)

res@lbLabelBarOn         = False
res@tmYRLabelsOn         = False              ; no right labels
res@tmYROn               = False              ; no right tickmarks
res@tmXTLabelsOn         = False              ; no top labels
res@tmXTOn               = False              ; no top tickmarksi
res@tmXBLabelFontHeightF = 0.025

res@tfDoNDCOverlay = True

res@pmTickMarkDisplayMode = "Always"
res@gsnMaximize           = True
res@tmXBMajorOutwardLengthF = 0

plot1 = new(5,graphic)
poly1 = new(5,graphic)

do i = 0,4

plot1(i) = gsn_csm_contour_map(wks, lwp(i,:,:), res)

;india state boundary
dir    ="gadm41_IND_shp/"
fname1 = "gadm41_IND_" + ispan(0,3,1) + ".shp"
lnres             = True
  lnres@gsLineColor = "black"
  poly1(i) = gsn_add_shapefile_polylines(wks,plot1(i),dir+fname1(1),lnres)

end do


  ;--------panel resousers

  pres                = True
  pres@gsnMaximize    = True
  pres@gsnFrame       = False
 ; pres@gsnPanelSave  = True   ; a secret resource telling NCL to leave the paneled plots in their resized f
  pres@gsnPanelFigureStrings = (/"6.30am","7.30am","8.30am" ,"9.30am","10.30am"/)
  pres@gsnPanelFigureStringsFontHeightF = 0.008
  pres@gsnPanelFigureStringsPerimOn     = False
  pres@amJust   = "TopRight"
  pres@gsnPanelFigureStringsBackgroundFillColor = 0
  pres@gsnPanelLabelBar = True                ; commom label bar
  pres@pmLabelBarOrthogonalPosF = -0.005
  pres@lbTitleString         = "LWP (g/m2)"
  pres@lbTitlePosition       = "Right"
  pres@lbTitleOffsetF        = -0.004
  pres@lbTitleFontHeightF    = 0.010
  pres@lbTitleDirection      = "across"
  pres@lbLabelFontHeightF    =  0.011  ;label bar font size
  pres@lbBoxSeparatorLinesOn = False
  gsn_panel(wks,plot1,(/1,5/),pres)


