This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
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