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/csm/contributed.ncl" 
begin 
fi = addfile("./20070708_06.nc", "r") 

wks = gsn_open_wks("eps", "20070708_06") 


lat = fi->lat
lon = fi->lon

rain  = fi->rain

printVarSummary(rain)

rain@lat2d = lat
rain@lon2d = lon

res = True
res@gsnAddCyclic = False 

;;;;;;;;;;;;;;set for the map
res@mpDataBaseVersion       = "HighRes"
res@mpDataSetName           = "Earth..4"
res@mpGridAndLimbOn = False                         ; Turn off lat/lon lines
res@mpOutlineOn           = True         ; Turn on map outlines
res@mpGeophysicalLineColor = "Black"      ; Change the outline line color
res@mpGeophysicalLineThicknessF= 2.      ; double the thickness of geophysical boundaries
res@mpOutlineSpecifiers    = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/) 
;;;;;;;;;;;;;;;ZOOM
res@mpMinLatF = 30 
res@mpMaxLatF = 34 
res@mpMinLonF = 114 
res@mpMaxLonF = 122 

res@pmTickMarkDisplayMode = "Always"    

; contours
res@cnFillOn = True 
res@cnLinesOn = False           ; Turn off contour lines 
res@cnLevelSelectionMode = "ExplicitLevels" 
res@cnLevels = (/1, 2, 6, 10, 15, 20, 30, 40, 50, 70, 90, 110, 130, 150, 200/)
; label bar
res@lbBoxLinesOn = False ; Turn off label bar box lines 
plot = gsn_csm_contour_map(wks, rain, res)


end