[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #QYG-706948]: gdfile- streamlines
- Subject: [GEMPAK #QYG-706948]: gdfile- streamlines
- Date: Wed, 09 Sep 2009 14:45:55 -0600
Keshab,
The GFS grid #225 is Hawaii 80km... you may not be ingesting/decoding it
anymore. If you have grid 211 available (CONUS 80km) you might try simply
using that one and see if this solves your problem.
-Michael
> We have a script to plot streamlines over US. It used to work fine but
> for some reasons it has stopped working. I suspect that something is
> missing or outdated. Here is an excerpt from the script and I am not
> able to locate gfs225.gem. There is a jump in the sequence from
> gfs213.gem to gfs254.gem. Could this be the reason?
>
> # Build all the variables that use times. Set the analysis to be that of
> # the ETA.
>
> set filestamp=`date -u "+%Y%m%d"`
> set gdfile="$GEMDATA/model/nam/${filestamp}${fhour}_icwf212.gem"
>
> set gdfile_ak="$GEMDATA/model/nam-ak/${filestamp}${fhour}_nam217.gem"
>
> set gdfile_hi="$GEMDATA/model/gfs/${filestamp}${fhour}_gfs225.gem"
>
> set datelabel=`date -u "+%h %d %Y"`
>
> #--------------------------------------------------------------------
> #
> # Set the output plot file
> #
> #--------------------------------------------------------------------
>
> # Set the GIF file names. Real-time and archive derived files go in
> # different places to avoid trampling on each other.
>
> set tmpgif=streamlines.gif.tmp
> set finalgif=$tranDIR/streamlines.gif
> set tmpgif1=streamlines_ak.gif.tmp
> set finalgif1=$tranDIR/streamlines_ak.gif
> set tmpgif2=streamlines_hi.gif.tmp
> set finalgif2=$tranDIR/streamlines_hi.gif
>
> #--------------------------------------------------------------------
> #
> # This section is boring. Skip reading it if you don't have to.
> #
> #--------------------------------------------------------------------
>
> # Copy in the white color table
>
> cp $utilDIR/coltbl.xwp coltbl.xwp
>
> # Remove the temporary gif file if it exists
>
> if ( -e $tmpgif ) rm $tmpgif
>
> #--------------------------------------------------------------------
> # Make sure the files we need exist, else we'll have problems
> #--------------------------------------------------------------------
>
> @ count = 0
> while (1)
> if (!(-e $gdfile) || !(-e $gdfile_ak) || !(-e $gdfile_hi)) then
> echo "$gdfile, $gdfile_ak or $gdfile_hi does not exist; waiting"
> sleep 600
> @ count = ($count + 1)
> else
> break
> endif
>
> if ($count >= 10) then
> echo "Files cannot be found, exiting"
> exit
> endif
> end
>
> #--------------------------------------------------------------------
> # Run the gempak program(s)
> #--------------------------------------------------------------------
> # First, plot the pretty streamlines, logo and box
> # over the good 'ol US of A
> #--------------------------------------------------------------------
>
> gdstream<< GDSTREAM_FLAG
> stnplt=
> satfil=
> radfil=
> lutfil=
> stream=.5
> gdattim=$gdattim
> glevel=10
> gvcord=hght
> gvect=wnd
> gdfile=$gdfile
> wind=////.8
> line=2/1/1/2/2
> title=32/-1/${hour}00Z $datelabel Surface Streamlines
> clear=yes
> panel=0
> device=gif|$tmpgif|800;600
> area=dset
> garea=23.0;-120.0;47.0;-65.0
> text=1/2/hw
> proj=lcc/25;-103;60
> map=32
> latlon=0
>
> run
>
> exit
>
> GDSTREAM_FLAG
>
> Thanks
> Keshab
>
>
>
Ticket Details
===================
Ticket ID: QYG-706948
Department: Support GEMPAK
Priority: High
Status: Open