Hi, I am using gsn_csm_contour_map_ce and am getting the following
warning messages per frame:
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransInitialize: error creating spline approximation for
trXCoordPoints; defaulting to linear
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransSetValues: error creating spline approximation for
trXCoordPoints; defaulting to linear
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransSetValues: error creating spline approximation for
trXCoordPoints; defaulting to linear
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransSetValues: error creating spline approximation for
trXCoordPoints; defaulting to linear
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransSetValues: error creating spline approximation for
trXCoordPoints; defaulting to linear
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for
X axis failed: consider adjusting trXTensionF value
warning:IrTransSetValues: error creating spline approximation for
trXCoordPoints; defaulting to linear
The plots look fine -- how can I suppress these messages, or fix what
they are warning about? My script is attached, as well as an ncdump
of the data file. (Also, I would like to make this same plot but without
any continental outlines) Thanks,
--Ben
-----------------------------------------------------------------------
Ben Foster High Altitude Observatory (HAO)
foster@ucar.edu phone: 303-497-1595 fax: 303-497-1589
Nat. Center for Atmos. Res. P.O. Box 3000 Boulder CO 80307 USA
-----------------------------------------------------------------------
;
; Plot ncep data. netcdf file made from /home/tgcm/mkncep.
;
; For dataproc:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;
; For dave (blackforest):
;load "$NCARG_ROOT/lib32/r4i4/ncarg/nclscripts/csm/gsn_code.ncl"
;load "$NCARG_ROOT/lib32/r4i4/ncarg/nclscripts/csm/gsn_csm.ncl"
;
;dimensions:
; day = UNLIMITED ; // (177 currently)
; lat = 36 ;
; lon = 73 ;
; nw = 12 ;
; lev = 3 ;
;variables:
; double lat(lat) ;
; lat:long_name = "geographic longitude (-west, +east)" ;
; lat:units = "degrees" ;
; double lon(lon) ;
; lon:units = "degrees" ;
; double lev(lev) ;
; lev:long_name = "pressure levels" ;
; lev:units = "mb" ;
; int day(day) ;
; day:long_name = "calendar day (yyddd)" ;
; day:units = "days" ;
; double ZM(day, lev, lat) ;
; ZM:long_name = "ZM" ;
; double ZS(day, lev, nw, lat) ;
; ZS:long_name = "ZS" ;
; double ZC(day, lev, nw, lat) ;
; ZC:long_name = "ZC" ;
; double HT(day, lat, lon) ;
; HT:long_name = "HT" ;
;
begin
filename = "2002001-2002181.nc"
ncfile = addfile(filename,"r")
zm = ncfile->ZM
zs = ncfile->ZS
zc = ncfile->ZC
ht = ncfile->HT
;
lat = ncfile->lat
nlat = dimsizes(lat)
;
lev = ncfile->lev
nlev = dimsizes(lev)
;
lon = ncfile->lon
nlon = dimsizes(lon)
;
days = ncfile->day
ndays = dimsizes(days)
;
print("ndays="+ndays+" nlat="+nlat+" nlev="+nlev+" nlon="+nlon)
;
; x11, ps, or ncgm:
;
x11 = 0
ps = 0
cgm = 1
wksflnm = "ncep"
if (x11.eq.1) then
type = "x11"
wks = gsn_open_wks(type,wksflnm)
print("Opened X11 workstation: file "+wksflnm)
else
if (ps.eq.1) then
type = "ps"
wks = gsn_open_wks(type,wksflnm)
print("Opened PS workstation: file "+wksflnm)
else
if (cgm.eq.1) then
type = "ncgm"
wks = gsn_open_wks(type,wksflnm)
print("Opened CGM workstation: file "+wksflnm)
end if
end if
end if
res = True
; res@mpFillOn = False
res@cnLevelSpacingF = 0.2
;
; do n=0,ndays-1
do n=0,10
; res@tiMainString = "Height (km) at 10 mb"+" YearDay = "+days(n)
title = "Height (km) at 10 mb"+" YearDay = "+days(n)
res@gsnCenterString = "YearDay = "+days(n)
print(""+title)
contour = gsn_csm_contour_map_ce(wks,ht(n,:,:),res)
end do
end
netcdf 2002001-2002181 {
dimensions:
day = UNLIMITED ; // (177 currently)
lat = 36 ;
lon = 73 ;
nw = 12 ;
lev = 3 ;
variables:
double lat(lat) ;
lat:long_name = "geographic latitude (-south, +north)" ;
lat:units = "degrees_north" ;
double lon(lon) ;
lon:long_name = "geographic longitude (-west, +east)" ;
lon:units = "degrees_east" ;
double lev(lev) ;
lev:long_name = "pressure levels" ;
lev:units = "mb" ;
int day(day) ;
day:long_name = "calendar day (yyyyddd)" ;
day:units = "days" ;
double ZM(day, lev, lat) ;
ZM:long_name = "ZM" ;
double ZS(day, lev, nw, lat) ;
ZS:long_name = "ZS" ;
double ZC(day, lev, nw, lat) ;
ZC:long_name = "ZC" ;
double HT(day, lat, lon) ;
HT:long_name = "Height at 10 mb" ;
HT:units = "km" ;
data:
lat = -87.5, -82.5, -77.5, -72.5, -67.5, -62.5, -57.5, -52.5, -47.5, -42.5,
-37.5, -32.5, -27.5, -22.5, -17.5, -12.5, -7.5, -2.5, 2.5, 7.5, 12.5,
17.5, 22.5, 27.5, 32.5, 37.5, 42.5, 47.5, 52.5, 57.5, 62.5, 67.5, 72.5,
77.5, 82.5, 87.5 ;
lon = -180, -175, -170, -165, -160, -155, -150, -145, -140, -135, -130,
-125, -120, -115, -110, -105, -100, -95, -90, -85, -80, -75, -70, -65,
-60, -55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0, 5, 10, 15,
20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105,
110, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175, 180 ;
lev = 30, 10, 5 ;
day = 2002005, 2002006, 2002007, 2002008, 2002009, 2002010, 2002011,
2002012, 2002013, 2002014, 2002015, 2002016, 2002017, 2002018, 2002019,
2002020, 2002021, 2002022, 2002023, 2002024, 2002025, 2002026, 2002027,
2002028, 2002029, 2002030, 2002031, 2002032, 2002033, 2002034, 2002035,
2002036, 2002037, 2002038, 2002039, 2002040, 2002041, 2002042, 2002043,
2002044, 2002045, 2002046, 2002047, 2002048, 2002049, 2002050, 2002051,
2002052, 2002053, 2002054, 2002055, 2002056, 2002057, 2002058, 2002059,
2002060, 2002061, 2002062, 2002063, 2002064, 2002065, 2002066, 2002067,
2002068, 2002069, 2002070, 2002071, 2002072, 2002073, 2002074, 2002075,
2002076, 2002077, 2002078, 2002079, 2002080, 2002081, 2002082, 2002083,
2002084, 2002085, 2002086, 2002087, 2002088, 2002089, 2002090, 2002091,
2002092, 2002093, 2002094, 2002095, 2002096, 2002097, 2002098, 2002099,
2002100, 2002101, 2002102, 2002103, 2002104, 2002105, 2002106, 2002107,
2002108, 2002109, 2002110, 2002111, 2002112, 2002113, 2002114, 2002115,
2002116, 2002117, 2002118, 2002119, 2002120, 2002121, 2002122, 2002123,
2002124, 2002125, 2002126, 2002127, 2002128, 2002129, 2002130, 2002131,
2002132, 2002133, 2002134, 2002135, 2002136, 2002137, 2002138, 2002139,
2002140, 2002141, 2002142, 2002143, 2002144, 2002145, 2002146, 2002147,
2002148, 2002149, 2002150, 2002151, 2002152, 2002153, 2002154, 2002155,
2002156, 2002157, 2002158, 2002159, 2002160, 2002161, 2002162, 2002163,
2002164, 2002165, 2002166, 2002167, 2002168, 2002169, 2002170, 2002171,
2002172, 2002173, 2002174, 2002175, 2002176, 2002177, 2002178, 2002179,
2002180, 2002181 ;
}
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Thu Aug 08 2002 - 11:10:27 MDT