GRIB Lambert proj errors

I could use some help on getting a GRIB to project over a Lambert Conformal.
I get these errors:
warning:Attempt to reference attribute (Latin1) which is undefined
warning:Attempt to reference attribute (Latin2) which is undefined
warning:GKS:GSVP: --RECTANGLE DEFINITION IS INVALID
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT/GKS ERROR REPORTED
FROM libhlu.a(Error.o)
warning:GKS:GSVP: --RECTANGLE DEFINITION IS INVALID
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT/GKS ERROR REPORTED
FROM libhlu.a(Error.o)
warning:["SetValues.c":474]:SetValuesChild never occurred on
map.PlotManager: Error in mapPlotClass
warning:tiMainString is not a valid resource in map at this time
Segmentation fault

I'm not sure what all these errors mean since I'm new to the language. Here
is the code:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin

  fili = "/home/jlahowet/THOR/test/wind/narr-a_221_20070314_0000_000.grb"
  f = addfile (fili, "r")

  names = getfilevarnames(f)
  print(names)

  atts = getfilevaratts(f,names(0))
  dims = getfilevardims(f,names(0))
  print(atts)
  print(dims)

  vpt = f->VSTM_221_HTGY
  lat2d = f->gridlat_221
  lon2d = f->gridlon_221

  wks = gsn_open_wks("ps","lcnative")
  gsn_define_colormap(wks,"gui_default")

  res = True
  res_at_tiMainString = "Native Lambert Conformal"

  res_at_mpLimitMode = "Corners"
  res_at_mpLeftCornerLatF = lat2d_at_corners(0)
  res_at_mpLeftCornerLonF = lon2d_at_corners(0)
  res_at_mpRightCornerLatF = lat2d_at_corners(2)
  res_at_mpRightCornerLonF = lon2d_at_corners(2)

  res_at_mpProjection = "LambertConformal"
  res_at_mpLambertParallel1F = lat2d_at_Latin1
  res_at_mpLambertParallel2F = lat2d_at_Latin2
  res_at_mpLambertMeridianF = lat2d_at_Lov

  res_at_pmTickMarkDisplayMode = "Always"
  res_at_mpFillOn = False
  res_at_mpOutlineDrawOrder = "PostDraw"
  res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"

  res_at_tfDoNDCOverlay = True

  res_at_cnFillOn = True
  res_at_cnLinesOn = False
  res_at_gsnSpreadColors = True
  res_at_gsnAddCyclic = False

  plot = gsn_csm_contour_map(wks,vpt,res)
end

Thanks

-- 
Jamie Ryan Lahowetz
University of Nebraska - Lincoln
Graduate Student - Geosciences
402.304.0766
jlahowe2_at_bigred.unl.edu

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 13 2008 - 08:40:12 MDT