variables undefined?

I am new to NCL but have a weird bug:
I am running examples form NCL website
on a Mac PowerBook G4, NCL version 4.2.0.a032
error message says some variables undefined in gsn_code.ncl
but when I look there they are defined.
Not sure what the problem is. On the example below, I manage to get by
with a x11 or pdf plot but not a ps (when defining wks). Other variables
undefined as well (polytype).

output below, as well as source ncl code and my calls:
------------------------------------------------------------------
fatal:syntax error: line 2054 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl before or near
reset_device_coordinates
procedure reset_device_coordinates
---------------------------------^

fatal:error in statement
fatal:Variable (wks) is undefined
fatal:Execute: Error occurred at or near line 2062 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl

fatal:syntax error: line 2923 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl before or near
gsn_primitive
procedure gsn_primitive
----------------------^

fatal:error in statement
fatal:syntax error: line 2926 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl before or near local
local
----^

fatal:error in statement
fatal:Variable (polytype) is undefined
fatal:Execute: Error occurred at or near line 2929 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl

----------------------------------------------------------------------

here is the source for those spots:

procedure reset_device_coordinates(wks:graphic)
begin
  setvalues wks
    "wkOrientation" : get_res_value(wks,"wkOrientation", 0)
    "wkDeviceLowerX" : get_res_value(wks,"wkDeviceLowerX", 36)
    "wkDeviceLowerY" : get_res_value(wks,"wkDeviceLowerY",126)
    "wkDeviceUpperX" : get_res_value(wks,"wkDeviceUpperX",576)
    "wkDeviceUpperY" : get_res_value(wks,"wkDeviceUpperY",666)
  end setvalues
end

AND

procedure gsn_primitive(wks:graphic,plotid:graphic,x[*]:numeric,\
                        y[*]:numeric,isndc:logical,polytype:string, \
                        resources:logical)
local res2, gsid, gsres, gmres, canvas, xf, yf, x2, y2
begin
  if(.not.any(polytype.eq.(/"polymarker","polygon","polyline"/)))
    print("Warning: gsn_primitive: Do not recognize primitive type '"+
polytype + "'.")
    return
  end if
.......... etc.........

----------------------- here are my calls -------------------------

;*************************************************
; adapted solution and initial solution
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
;************************************************
; open file and read in variable
;************************************************
x1=asciiread("fort.16",(/13/),"float")
x2=asciiread("fort.36",(/49/),"float")
y1=asciiread("fort.26",(/13/),"float")
y2=asciiread("fort.46",(/49/),"float")
;************************************************
; plot the correlations
;************************************************
  wkspdf = gsn_open_wks("pdf","adapted") ; open a pdf plot
  wksx11 = gsn_open_wks("x11","adapted") ; open a x11 plot
  wksps1 = gsn_open_wks("ps","adapted1") ; open a ps plot
  wksps2 = gsn_open_ps("ps","adapted2") ; open a ps plot
  res1 = True ; make plot mods

  res1_at_xyLineColor = "blue" ; change line color
  res1_at_xyLineThicknesses = 2. ; 1 is default

  res2 = True ; make plot mods
  res2_at_tiMainString = "Initial Condition (blue) and Solution (red) at
t" ; title
  res1_at_tiXAxisString = "X" ; x-axis label
  res1_at_tiYAxisString = "U" ; y-axis label

  res2_at_xyLineColor = "red" ; change line color
  res2_at_xyLineThicknesses = 3. ; 1 is default

  plot =
gsn_csm_x2y2(wksps1,x1,x2,y1,y2,res1,res2) ; plot
correlation
  plot =
gsn_csm_x2y2(wksps2,x1,x2,y1,y2,res1,res2) ; plot
correlation
  plot =
gsn_csm_x2y2(wkspdf,x1,x2,y1,y2,res1,res2) ; plot
correlation
  plot =
gsn_csm_x2y2(wksx11,x1,x2,y1,y2,res1,res2) ; plot
correlation

;************************************************
end

any help appreciated
Thanks

-- 
Catherine Mavriplis, PhD
Research Scientist, 
University of Oklahoma 
Cooperative Institute for Mesoscale Meteorological Studies (CIMMS)
NOAA National Severe Storms Laboratory
1313 Halley Circle,
Norman, OK 73069 
(405) 366-0534 Voice
(405) 366-0472 FAX
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk