[Met_help] [Fwd: MET point-stat]

Fri Mar 12 11:23:54 MST 2010

Jianting,

MET does not currently support data on a rotated lat/lon projection.  If there's a need from WRF-ARW users for us to support this directly, we could add code to do so.

However, I was going to recommend that you run your data through the copygb tool to interpolate it to one of the projections that MET does currently support: lat/lon, lambert conformal, polar
stereographic, or mercator.  However, when I tried running the file you sent through copygb, it just ran silently and produced an output file of zero length.

I'll ask the WPP support people what's going on to find out if there is a version of copygb that can read these files.

Thanks,
John

Jianting Chu wrote:
> Dear John,
>    I am sorry to bother you again.
>    After testing the WPP-beta and some discussion with wrfhelp. Now,
> this tool can be used to rotated lat-lon project. But, the ./point_stat
> still can not work. Whether I use Z0 or Z2 in the PointStatConfig. I
> think the TMP is correct in the GRIB file.
>    When I run: ./point_stat WRFPRS_d01.000 wess_198901.nc
> PointStatConfig_t2m, it gave the following message:
> ---------------------------------------------------------------
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1125292619
> Forecast File: WRFPRS_d01.000
> Climatology File: none
> Configuration File: PointStatConfig_t2m
> Observation File: wess_198901.nc
> 
> 
> ERROR: read_gds() -> Grid type 205 not currently supported.
> -----------------------------------------------------------
>    Please check the GRIB file and PointStatConfig in the attached file.
>    Thank you very much!
> 
> Best wishes,
> Jianting
> 
> 
> John Halley Gotway wrote:
>> Jianting,
>>
>> Your question about the MET Point-Stat tool was forwarded to us here
>> at met_help at ucar.edu.  Please send future MET-related questions
>> directly to "met_help at ucar.edu".
>> And any questions you have about WRF or the WRF Post-Processor should
>> be sent to "wrfhelp at ucar.edu".
>>
>> OK, there's two ways we can go with this...
>>
>> (1) First, I took a look at the GRIB file you sent.  When you run it
>> through the "wgrib" tool, you'll see the following output:
>>
>> wgrib COPS_T2M_1989-01-15.grb
>>   
>> 1:0:d=89011500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 2:3248:d=89011503:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 3:6496:d=89011506:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 4:9744:d=89011509:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 5:12992:d=89011512:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 6:16240:d=89011515:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 7:19488:d=89011518:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>   
>> 8:22736:d=89011521:TMP:kpds5=11:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
>>
>>
>> It looks like this file contains 8 records of temperature.  However,
>> the level value for all of these records is set to 0.  For 2-meter
>> temperature, the level should be set to a value of 2.  My guess
>> is that these should really be levels like 2-meters, 1000mb, 800mb,
>> and so on.  In your configuration file, I tried changing "TMP/Z2" to
>> "TMP/Z0" and MET was able read the data.  However, you really
>> should fix how the level value is being written in the GRIB file.
>>
>> (2) Second, there is a new beta version of the WRF PostProcessor that
>> DOES support the ARW rotated lat-long grid.  We're currently in the
>> testing phase with these changes.  And we'd like to find users
>> to try the code out.  I've attached this WPP beta version to this
>> message.  If you're willing to test it out for us, we'd appreciate
>> it.  You can send any questions or comments that come up directly
>> to wrfhelp at ucar.edu.
>>
>> Hope that helps,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>> wrfhelp wrote:
>>  
>>> John,
>>>
>>> Would you please help take a look at this email? Thanks!
>>>
>>> Ming
>>>
>>>
>>> ------------------------------------------------------------------------
>>>
>>> Subject:
>>> MET point-stat
>>> From:
>>> Jianting Chu <Jianting.Chu at uni-hohenheim.de>
>>> Date:
>>> Sun, 28 Feb 2010 18:04:41 +0100
>>> To:
>>> wrfhelp <wrfhelp at ucar.edu>
>>>
>>> To:
>>> wrfhelp <wrfhelp at ucar.edu>
>>>
>>>
>>> Dear Sir,
>>>    I have some trouble with point-stat of METv2.0. Would you please
>>> giving me some advices?
>>>    Because I use rotated lat-lon projection, I can not use WPP to get
>>> GRIB files from wrfout*.nc. In stead, I transfer wrfout*.nc to CF.nc
>>> with NCL firstly, then, transfer it to GRIB file with CDO.
>>>    Please check the attenched file for my GRIB file and PointStatConfig.
>>> However, it always said that:
>>> ----------------------------------------------------
>>>
>>> ERROR: process_grib_codes() -> no records matching GRIB code 11 with
>>> level indicator of Z2 found in GRIB file: COPS_T2M_1989-12-30.grb
>>>
>>> ----------------------------------------------------
>>>    I have changed the variable name from T_2m to var11.
>>>    Are there some special name or rule for GRIB files to METool? Would
>>> you please provide a sample GRIB file from WPP?
>>>    Thank you very much in advance.
>>>
>>> Best regards,
>>> Jianting
>>>
>>>
>>> ------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> Met_help mailing list
>>> Met_help at mailman.ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/met_help