Hi Jonathan,
You are correct in using GDRADR to create a radar composite. There are some
good examples of using GDRADR at
http://www.unidata.ucar.edu/software/gempak/examples/gdradr/
. Instead of trying to create a COMPOSITE image myself, I've found that it's
easier to use the pre-generated radar composite fed through the LDM feed.
This pre-generated data uses GDRADR which saves us a step. If you don't
have access to this data it can be accessed at:
http://motherlode.ucar.edu/cgi-bin/ldm/genweb?decoded/gempak/radar
http://motherlode.ucar.edu/cgi-bin/ldm/genweb?decoded/gempak/radar_idd
You can use GDPLOT2 (I think) and GDPLOT3 (definitely) to view these images.
You will have to tweak the area and projection to fit your WRF model domain.
I'm not sure of an easy way to do that, brute force trial and error will
eventually get what you need. In my case, I had to create a composite radar
image of the U.S. which matched exactly to an image generated elsewhere.
After a lot of trial and error tweaking "GAREA and PROJ" I was able to match
them up. Here's a sample of a simple shell script which generates a U.S.
composite radar image:
-----
#!/bin/csh
source /data/gempak/Gemenviron
gdplot3 << GDPLOT3_EOF
GDFILE=/home/ldm/data/gempak/radar/2010011215_rad1.gem
GDATTIM=100112/1507F000
GLEVEL=0
GVCORD=none
PANEL=0
SKIP=0
SCALE=0
GDPFUN=sge(n0r,15)
TYPE=z
CONTUR=0
CINT=0
LINE=0
FINT=15;20;25;30;35;40;45;50;55;60;65;70;75
FLINE=3;4;5;6;7;8;9;10;11;12;13;14;15;16
HILO=0
HLSYM=
CLRBAR=0
WIND=
REFVEC=0
TITLE=0
TEXT=
CLEAR=yes
GAREA=14.2;-120.3;47.3;-40
IJSKIP=0
PROJ=str/90;-96;0
MAP=0
MSCALE=
LATLON=
DEVICE=gif | /home/ldm/images/current/radar/us_201001121507_radar.gif |
960 ; 720
STNPLT=
SATFIL=
RADFIL=
IMCBAR=
LUTFIL=none
STREAM=
POSN=
COLORS=0
MARKER=
GRDLBL=
FILTER=
BOXLIN=
REGION=
TXTCOL=
TXTFIL=
TXTLOC=
COLUMN=
SHAPE=
INFO=
LOCI=
ANOTLN=
ANOTYP=
r
e
GDPLOT3_EOF
gpend
-----
Hope this helps you out...
Evan
Business Applications Meteorologist
Weather Trends International, Inc.
Phone 610-807-3197
http://www.wxtrends.com
This communication is privileged and may contain confidential information.
It's intended only for the use of the person or entity named above. If you
are not the intended recipient, do not distribute or copy this
communication. If you have received this communication in error, please
notify the sender immediately and return the original to the email address
above. © Copyright 2010 Weather Trends International, Inc.
-----Original Message-----
From: gembud-bounces@xxxxxxxxxxxxxxxx
[mailto:gembud-bounces@xxxxxxxxxxxxxxxx] On Behalf Of Case.Jon
Sent: Tuesday, January 12, 2010 10:03 AM
To: gembud@xxxxxxxxxxxxxxxx
Cc: Jonathan.Case-1@xxxxxxxx
Subject: [gembud] GDRADR Segmentation fault
Dear GEMPAK Users,
I'm attempting to create a radar mosaic/composite on my WRF model forecast
domain, and I have not yet been successful with this undertaking. At first,
I considered using nex2img or nex2gini to take Level III radar products from
NCDC and mosaic them. However, I would like to be able to create a
multi-panel display of forecast reflectivitiy and observed reflectivity in a
single GIF image.
So, I switched to GDRADR, which looks to be the program I want to use. So
far, I have 2 radar sites and it appears that I have set up all the
environmental variables correctly. I have the archived base reflectivity
data with the following naming convention:
$RAD/NIDS/%SID%/%PROD%/%PROD%_YYYYMMDD_HHMM
where %SID% is the station ID of the radar site, and %PROD% is the product
name, in my case "N0R".
When I try to run GDRADR, I consistently receive a Segmentation fault. I've
tried both setting up the grid parameters and setting a CPYFIL -- both the
same result. Does anyone have a recommendation for me?
My parameter settings are given below.
Many thanks,
Jonathan
lc2:/raid1/casejl/DATA/RADAR > gdradr
GRDAREA Area covered by grid 39.62;-91.83;49.57;-72.50
PROJ Map projection/angles/margins|dr lcc//6;1;1;2
KXKY Number of grid points in x;y 380;280
GDPFUN Scalar grid or vector grid funct N0R
GDFILE Grid file rad.gem
RADTIM Radar composite current/dattim 091209/1205
RADDUR Radar time window (minutes prior 10
RADFRQ Update Frequency
CPYFIL Grid file whose navigation is to
STNFIL Station information file nexrad.tbl
MAXGRD Maximum number of grids 9999
RADMODE Radar operational mode pc
NDVAL Replacement value for "ND" (defa
Parameters requested: GRDAREA,PROJ,KXKY,GDPFUN,GDFILE,RADTIM,RADDUR,RADFRQ,
CPYFIL,STNFIL,MAXGRD,RADMODE,NDVAL.
GEMPAK-GDRADR>r
Segmentation fault
________________________________
The information contained in this email message is intended only for the use
of the individual(s) to whom it is addressed and may contain information
that is privileged and sensitive. If you are not the intended recipient, or
otherwise have received this communication in error, please notify the
sender immediately by email at the above referenced address and note that
any further dissemination, distribution or copying of this communication is
strictly prohibited.
The U.S. Export Control Laws regulate the export and re-export of technology
originating in the United States. This includes the electronic transmission
of information and software to foreign countries and to certain foreign
nationals. Recipient agrees to abide by these laws and their regulations --
including the U.S. Department of Commerce Export Administration Regulations
and the U.S. Department of State International Traffic in Arms Regulations
-- and not to transfer, by electronic transmission or otherwise, any content
derived from this email to either a foreign national or a foreign
destination in violation of such laws.
_______________________________________________
gembud mailing list
gembud@xxxxxxxxxxxxxxxx
For list information or to unsubscribe, visit:
http://www.unidata.ucar.edu/mailing_lists/
