<div class="gmail_quote">On Fri, Jun 20, 2008 at 9:32 PM, Mike Johnson <<a href="mailto:mike.w.johnson@noaa.gov">mike.w.johnson@noaa.gov</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Users, Can anyone help with what I suspect is either a regrid problem<br>
(using re.gex) - or - in the methodology I'm using for output fwrite - or -<br>
an issue with my understanding of how options yrev works. </blockquote><div><br></div><div>This shouldn't be an issue in this case. Usually, when you misuse yrev you end up with an upside down plot, something easy to diagnose.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">I am attempting<br>
to regrid a subdomain of a large regional netcdf topo file<br>
(US_Rockies_GLOBE_30sec_topo.nc resolution - 30 Sec resolution) and write<br>
that regridded subdomain to a binary file (Idaho.dat) which I then want to<br>
overlay with high resolution surface wind field from CALMET (0.25km<br>
resolution).<br>
<br>
My basic procedure is to (1) convert grid spacing for control files and re<br>
() to DeltaX and DeltaY in degrees; (2) select a subdomain of the large<br>
NetCDF TOPO file (US_Rockies_GLOBE_30sec_topo.nc) which has Y dims written<br>
from north-to-south and (3) to regrid this smaller domain using re() in<br>
order to match domain size and dimensions of CALMET surface wind field; and<br>
(4) write a binary file (Idaho.dat) which I'm hoping to overlay with the<br>
CALMET wind data (and future runs of this model over the same domain).<br>
<br>
I am able to open large domain TOPO file (US_Rockies_GLOBE_30sec_topo.nc),<br>
regrid to higher resolution and plot subdomain. However...after writing<br>
this re-gridded subdomain file to binary output (Idaho.dat), when I try to<br>
plot it using (Idaho.ctl) the file appears to be garbage with GrADS<br>
echoing "Contouring: 1e+36 to 9e+36" whereas in reality values should be<br>
something like "1000 to 2800".<br></blockquote><div><br></div><div>Hmm... In your ctl you have</div><div><br></div><div><span class="Apple-style-span" style="border-collapse: collapse; ">UNDEF -9999.0</span><br>
</div><div><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">are you sure of that? Do you have any UNDEF in your dataset? Make sure re() is not adding undef at the borders. Here is a couple of things you could do to further diagnose your problem:</span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">1) Specify the clevs when contouring your data:</span></div><div>
<span class="Apple-style-span" style="border-collapse: collapse;"><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">ga-> set clevs 1000 1500 2000 2500 3000</span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">ga-> d tpo</span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">2) Examine the min/max of your data:</span></div><div><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">ga-> set gxout stat</span></div><div><span class="Apple-style-span" style="border-collapse: collapse;">ga-> d tpo</span></div><div>
<span class="Apple-style-span" style="border-collapse: collapse;"><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse;"> Arlindo</span></div><div><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></div><div><br></div></div><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>