John,<br> Thanks for your explanation! I will apply the similar method to my nested domain and let you know the result next week. <br>Have a nice weekend!<br><br>Hao<br><br><div class="gmail_quote">On Fri, May 7, 2010 at 3:32 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;">Hao,<br>
<br>
You are welcome.<br>
<br>
The (24.828,100.01) is the lower left corner of your CMAQ grid. I simply dumped the x,y data, sorted it, and found the point.<br>
<br>
95.969 is the longitude of your upper left corner and 131.2617 is the longitude of your upper right corner. Take the difference and divide by 94 to get the x delta.<br>
<br>
Similarly, take the difference between your lower and upper latitude and divide by 67 to get the y delta.<br>
<div class="im"><br>
John Huddleston, PhD<br>
<br>
From: Hao He [mailto:<a href="mailto:haohe.umcp@gmail.com">haohe.umcp@gmail.com</a>]<br>
</div>Sent: Friday, May 07, 2010 1:24 PM<br>
<div><div></div><div class="h5">To: Huddleston, John<br>
Subject: Re: [gradsusr] data not in GRIB format<br>
<br>
John,<br>
Thanks a lot! It looks very nice. Can you tell me how to calculate the latref (24.828N), lonref(100.01), xdef and ydef? I have another nested domain, and I want to go through the process.<br>
<br>
Best<br>
<br>
Hao<br>
On Fri, May 7, 2010 at 3:19 PM, Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>> wrote:<br>
Hao, attached is an image from your data file using<br>
<br>
d ave(so2,t=1,t=25))<br>
<br>
where the control file is<br>
<br>
dset ^SO2_CHINA_D1.ncf<br>
dtype netcdf<br>
title Vertically integrated atmos chem/aerosol vbles<br>
pdef 94 67 lcc 24.828 100.01 1 1 30 60 113.7 30000 30000<br>
undef -9999<br>
xdef 94 linear 95.769 0.3659<br>
ydef 67 linear 24.828 0.2644<br>
zdef 34 linear 1 1<br>
tdef 25 linear 00:00Z20apr2006 1hr<br>
vars 3<br>
LAT=>lat 1 t,z,y,x Latitude<br>
LON=>lon 1 t,z,y,x Longitude<br>
SO2=>so2 1 t,z,y,x Sulfur Dioxide<br>
endvars<br>
<br>
I don't know your data as well as you do; however, it appears to plot OK.<br>
<br>
What do you think is the problem?<br>
<br>
John Huddleston, PhD<br>
<br>
From: Hao He [mailto:<a href="mailto:haohe.umcp@gmail.com">haohe.umcp@gmail.com</a>]<br>
Sent: Friday, May 07, 2010 12:26 PM<br>
To: Huddleston, John<br>
<br>
Subject: Re: [gradsusr] data not in GRIB format<br>
<br>
John,<br>
I created a file only with SO2 data. It contains 34 layers and 25 time frames. Please see attached file. I also attach one snapshot by ncview.<br>
<br>
Hao<br>
On Fri, May 7, 2010 at 2:01 PM, Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>> wrote:<br>
Hao<br>
<br>
Put out a dataset for me to download. All I have is that one latitude longitude file.<br>
<br>
What doesn't happen? Where is the data plotted?<br>
<br>
John<br>
<br>
John Huddleston, PhD<br>
<br>
From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Hao He<br>
Sent: Friday, May 07, 2010 11:11 AM<br>
<br>
To: GrADS Users Forum<br>
Subject: Re: [gradsusr] data not in GRIB format<br>
<br>
John,<br>
I use your template to open the file. But it is not right. Please have a look at my .ctl file.<br>
<br>
dset ^EMIS_CHINA_D1.ncf<br>
title INTEX-B SO2 emission<br>
undef -9999<br>
pdef 94 67 lcc 34.8 113.7 1 1 30.0 60.0 113.7 30000 30000<br>
xdef 94 linear 95.769 0.3659<br>
ydef 67 linear 24.828 0.2644<br>
tdef 25 linear 00:00Z03APR2008 1HR<br>
zdef 34 linear 1 1<br>
dtype netcdf<br>
vars 1<br>
SO2=>so2 0 t,z,y,x Sulfur Dioxide<br>
endvars<br>
<br>
Best<br>
<br>
Hao<br>
On Fri, May 7, 2010 at 12:06 PM, Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>> wrote:<br>
Hao,<br>
<br>
Based on the information you sent, here is a representative control file to plot your (COL,ROW)=(94,67) grid.<br>
<br>
dset ^YourFile.ncf<br>
dtype netcdf<br>
title Vertically integrated atmos chem/aerosol vbles<br>
pdef 94 67 lcc 24.828 100.01 1 1 30 60 113.7 30000 30000<br>
undef -9999<br>
xdef 94 linear 95.769 0.3659<br>
ydef 67 linear 24.828 0.2644<br>
.<br>
<br>
<br>
John Huddleston, PhD, PE<br>
<br>
<br>
From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Hao He<br>
Sent: Friday, May 07, 2010 8:43 AM<br>
<br>
To: GrADS Users Forum<br>
Subject: Re: [gradsusr] data not in GRIB format<br>
<br>
John,<br>
<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>
_______________________________________________<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>
<br>
<br>
<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>
</div></div></blockquote></div><br>