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.
>From: Brett Newkirk <address@hidden> >Organization: . >Keywords: 199903162219.PAA23320 >Hi- > I was wondering if it was possible to do subtraction of two fields >in 2 different GEMPAK files using GDCROSS. I have a cross section defined >and two separate model runs a.gem and b.gem. I would like to subtract the >potential temperature in the cross sectional plane from each other. I have >tried using GDFILE=a.gem + b.gem and doing GFUNC=thta//one (save it in a >grid called one) then doing GFUNC=SUB(one,thta) with GDFILE=b.gem + a.gem >and this did not work. Is there a simpler way of doing this task from >within GEMPAK? > Thanks for your help-- >Brett Newkirk > >Brett Newkirk E-MAIL: address@hidden >Office: ATG 424 Atmospheric Science/Geophysics Building >Mailing Address: Department of Atmospheric Sciences, University of Washington > Box 351640 > Seattle, WA 98195-1640 >Office Phone: (206) 685-2183 > Brett, To subtract two fields from different grid files, use the + operator to list both files in gdfile. Then, assume I want to compare the difference in the F012 tmpc forecast for the 00Z eta run and the F000 tmpc forecast for the 12Z run. GDFILE = $HDS/99031612_eta211.gem+$HDS/99031600_eta211.gem GDATTIM = F000 GFUNC = sub(tmpc,tmpc+2^f012) in the gfunc, tmpc will use the f000 forecast time from the first file. In order to get the 12 hour forecast from the second file, you have to use the inline modifiers to specify F012 from the second file. Steve Chiswell