I'm using gsn_csm_contour_map_polar to contour some data
over the Antarctic and I only get 5 of the expected 12
longitude labels. Namely, with 0 degrees at the top,
I only see those in the lower half of the plot, i.e.,
180 deg and two more to either side (@30 deg intervals).
I'm also getting the following set of error messages,
once for each missing label:
fatal:NhlCvtStringToEnum: Unable to convert string "missing"
to requested type
warning:Error retrieving resource amJust from args - Ignoring Arg
fatal:NhlCvtStringToEnum: Unable to convert string "missing" to
requested type
warning:Error retrieving resource amSide from args - Ignoring Arg
Thus this appears to be something to do with amJust and amSide.
I tried setting amTrackData (somewhat blindly) but get another
error message (warning:amTrackData is not a valid resource in
ezplot1p_map at this time). Rather than continuing to thrash
(I'm very new to ncl), I'm hoping someone else has seen this
or otherwise knows how to fix it.
We are using ncl 4.2.0.a024.
The complete script is below.
thanks,
dave reusch
dbr AT geosc.psu.edu
--------------
load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
Z700=asciiread("Zsub.txt",(/721/),"float")
ll=asciiread("EASEsub.txt",(/721,2/),"float")
d0 = ll(:,0)
d1 = ll(:,1)
znew = idsfft( d0, d1, Z700, (/30,30/) )
znew!0="lat"
znew!1="lon"
znew&lat@units="degrees_north"
znew&lon@units="degrees_east"
xwks = gsn_open_wks( "x11","test") ; Open an X11 workstation.
gsn_define_colormap(xwks,"temp1") ; choose colormap
res = True
res@cnFillOn = True ; turn on color
res@gsnPolar = "SH"
res@mpFillOn = 1
res@mpOutlineOn = "True"
res@vpYF = 0.9
res@lbLabelAngleF = 45
res@lbTitleOn = True ; turn on title
res@lbTitleString = "700 mb Geopotential Height"
res@lbTitleFontHeightF= .015
res@gsnSpreadColors = True ; use full colormap
res@mpMaxLatF = max(znew&lat) ; zoom in on region
res@mpMinLatF = min(znew&lat)
res@mpMinLonF = min(znew&lon)
res@mpMaxLonF = max(znew&lon)
res@tiMainString = "Z700 Annual Average"
plot = gsn_csm_contour_map_polar(xwks,znew,res)
end
_______________________________________________
ncl-talk mailing list
ncl-talk AT ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Mon Aug 05 2002 - 17:27:15 MDT