Hi Helen,
the warning tells you that it cannot interpolate in the y-axis. Means the numbers must be rubbish.
have you checked, if you produce meaningfull numbers?
I think your averaging is somehow wrong. Better than reordering would be to use dim_avg_n_Wrap where you can chose the dimension(s) to be averaged.
And than check them with printVarSummary before summing up everything.
Only, if you are sure you have something reasonable in the longitude and time, try to plot.
Good luck, Kristina
"Parish, Helen" <hparish_at_epss.ucla.edu> schrieb am 12:37 Freitag, 24.Januar 2014:
I am getting the following message when creating a Hovmuller plot, but I don't know what this means :
" warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for Y axis failed: consider adjusting trYTensionF value
warning:IrTransInitialize: error creating spline approximation for trYCoordPoints; defaulting to linear "
Can anyone tell me what this warning means, whether I need to be concerned, and whether I need to adjust something to produce a meaningful plot?.
I enclose a copy of my script below:
Thanks,
Helen.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
; ********************
wks = gsn_open_wks ("pdf", "hovdailyzon350days4" ) ; open workstation
gsn_define_colormap(wks,"rainbow") ; choose colormap
; ************************
diri = "./"
fili = systemfunc("cd "+diri+" ; ls cesm*cam2.h0*nc")
print(fili)
nfili = dimsizes(fili)
; define pi, universal gas constant (rgas), Venus planetary rotation rate (omegarot) etc.
pi = 3.1415926
twopi = 2.0 * 3.1415926536
rgas = 8.314468 ; gas constant (J/K/mole)
radi = 6.05180e6 ; radius of Venus in meters
omegarot = 2.992677e-7 ; rotation rate of Venus
mwt = 43.45 ; mean molecular weight (kg/mole)
grav = 8.87 ; gravitataional acceleration (m/s^2)
cp = 8.5e2 ; specific heat at constant pressure (J/Kg/K or m^2/s^2/K)
do nf = 0, nfili-1
nf100 = nf + 100
fi = addfile (diri+fili(nf), "r")
diro = "./"
filo = "uhelen"+nf100+".nc"
print(filo)
system ("/bin/rm -f "+diro+filo) ; remove any pre-exist file
fo = addfile (diro+filo, "c")
hyam = fi->hyam ; read hybrid info
hybm = fi->hybm
hyai = fi->hyai ; read hybrid info
hybi = fi->hybi
U = fi->U
lat = fi->lat
lev = fi->lev
lon = fi->lon
fo->U = U ; write to netCDF file
fo->lat = lat ; write to netCDF file
fo->lon = lon ; write to netCDF file
fo->lev = lev ; write to netCDF file
end do
; ******************************************
;
; *******************************************
diri = "./"
fili2 = systemfunc("cd "+diri+" ; ls uhelen*nc")
print(fili2)
nfili = dimsizes(fili2)
f2 = addfiles (fili2,"r")
U2= addfiles_GetVar(f2,fili2,"U")
printVarSummary( U2 )
dimt4 = dimsizes( U2 )
ntim4 = dimt4(0)
klvl4 = dimt4(1)
nlat4 = dimt4(2)
mlon4 = dimt4(3)
ureorder = U2(lev | :, lat | :, lon |:, time | :) ; (lev,lat,lon,time)
printVarSummary( ureorder )
uzonav = dim_avg_Wrap(U2) ; (time,lev,lat)
printVarSummary( uzonav )
utimeav = dim_avg_Wrap(ureorder) ; (lev,lat,lon)
printVarSummary( utimeav )
UTIMEAV = U2 ; copy U2 coordinate variables to UZONAV
printVarSummary( UTIMEAV )
UTIMEAV = conform(U2,utimeav,(/1,2,3/)) ; (time,lev,lat,lon)
udiff = U2 ; copy U2 coordinate variables to udiff
udiff = U2 - UTIMEAV ; (time,lev,lat,lon)
printVarSummary( udiff )
;***********************
; Create Plot
;***********************
res = True
res_at_gsnLeftString = ""
res_at_gsnCenterString = "Hovmuller Plot U-Utime-mean - Equator"
res_at_gsnRightString = ""
res_at_cnFillOn = True
res_at_lbLabelAutoStride = True
res_at_gsnMaximize = True ; if [ps, eps, pdf] make large
res_at_gsnSpreadColors = True ; span color map
res_at_cnLineLabelsOn = False ; Turn off contour line labels
res_at_cnInfoLabelOn = False ; Turn off informational label
res_at_tiXAxisString = "Longitude (degrees)"
res_at_tiYAxisString = "Time (Earth days)"
res_at_pmTitleZone = 3
res_at_pmLabelBarOrthogonalPosF = 0.05
res_at_lbTitleOn = True
res_at_lbTitleDirection = "Across" ; title direction
res_at_lbTitlePosition = "Right"
res_at_lbTitleString = "(m s~S~-1~N~)"
res_at_lbTitleFontHeightF = 0.02
res_at_lbTitleOffsetF = 0.0
plot = gsn_csm_hov(wks, udiff(:,19,95,:), res ) ; (time,lon)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jan 24 05:07:19 2014
This archive was generated by hypermail 2.1.8 : Fri Feb 07 2014 - 16:39:11 MST