[Wrf-users] ERA-Interim and calc_ecmwf_p

[Bart Brashers]

I'm using WRFv3.4.1, trying to initialize from ERA-Interim data.  I've tried both ds627.0 and directly from http://data-portal.ecmwf.int/data/d/dataserver.  For this discussion, let's assume it's the more commonly used ds627.0.



I ran ungrib.exe using Vtable.ERA-interim.ml using the types of files specified at the end of the Vtable file.  Using



# rd_intermediate.exe FILE:2005-12-31_12 | grep LEVEL



I can see the 200100 level, and levels labeled 60.000000 down to 1.000000.  I extracted the 60-layer table from

http://www.ecmwf.int/products/data/technical/model_levels/model_def_60.html and put it in a file ecmwf_coeffs.



I used the same namelist.wps and Vtable to run calc_ecmwf_p.exe.  I get PRES:* files, but they still contain the same 200100 level and levels from 1.00000 to 60.0000.  Wasn't I supposed to get values like 100000 and 85000?  My namelist DOES have the block



&mod_levs

press_pa = 201300 , 200100 , 100000 ,

             95000 ,  90000 ,

             85000 ,  80000 ,

             75000 ,  70000 ,

             65000 ,  60000 ,

             55000 ,  50000 ,

             45000 ,  40000 ,

             35000 ,  30000 ,

             25000 ,  20000 ,

             15000 ,  10000 ,

              5000 ,   1000



included.



For future searches of this list, the error in the next step (running obsgrid.exe) is found its standard output:



Error finding the 850, 700, 500 hPa levels.

pressure=    1001.000        1012.000        1009.000        1005.000

    999.0000        990.0000        979.0000        966.0000

    950.0000        931.0000        909.0000        884.0000

    858.0000        828.0000        797.0000        763.0000

    728.0000        692.0000        654.0000        616.0000

    577.0000        539.0000        500.0000        462.0000

    425.0000        388.0000        353.0000        320.0000

    288.0000        257.0000        229.0000        202.0000

    177.0000        154.0000        133.0000        113.0000

    96.00000        80.00000        67.00000        55.00000

    44.00000        36.00000        29.00000        23.00000

    19.00000        15.00000        12.00000        10.00000

    8.000000        6.000000        5.000000        4.000000

    3.000000        3.000000        2.000000        2.000000

    1.000000        1.000000        1.000000        0.000000

    0.000000



Any suggestions?  Is there another version of calc_ecmwf_p.f90 that is known to work with ds627.0 data?



Thanks,



Bart Brashers















________________________________
This message contains information that may be confidential, privileged or otherwise protected by law from disclosure. It is intended for the exclusive use of the Addressee(s). Unless you are the addressee or authorized agent of the addressee, you may not review, copy, distribute or disclose to anyone the message or any information contained within. If you have received this message in error, please contact the sender by electronic reply to email at environcorp.com and immediately delete all copies of the message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20120910/55390da6/attachment.html