<HTML>
<HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META content="OPENWEBMAIL" name=GENERATOR>
</HEAD>
<BODY bgColor=#ffffff>
<P><FONT size=2>Dan,</FONT></P>
<P><FONT size=2>"d hgtprs(lev=700)"</FONT></P>
<P><FONT size=2> should work, I guess.</FONT></P>
<P><FONT size=2>Ricardo</FONT></P>
<P><FONT size=2><STRONG>On Sat, 23 Jun 2007 08:18:51 -0500, Dan Leins wrote</STRONG>
<BR>> All,
<BR>>
<BR>> I have a series of scripts I use to plot fields at numerous model levels, (i.e. heights at 925mb, 850mb, 700mb, etc...) from a grib file. I dump the grib files out to a control file using grib2ctl and go from there. To do this, I simply "d hgtprs(z=7)" or "d hgtprs(z=11)" or whatever z level I'm interested in. However, is there a way to plot a specific level without referring to Z? In short, I wrote my scripts to use on one specific set of grib files, but occasionally I will need to use my scripts on a different set of grib files which do not contain the same levels. So usually Z=7 would be equivalent to 850mb, but in my other set of grib files z=7 is actually equivalent to 250mb which is a big problem. Any workarounds or other ways to refer to the model levels or will I need a completely different set of scripts for my other grib files?
<BR>>
<BR>> Thanks,
<BR>> Dan
<BR>
<BR>
<BR>Hallak
<BR>
<BR>Laboratório de Hidrometeorologia (LABHIDRO)
<BR>Departamento de Ciências Atmosféricas (DCA)
<BR>Instituto de Astronomia, Geofísica e Ciências Atmosféricas (IAG)
<BR>Universidade de São Paulo (USP)
<BR>Rua do Matão, 1226 - Cidade Universitária -São Paulo SP -05508-900
<BR>Fone.: +55 (0-XX-11) 3091-4735
<BR>Fax: +55 (0-XX-11) 3091-4714
<BR>
<BR></P>
</FONT>
</BODY>
</HTML>