[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 -- address@hidden
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