20050930: Gddiag - no joy

John,

To reference the variable gfile in your shell script,
you need to use "$gfile" within the script after you have set it.
Eg, the lines in particular:
   GDFILE = $gfile
and
   CPYFIL = $gfile

See if that helps.....

Chiz


----------------------------------------
Steve Chiswell
Unidata User Support




>From: John Merrill <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200509302055.j8UKt4G7014713

>Hello, Steve:
>   I finally got to trying the setup you crafted for making
>a tropopause diagnostic grid. I must be missing something,
>because no output file is created. Here's the script I
>used:
>
>#!/bin/csh
>#
>#  This is used to plot a cross section with PV, Theta and the
>#  tropopause height.
>#
>#  Argument 1 is the date and time, mmddhh (as in 042418).
>#
>source /ljuka/gemadm/gempak/Gemenviron
>#
>set gfile=/chili/jdubois/gempak.d/conduit.{$1}.gem
>#
>#  First, prepare file with tropopause height diagnostic.
>#
>@ LEV = 1000
>while ($LEV >= 50)
>   gddiag << EOF
>   GDFILE = gfile
>   GDOUTF = ./tropo.ruc
>   GFUNC = quo(pres@0%trop,pres)
>   GDATTIM = f000
>   GLEVEL = ${LEV}
>   GVCORD = pres
>   GRDNAM = TIDX@${LEV}%pres
>   GRDTYP = S
>   GPACK =
>   GRDHDR =
>   PROJ = grid
>   GRDAREA =
>   KXKY =
>   MAXGRD = 5000
>   CPYFIL = gfile
>   ANLYSS =
>r
>
>e
>EOF
>@ LEV = $LEV - 25
>end
>exit(0)
>
>Gddiag complains a lot that "only one output file is allowed," but I
>hoped this was only a warning. The script runs through the levels
>and exits, but no output fils is created. Thanks for suggestions.  John
>
--
NOTE: All email exchanges with Unidata User Support are recorded in the
Unidata inquiry tracking system and then made publicly available
through the web.  If you do not want to have your interactions made
available in this way, you must let us know in each email you send to us.