[
Date Prev][
Date Next][
Thread Prev][
Thread Next][
Date Index][
Thread Index]
Re: [GEMPAK #YCE-717567]: Question on GDPLOT (well, 2)
Chiz,
Thanks for the info. I'll add _gif support (and improve our processing
somewhat) tomorrow!
We didn't include the map in that display because we're using a
web-mapping package (mapserver) to display the result. This will be an
intermediate product 'til I'm happy with the netcdf data I'm producing
from the MM5 runs, and start using the OGC Web Coverage Service for
similar work. We've got maps of the region in the mapserver dataset, so
we would prefer to use that as the two map datasets don't always
correspond when we've tried something like this in the past.
Our grid is an LCC 60/-101/30 grid with corners (via GAREA) as:
GAREA 23.84;-106.81;37.36;-84.1
I'd have thought that we'd have some blank (or black) data in this
region but the image we're seeing isn't consistent with the reprojection
geometry I'd expected for a CED result, so either I'm coffee-deprived or
confused. I'll look at interpolating using GDDIAG. Thanks.
The scripts we're using and the .nts files are theoretically attached to
this one. If they don't come thru, look on http://mesonet.tamu.edu/Ver/
for the files.
Thanks!
See you at AMS.
gerry
Unidata GEMPAK Support wrote:
Gerry,
I didn't get your .nts file attachments.
In answer to your question 1), there isn't a target in the gdplot2 Makefile to
build a gdplot2_gif,
but, it is very simple to add this into the Makefile. I have attached an update
for
$GEMPAK/source/prograns/gd/gdplot2/Makefile
Replace the current gdplot2 Makefile, and then build with:
cd $GEMPAK/source/prograns/gd/gdplot2/Makefile
make programs_gif
make clean
2) Your graphic doesn't have a map, so I'm not sure what you are referring to.
To define the display of a gridded data set using the CED projection, set the parameter
PROJ=ced instead of DEF or whatever other projection you are using.
If you want to interpolate your grid from a different projection to CED points,
you can do that using GDDIAG, such that a CED display of the grid will have the
interpolated
point locations on straight rows/columns.
Note that anytime you display a grid in something other than its native
projection, the
rows and columns will not necessarily look straight!
If you want to send me your .nts input, I'll see what you are doing,
Steve Chiswell
Unidata User Support
Ticket Details
===================
Ticket ID: YCE-717567
Department: Support GEMPAK
Priority: Normal
Status: On Hold
--
Gerry Creager -- gerry.creager@xxxxxxxx
Texas Mesonet -- AATLT, Texas A&M University
Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
GDFILE grid2_ps.gem
GDATTIM LAST
GLEVEL 0
GVCORD none
PANEL 0
SKIP /3;3
SCALE 0
GFUNC
CTYPE
CONTUR
CINT
LINE
FINT
FLINE
HILO
HLSYM
CLRBAR
GVECT vecr(u10,v10)
WIND bk32/0.8/1.5/0212
REFVEC
TITLE
TEXT
CLEAR no
GAREA 23.84;-106.81;37.36;-84.1
IJSKIP
PROJ ced
MAP
LATLON
DEVICE
STNPLT
GDFILE grid2_ps.gem
GDATTIM LAST
GLEVEL 0
GVCORD none
PANEL 0
SKIP /3;3
SCALE 0
GFUNC dbzx
CTYPE f
CONTUR 1/1
CINT
LINE
FINT 15;20;25;30;35;40;45;50;55;60;65;70;75
FLINE 0;4;3;22;23;5;19;17;2;15;14;29;30;1
HILO
HLSYM
CLRBAR
GVECT
WIND
REFVEC
TITLE
TEXT
CLEAR yes
GAREA 23.84;-106.81;37.36;-84.1
IJSKIP
PROJ ced
MAP 0
$MAPFIL
LATLON 0
DEVICE
STNPLT
GDFILE grid2_ps.gem
GDATTIM LAST
GLEVEL 0
GVCORD none
PANEL 0
SKIP /3;3
SCALE 0
GFUNC add(32,mul(1.8,sub(t2,273.15)))
CTYPE f
CONTUR 1/1
CINT
LINE
FINT -10;-5;0;5;10;15;20;25;30;35;40;45;50;55;60;65;70;75;80;85;90;95;100;105
FLINE 11;12;30;29;28;27;26;25;4;24;23;22;3;5;20;19;10;9;2;16;15;14;30;29;1
HILO
HLSYM
CLRBAR
GVECT vecr(u10,v10)
WIND bk32/0.8/1.5/0212
REFVEC
TITLE
TEXT
CLEAR yes
GAREA 23.84;-106.81;37.36;-84.1
IJSKIP
PROJ ced
MAP 0
$MAPFIL
LATLON 0
DEVICE
STNPLT
Attachment:
ver.csh
Description: C-Shell script
#! /bin/bash
#
#Can be run from cron
#
#If run from cron, make sure an X server is available
#source /data/gempak/hds/
#
file=$1
hour=$2
gdir=$3
# \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
# SOUTHEAST US PRECIPITATION MAP SCRIPT
# ///////////////////////////////////////////////////////
if [ "$file" = "eta" ] ; then
mod=na
restore_var1=southeast_us_prcp
restore_var2=southeast_us_slp
size="900;750"
fi
if [ "$file" = "avn211" ] ; then
mod=gf
restore_var1=southeast_us_prcp
restore_var2=southeast_us_slp
size="900;750"
fi
if [ "$file" = "ruc" ] ; then
mod=ru
restore_var1=southeast_us_prcp
restore_var2=southeast_us_slp
size="900;750"
fi
if [ "$file" = "grid2_ps.gem" ] ; then
mod=ta
restore_var1=southeast_us_prcp2_ver
size="900;800"
fi
rm -f $gdir/ver_${mod}sepr.$hour.gif
gdplot << GDPLOT_EOF
restore $restore_var1
GDFILE = $file
CONTUR = 1/1
GDATTIM =last
DEVICE = gif|$gdir/ver_${mod}sepr.$hour.gif|$size
run
exit
GDPLOT_EOF
gpend
#! /bin/bash
#
#Can be run from cron
#
#If run from cron, make sure an X server is available
#source /data/gempak/hds/
#
file=$1
hour=$2
gdir=$3
# \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
# SOUTHEAST US SURFACE TEMPERATURE/DEWPOINT MAP SCRIPT
# ////////////////////////////////////////////////////
if [ "$file" = "eta" ] ; then
mod=na
restore_var1=southeast_us_tmpf
restore_var2=southeast_us_dwpf
size="900;750"
fi
if [ "$file" = "avn211" ] ; then
mod=gf
restore_var1=southeast_us_tmpf
restore_var2=southeast_us_dwpf
size="900;750"
fi
if [ "$file" = "ruc" ] ; then
mod=ru
restore_var1=southeast_us_tmpf
restore_var2=southeast_us_dwpf
size="900;750"
fi
if [ "$file" = "grid2_ps.gem" ] ; then
mod=ta
restore_var1=southeast_us_tmpf2_ver
restore_var2=southeast_us_dwpf2_ver
size="900;800"
fi
rm -f $gdir/ver_${mod}sest.$hour.gif
gdplot << GDPLOT_EOF
restore $restore_var1
GDFILE = $file
CONTUR = 1/1
GDATTIM =last
DEVICE = gif|$gdir/ver_${mod}sest.$hour.gif|$size
run
exit
GDPLOT_EOF
gpend
