Hovmuller plot

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@gsnLeftString = ""
    res@gsnCenterString = "Hovmuller Plot U-Utime-mean - Equator"
    res@gsnRightString = ""

 res@cnFillOn = True
 res@lbLabelAutoStride = True
 res@gsnMaximize = True ; if [ps, eps, pdf] make large
 res@gsnSpreadColors = True ; span color map
 res@cnLineLabelsOn = False ; Turn off contour line labels
 res@cnInfoLabelOn = False ; Turn off informational label

 res@tiXAxisString = "Longitude (degrees)"
    res@tiYAxisString = "Time (Earth days)"

   res@pmTitleZone = 3
   res@pmLabelBarOrthogonalPosF = 0.05

       res@lbTitleOn = True
       res@lbTitleDirection = "Across" ; title direction
       res@lbTitlePosition = "Right"
       res@lbTitleString = "(m s~S~-1~N~)"
       res@lbTitleFontHeightF = 0.02
       res@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
Received on Fri Jan 24 04:37:09 2014