Error with plotting Hovmueller gsn_csm_hov

I used the 1st Hovmueller example on at http://www.ncl.ucar.edu/Applications/time_lon.shtml
to try and create a plot of WRF model data.

The data is reading fine but I am getting a plot error that I do not
understand. I am getting an error due to the gsn_csm_hov function.
(All readout and errors are attached to this email)

May I have some help with this?
-Erik

0) Working on time: 2006-09-01_00:00:00
.
. ( I cut for space in email)
.
(0) Working on time: 2006-09-30_00:00:00

Variable: v_plane
Type: float
Total Size: 60604 bytes
             15151 values
Number of Dimensions: 2
Dimensions and sizes: [109] x [139]
Coordinates:
Number Of Attributes: 6
   PlotLevelID : 700 hPa
   units : m/s
   description : v Velocity
   lon2d : <ARRAY of 15151 elements>
   lat2d : <ARRAY of 15151 elements>
   _FillValue : -999999

warning:tmYLMinorLengthF is not a valid resource in
hov_contour.PlotManager at this time
warning:NhlGetValues:Error retrieving tmYLMinorLengthF
fatal:Execute: Error occurred at or near line 2293 in file $NCARG_ROOT/
lib/ncarg/nclscripts/csm/gsn_code.ncl

fatal:Execute: Error occurred at or near line 7875 in file $NCARG_ROOT/
lib/ncarg/nclscripts/csm/gsn_code.ncl

fatal:Execute: Error occurred at or near line 10500 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl

fatal:Execute: Error occurred at or near line 10632 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl

fatal:Execute: Error occurred at or near line 36 in file hov_WRFV700.ncl

; ===========================================
; hov_1.ncl
; ===========================================
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
;WRF Out file
diri1 ="/Volumes/Data_and_Models/Model-Output/Athena/"
ifile1 ="WRF-SOP3_Athena_3_1_2.nc"
a=addfile(diri1+ifile1,"r")
; Which times and how many time steps are in this data set?
      times = wrf_user_list_times(a) ; get times in the file
      ntimes = dimsizes(times) ; number of times in the file
;;;;;;;Get Variables;;;;;;;;;;
times = wrf_user_list_times(a)
do it = 248,ntimes-1, 8 ;;;;;Time step you
want!!!!
        print("Working on time: " + times(it) )
  pressure = 700.
  v = wrf_user_getvar(a,"va", it) ; v averaged to mass points
  p = wrf_user_getvar(a, "pressure", it) ; pressure
  v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
  v_plane_at_lat2d = a->XLAT(0,:,:)
  v_plane_at_lon2d = a->XLONG(0,:,:)
end do
printVarSummary(v_plane)
;==============================================
; creat plot
;==============================================
  wks = gsn_open_wks ("x11", "hov" ) ; open workstation
  res = True ; plot mods desired
  res_at_tiMainString = "Default Hovmueller" ; title
  plot = gsn_csm_hov(wks, v_plane, res ) ; default hovmueller
end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Nov 09 2007 - 13:45:57 MST