PLCHHQ error

From: 감은아 <fresia01_at_nyahnyahspammersnyahnyah>
Date: Tue, 01 Sep 2009 21:47:32 +0900

Dear NCL users,
when I run ncl, i am getting the following error,
 PLCHHQ - CHARACTER NUMBER 34 (t) IS NOT A LEGAL FUNCTION CODE PLCHHQ - CHARACTER NUMBER 35 (o) IS NOT A LEGAL FUNCTION CODE PLCHHQ - CHARACTER NUMBER 42 (b) IS NOT A LEGAL FUNCTION CODE PLCHHQ - CHARACTER NUMBER 43 (y) IS NOT A LEGAL FUNCTION CODE
I would appreciate if anybody could help me how to fix it.
hear is my code:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
  FILES = systemfunc( "ls ../geo_em.*.nc" )  numFILES = dimsizes(FILES)
  t = numFILES
   type = "x11";   type = "png";   type = "ncgm";   type = "pdf";   type = "ps"   wks = gsn_open_wks(type,"geo_em_test")
   setvalues wks     "wkBackgroundColor" : "white"     "wkForegroundColor" : "black"   end setvalues
   gsn_define_colormap(wks,"gui_default")
  opts = True  opts_at_InitTime = False  opts_at_Footer = False
  cnres = opts  cnres_at_gsnDraw = True  cnres_at_gsnFrame = True
  cnres_at_cnInfoLabelOn = False
  lnres = True  lnres_at_gsLineThicknessF = 1.8  lnres_at_gsLineColor = "white"  lnres_at_gsLineLabelString = "d02"
  a = addfile(FILES(0),"r")
  ter = wrf_user_getvar(a,"HGT_M",0)  cnres_at_cnFillOn = True  cnres_at_ContourParameters =(/0,6000,50/)
  res = cnres
  pltres = True  pltres_at_FramePlot = False
  mpres = True  mpres_at_mpGeophysicalLineColor = "black"  mpres_at_mpGridLineColor = "black"  mpres_at_gsnMaximize = True
  contour=wrf_contour(a,wks,ter,res)  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
  if ( t .eq. 3  ) then
     overlay(plot, contour)
     b = addfile(FILES(1),"r")
     bter = wrf_user_getvar(b,"HGT_M",0)
     btlat = b->XLAT_M(0,:,:)     btlon = b->XLONG_M(0,:,:)    res2 = cnres     res2_at_lbLabelBarOn = False     res2_at_gsnMaximize = False
     dims = dimsizes(bter)     r2 = dimsizes(dims)     r21 = dims(r2-1)     r22 = dims(r2-2)
     xbox = (/btlon(0,0), btlon(0,r21-1), btlon(r22-1, r21-1), \              btlon(r22-1,0), btlon(0,0)/)     ybox = (/btlat(0,0), btlat(0,r21-1), btlat(r22-1, r21-1), \              btlat(r22-1,0), btlat(0,0)/)
     x_out = new(dimsizes(xbox), typeof(xbox))     y_out = new(dimsizes(ybox), typeof(ybox))
     datatondc(plot, xbox, ybox, x_out, y_out)     gsn_polyline_ndc(wks, x_out, y_out, lnres)
     c = addfile(FILES(2),"r")
     cter = wrf_user_getvar(c,"HGT_M",0)     ctlat = c->XLAT_M(0,:,:)     ctlon = c->XLONG_M(0,:,:)
     xbox = (/ctlon(0,0), ctlon(0,r21-1), ctlon(r22-1, r21-1), \              ctlon(r22-1,0), ctlon(0,0)/)     ybox = (/ctlat(0,0), ctlat(0,r21-1), ctlat(r22-1, r21-1), \              ctlat(r22-1,0), ctlat(0,0)/)
     x_out = new(dimsizes(xbox), typeof(xbox))     y_out = new(dimsizes(ybox), typeof(ybox))
     datatondc(plot, xbox, ybox, x_out, y_out)     gsn_polyline_ndc(wks, x_out, y_out, lnres)
        frame(wks)     end if
  if ( t .eq. 2  ) then
     overlay(plot, contour)
     b = addfile(FILES(1),"r")
     bter = wrf_user_getvar(b,"HGT_M",0)
     btlat = b->XLAT_M(0,:,:)     btlon = b->XLONG_M(0,:,:)
     res2 = cnres     res2_at_lbLabelBarOn = False     res2_at_gsnMaximize = False
     dims = dimsizes(bter)     r2 = dimsizes(dims)     r21 = dims(r2-1)     r22 = dims(r2-2)
     xbox = (/btlon(0,0), btlon(0,r21-1), btlon(r22-1, r21-1), \              btlon(r22-1,0), btlon(0,0)/)     ybox = (/btlat(0,0), btlat(0,r21-1), btlat(r22-1, r21-1), \              btlat(r22-1,0), btlat(0,0)/)
     x_out = new(dimsizes(xbox), typeof(xbox))     y_out = new(dimsizes(ybox), typeof(ybox))
     datatondc(plot, xbox, ybox, x_out, y_out)     gsn_polyline_ndc(wks, x_out, y_out, lnres)
     frame(wks)     end if
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 01 2009 - 06:47:32 MDT

This archive was generated by hypermail 2.2.0 : Thu Sep 03 2009 - 07:55:08 MDT