[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
19990304: how to store potential vorticity to an existing gempak file using gddiag or other gempak subroutines
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.
- Subject: 19990304: how to store potential vorticity to an existing gempak file using gddiag or other gempak subroutines
- Date: Thu, 04 Mar 1999 13:51:53 -0700
Shangyao, Your script is missing a blank line following the first "run" line. This might be leading to some confusion in the grdnam variable. There could also be a problem is the csh interpretation if there are any hidden control characters floating around in the script. Steve Chiswell Unidat User Support >From: Shangyao Nong <address@hidden> >Organization: . >Keywords: 199902081551.IAA13205 >To whom it may concern, > >I tried to use gddiag to add potential vorticity to one gempak data file. I us > ed the following >script >***************************** >#!/bin/csh ># SCRIPT TO create new grid file. This script adds summed up surface ># u-direction wind to an existed grid file. > >$GEMEXE/gddiag << endgddiag > > gdfile=95080800.rnl_gem > gdoutf=95080800.rnl_gem > gdattim=08/0000f00 > glevel=925:1000 > gvcord=pres > gfunc=mul(-9.8e4,mul(avor(vlav(obs), quo(ldf(thta),ldf(pres))) > grdnam=q > run > gdfile=95080800.rnl_gem > gdoutf=95080800.rnl_gem > glevel=850:925 > run > > exit > >endgddiag > >#$GEMEXE/gpend > ># gpend >**************************************** >There was no error message after running the above script. But when I used gdi > nfo to check >whether I already put Q in the data file, I found >*********************************** > 174 950808/0000F000 925 PRES VREL > 175 950808/0000F000 925 PRES HGHT >Enter <CR> to page or type EXIT: > 176 950808/0000F000 925 1000 PRES GDFILE > =9508 > 177 950808/0000F000 850 PRES TMPK > 178 950808/0000F000 850 PRES UREL > 179 950808/0000F000 850 PRES VREL > 180 950808/0000F000 850 PRES HGHT > 181 950808/0000F000 850 925 PRES Q > 182 950808/0000F000 700 PRES TMPK > 183 950808/0000F000 700 PRES UREL >*********************************** >I don't understand the above result, especially grid 176. Why does the variabl > e name become >GDFILE=9508 instead of Q? But grid 181 has correct variable name. > >When I tried to use gdcntr to plot grid 181, I couldn't. Why? I can calculate > potential >vorticity between 850 and 925. I need to know how to store it to an existing g > empak file. > >Thank you. > >Shangyao Nong >address@hidden > >From address@hidden Mon Feb 8 18:59:25 1999 >Received: from MIT.EDU (SOUTH-STATION-ANNEX.MIT.EDU [18.72.1.2]) > by unidata.ucar.edu (8.8.8/8.8.8) with SMTP id SAA00259 > for <address@hidden>; Mon, 8 Feb 1999 18:59:24 -0700 (MST) >Organization: . >Keywords: 199902090159.SAA00259 >Received: from ALL-NIGHT-TOOL.MIT.EDU by MIT.EDU with SMTP > id AA24752; Mon, 8 Feb 99 20:59:21 EST >Received: by all-night-tool.mit.edu (8.8.7/4.7) id UAA10153; Mon, 8 Feb 1999 2 > 0:59:23 -0500 (EST) >Message-Id: <address@hidden> >To: address@hidden >Subject: a problem with using gddiag >Date: Mon, 08 Feb 1999 20:59:23 EST >From: Shangyao Nong <address@hidden> > >To whom it may concern, > >I am using gddiag to add potential vorticity my gempak files. I used the >following script >*************************************************** >#!/bin/csh ># SCRIPT TO create new grid file. This script adds pvor to an existing ># gempak file. > >$GEMEXE/gddiag << endgddiag > > gdfile=95080800.rnl_gem > gdoutf=95080800.rnl_gem > gdattim=08/0000f00 > glevel=925:1000 > gvcord=pres > gfunc=mul(9.8e6,pvor(thta,obs)) > grdnam=PVOR@950 > run > > exit > >endgddiag > >#$GEMEXE/gpend > > gpend >************************************************* >But when I used gdcntr to plot pvor@950, what I got was just a constant. Here >is what I had >************************************************* >Creating process: xw for queue 1602 > > Grid file: 95080800.rnl_gem > > GRID IDENTIFIER: > TIME1 TIME2 LEVL1 LEVL2 VCORD PARM >950808/0000F000 950 PRES PVOR > > GAREA: 2.5;-140;60;-32.5 > SCALE: 0 > > MINIMUM AND MAXIMUM VALUES -28.25 -28.25 > > LINE CONTOURS: > > LEVELS: -28.25 > COLORS: 1 > LINTYP: 2 > LINWID: 1 > LABEL: 1 >Enter <cr> to accept parameters or type EXIT: > Parameters requested: GDATTIM,GLEVEL,GVCORD,GFUNC,GDFILE,CINT,LINE,MAP, > TITLE,DEVICE,SATFIL,RADFIL,PROJ,GAREA,CLEAR,PANEL,TEXT,SCALE,LATLON,HILO, > HLSYM,CLRBAR,CONTUR,SKIP,FINT,FLINE,CTYPE,LUTFIL,STNPLT. >*************************************************************** > >However, when I used gdcntr to plot pvor between 925 and 1000 directly, >I had the correct pvor distribution. Can you tell me what is wrong with my >script? It seems to me my script is fine. > >Thanks. > >Shangyao Nong >address@hidden > >From address@hidden Tue Feb 9 10:11:34 1999 >Received: from MIT.EDU (PACIFIC-CARRIER-ANNEX.MIT.EDU [18.69.0.28]) > by unidata.ucar.edu (8.8.8/8.8.8) with SMTP id KAA20064 > for <address@hidden>; Tue, 9 Feb 1999 10:11:34 -0700 (MST) >Organization: . >Keywords: 199902091711.KAA20064 >Received: from ALL-NIGHT-TOOL.MIT.EDU by MIT.EDU with SMTP > id AA24514; Tue, 9 Feb 99 12:11:57 EST >Received: by all-night-tool.mit.edu (8.8.7/4.7) id MAA11480; Tue, 9 Feb 1999 1 > 2:11:32 -0500 (EST) >Message-Id: <address@hidden> >To: address@hidden >Subject: a problem with using gddiag >Date: Tue, 09 Feb 1999 12:11:32 EST >From: Shangyao Nong <address@hidden> > >To whom it may concern, > >I am using gddiag to add potential vorticity my gempak files. I used the >following script >*************************************************** >#!/bin/csh ># SCRIPT TO create new grid file. This script adds pvor to an existing ># gempak file. > >$GEMEXE/gddiag << endgddiag > > gdfile=95080800.rnl_gem > gdoutf=95080800.rnl_gem > gdattim=08/0000f00 > glevel=925:1000 > gvcord=pres > gfunc=mul(9.8e6,pvor(thta,obs)) > grdnam=PVOR@950 > run > > exit > >endgddiag > >#$GEMEXE/gpend > > gpend >************************************************* >But when I used gdcntr to plot pvor@950, what I got was just a constant. Here >is what I had >************************************************* >Creating process: xw for queue 1602 > > Grid file: 95080800.rnl_gem > > GRID IDENTIFIER: > TIME1 TIME2 LEVL1 LEVL2 VCORD PARM >950808/0000F000 950 PRES PVOR > > GAREA: 2.5;-140;60;-32.5 > SCALE: 0 > > MINIMUM AND MAXIMUM VALUES -28.25 -28.25 > > LINE CONTOURS: > > LEVELS: -28.25 > COLORS: 1 > LINTYP: 2 > LINWID: 1 > LABEL: 1 >Enter <cr> to accept parameters or type EXIT: > Parameters requested: GDATTIM,GLEVEL,GVCORD,GFUNC,GDFILE,CINT,LINE,MAP, > TITLE,DEVICE,SATFIL,RADFIL,PROJ,GAREA,CLEAR,PANEL,TEXT,SCALE,LATLON,HILO, > HLSYM,CLRBAR,CONTUR,SKIP,FINT,FLINE,CTYPE,LUTFIL,STNPLT. >*************************************************************** > >However, when I used gdcntr to plot pvor between 925 and 1000 directly, >I had the correct pvor distribution. Can you tell me what is wrong with my >script? It seems to me my script is fine. > >Thanks. > >Shangyao Nong >address@hidden > **************************************************************************** < Unidata User Support UCAR Unidata Program < (303)497-8644 P.O. Box 3000 < address@hidden Boulder, CO 80307 < ---------------------------------------------------------------------------- < Unidata WWW Service http://www.unidata.ucar.edu/ < **************************************************************************** <
- Prev by Date: 19990303: Garp
- Next by Date: 19990305: boundary layer winds and GARP
- Previous by thread: 19990303: Garp
- Next by thread: 19990305: boundary layer winds and GARP
- Index(es):