Doug,
One idea that comes to mind is to override the times with some fixed
placeholder in your GRDNAM settings. So, instead of
GRDNAM = TMP3
have
GRDNAM = TMP3^200101/0000F000
where the date/time is arbitrary but the same for each intermediate result.
No guarantees, but it is worth a shot.
Sincerely,
Steve
On 3/11/22 2:37 PM, Douglas Miller wrote:
Thank you to all who helped steer me in a fruitful direction
with regards to building ERA5-based GEMPAK format files!
I'm now trying to create horizontal and vertical sections of
ERA5-based composites and am close, but yet so far, to being able to
display them. Each composite consists of seven members so, again, I
use GDDIAG to add three time periods at a time (ending with the
addition of the seventh time period) to create three sets of summed
grid fields.
==> Challenge #1
Each time I add fields across three different time periods, I get
a non-blank time in the Time2 designator (find attachment
"...twotimes2.png"). The problem/challenge happens later in the script
when I try to add the newly created fields at the end. I haven't
figured out how to properly designate the time in the script where I'm
adding the sum of the first three times (TMP3), the sum of the second
three times (TMP6), and the TMPK at the seventh time period. Manually,
I can plot the isobaric summed temperature via "garp" using
"TMP3^030405/0900F000:131222/0700F000" in the command line (attachment
"...twotimes2.png"), but defining a time in this fashion via the
script gives the error message...
[DG -7] Input grid TMP3^030405/0900F000:131222/0700F000 @1000 %PRES
in ADD cannot be found.
==> Challenge #2
As mentioned above, I can plot the isobaric summed temperature
maps with no problem. However, when I try to manually plot vertical
cross sections via "garp," I get the GDCROSS "No levels..." error
message (find attachment "...twotimes3.png"). I have successfully
plotted vertical cross sections pointing to each of the individual
ERA5 GEMPAK files, so the problem seems to happen when I create these
composite (summed) fields, even though I can see the fields exist on
their isobaric levels via the FileBrowser utility in "garp."
Any insights on how to generate composites would be most appreciated!
Regards,
Doug
##################################################
A sample block of the script...
set grido=CAUSo_era5.gem
set gdat1 = "030405/0900F000"
set gdat2 = "080304/2200F000"
set gdat3 = "131222/0700F000"
.
set gdat8 = "010101/0100F000"
.
.
.
#--------------------------------------------------------------------------
echo "creating composite files based on ERA5"
$GEMEXE/gddiag <<EOF>> $LOGFILE
GDATTIM = $gdat8
LEVELS = ALL
MAXGRD = 400
GVCORD = pres
GRDAREA = 24;-130;51;-65
GDOUTF = $grido
KXKY = 261;109
GDFILE = 2003040509_era5.gem + 2008030422_era5.gem + 2013122207_era5.gem
GLEVEL = 1000
GFUNC = ADD( ADD(TMPK+1^$gdat1,TMPK+2^$gdat2), TMPK+3^$gdat3)
GRDNAM = TMP3
r
GFUNC = ADD( ADD(RELH+1^$gdat1,RELH+2^$gdat2), RELH+3^$gdat3)
GRDNAM = REL3
r
GFUNC = ADD( ADD(HGHT+1^$gdat1,HGHT+2^$gdat2), HGHT+3^$gdat3)
GRDNAM = HGH3
r
##################################################
Douglas K. Miller
Professor
Atmospheric Sciences Dept.
UNC-Asheville
RBH 236A, ATMS
One University Heights
Asheville, NC 28804
Phone: (828) 232-5158
FAX: (828) 232-5046
_______________________________________________
NOTE: All exchanges posted to Unidata maintained email lists are
recorded in the Unidata inquiry tracking system and made publicly
available through the web. Users who post to any of the lists we
maintain are reminded to remove any personal information that they
do not want to be made public.
gembud mailing list
gembud@xxxxxxxxxxxxxxxx
For list information or to unsubscribe, visit:
https://www.unidata.ucar.edu/mailing_lists/
--
Steve Decker, Associate Teaching Professor
Director, Meteorology Undergraduate Program
Department of Environmental Sciences Phone: (848) 932-5750
Rutgers University Fax: (732) 932-8644
14 College Farm Rd Email: decker@xxxxxxxxxxxxxxxxxx
New Brunswick, NJ 08901-8551
