begin

;::::::::::::::::::::::::::::::::::
fils=systemfunc("ls ./xMOD06_L2*.hdf")
f=addfiles(fils,"r")
;::::::::::::::::::::::::::::::::::

ListSetType(f,"cat")
cl = f[:]->Cloud_Top_Pressure

;::::::::::::::::::::::::::::::::::

printVarSummary(cl)

cl@lat2d = f[:]->Latitude
cl@lon2d = f[:]->Longitude 

wks=gsn_open_wks("png","1_modis_ctp_349")          
;:::::::::::::::::::::::::::::::::::::::::::

res=True						; plot mods desired
res@gsnMaximize		=True				; make plot large
res@cnLinesOn		=False				; turn off contour lines
;res@cnLineLabelsOn	=True				; turn off contour line labels
res@gsnSpreadColors	=True 				; use the entire color spectrum
res@cnFillOn		=True				; enable contour fill
res@gsnMaximize		=True				; make plot large
res@mpGeophysicalLineThicknessF= 2. 			; double the thickness of geophysical boundaries
res@cnMissingValFillPattern = 0

;::::::ZOOM:::::::::::;
res@mpMinLatF = 20
res@mpMaxLatF = 32
res@mpMinLonF = 68
res@mpMaxLonF = 95

;:::::::::::::Title&Plotting::::::::::::::::
  res@tiMainString="Cloud Top Pressure"			; create title
  plot=gsn_csm_contour_map_ce(wks,cl,res) 		; plot on world map first
  delete(res@tiMainString) 				; delete title for subsequent plots

;:::::::::::::Clean up resources used:::::::::::::::::::::
  delete(plot) 
  delete(wks)
  delete(res)
end