[gradsusr] any grads users know why ARWpost interp_method=1 does not work for WRF v3.3.1

[Heath, Nicholas]

Hey Jonathan,

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.

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:

Original:

     IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 )     iprogram = 8 !! wrf.exe output
     IF ( iprogram == 0 ) THEN
         print*," "
         print*,"  WARNING --- I do not recognize this data."
         print*,"             ",trim(title)
         print*,"              Will make an attempt to read it."
         print*," "
      END IF

Add the following line underneath line 193:

      IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 )     iprogram = 8 !! wrf.exe output
      IF ( INDEX(title,'OUTPUT FROM *') /=0 )        iprogram = 8 !! wrf.exe output
      IF ( iprogram == 0 ) THEN
         print*," "
         print*,"  WARNING --- I do not recognize this data."
         print*,"             ",trim(title)
         print*,"              Will make an attempt to read it."
         print*," "
      END IF

After you make this change, do a ./clean -a, then re-compile.  After that, everything should work fine!  Good luck!

-Nick Heath

________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on behalf of SMITH, JONATHAN WYNN. (GSFC-6140) [jonathan.smith at nasa.gov]
Sent: Tuesday, June 19, 2012 2:10 PM
To: gradsusr at gradsusr.org
Subject: [gradsusr] any grads users know why ARWpost interp_method=1 does not work for WRF v3.3.1

Hello:

I have submitted this issue to wrfhelp at ucar.edu<mailto:wrfhelp at ucar.edu>.   But I thought I would try here to see if any grads users have run into the problem below using ARWpost.

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.


For example even with interp_method = 1 with WRF/Chem v3.3.1 output:

I get this message first

START PROCESSING DATA

   WARNING --- I do not recognize this data.
               OUTPUT FROM *             PROGRAM:WRF/CHEM V3.3.1 MODEL
               Will make an attempt to read it.

Then it just prints out the model levels

   16.19018
   16.38956
   16.58874
   16.78770
zdef   24 linear 1 1
tdef   13 linear 00Z20MAY2006      60MN
VARS    5
no2           24  0  NO2 concentration (ppmv)
no            24  0  NO concentration (ppmv)
o3            24  0  O3 concentration (ppmv)
hno3          24  0  HNO3 concentration (ppmv)
co            24  0  CO concentration (ppmv)

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:

 225.00000
 200.00000
 175.00000
 150.00000
 125.00000
 100.00000
tdef   13 linear 01Z30MAY2006      60MN
VARS    8
no2           36  0  NO2 concentration (ppmv)
no            36  0  NO concentration (ppmv)
o3            36  0  O3 concentration (ppmv)
hno3          36  0  HNO3 concentration (ppmv)
co            36  0  CO concentration (ppmv)
tc            36  0  Temperature (C)
dbz           36  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)


Thanks

Jonathan Wynn Smith, Ph. D.
jonathan.smith at nasa.gov
NASA-Goddard Space Flight Center
Code 614-Atmospheric Chemistry and Dynamics Branch
Greenbelt, MD 20771
301-614-5956

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120619/882047d5/attachment-0003.html