if you can plot via gempak but not garp that likely means that the
fields or not being written properly to file, or that garp can not read
them... if you find yourself stuck, feel free to contact me off list
again, and i will walk you through a fix.
cheers,
--patrick
------ Original Message ------
From: "Douglas Miller" <dmiller@xxxxxxxx>
To: gembud@xxxxxxxxxxxxxxxx
Sent: 3/11/2022 2:37:05 PM
Subject: [gembud] creating ERA5-based moisture maps, Part II: Composites
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
