[Met_help] verifying 2m temperaures with PointStat

Joe Olson Joseph.B.Olson at noaa.gov
Thu Dec 4 17:25:12 MST 2008


I'm trying to verify 2m temps output from WRF-ARW v3.0 and using  
METv1.1 (PointStat).

First, I have performed the following to prepare the prepbufr:

1) cwordsh block prepbufr_file prepbufr_file.blk
2) Converted Prepbufr to netcdf:
         /METdir/bin/pb2nc prepbufr_file.blk prepbufr.ncdf / 
METconfigdir/PB2NCConfig

So far, this seems to work fine, since I have successfully used this  
file to verify upper-level temperatures...

Then, convert wrfouts (netcdf) to an A-grid (grib)
3) run_wrfpost
     This outputs my wrfprs_d01.XX file.
      In the wrf_cntrl.parm configuration file, I made sure to specify  
various forms of temperatures:

  (SHELTER TEMPERATURE ) SCAL=( 5.0)
  L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000  
00000 00000 00000)
  (TEMP ON PRESS SFCS  ) SCAL=( 4.0)
  L=(11111 11111 11111 11111 11111 11111 11111 11111 11111 10000 00000  
00000 00000 00000)
etc...

Then, when running PointStat with the output from (2) and (3) and  
using the following edits to the config file:

vx_grib_code[] = [ "TMP/Z2"];
thresholds[] = [ "gt0"];
message_type[] = [ "ADPSFC", "SFCSHP" ];

I get the requested output, but looking at the mpr file, the model  
forecast is way too cold (229-233 K) while the observations seem  
correct (270-278 K)
:

     TOTAL INDEX OBS_LAT OBS_LON OBS_LVL    OBS_ELV     FCST    OBS  
CLIMO
1     286     1   58.68 -156.65   989.3   14.98281 232.3029 277.05 -9999
2     286     2   59.75 -154.92   995.2   49.12878 233.0589 275.95 -9999
3     286     3   57.75 -152.50   997.9   33.93378 232.8184 278.15 -9999
4     286     4   59.52 -139.67  1026.3    9.02245 230.9946 272.55 -9999
5     286     5   58.37 -134.58  1029.5    6.98383 230.5191 270.95 -9999
6     286     6   62.30 -150.10  1004.5  109.13992 233.0715 269.25 -9999
7     286     7   55.13 -131.58  1031.8    0.00000 230.0803 277.15 -9999
8     286     8   55.13 -131.58  1031.8    0.00000 230.0803 276.15 -9999
9     286     9   55.35 -131.70  1027.3   28.94379 230.0819 275.35 -9999
10    286    10   56.97 -133.95  1031.5    0.00000 229.8285 274.15 -9999
11    286    11   56.97 -133.95  1031.2    0.00000 229.8285 273.15 -9999
etc....

Has anyone seen a problem like this? Are there any other steps I need  
to consider?

-joe


More information about the Met_help mailing list