Re: [gembud] creating ERA5-based moisture maps, Part II: Composites

  • To: gembud@xxxxxxxxxxxxxxxx
  • Subject: Re: [gembud] creating ERA5-based moisture maps, Part II: Composites
  • From: Steve Decker <sgdecker@xxxxxxxxxxxxxxxxxx>
  • Date: Fri, 11 Mar 2022 15:28:51 -0500
  • Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=envsci.rutgers.edu; dmarc=pass action=none header.from=envsci.rutgers.edu; dkim=pass header.d=envsci.rutgers.edu; arc=none
  • Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector9901; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=775vxQxVcvkAm8TBx1uNKu1d/DZTNxyzxAakwaM2dz0=; b=m7FYk2awGBspSx+4upWWKQOrEc9Pla9+cWxHpEE4YRfE46GIGoaQ5/iBNffg5SGorx5ZeJ+j6DBSsQwSYDosn3eZ4tBc3eKkCjDlYVi8PH3q/+JfU3oZx9H+EGqxhquBINxTwMrWcSK+nO2G0kFRlnEs9c3bHirPTp0FryxDJ7ocWZ/q6Mnt4P/+XxNjmpQX+Fqb0bY28/qYUzdtLcCNetN0F/VwizAvo8sAcCTvCgjF+ZmuQRJ6sWe/CJLfv9ZzGexwTTQpjEUWL+VINn3FttCSCor+q567o6gTspi08126GyUVRJ76TfATLHKbixJptlTDsSWXCW6RtOWDAhOB3A==
  • Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=X6MZhqYoA8C7y5kIhWHWSTCwaj7b6nSioXu91K//HknVou2mwrSdfgXO/bGDaObiqCfzv8x/UNyO00/WUE+BgCq0/OLGyYZr3Nt7x5KFZbbVSHEa7ouF7PcgAb8vk0nQsFvBfM8Cnd3YR4RBbFeioDRjaqdzfrJTrrdN4JKqrDRYfgmqT4h14hHqOCZ/lLFWHxX9HzGoiEApMJ0XMen33z3wIb+72fehFNGX8H3nPhiyDDBXkzBQ9E4pXhGGfp4YQ0pzhHXqWnW9bf/Ya2YXHBqzJlJ7hEHxbuWzqGZ02Yam0Ur30h2PBmP1K8l7NeFOoAhXpbPjd1pevyMrd8HB7w==
  • Authentication-results: dkim=none (message not signed) header.d=none;dmarc=none action=none header.from=envsci.rutgers.edu;
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


  • 2022 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the gembud archives: