Re: GRIB Lambert proj errors

On Wed, 13 Aug 2008, Jamie Lahowetz wrote:

> 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:

Jamie,

The main problem is that you are trying to reference a couple of
attributes that don't exist:

   res_at_mpLambertParallel1F = lat2d_at_Latin1
   res_at_mpLambertParallel2F = lat2d_at_Latin2

This code is trying to reference attributes "Latin1" and "Latin2"
attached to the variable "lat2d".

I think this was an original example where the "lat2d" variable did
indeed have these attributes, because they were originally on the
file. However, you are using your own file, and it looks like
f->gridlat_221 does not have "Latin1" or "Latin2" attached to it.

What you need to do is first print the contents of your file, either
by adding a "print(f)" after the "addfile" call, or by doing an
ncl_filedump on the UNIX command line:

   ncl_filedump /home/jlahowet/THOR/test/wind/narr-a_221_20070314_0000_000.grb

Check the output closely to see if there's anything that looks like it
might be related to the map projection of your data, and especially,
that might give you the two lambert parallel values that you need.

For example, in one NARR example that we have, the two parallels are
attached to a variable called "Lambert_Conformal" and are stored in an
array called "standard_parallel" that has two elements:

   res_at_mpLambertParallel1F = f->Lambert_Conformal_at_standard_parallel(0)
   res_at_mpLambertParallel2F = f->Lambert_Conformal_at_standard_parallel(1)

Hopefully this will help you get started.

--Mary

> 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 - 09:20:20 MDT