[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[GEMPAK #QYG-706948]: gdfile- streamlines

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