I am not sure what I am doing wrong. I am trying to plot one variable
(sea-level pressure) from a gridded GFS file in .nc format.
2 questions:
Do I really need to specify a sub-region tp plot the contours?
The data only pertains to a region; the file itself only has data from
Lon35E, 35W: Lat:35 N to -25 S.
Why won't my script plot my conour intervals?
The contour plot I am resulting does not draw all of the
contourintervals that I want it to
Any help would be greatly appreciated.
Sincerely,
Erik
a = addfile("met_em.d01.2006-08-07_00:00:00.nc","r")
;************************************************
;************************************************
;do it = 0,ntimes-1,8 ; TIME LOOP
; print("Working on time: " + times(it) )
; print("Working on file: " + numFILES(ifil) )
slp = a->PMSL(0,:,:)
; slp = slp/100 ; pa to mb
;************************************************
; create default plot
;************************************************
wks = gsn_open_wks("x11","ce") ; open a ps file
; Set some Basic Plot options
inputres = True
inputres_at_MainTitle = "FNL Boundary Condition"
mpres= True
mpres= inputres
mpres_at_mpMinLonF = -35. ; select a subregion
mpres_at_mpMaxLonF = 35
mpres_at_mpMinLatF = -25.
mpres_at_mpMaxLatF = 35
mpres_at_cnLineColor = "NavyBlue"
mpres_at_cnHighLabelsOn = True
mpres_at_cnLowLabelsOn = True
mpres_at_cnLevelSelectionMode = "ManualLevels" ; set the
contour levels
mpres_at_cnMinLevelValF = 90000. ; set the
minimum contour level
mpres_at_cnMaxLevelValF = 104000. ; set the
maximum contour level
mpres_at_cnLevelSpacingF = 400. ; set the
interval between contours
mpres_at_gsnContourLineThicknessesScale = 2.0
plot = gsn_csm_contour_map_ce(wks,slp,mpres) ; create a default plot
end
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 23 2007 - 11:09:02 MDT
This archive was generated by hypermail 2.2.0 : Mon Aug 27 2007 - 14:59:53 MDT