<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style>
<!--
@font-face
{font-family:Calibri}
@font-face
{font-family:Tahoma}
@font-face
{font-family:"Comic Sans MS"}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif"}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline}
p.Computer, li.Computer, div.Computer
{margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Courier New";
color:black}
span.EmailStyle18
{font-family:"Calibri","sans-serif";
color:#1F497D}
.MsoChpDefault
{font-family:"Calibri","sans-serif"}
@page WordSection1
{margin:1.0in 1.0in 1.0in 1.0in}
div.WordSection1
{}
-->
</style>
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif"; color:#1F497D">If you mean you’d like to plot the output from OBSGRID (as opposed to WRFDA) then this script will make a file with location extracted from a qc_obs_used*
file that can be written by OBSGRID. You could plot it with a symbol scaled to your radius. It was re-projected from lon,lat to LCC coordinates using cs2cs from the PROJ package.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif"; color:#1F497D"> </span></p>
<p class="Computer"><span lang="FR">$ cat obsgrid_locations.csh</span></p>
<p class="Computer"><span lang="FR">#!/bin/csh -f</span></p>
<p class="Computer"><span lang="FR">if ($#argv == 0 | "$1" =~ -*h*) then</span></p>
<p class="Computer"><span lang="FR"> echo "Usage: $0:t path/to/geo_em_d03.nc path/to/qc_obs_used_file"</span></p>
<p class="Computer"><span lang="FR"> exit</span></p>
<p class="Computer"><span lang="FR">endif</span></p>
<p class="Computer"><span lang="FR"> </span></p>
<p class="Computer"><span lang="FR">set domain = $1</span></p>
<p class="Computer"><span lang="FR">set inp = $2</span></p>
<p class="Computer"><span lang="FR"> </span></p>
<p class="Computer"><span lang="FR">set sites = (`grep "[A-Za-z]" $inp | cut -c41-46 | sort | uniq`)</span></p>
<p class="Computer"><span lang="FR"> </span></p>
<p class="Computer"><span lang="FR">set ref_lon = `ncdump -h $domain | grep STAND_LON | cut -d= -f2 | cut -df -f1`</span></p>
<p class="Computer"><span lang="FR">set ref_lat = `ncdump -h $domain | grep MOAD_CEN_LAT | cut -d= -f2 | cut -df -f1`</span></p>
<p class="Computer"><span lang="FR">set truelat1 = `ncdump -h $domain | grep TRUELAT1 | cut -d= -f2 | cut -df -f1`</span></p>
<p class="Computer"><span lang="FR">set truelat2 = `ncdump -h $domain | grep TRUELAT2 | cut -d= -f2 | cut -df -f1`</span></p>
<p class="Computer"><span lang="FR"> </span></p>
<p class="Computer"><span lang="FR">echo "Lon, Lat, Xlcc, Ylcc, Site"</span></p>
<p class="Computer"><span lang="FR">foreach site ($sites)</span></p>
<p class="Computer"><span lang="FR"> set lat = `grep $site $inp | head -1 | cut -c1-20`</span></p>
<p class="Computer"><span lang="FR"> set lon = `grep $site $inp | head -1 | cut -c21-40`</span></p>
<p class="Computer"><span lang="FR"> set tmp = (`echo "$lat $lon" | proj -r -f '"%.5f"' +proj=lcc +lat_1=$truelat1 +lat_2=$truelat2 +lat_0=$ref_lat +lon_0=$ref_lon +ellip=sphere +a=6370000 +b=6370000 +units="km" | tr -d '"'`)</span></p>
<p class="Computer"><span lang="FR"> set Xlcc = $tmp[1]</span></p>
<p class="Computer"><span lang="FR"> set Ylcc = $tmp[2] </span></p>
<p class="Computer"><span lang="FR"> </span></p>
<p class="Computer"><span lang="FR"> echo "$lon, $lat, $Xlcc, $Ylcc, $site" </span>
</p>
<p class="Computer"><span lang="FR">end</span><span lang="FR"></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif"; color:#1F497D"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif"; color:#1F497D">Bart Brashers</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif"; color:#1F497D"> </span></p>
<div style="border:none; border-left:solid blue 1.5pt; padding:0in 0in 0in 4.0pt">
<div>
<div style="border:none; border-top:solid #B5C4DF 1.0pt; padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt; font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt; font-family:"Tahoma","sans-serif""> wrf-users-bounces@ucar.edu [mailto:wrf-users-bounces@ucar.edu]
<b>On Behalf Of </b>Lyndon Mark Olaguera<br>
<b>Sent:</b> Wednesday, October 16, 2013 7:15 AM<br>
<b>To:</b> wrfhelp; wrf-users<br>
<b>Subject:</b> [Wrf-users] Plotting a gridpoint with a the radius of influence</span></p>
</div>
</div>
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal">Hi all,</p>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">I'd like to know how to plot the radius of influence for a grid point for the objective analysis in WRF. Any suggestion on how to do this?</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Ill appreciate any help.Thanks</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal"><br clear="all">
</p>
<div>
<div>
<p class="MsoNormal"><b><i><u><span style="font-family:"Comic Sans MS"">Lyndon Mark P. Olaguera</span></u></i></b></p>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<hr size="1" color="#004250">
This message contains information that may be confidential, privileged or otherwise protected by law from disclosure. It is intended for the exclusive use of the Addressee(s). Unless you are the addressee or authorized agent of the addressee, you may not review,
copy, distribute or disclose to anyone the message or any information contained within. If you have received this message in error, please contact the sender by electronic reply to email@environcorp.com and immediately delete all copies of the message.
</body>
</html>