NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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/
gembud
archives: