NCL error message

Hi all,

Has anyone encountered this error message using NCL? We are running
quite a lot of NCL jobs every day and this error message occurrs
sporadically und unreproducibly...

(0) jmb_getvar: reading HSURF named HH_GDS10_SFC in file
fatal:ContourPlotDraw: ARSCAM/ARPRAM - ALGORITHM FAILURE
fatal:ContourPlotDraw: draw error
warning:WorkstationDeactivate: workstation not active or not opened
(0) jmb_getvar: reading TOT_PREC named TOT_PREC_GDS10_SFC_acc3h in
file

Cheers,
Oli

The script where the error occurs is the following...

  ; setup plot
  name = "${model}_${NAME}_${start}.gin"
  format = "png"
  xsize = 500.0
  ysize = 400.0
  shift =392

  ; set plot domain
  xs = 400.0
  xe = 900.0
  ys = 0.0
  ye = 400.0

  ; open files
  c = addfile("lfff00000000c.grb", "r")
  f = addfile("_NCLINP_${start}_delta03h_c.grb", "r")

  ; read data
  d = jmb_getvar(f, "TOT_PREC") ; reads a variable from the grib file
  
  ; geo-reference
  jmb_getgrid(f, c, d) ; attaches geo-referencing data to
the field d

  ; compute swissgrid coordinates
  x = latlon2swissx(d@lat2d, d@lon2d)
  y = latlon2swissy(d@lat2d, d@lon2d)

  ; remove geo-referencing
  jmb_rmgrid(d) ; removed geo-referencing data from
field d

  ; set plot size in inch
  dpi = ceil(xsize/8.5/72.0)*72.0
  xinch = xsize/dpi
  yinch = ysize/dpi

  ; open graphic port
  r = "ps"
  r@wkDeviceLowerX = 0.0
  r@wkDeviceUpperX = xinch*72.0
  r@wkDeviceLowerY = 0.0
  r@wkDeviceUpperY = yinch*72.0
  wks = gsn_open_wks(r, name)
  delete(r)
  rct = jmb_set_ct(wks,"rr_ninjo_21lev",False)

  ; setup plot
  r = rct
  r@gsnMaximize = False
  r@gsnTickMarksOn = False
  r@vpWidthF = 1.0
  r@vpHeightF = 1.0
  r@vpXF = 0.0
  r@vpYF = 1.0
  r@vpKeepAspect = True
  r@sfXArray = x
  r@sfYArray = y
  r@trGridType = "TriangularMesh"
  r@trXMinF = xs
  r@trYMinF = ys
  r@trXMaxF = xe
  r@trYMaxF = ye
  r@cnFillOn = True
  r@cnLinesOn = False
  r@cnLineLabelsOn = False
  r@cnInfoLabelOn = False

  ; make plot
  cn = gsn_contour(wks, d, r)
  delete(r)

  ; close graphic port
  delete(wks)

 
________________________________________

Oliver Fuhrer
Numerical Models

Federal Departement of Home Affairs FDHA
Federal Office of Meteorology and Climatology MeteoSwiss

Kraehbuehlstrasse 58, P.O. Box 514, CH-8044 Zurich, Switzerland

Tel. +41 44 256 93 59
Fax +41 44 256 92 78
oliver.fuhrer@meteoswiss.ch
www.meteoswiss.ch - First-hand information
  
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Mar 26 05:21:48 2010