Function wrf_contour.

Hi dear WRF Help,

I am currently encountering a problem while trying to use the function
wrf_contour.
http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_contour.shtml
Using it as in the script example, I found a strange output such as
PLCHHQ - CHARACTER NUMBER 33 (t) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 34 (o) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 41 (b) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 42 (y) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 13 (9) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 17 (-) IS NOT A LEGAL FUNCTION CODE
PLCHHQ - CHARACTER NUMBER 20 (_) IS NOT A LEGAL FUNCTION CODE

The only thing that may be of concern is that the projection of the
netcdf file is polar (AMP domain)
Here is the script

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load
"$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;-------------------------------------------------------------------------------;

begin
  
;-------------------------------------------------------------------------------;

  ; open file and read in data
  
;-------------------------------------------------------------------------------;

 works="ps" ; open a ps file
 region = "ampsrt"
 ;
---------------------------------------------------------------------------
 res = True
 res_at_gsnDraw = True ; Forces the plot to be drawn
 res_at_gsnFrame = True ; Frame advance
 res_at_cnFillOn = True

 filename_inc = "/ptmp/yann/data/ampsrt/rc/2009052400/wrfinput_d01.nc"
 ff1 = addfile(filename_inc,"r")
 hgt1 = wrf_user_getvar(ff1,"HGT",0) ; Get terrain height for time 0

 wks = gsn_open_wks("ps","height_"+region) ; open a ps file
 plts = wrf_contour(ff1,wks,hgt1,res)

end

NCARG_ROOT is /users/local on bluefire
Many thanks for your help!

Yann MICHEL

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 29 2009 - 09:42:15 MDT