Hi John,<br>I run point_stat tool and did two test with different fcst_thresh[] in config file. I expect the different MPR output, but they create same MPR. I could not figure out why. In my understanding, if fcst_thresh eq 0, fcst value in MPR should be 0. if  fcst_thresh le 100 fcst value in MPR should be le 100, right?<br>
<br>Here are two fcst_thresh:  <br><br>fcst_field[] = [&quot;UGRD/Z10&quot;,&quot;UGRD/P1000&quot;,&quot;UGRD/P925&quot;,&quot;UGRD/P850&quot;,&quot;UGRD/P700&quot;,&quot;UGRD/P500&quot;,&quot;UGRD/P400&quot;,&quot;UGRD/P300&quot;,&quot;UGRD/P250&quot;,&quot;UGRD/P200&quot;,&quot;UGRD/P150&quot;,&quot;UGRD/P100&quot;,&quot;UGRD/P70&quot;,&quot;UGRD/P50&quot;,&quot;UGRD/P30&quot;,&quot;UGRD/P20&quot;,&quot;UGRD/P10&quot; ];<br>
<br>//fcst_thresh[] = [ &quot;le100&quot;,
&quot;le10&quot;,&quot;le10&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot;,
&quot;le100&quot;,&quot;le100&quot;, &quot;le100&quot;,&quot;le100&quot;,&quot;le100&quot; ];<br>
fcst_thresh[] = [ &quot;eq0&quot;, &quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;le10&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot;, &quot;eq0&quot;,&quot;eq0&quot;, &quot;eq0&quot;,&quot;eq0&quot;,&quot;eq0&quot; ];<br>
<br>Here is MPR file:<br><br>/ptmp/lir/test/2009_t8_15km/verify_met2.0/test% ls -l /ptmp/lir/data/2009_t8_15km/verify_met/gsi/point_stat/UGRD/2007081512/<br>
total 14208<br>
-rw-r--r--    1 lir      ncar         786984 Aug 20 16:05 point_stat_000000L_20070815_120000V.stat<br>
-rw-r--r--    1 lir      ncar         652568 Aug 20 16:05 point_stat_000000L_20070815_120000V_mpr.txt<br>
-rw-r--r--    1 lir      ncar         698412 Aug 20 16:05 point_stat_120000L_20070816_000000V.stat<br>
-rw-r--r--    1 lir      ncar         579124 Aug 20 16:05 point_stat_120000L_20070816_000000V_mpr.txt<br>
-rw-r--r--    1 lir      ncar         761574 Aug 20 16:05 point_stat_240000L_20070816_120000V.stat<br>
-rw-r--r--    1 lir      ncar         631498 Aug 20 16:05 point_stat_240000L_20070816_120000V_mpr.txt<br>
-rw-r--r--    1 lir      ncar         740520 Aug 20 16:05 point_stat_360000L_20070817_000000V.stat<br>
-rw-r--r--    1 lir      ncar         614040 Aug 20 16:05 point_stat_360000L_20070817_000000V_mpr.txt<br>
-rw-r--r--    1 lir      ncar         791340 Aug 20 16:05 point_stat_480000L_20070817_120000V.stat<br>
-rw-r--r--    1 lir      ncar         656180 Aug 20 16:05 point_stat_480000L_20070817_120000V_mpr.txt<br>
<br>Script:<br>/ptmp/lir/test/2009_t8_15km/verify_met2.0/test% point_stat_gsi.ksh <br>*** Running POINT_STAT on WRF OUTPUT ***<br>***          UGRD                ***<br>GSL_RNG_TYPE=mt19937<br>GSL_RNG_SEED=2412285258<br>
Forecast File: /ptmp/lir/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000<br>Climatology File: none<br>Configuration File: ./config/PointStatConfig_GSI.UGRD<br>Observation File: /ptmp/lir/data/2009_t8_15km/verify_met2.0/ob_nc/2007081512/ob.nc.qc2<br>
<br>----------------------------------------<br><br>Reading records for UGRD/Z10.<br>For UGRD/Z10 found 1 forecast levels and 0 climatology levels.<br><br clear="all"><br>Thanks for your support,<br>-- <br>Ruifang <br><br>
Mesoscale &amp; Microscale Meteorology Division<br>Phone: 303-497-8938<br>Office:  FL3-3085<br>Email:  <a href="mailto:lir@ucar.edu">lir@ucar.edu</a><br>