Re: Hovmuller plot

From: Kristina Fröhlich <kristina_froehlich_at_nyahnyahspammersnyahnyah>
Date: Fri Jan 24 2014 - 05:07:06 MST

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