[Met_help] verifying upper-level wind speed with PointStat
John Halley Gotway
johnhg at rap.ucar.edu
Mon Jan 5 15:07:00 MST 2009
Joe,
I believe this is being caused by a bug that's posted to the known issues page on the MET website.
On the known issues page (http://www.dtcenter.org/met/users/support/known_issues/METv1.1/index.php), please search for the following bug report:
Bug in the vx_met_util and the vx_data_grids libraries when rotating winds from grid-relative to earth-relative.
Posted 10/14/2008
Please try grabbing the fix, rebuilding MET, and see if the problem is corrected. To make things a bit easier, you could download a tarball containing all of the recent fixes for METv1.1 bundled
together. You can download the tarball of patches from:
http://www.dtcenter.org/met/users/support/known_issues/METv1.1/patches/METv1.1_patches_20081222.tar.gz
Once you download it, move that tarball file into the METv1.1 directory. And then, untar it. This will overwrite several of the original source files with their updated versions. Then do a "make
clean" followed by a "make".
Please let me know if you're still seeing problems after implementing this fix.
Thanks,
John
Joe Olson wrote:
> Hi John,
>
> I've uploaded some files I used to verify the WRF-based RapidRefresh
> System for 20090105 02Z cycle. The verification plots for this 12 hour
> forecast can be found at:
>
> http://rapidrefresh.noaa.gov/AKverify/all.cgi?R1=RAPREFCYC&R2=NONE&run_time_27=05%20Jan%202009%20-%2002Z&D1=ph1
>
>
> Forecast hour 10 (valid at 12Z) has both surface and upper-air
> verification. By clicking on the "Mon 12" column for "sfc wsp" and "ua
> wsp", you can see that the model surface biases seem reasonable, but the
> upper-air biases are too large. Note: I only use the mpr output from
> pointstat to make these figures.
>
> The files I used to make these plots are:
>
>> sample forecast file (GRIB output of WPP): wrfgrbALL_d01_2009010502.10
>> observation file (output of PB2NC): PBFR_2009010512.nc
>> Point-Stat configuration: PointStatConfig_UA_WSP2
>
> I had permission errors while trying to upload these files, so I posted
> them at:
> http://yoda.msrc.sunysb.edu/joe/MET/
>
> Thanks in advance,
> Joe
>
> PS: Its not important but in case you are interested, here are the
> errors I had while attempting the upload:
>
> [joe at yoda ~]$ ftp ftp.rap.ucar.edu
> Connected to muffet.rap.ucar.edu.
> 220-
> 220- *** National Center for Atmospheric Research ***
> 220- *** Research Applications Program ***
> 220-
> 220- *** Use of this system for unauthorized purposes is a ***
> 220- *** violation of federal law. See Public Law 99-474. ***
> 220- *** Any activity may be logged and monitored. By ***
> 220- *** using this system, you consent to such logging and ***
> 220- *** monitoring. ***
> 220-
> 220- The local time is: Mon Jan 5 11:39:59 2009
> 220-
> 220- If you have problems logging in, please try using a dash (-) as the
> 220- first character of your password--this will turn off the continuation
> 220- messages that may be confusing your ftp client.
> 220-
> 220- If you have unusual problems, please contact: ftpadmin at rap.ucar.edu
> 220-
> 220 muffet.rap.ucar.edu FTP server (Version wu-2.6.2(1) Mon Dec 11
> 13:26:04 UTC 2006) ready.
> 530 Please login with USER and PASS.
> 530 Please login with USER and PASS.
> KERBEROS_V4 rejected as an authentication type
> Name (ftp.rap.ucar.edu:joe): anonymous
> 331 Guest login ok, send your complete e-mail address as password.
> Password:
> 230 Guest login ok, access restrictions apply.
> Remote system type is UNIX.
> Using binary mode to transfer files.
> ftp> put PBFR_2009010512.nc
> local: PBFR_2009010512.nc remote: PBFR_2009010512.nc
> 227 Entering Passive Mode (128,117,192,211,192,4)
> 553 PBFR_2009010512.nc: Permission denied on server. (Upload)
>
>
> On Jan 5, 2009, at 10:07 AM, John Halley Gotway wrote:
>
>> Joe,
>>
>> Sorry for the delay in getting back to you. If you're still
>> experiencing these issues, please send me a sample forecast file (GRIB
>> output of WPP), observation file (output of PB2NC), and the
>> Point-Stat configuration file you're using.
>>
>> I'll run it through the debugger and see if I can figure out what's
>> going on. You can post these file to our anonymous FTP site using the
>> following commands:
>> ftp ftp.rap.ucar.edu
>> username = anonymous
>> password = "your email address"
>> cd incoming/irap/johnhg
>> put "your files"
>> bye
>>
>> FYI, when you request WIND speed in MET, it'll look in the GRIB file
>> for a WIND variable at the requested level. If MET can't find one,
>> it'll look for the U and V components at the level and derive a
>> wind speed from them.
>>
>> Thanks,
>> John
>>
>>
>> Joe Olson wrote:
>>> I'm trying to verify 300, 500, 700, 850, and 925 mb wind speeds from
>>> WRF-ARW v3.0 and using METv1.1 (PointStat).
>>>
>>> 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
>>>
>>> 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
>>> winds on pressure surfaces:
>>>
>>> (U WIND ON PRESS SFCS) SCAL=( 4.0)
>>> L=(11111 11111 11111 11111 11111 11111 11111 11111 11111 10000 00000
>>> 00000 00000 00000)
>>> (V WIND ON PRESS SFCS) SCAL=( 4.0)
>>> L=(11111 11111 11111 11111 11111 11111 11111 11111 11111 10000 00000
>>> 00000 00000 00000)
>>>
>>> (I do not see any wind speed variables...)
>>>
>>> Then, when running PointStat with the output from (2) and (3) and using
>>> the following edits to the config file:
>>>
>>> vx_grib_code[] = [
>>> "WIND/P300","WIND/P500","WIND/P700","WIND/P850","WIND/P925"];
>>> thresholds[] = [ "gt0", "gt0", "gt0", "gt0", "gt0"];
>>> message_type[] = [ "ADPUPA" ];
>>>
>>> I get the requested mpr output file, however, the model winds do not
>>> compare well with that seen in the wrfouts (via ncview). The
>>> observations, on the other hand, seem reasonable (soundings had wind
>>> speeds ~25 m/s). Here's a sample output for a single station in alaska:
>>>
>>> VERSION MODEL FCST_LEAD FCST_VALID OBTYPE VX_MASK
>>> TYPE VAR LEVEL = TOTAL
>>> INDEX OBS_LAT OBS_LON OBS_LVL
>>> OBS_ELV FCST OBS CLIMO
>>> INTERP_MTHD INTERP_PNTS
>>> V1.1 WRF 000000 20081229_120000 ADPUPA FULL
>>> MPR WIND P300-300 = 13
>>> 9 60.79000 -161.84000 300.00000
>>> -9999 1.80940 25.65638 -9999
>>> DW_MEAN 4
>>> V1.1 WRF 000000 20081229_120000 ADPUPA FULL
>>> MPR WIND P500-500 = 13
>>> 9 60.79000 -161.84000 500.00000
>>> -9999 2.47792 23.65354 -9999
>>> DW_MEAN 4
>>> V1.1 WRF 000000 20081229_120000 ADPUPA FULL
>>> MPR WIND P700-700 = 13
>>> 9 60.79000 -161.84000 700.00000
>>> -9999 1.22441 24.68299 -9999
>>> DW_MEAN 4
>>> V1.1 WRF 000000 20081229_120000 ADPUPA FULL
>>> MPR WIND P850-850 = 13
>>> 9 60.79000 -161.84000 850.00000
>>> -9999 2.20001 26.26119 -9999
>>> DW_MEAN 4
>>> V1.1 WRF 000000 20081229_120000 ADPUPA FULL
>>> MPR WIND P925-925 = 13
>>> 9 60.79000 -161.84000 925.00000
>>> -9999 1.19297 24.67387 -9999
>>> DW_MEAN 4
>>>
>>> using the "-v 3" option with point_stat shows that the "Grib Record
>>> Index" is never ==32, which I *think* it should be if wind speed is to
>>> be used. My guess is that I am missing the wind speed variable in my
>>> wrf_cntrl.parm configuration file...(?) Also, if I use copygb to strip
>>> down the grib file [output from (3)] such that it only contains pressure
>>> level data (i.e. copygb.exe -k'5*-1 100' -x WRFgrib_ALL WRFgrib_PRS), I
>>> still have the same problem, so I don't think its using the wrong
>>> vertical scale (like meters above ground level).
>>>
>>>
>>> thanks,
>>> -joe
>>> _______________________________________________
>>> Met_help mailing list
>>> Met_help at mailman.ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/met_help
More information about the Met_help
mailing list