<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'>Hi Robert,<div><br></div><div>You can define a variable on the 25 km grid that contains the lapse rate you want for each grid point. Since most points will have the 6.5 lapse rate, define a constant variable:</div><div><br></div><div>'lapse=const(var,0.0065,-a)'</div><div><br></div><div>where 'var' is a variable already defined on the 25 km grid. Then use 'set defval' to change the value of lapse at selected points to 0.0045 (see <a href="http://cola.gmu.edu/grads/gadoc/gradcomdsetdefval.html">defval</a> documentation). You'll have to do this one grid point at a time, since it's a custom mask that can't be calculated from any other variables you have. Once this is done, you can use the lapse variable to calculate the temperature adjustment, something like</div><div><br></div><div>'adjust=lapse*(hs2-hs1)'</div><div><br></div><div>where hs1 and hs2 are the surface topography of the two datasets, in meters, interpolated to the 25 km grid. The way you've calculated the temperature difference, hs1 should be the NARR topography and hs2 the HadRM topography so the adjustment has the correct sign.</div><div><br></div><div>Best regards,</div><div><br><br><div><span name="x"></span>Eric L. Altshuler<br>Research Scientist<br>Center for Ocean-Land-Atmosphere Studies<br>253 Research Hall, Mail Stop 6C5<div>George Mason University</div><div>4400 University Drive<br>Fairfax, VA 22030 USA<br><br>E-mail: ela@cola.iges.org<br>Phone: (703) 993-5725<br>Fax: (703) 993-5770</div><span name="x"></span><br></div><br><hr id="zwchr"><b>From: </b>"Roberto Mera" <RMera@ucsusa.org><br><b>To: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br><b>Sent: </b>Tuesday, February 18, 2014 1:15:20 PM<br><b>Subject: </b>Re: [gradsusr] Masking out grid boxes for elevation correction<br><br>
<style><!--
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
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
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Sure thing.</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">I’m putting output from HadRM and NARR onto a common 25 km grid.</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">I will only be adjusting temperature such that once the adjustment has been calculated, I can do the following:
</span><span style="font-size:11.0pt;font-family:"Courier New"">Tavg_dif = Tavg_model - Tavg_narr + adjust</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">So, what I need is to be able to make a mask with the adjustment such that it gets computed using 6.5 c/km lapse rate except for over the grid cells over the Cascades and Sierra as
I specified below. Over the mountains the lapse rate is 4.5 c/km.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">I want to know how I can mask out a couple of boxes of grid points over the mountains. In the past I have been able to mask out the ocean, for example. But in this work it’s not necessarily
a mask over a certain height (say 2000 meters). Doing so misses most of the Cascades and adds other areas in the high desert that we prefer to leave at 6.5 c/km. Instead we want to do it specifically for the mountain ranges.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">Thanks!</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">Robert</span><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"> </span></p>
<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""> gradsusr-bounces@gradsusr.org [mailto:gradsusr-bounces@gradsusr.org]
<b>On Behalf Of </b>Jeff Duda<br>
<b>Sent:</b> Tuesday, February 18, 2014 12:22 PM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] Masking out grid boxes for elevation correction</span></p>
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal">Robert,</p>
<div>
<p class="MsoNormal">Can you be a little more precise in what you want to do? Which two grids are you putting onto a common 25 km grid? Are you trying to adjust just the temperature?</p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Jeff Duda</p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"> </p>
<div>
<p class="MsoNormal">On Tue, Feb 18, 2014 at 8:20 AM, Roberto Mera <<a href="mailto:RMera@ucsusa.org" target="_blank">RMera@ucsusa.org</a>> wrote:</p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">All,</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p>I’m trying to apply an elevation correction to model data by adjusting the lapse rate based on NARR. Basically, I’m looking to
<span style="font-family:"Courier New"">use a rate of 4.5 C/km for mountains, otherwise 6.5. where the adjustment is the rate the times the elevation difference.
</span></p>
<p><span style="font-family:"Courier New""> </span></p>
<p><span style="font-family:"Courier New"">The location where I’d like to apply the masks are</span> lon -122.5 -121 lat 43.5 45.2 for the Oregon Cascades and lat 46 49 for Washington Cascades and lon -121 -117.9, lat 35.5 - 40.5 for the Sierra.</p>
<p> </p>
<p>I’m placing the two grids on a common 25km grid and taking the difference between elevation.</p>
<p> </p>
<p>I appreciate the help.</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Robert</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a></p>
</div>
<p class="MsoNormal"><br>
<br clear="all">
</p>
<div>
<p class="MsoNormal"> </p>
</div>
<p class="MsoNormal">-- <br>
Jeff Duda<br>
Graduate research assistant<br>
University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms</p>
</div>
</div>
<br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></div></body></html>