[
Date Prev][
Date Next][
Thread Prev][
Thread Next][
Date Index][
Thread Index]
20000731: no contours for lifted index
- Subject: 20000731: no contours for lifted index
- Date: Mon, 31 Jul 2000 13:01:50 -0600
Tom,
I mentioned gdinfo since it is a program to let you see all the grids in your
grid file. You do not need it in you script, but it is educational to see
the glevel and gvcord of the parameters that you will want to plot.
Presumably, you have run oabsnd to grid your upperair observations to
a gempak gridded data file. The STNDEX parameter in oabsnd specified
which stability indicies you will compute and grid into your data set.
The stability indicies (or whole sounding variables) are generally written
using GLEVEL=0 and GVCORD=none since they are for the entire sounding, and
not a specific level.
Running GDINFO on your gridded data set will show you all of the parameters,
levels and times you have in your grid file. If you want to run gdcntr with
cape, swet, totl etc, then you need to have gridded those with oabsnd.
If you aren't getting the output you expect with gdcntr, then look to make sure
that the parameters do exist in your grid file. You may need to go back
and look at your script which grids the data to make sure you are computing
all the indicies you want.
Steve Chiswell
>From: Tom Priddy <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200007311728.e6VHSQT08429
>From: Tom Priddy
>
>OK....I made those changes....but a little confused. I wasn't using GDINFO at
> all
>and
>not sure how to incorporate it in my script. Maybe I shouldn't be using snmap
> at
>all.
>I can get the lifted index to plot and the numbers look realistic....but after
> I
>made
>the changes I still only got a plot of LI. Do you have an example script
>available
>for LI, cape, etc...and if not what am I doing wrong. Maybe I'm missing somet
> hing
>
>concerning the whole script./ktp
>....new script below....
>#!/bin/csh -f
>#
># This is the program used to produce the lifted index contours. Parameters th
> at
># must be passed to it are the hour.
>#
># Approximate run time: 1 minute
>#
>cd /export/home/gempak/scripts/LDM
>source ~gempak/.cshrc
>source ~gempak/GEMsetup
>
>if ( $#argv == 1 ) then
> set hh = $argv[1]
> set level = 100
> set yy = `date -u +%y`
> set mm = `date -u +%m`
> set dd = `date -u +%d`
> set date1 = ${yy}${mm}${dd}
> set datetime = ${date1}/${hh}
>else
> echo ""
> echo USAGE: uppermap hour level
> echo ""
> exit
>endif
>
>set ldmdata = "/ldmdata"
>
># Path and file name convention for GEMPAK upper air data
>set upperpath = ${ldmdata}/gempak/upperair
>set snfile = ${upperpath}/${date1}_upa.gem
>
># Path and file name convention for GEMPAK upper air obj analysis grid file
>set gridpath = ${ldmdata}/gempak/grids
>set uprgrid = ${gridpath}/${date1}_${hh}Z_ua.gem
>
>set fname = "${date1}${hh}Z_${level}"
>set output = "${fname}.gif"
>set device_name = "gfw|${output}|792;612"
>
>rm -f *.nts ${output}
>
>snmap << EOF0
> AREA = garea
> GAREA = 15;-118;58;-55
> SNPARM = lift
> DATTIM = ${datetime}
> LEVELS = 0
> VCOORD = none
> SNFILE = ${snfile}
> COLORS = 4
> MAP = 32
> LATLON = 0
> MARKER = 0
> TITLE = 32/-3/U.K Ag. Weather Center
> CLEAR = YES
> PANEL =
> DEVICE = ${device_name}
> PROJ = lcc
> FILTER = YES
> TEXT = 1/21//hw
> RUN
>
> EXIT
>EOF0
>
>gdcntr << EOF1
> GDATTIM = ${datetime}
> GLEVEL = 0
> GVCORD = none
> GFUNC = lift
> GDFILE = ${uprgrid}
> CINT = -6;-4;-2;0;2;4;6
> LINE = 32/1/1/1
> MAP = 32
> TITLE = 32/-2/ Lifted Index Analysis
> DEVICE = ${device_name}
> PROJ = lcc
> GAREA = 15;-118;58;-55
> CLEAR = NO
> PANEL =
> TEXT = 1/21//hw
> SCALE =
> LATLON = 0
> CONTUR =
> SKIP =
> FINT = -6;-4;-2;0;2;4;6
> FLINE = 11;12;13;13;21;22;23;24
> CTYPE = f/c
> HILO =
> HLSYM =
> RUN
>
> EXIT
>EOF1
>
>gpend
>
>echo " "
>echo " "
>
>/export/home/gempak/scripts/add_logo /export/home/gempak/scripts/LDM/${fname}
>latestlifted logo.gif
>
>rm -f *.nts ${output}
>dalton%
>
>
>Unidata Support wrote:
>
>> Tom,
>>
>> When you grid LIFT from your sounding data, you will see
>> in GDINFO that the gvcord=NONE and glevel=0.
>> In the script fragment below, I don't see level defined, but
>> for gdcntr it should be 0 and you need to change the coordinate
>> from pres to none.
>>
>> This is because LIFT is a stability index (STNDEX), rather than a
>> variable for a particular level.
>>
>> Steve Chiswell
>> Unidata User Support
>>
>> >From: Tom Priddy <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200007311531.e6VFVBT04439
>>
>> >From: Tom Priddy/Erin Maxwell
>> >
>> >Steve....we can get the lifted index to plot but we can't seem to
>> >figure out how to get the data contoured and color-filled. Could
>> >you please take a look at the below listed gempak script and
>> >tell us what we are doing wrong? Thanks!
>> >
>> >
>> >
>> >
>> >
>> >set ldmdata = "/ldmdata"
>> >
>> ># Path and file name convention for GEMPAK upper air data
>> >set upperpath = ${ldmdata}/gempak/upperair
>> >set snfile = ${upperpath}/${date1}_upa.gem
>> >
>> ># Path and file name convention for GEMPAK upper air obj analysis grid
>> >file
>> >set gridpath = ${ldmdata}/gempak/grids
>> >set uprgrid = ${gridpath}/${date1}_${hh}Z_ua.gem
>> >
>> >set fname = "${date1}${hh}Z_${level}"
>> >set output = "${fname}.gif"
>> >set device_name = "gfw|${output}|792;612"
>> >
>> >rm -f *.nts ${output}
>> >
>> >snmap << EOF0
>> > AREA = garea
>> > GAREA = 15;-118;58;-55
>> > SNPARM = lift
>> > DATTIM = ${datetime}
>> > LEVELS = ${level}
>> > VCOORD = pres
>> > SNFILE = ${snfile}
>> > COLORS = 4
>> > MAP = 32
>> > LATLON = 0
>> > MARKER = 0
>> > TITLE = 32/-3/U.K Ag. Weather Center
>> > CLEAR = YES
>> > PANEL =
>> > DEVICE = ${device_name}
>> > PROJ = lcc
>> > FILTER = YES
>> > TEXT = 1/21//hw
>> > RUN
>> >
>> > EXIT
>> >EOF0
>> >
>> >gdcntr << EOF1
>> > GDATTIM = ${datetime}
>> > GLEVEL = ${level}
>> > GVCORD = pres
>> > GFUNC = lift
>> > GDFILE = ${uprgrid}
>> > CINT = -6;-4;-2;0;2;4;6
>> > LINE = 32/1/1/1
>> > MAP = 32
>> > TITLE = 32/-2/ Lifted Index Analysis
>> > DEVICE = ${device_name}
>> > PROJ = lcc
>> > GAREA = 15;-118;58;-55
>> > CLEAR = NO
>> > PANEL =
>> > TEXT = 1/21//hw
>> > SCALE =
>> > LATLON = 0
>> > CONTUR =
>> > SKIP =
>> > FINT = -6;-4;-2;0;2;4;6
>> > FLINE = 11;12;13;13;21;22;23;24
>> > CTYPE = f/c
>> > HILO =
>> > HLSYM =
>> > RUN
>> >
>> > EXIT
>> >EOF1
>> >
>> >gpend
>> >
>> >echo " "
>> >echo " "
>> >
>> >/export/home/gempak/scripts/add_logo
>> >/export/home/gempak/scripts/LDM/${fname} latestlifted logo.gif
>> >
>> >rm -f *.nts ${output}
>> >
>> >
>> >--
>> >===========================================================================
>> >
>> >Tom Priddy, Extension Ag. Meteorologist.....email: address@hidden
>> >Ag. Weather Center.............................fax: 606-257-5671
>> >243 Ag. Engineering Bldg.....................phone: 606-257-8803 ext 245
>> >
>> >Dept. of Biosystems and Ag. Engr.
>> >University of Kentucky
>> >Lexington, KY 40546-0276
>> >HOMEPAGE......................................http://wwwagwx.ca.uky.edu/
>> >
>> >===========================================================================
>> >
>> >
>> >
>>
>> ****************************************************************************
>> Unidata User Support UCAR Unidata Program
>> (303)497-8644 P.O. Box 3000
>> address@hidden Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service http://www.unidata.ucar.edu/
>> ****************************************************************************
>
>--
>===========================================================================
>Tom Priddy, Extension Ag. Meteorologist.....email: address@hidden
>Ag. Weather Center.............................fax: 606-257-5671
>243 Ag. Engineering Bldg.....................phone: 606-257-8803 ext 245
>Dept. of Biosystems and Ag. Engr.
>University of Kentucky
>Lexington, KY 40546-0276
>HOMEPAGE......................................http://wwwagwx.ca.uky.edu/
>===========================================================================
>
>
