[Met_help] PWAT verification with MET Tools

[Thomas Schwitalla]

John,

when I try your suggestion, I got the following message:

ERROR: process_grib_codes() -> no records matching GRIB code 54 with 
level indicator of Z0 found in GRIB file: 
/taifun/20070720_cops_WRF_sensitivity/assim/wrf_files/postprd/WRFPRS_2007072000_d01.002

"wgrib -v" on this file gives:

171:80565400:D=2007072001:PWAT:atmos col:kpds=54,200,0:1hr fcst:winds in 
grid direction:"Precipitable water [kg/m^2]

A sample line of my ASCII file reads:

ADPSFC  0007  20070720_233700     52.9480     12.3922     60.9550  54  
-9999     60.9550     31.3000

Now I am a bit confused...

Thomas

John Halley Gotway schrieb:
> Thomas,
>
> Just to be clear...
> (1) Your forecast from WPP of precipitable water is a value that's integrated over the whole column of air above each grid point.
> (2) Your point observations of precipitable water are the same - a value for the whole column of air above the observation point at the surface.
>
> If that's true, I'd set up the PointStatConfig file as follows:
> fcst_field[] = [ "PWAT/Z0" ];
> message_type[] = [ "ADPSFC" ];
>
> That means "give me the precipitable water field with a level value of 0", which is how's it's stored in your GRIB file.  And verify against "surface" observations.  When you call ASCII2NC, be sure to
> encode these PWAT observations using the "ADPSFC" message type.
>
> Give that a try and let me know if you have any problems.
>
> Thanks,
> John
>
> Thomas Schwitalla wrote:
>   
>> John,
>>
>> I want to verify WRF forecasts (postprocessed with WPP)  with PWAT data
>> from GPS stations. I already wrote it to NETCDF with ascii2nc but I
>> haven't figured out how to edit PointStatConfig that I can use PWAT. My
>> problem is: What value to set for the level type/code as PWAT is an
>> integrated value?
>>
>> wgrib output looks as follows:
>>
>> PWAT kpds5=54 kpds6=200 kpds7=0 levels=(0,0) grid=255 atmos col 1hr fcst:
>>
>>
>> Best regards,
>> Thomas
>>
>>     


-- 
*****************************************
*Thomas Schwitalla                      *
*Institute of Physics and Meteorology   *
*University of Hohenheim                *
*Garbenstrasse 30                       *
*70599 Stuttgart                        *
*Germany                                *
*Tel.: +49 711 459 22154                *
*****************************************