<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style>
<!--
@font-face
{font-family:Calibri}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif"}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline}
span.EmailStyle17
{font-family:"Calibri","sans-serif";
color:windowtext}
@page WordSection1
{margin:1.0in 1.0in 1.0in 1.0in}
-->
</style><style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1" lang="EN-US" link="blue" vlink="purple">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hey Jonathan,<br>
<br>
You're in luck...there is a fairly easy fix to get this working! The problem arises with the header of the netcdf files that the newer versions of WRF-Chem produce---that "Output from *" line.<br>
<br>
In your ARWpost folder, go into the 'src' directory. From there, open the 'input_module.f90'. Go to line 193, and make the following change:<br>
<br>
<span style="font-weight: bold;">Original:<br>
</span><br>
IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 ) iprogram = 8 !! wrf.exe output<br>
IF ( iprogram == 0 ) THEN<br>
print*," "<br>
print*," WARNING --- I do not recognize this data."<br>
print*," ",trim(title)<br>
print*," Will make an attempt to read it."<br>
print*," "<br>
END IF<span style="font-weight: bold;"><br>
<br>
</span>Add the following line underneath line 193:<br>
<br>
IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 ) iprogram = 8 !! wrf.exe output<br>
<span style="font-style: italic;"> </span><span style="font-weight: bold; font-style: italic;">IF ( INDEX(title,'OUTPUT FROM *') /=0 ) iprogram = 8 !! wrf.exe output</span><br>
IF ( iprogram == 0 ) THEN<br>
print*," "<br>
print*," WARNING --- I do not recognize this data."<br>
print*," ",trim(title)<br>
print*," Will make an attempt to read it."<br>
print*," "<br>
END IF<span style="font-weight: bold;"><span style="font-weight: bold;"><br>
<br>
</span></span>After you make this change, do a ./clean -a, then re-compile. After that, everything should work fine! Good luck!<br>
<br>
-Nick Heath<span style="font-weight: bold;"><span style="font-weight: bold;"><br>
<br>
</span></span>
<div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF998664"><font color="#000000" face="Tahoma" size="2"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of SMITH, JONATHAN WYNN. (GSFC-6140) [jonathan.smith@nasa.gov]<br>
<b>Sent:</b> Tuesday, June 19, 2012 2:10 PM<br>
<b>To:</b> gradsusr@gradsusr.org<br>
<b>Subject:</b> [gradsusr] any grads users know why ARWpost interp_method=1 does not work for WRF v3.3.1<br>
</font><br>
</div>
<div></div>
<div>
<div class="WordSection1">
<p class="MsoNormal">Hello:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I have submitted this issue to <a href="mailto:wrfhelp@ucar.edu" target="_blank">
wrfhelp@ucar.edu</a>. But I thought I would try here to see if any grads users have run into the problem below using ARWpost.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">For both ARWpost v2.2 and 3.1, the interp_method = 1 is not produce .ctl and .dat files on the specified pressure levels for WRF/Chem v3.3.1 output. However interp_method = 1 prints the .ctl and .dat on the specified pressure levels
with WRF/Chem v3.2.1 data. I have no choice but to use WRF/Chem v 3.3.1 for my research because the pre-processing emissions software for WRF/Chem includes a mercator projection that only works with v3.3.1 and not wit v3.2.<br>
<br>
<i><br>
For example even with interp_method = 1 with WRF/Chem v3.3.1 output:</i><br>
<br>
I get this message first<br>
<br>
START PROCESSING DATA<br>
<br>
WARNING --- I do not recognize this data.<br>
OUTPUT FROM * PROGRAM:WRF/CHEM V3.3.1 MODEL<br>
Will make an attempt to read it.<br>
<br>
Then it just prints out the model levels<br>
<br>
16.19018<br>
16.38956<br>
16.58874<br>
16.78770<br>
zdef 24 linear 1 1<br>
tdef 13 linear 00Z20MAY2006 60MN<br>
VARS 5<br>
no2 24 0 NO2 concentration (ppmv)<br>
no 24 0 NO concentration (ppmv)<br>
o3 24 0 O3 concentration (ppmv)<br>
hno3 24 0 HNO3 concentration (ppmv)<br>
co 24 0 CO concentration (ppmv)<br>
<br>
<i>However with the v3.2 or v3.2.1, I do not get the messages and the data is printed on the specified pressure levels:</i><br>
<br>
225.00000<br>
200.00000<br>
175.00000<br>
150.00000<br>
125.00000<br>
100.00000<br>
tdef 13 linear 01Z30MAY2006 60MN<br>
VARS 8<br>
no2 36 0 NO2 concentration (ppmv)<br>
no 36 0 NO concentration (ppmv)<br>
o3 36 0 O3 concentration (ppmv)<br>
hno3 36 0 HNO3 concentration (ppmv)<br>
co 36 0 CO concentration (ppmv)<br>
tc 36 0 Temperature (C)<br>
dbz 36 0 Reflectivity (-)<br>
max_dbz 1 0 Max Reflectivity (-)<br>
<br>
<br>
Thanks</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><b><span style="font-size: 12pt;">Jonathan Wynn Smith, Ph. D.</span></b></p>
<p class="MsoNormal"><u>jonathan.smith@nasa.gov</u></p>
<p class="MsoNormal">NASA-Goddard Space Flight Center</p>
<p class="MsoNormal">Code 614-Atmospheric Chemistry and Dynamics Branch </p>
<p class="MsoNormal">Greenbelt, MD 20771</p>
<p class="MsoNormal">301-614-5956</p>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
</div>
</body>
</html>