[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[GEMPAK #APZ-498322]: Resend of gempak question with script attached



Hi Bob,

I had to make a number of edits to your script to get it to work.  Some of 
these may be peculiarities of your scripting environment and WRF grids, but to 
me they look like input errors.

1) $RESPOND is a hidden variable in GEMPAK and is not listed in programs like 
GDPROF.  If you set \$RESPOND = no and run a script, it is saved as "no" even 
after exiting the program and calling GPEND.  running another GEMPAK program 
will use this definition.  Therefore, to unset $RESPOND = no you must enter 
\$RESPOND=yes at least once and run it.  I have confirmed that this will is the 
case with your script.  That's the big problem solved, but there are other 
smaller issues.

2) device  = gif | ${date}${cycle}-${fh}.gif | 6;10
This directs GPLT to create a gif image size 6x10 pixels.  I could not get GPLT 
to create this image likely because the dimensions are below some minimum.  I 
had to remove size definition (DEVICE = gif | filename.gif) but also checked 
that (DEVICE = gif | filename.gif | 600;1000 ) worked.

3) I had to change XAXIS to "-20/40/5" using a recent GFS run to confirm both 
TEMP and DEWP are plotted vertically, which they are.

I have attached my example script.  Let me know if you have more questions.

Michael James
Unidata


> Michael:
> Thanks for the suggestion.  Unfortunately, it didn't
> make a difference removing the RESPOND=no.  I
> still only get tmpc and winds plotted.  We usually
> run WRF out to 36 hours.  I will try to send a much
> smaller gempak file (0 to 3 hours) and a sample .gif
> image from my script sometime nest week.
> Thanks for looking at my problem with the script.
> Bob
> 
> 
> address@hidden> wrote:
> 
> > Hi Bob, I have a guess that \$RESPOND=no in your script is the reason.
> > After the first "run" with temperature GDPROF exits.
> >
> > Michael
> >
> >
> > > Support:
> > >
> > > Sorry, I forgot to attach the script in previous message.
> > > I am trying to develop a script that generates a series
> > > of soundings from a WRF file in gempak format at a
> > > user-defined location (specified by gpoint). When I
> > > run gdprof in the script, it plots wind flags and tmpc
> > > correctly, but I cannot figure out how to get the script
> > > to overlay dwpc on the plot. I know how to do this when
> > > I run gdprof in command mode, but it doesn't work in
> > > this script. I have attached the script that I am using.
> > >
> > > I run script called sndg by typing:
> > > ./sndg 20101216 00 43.35 -76.39
> > > where 20101216 is the date, 00 the cycle and 43.35
> > > and -76.39 are the coordinates of FZY (Fulton, NY).
> > > Any suggestions will be greatly appreciated.
> > > Thanks,
> > > Bob
> > >
> > >
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: APZ-498322
> > Department: Support GEMPAK
> > Priority: Normal
> > Status: Open
> >
> >
> 
> 


Ticket Details
===================
Ticket ID: APZ-498322
Department: Support GEMPAK
Priority: Normal
Status: Open
#!/bin/csh
source /machine/gempak/NAWIPS/Gemenviron

$OS_BIN/gdprof << EOF
\$RESPOND=yes
gpoint  = 50;-100
gdattim = f012
gvcord  = pres
gfunc   = tmpc
gvect   = wnd
gdfile  = $MODEL/gfs/2013082706_gfs003.gem
line    = 5/1/2/1
scale   =
border  = 1
ptype   = skewt
xaxis   = -20/40/5
yaxis   = 1000/500/50
wind    = brbk31/0.8/0.9
winpos  = 1
filter  = no
clear   = no
text    = 0.9
title   = 1/-2/
device  = gif | gdprof.gif
thtaln  = 4/2/1/260/350/10
thteln  = 5/3/1/250/350/10
mixrln  = 3/10/1/0.5/5/0.5
run

gfunc   = dwpc
line    = 1/2/2/1
run

exit
EOF

$OS_BIN/gpend