[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
- 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/ <
**************************************************************************** <