Bill,<br><br>Thanks for the suggestion. I attempted to set up a mask using the difference between hgtprs and hgtsfc. Unfortunately this produced similar "boxy" results like the images I attached initially. <br><br>
What routine does maskout use to "mask out" the negative areas? Is there a way for it to use the same smooth routine used when contouring?<br><br>Dan<br><br><div class="gmail_quote">On Wed, Jul 2, 2008 at 11:18 AM, Bill Bua <<a href="mailto:Bill.Bua@noaa.gov">Bill.Bua@noaa.gov</a>> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dan --<br>
<br>
I suspect it's because you're using lev, which is discretized into<br>
layers, and shading, resulting in 'boxy' plots to the nearest level set<br>
up in your control file and "masking out" anything that has a value of<br>
your variable t1pm1-lev < 0. The plot of t1pm1-lev is smooth because of<br>
the plotting routine used for contours.<br>
<br>
I'm not sure if this will work better, but perhaps you can use the model<br>
topographic height if it's available and subtract the height of the<br>
'lev' pressure surfaces as an alternative. Your mask might look more<br>
like the topography than the levels if you do so.<br>
<br>
Hope this helps.<br>
<br>
Bill Bua<br>
<br>
Dan Leins wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="Wj3C7c">
All,<br>
<br>
I am attempting to draw a simple cross section showing temperature<br>
across a model domain and I would like to maskout the area below<br>
ground. I have seen some traffic on this issue. I see that I should be<br>
able to define the terrain by subtracting pressfc/100 from lev, with<br>
all negative values being areas that are underground. Then using<br>
maskout, I can do the following:<br>
<br>
d maskout(tmpprs,t1pm1-lev)<br>
<br>
This gives me something in the ballpark of what I'm looking for<br>
(image1.png), but the shaded area looks very gridded in nature. When I<br>
simply plot a contoured plot of 't1pm1-lev', it doesn't look gridded<br>
at all and appears nice and smooth (image2.png). The model grid<br>
spacing is 4km and I'm covering about 500km in my cross section so I'm<br>
guessing things should appear relatively smooth.<br>
<br>
Why does maskout produce results that look so gridded? Is there any<br>
other way I can maskout the underground areas but still have a nice<br>
crisp outline of the terrain?<br>
<br>
Thanks!<br>
Dan Leins<br>
<br>
<br></div></div>
------------------------------------------------------------------------<br>
<br>
<br>
------------------------------------------------------------------------<br>
<br>
</blockquote>
</blockquote></div><br>