data=addfile("Phailin_surface.grib","r")
f=data->MSL_GDS0_SFC(32,:,:)
time=data->initial_time0_hours
lat =data->g0_lat_1
lon =data->g0_lon_2
minlat = 0 
maxlat = 40
minlon = 60
maxlon = 110
f_sub = f({minlat:maxlat},{minlon:maxlon})  ;subsetting from slp and put the values in slp_sub
print(time)
printVarSummary(f)
printVarSummary(f_sub) 
printMinMax(f_sub,0)              
wks =gsn_open_wks("png","sst")
res =True
res@gsnDraw =False
res@gsnFrame =False
res@cnFillOn =True
res@cnFillPalette ="WhBlGrYeRe" 
res@mpMinLatF =0
res@mpMaxLatF =40
res@mpMinLonF =60
res@mpMaxLonF =110
plot	=gsn_csm_contour_map(wks,f,res)
draw(plot)
frame(wks)