<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content=text/html;charset=iso-8859-1>
<META content="MSHTML 6.00.6001.18226" name=GENERATOR></HEAD>
<BODY id=MailContainerBody
style="PADDING-RIGHT: 10px; PADDING-LEFT: 10px; PADDING-TOP: 15px" leftMargin=0
topMargin=0 CanvasTabStop="true" name="Compose message area">
<DIV><FONT face=Arial size=2>Hi all.</FONT></DIV>
<DIV><FONT face=Arial size=2>I would like to use MET, in particular GRID-STAT,
to compare WRF_NMM at a resolution of 8km with the ECMWF 0.25° model. I use
ggrid to extract the domain of interest (Italy) and copygb to regrid WRF_NMM
fields over ECMWF grid.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I noted that the two grib files use different
parameter tables (in particular kpds5).</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV>
<DIV><FONT face=Arial size=2>ECMWF</FONT></DIV>
<DIV><FONT face=Arial
size=2>14596:100267650:d=07103112:<STRONG>T:kpds5=130</STRONG>:kpds6=100:kpds7=500:TR=1:P1=0:P2=0:TimeU=1:500
mb:anl:NAve=0<BR>14597:100274520:d=07103112:<STRONG>Z:kpds5=129</STRONG>:kpds6=100:kpds7=500:TR=1:P1=0:P2=0:TimeU=1:500
mb:anl:NAve=0<BR>14598:100281390:d=07103112:<STRONG>R:kpds5=157</STRONG>:kpds6=100:kpds7=500:TR=1:P1=0:P2=0:TimeU=1:500
mb:anl:NAve=0<BR>14599:100288260:d=07103112:<STRONG>U:kpds5=131</STRONG>:kpds6=100:kpds7=500:TR=1:P1=0:P2=0:TimeU=1:500
mb:anl:NAve=0<BR>14600:100295130:d=07103112:<STRONG>V:kpds5=132</STRONG>:kpds6=100:kpds7=500:TR=1:P1=0:P2=0:TimeU=1:500
mb:anl:NAve=0</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>WRF_NMM</FONT></DIV>
<DIV><FONT face=Arial
size=2>7297:42763434:d=07123100:<STRONG>TMP:kpds5=11</STRONG>:kpds6=100:kpds7=500:TR=0:P1=12:P2=0:TimeU=1:500
mb:12hr fcst:NAve=0</FONT></DIV>
<DIV><FONT face=Arial
size=2>7296:42756296:d=07123100:<STRONG>HGT:kpds5=7</STRONG>:kpds6=100:kpds7=500:TR=0:P1=12:P2=0:TimeU=1:500
mb:12hr
fcst:NAve=0<BR>7298:42768502:d=07123100:<STRONG>RH:kpds5=52</STRONG>:kpds6=100:kpds7=500:TR=0:P1=12:P2=0:TimeU=1:500
mb:12hr
fcst:NAve=0<BR>7299:42774812:d=07123100:<STRONG>UGRD:kpds5=33</STRONG>:kpds6=100:kpds7=500:TR=0:P1=12:P2=0:TimeU=1:500
mb:12hr
fcst:NAve=0<BR>7300:42780708:d=07123100:<STRONG>VGRD:kpds5=34</STRONG>:kpds6=100:kpds7=500:TR=0:P1=12:P2=0:TimeU=1:500
mb:12hr fcst:NAve=0</FONT></DIV></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Running GRID-STAT with a file CONFIG with this
indication:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>fcst_field[] = [ "11/P850" ];<BR>obs_field[]
= [ "130/P850" ];</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I obtain an error like:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>***WARNING***: process_scores() -> TMP/P850 not
found in GRIB file: 2007_ecm_cut.grb.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I imagine the problem was on the parameter table
used by ECMWF. </FONT><FONT face=Arial size=2>Is there a way to rewrite the GRIB
file using the usual table (that defined by WMO)? I tried to use copygb with the
option "-k" but with any results.</FONT></DIV>
<DIV><FONT face=Arial size=2>Any suggestions?</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Francesco Piani</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Dott. Ing. Francesco Piani<BR>Weather
Forecaster<BR>LaMMA Consortium - Hydrological Service Tuscany Region<BR>Via
Madonna del Piano 10<BR>Sesto Fiorentino (FI)<BR><A
href="mailto:piani@lamma.rete.toscana.it">piani@lamma.rete.toscana.it</A><BR><A
href="mailto:francesco.piani@regione.toscana.it">francesco.piani@regione.toscana.it</A></FONT></DIV></BODY></HTML>