John,<br> I get the WRF WPS projection information as following:<br><br>e_we = 98, 118,<br> e_sn = 70, 103,<br> geog_data_res = '10m','2m', <br> dx = 30000,<br>
dy = 30000,<br> map_proj = 'lambert',<br> ref_lat = 34.8,<br> ref_lon = 113.7,<br> truelat1 = 30.0,<br> truelat2 = 60.0,<br> stand_lon = 113.7,<br><br>My key point here is : CMAQ domain is a little bit smaller than WRF domain because we eliminate the outside columns and rows to diminish lateral effects. So my WRF domain has 98*70 grids while my CMAQ domain has 94*67 grids.<br>
<br>so the question is: how to set (isize, jsize) and (iref, jref) in pdef variable? I am confused.<br><br>Best<br><br>Hao<br><br><br><div class="gmail_quote">On Thu, May 6, 2010 at 2:50 PM, Huddleston, John <span dir="ltr"><<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div link="blue" vlink="purple" lang="EN-US">
<div>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Hao</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">CMAQ files usually are ioapi or WRF compliant. If they are
ioapi, run the ‘ncdump –h file.ncf > h’ and look at the
global characteristics. It will have everything you need to compose the
control file.</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">If they are WRF compliant, you will have to go back to the modeler
who created them.</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Please see <a href="http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html#GRID138" target="_blank">http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html#GRID138</a>
as an example of the information you need.</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 10pt; color: rgb(31, 73, 125);">John Huddleston, PhD</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color; padding: 3pt 0in 0in;">
<p class="MsoNormal"><b><span style="font-size: 10pt;">From:</span></b><span style="font-size: 10pt;">
<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] <b>On
Behalf Of </b>Hao He<br>
<b>Sent:</b> Thursday, May 06, 2010 12:15 PM<div><div></div><div class="h5"><br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] data not in GRIB format</div></div></span></p>
</div><div><div></div><div class="h5">
<p class="MsoNormal"> </p>
<p class="MsoNormal" style="margin-bottom: 12pt;">John,<br>
I start to work on the .ctl file, and I have a question about the
xdef and ydef. As you plotted, the left corner is not the minimal longitude and
latitude of the domain. Which value should I use?<br>
<br>
Best<br>
<br>
Hao</p>
<div>
<p class="MsoNormal">On Wed, May 5, 2010 at 2:16 PM, Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu" target="_blank">Huddleston@cira.colostate.edu</a>>
wrote:</p>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">Hao, this might be closer if it is 30KM. Change the ytime and
year accordingly. John</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">dset ^YourFile.ncf<br>
dtype netcdf<br>
title Vertically integrated atmos chem/aerosol vbles</span></p>
</div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">pdef 94 67 lcc 24.87 -127.10 1 1 33 45 -114.00
30000 30000<br>
undef -9999<br>
xdef 94 linear -131.336 0.48<br>
ydef 67 linear 24.87 0.300</span></p>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"><br>
zdef 1 linear 1 1</span></p>
</div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">tdef 24 linear 00:00Z20apr2006 1hr</span></p>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"><br>
vars 8<br>
NO2=>no2 1 t,z,y,x NO2<br>
NO=>no 1 t,z,y,x NO<br>
O=>o 1 t,z,y,x O<br>
O3=>o3 1 t,z,y,x O3<br>
NO3=>no3 1 t,z,y,x NO3<br>
CO=>co 1 t,z,y,x CO<br>
CXO3=>cxo3 1 t,z,y,x cxo3<br>
NH3=>nh3 1 t,z,y,x NH3<br>
endvars</span></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
</div>
</div>
<div>
<div class="MsoNormal" style="text-align: center;" align="center"><span style="font-size: 10pt; color: black;">
<hr width="100%" align="center" size="2">
</span></div>
<div>
<p class="MsoNormal"><b><span style="font-size: 10pt; color: black;">From:</span></b><span style="font-size: 10pt; color: black;"> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>
[<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>]
On Behalf Of Hao He [<a href="mailto:haohe.umcp@gmail.com" target="_blank">haohe.umcp@gmail.com</a>]</span></p>
</div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><b><span style="font-size: 10pt; color: black;">Sent:</span></b><span style="font-size: 10pt; color: black;">
Wednesday, May 05, 2010 11:59 AM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] data not in GRIB format</span></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><span style="font-size: 10pt; color: black;">That is great! My grid size is
30 km instead of 36 km, and I think the difference is trivial. I will read the
format and modified it to my own grid. <br>
Thanks a lot!<br>
<br>
Hao</span></p>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">On Wed, May 5, 2010 at 1:52 PM, Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu" target="_blank">Huddleston@cira.colostate.edu</a>>
wrote:</span></p>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">Hao</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">If it is 36 kilometer CMAQ netCDF you can probably
start with this grid control definition. </span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">dset ^YourFile.ncf<br>
dtype netcdf<br>
title Vertically integrated atmos chem/aerosol vbles<br>
pdef 148 112 lcc 16.77 -121.51 1 1 33 45 -97.00 36000 36000<br>
undef -9999<br>
xdef 148 linear -139.285 0.57565<br>
ydef 112 linear 16.77 0.3600<br>
zdef 1 linear 1 1<br>
tdef 72 linear 00:00Z20apr2006 1hr<br>
vars 8<br>
NO2=>no2 1 t,z,y,x NO2<br>
NO=>no 1 t,z,y,x NO<br>
O=>o 1 t,z,y,x O<br>
O3=>o3 1 t,z,y,x O3<br>
NO3=>no3 1 t,z,y,x NO3<br>
CO=>co 1 t,z,y,x CO<br>
CXO3=>cxo3 1 t,z,y,x cxo3<br>
NH3=>nh3 1 t,z,y,x NH3<br>
endvars</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">john Huddleston, PhD, PE</span></p>
</div>
<div>
<div class="MsoNormal" style="text-align: center;" align="center"><span style="font-size: 10pt; color: black;">
<hr width="100%" align="center" size="2">
</span></div>
<div>
<p class="MsoNormal"><b><span style="font-size: 10pt; color: black;">From:</span></b><span style="font-size: 10pt; color: black;"> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>
[<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>]
On Behalf Of Hao He [<a href="mailto:haohe.umcp@gmail.com" target="_blank">haohe.umcp@gmail.com</a>]</span></p>
</div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><b><span style="font-size: 10pt; color: black;">Sent:</span></b><span style="font-size: 10pt; color: black;">
Tuesday, May 04, 2010 1:58 PM<br>
<b>To:</b> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<b>Subject:</b> [gradsusr] data not in GRIB format</span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;">Hi,<br>
I have a 2D data set with uneven Lon and Lat distribution. I
notice that the GRIB data for GrADS has even grids, so I am wondering whether
it is possible to display this data set on GrADS? Thanks for help!<br>
<br>
Hao</span></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><span style="font-size: 10pt; color: black;"><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></span></p>
</div>
<p class="MsoNormal"><span style="font-size: 10pt; color: black;"> </span></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><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"> </p>
</div></div></div>
</div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br>