NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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 #1Each 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 #2As 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
gembud archives: