plot of contours not doing what I tell it to

From: Erik Noble <enoble_at_nyahnyahspammersnyahnyah>
Date: Thu, 23 Aug 2007 13:09:02 -0400

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