Re: plot of contours not doing what I tell it to

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Thu, 23 Aug 2007 11:35:39 -0600

Hi Erik,

I just read in a global slp data file I had here, used your exact
resource list, and am getting contours.

When you say that the data goes from 35E, 35W, do you mean that there is
data from 35W eastwards towards 35E? Or, is there data from 35E->325E?

Eliminate your mpMin/mpMax resources and plot the entire globe, and see
if you get any contours. If not, make sure your data has values that
fall within your specified min/max contour values...

Good luck,
Adam

Erik Noble wrote:
> 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

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 23 2007 - 11:35:39 MDT

This archive was generated by hypermail 2.2.0 : Mon Aug 27 2007 - 14:59:53 MDT