Verify WRF Forecast using MET Ruifang Li MMM

Verify WRF Forecast using MET Ruifang Li MMM Lir@ucar. edu

METv 1. 1 Installation n n /ptmp/lir/wrf/METv 1. 1 (bluefire) /users/lir/wrf/met/METv 1. 1 -gnu (Linux)

MET Components

MET Tools used n PB 2 NC n n Point-Stat n n The PB 2 NC tool is used to stratify the contents of an input Prep. Bufr point observation file and reformat it into Net. CDF format for use by the Point-Stat tool The Point-Stat tool provides verification statistics forecasts at observation points, as opposed to over gridded analyses. VSDB-Analysic n The VSDB-Analysis tool ties together results from the Point-Stat and Grid-Stat tools by providing summary statistical information and a way to filter their VSDB output files.

Scripts Path n n /ptmp/lir/wrf_test/wrf_em/verification/nu dge/scripts I used some of the scripts which were wroten by our visitor Xiaqun Zhou. Thanks Xiaoqun

WRF Forcast and Obs Data n Obs Data (Prebufr Format) n n n � /ptmp/liuz/t 8_2008_45 km/ob/2007081500 /ptmp/liuz/t 8_2008_45 km/ob/2007081506. . . /ptmp/liuz/t 8_2008_45 km/ob/2007081700

WRF Forcast and Obs Data n grid nudging experiment n n n wrfout_d 01_2007 -08 -14_12: 00 wrfout_d 01_2007 -08 -15_00: 00. wrfout_d 01_2007 -08 -17_00: 00 control experiment (without grid nudging) n n n wrfout_d 01_2007 -08 -15_00: 00. wrfout_d 01_2007 -08 -17_00: 00

WRF Postprocessing n run_wrfpost_cyc. ksh n Change WRF output in netcdf format to grib format.

Running Process Obs in prebufr point_stat_6000 00 L_2007081�� 5_000000 V. vsdb Obs in prebufr Pb 2 nc. ksh � point_stat_600000 L _2007081�� 5_0000 06 V. vsdb point_stat. k sh fcst in netcdf fcst in grib run_wrfpost _cyc. ksh� . . point_stat_6000 00 L_2007081�� 7 _000000 V. vsdb cnt_point_stat_TMP. vsdb cnt_point_stat_SPFH. vsdb_anal ysis_point. ksh cnt_point_stat_UGRD. vsdb cnt_point_stat_VGRD. vsdb cnt_point_stat_WIND. vsdb

One record in cnt_point_stat_TMP. vsdb V 1. 1 WRF 180000 20070815_060000 ADPSFC FULL CNT/0. 050 TMP Z 02 = 1221 300. 77947 300. 55386 301. 00507 300. 55137 300. 93672 4. 02209 3. 86865 4. 18829 3. 85698 4. 23203 297. 25917 297. 06136 297. 45698 297. 04889 297. 42627 3. 52660 3. 39207 3. 67234 3. 35312 3. 70185 0. 77959 0. 75658 0. 80067 0. 74965 0. 79329 0. 68545 0. 51188 1221 717 1134 3. 52030 3. 37731 3. 66328 3. 43463 3. 64954 2. 54916 2. 45191 2. 65450 2. 46988 2. 66370 1. 01184 1. 01155 1. 01228 3. 70852 3. 62637 3. 83187 18. 88538 18. 48708 19. 92617 6. 49289 6. 09545 7. 08971 4. 34573 4. 29965 4. 46379 0. 19351 -0. 00492 0. 40473 1. 81395 1. 62959 1. 98116 3. 47373 3. 37264 3. 66670 5. 37370 5. 30947 5. 54463 6. 71369 6. 68062 6. 88146 UW_MEAN 1 Code to pick up statistic values: pick 2 ncl_control. f 90 pick 2 ncl_nudge. f 90 Get cnt_point_stat_TMP. txt for reading And cnt_point_stat_TMP. bin for plotting (NCL)

Final Result n vi cnt_point_stat_TMP. txt n . ADPUPA (Up air observations) RMSE n n n P 1000 -1000 P 925 -875 P 875 -825 P 725 -675 P 625 -575 P 525 -475 P 425 -375 P 325 -275 P 225 -175 P 125 -75 number 00 00 1. 69 19. 0. 84 52. 0. 75 63. 0. 56 93. 1. 16 105. 0. 69 83. 0. 67 96. 0. 84 89. 1. 31 77. 4. 96 224. RMSE number 06 06 ……………… 48 48 -999. 0. 67 2. 0. 36 1. 0. 23 2. 1. 06 3. 0. 31 1. 0. 89 3. 0. 66 2. 0. 64 2. 5. 47 5.

Final Result n vi cnt_point_stat_TMP. txt n . ADPSFC (surface observations: in land) n RMSE n 00 00 06 06 ……………… 48 48 2. 86 1284. 3. 48 1221. 1. 95 4. 0. 03 1. 2. 34 73. 5. 08 73. 2. 05 36. 1. 98 26. -999. 00 -999. n n n Z 02 1000 -1000 P 925 -875 P 875 -825 P 725 -675 P 625 -575 P 525 -475 P 425 -375 P 325 -275 P 225 -175 P 125 -75 number RMSE number

Final Result n vi cnt_point_stat_TMP. txt n . SFCSHP (surface observations: ship) RMSE n 00 00 06 Z 02 1. 36 148. 1. 27 1000 -1000 P 925 -875 P 875 -825 P 725 -675 P 625 -575 P 525 -475 P 425 -375 P 325 -275 P 225 -175 P 125 -75 -999. 00 -999. 00 n n n number RMSE 06 number ……………… 48 48 152. -999. 00 -999. 00 -999.

Final Result n vi cnt_point_stat_TMP. txt n AIRCFT(aircraft observation) n RMSE n 00 00 12 12 ……………… 48 48 -999. 00 -999. -999. 00 -999. 0. 48 1. -999. 00 -999. n n n Z 02 1000 -1000 P 925 -875 P 875 -825 P 725 -675 P 625 -575 P 525 -475 P 425 -375 P 325 -275 P 225 -175 P 125 -75 number RMSE number

Specific Humidity

Temperature

U component

V component

Wind

Appendix A n The steps to run on bluefire 1 . Run wpp to change wrf output. in netcdf format to grib format. n n 2 . Run point-stat to get matched fcst and obs pairs n n 3 n /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/vsdb_analysis_point_control. ksh /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/vsdb_analysis_point_nudge. ksh . Get RMSE and BIAS info of matched fcst and obs pairs n n 5 /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/point_stat_control. ksh /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/point_stat_nudge. ksh . Run vsdb tool to get summary info of matched fcst and obs pairs n 4 /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/run_wrfpost_cyc_control. ksh /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/run_wrfpost_cyc_nudge. ksh /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/pick 2 ncl_control. f 90 . Plot figures n /ptmp/lir/wrf_test/wrf_em/verification/nudge/scripts/plot/

� n Questions? n Thanks!