[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #APZ-498322]: Resend of gempak question with script attached
- Subject: [GEMPAK #APZ-498322]: Resend of gempak question with script attached
- Date: Tue, 27 Aug 2013 09:49:28 -0600
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