
Fuzzy Verification of Precipitation Using VAST 3.0 Software Package
Discover how VAST 3.0, a software package composed of Bash, Fortran 90, and R scripts, facilitates fuzzy verification of precipitation through GRIB files. Learn about the required software, input file specifications, and the purpose of VAST 3.0 in generating heatmap plots for data analysis.
Download Presentation

Please find below an Image/Link to download the presentation.
The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author. If you encounter any issues during the download, it is possible that the publisher has removed the file from their server.
You are allowed to download the files provided on this website for personal or commercial use, subject to the condition that they are used lawfully. All files are the property of their respective owners.
The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author.
E N D
Presentation Transcript
VAST 3.0 VAST 3.0 QUICK OVERVIEW QUICK OVERVIEW Naima Vela1 1Arpa Piemonte, Torino, Italy September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
WHAT IS VAST 3.0 WHAT IS VAST 3.0 VAST 3.0 VAST 3.0 is is a software package a software package composed composed by Bash, Fortran 90 and R scripts by Bash, Fortran 90 and R scripts Main Components: Main Components: 3 bash scripts, 3 Fortran 90 scripts (and 1 README pdf) 3 bash scripts, 3 Fortran 90 scripts (and 1 README pdf) Purpose Purpose: : Fuzzy Fuzzy verification verification of of precipitation precipitation INPUT: INPUT: GRIB GRIB files (forecast and observation) files (forecast and observation) OUTPUT: OUTPUT: heatmap plots heatmap plots September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
REQUIRED SOFTWARE REQUIRED SOFTWARE The grib files are elaborated using CDO (https://code.mpimet.mpg.de/projects/cdo/wiki/Cdo#Documentation) CDO functions The ECMWF package grib_api (https://www.ecmwf.int/en/elibrary/80529-grib-api , https://confluence.ecmwf.int/display/ECC ) grib_api or ecCodes ecCodes is needed by the user to elaborate the grib files Part of the fuzzy elaboration is run through Fortran 90 scripts compiled using the gfortran (https://gcc.gnu.org/onlinedocs/gfortran/) gfortran compiler The plots are created by VAST using R R (https://www.r-project.org/) A version 3.0 3.0 or higher is required. September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
INPUT FILES INPUT FILES VAST 3.0 works with GRIB files VAST 3.0 works with GRIB files The GRIB files must contain only the precipitation variable The GRIB files must contain only the precipitation variable One file per day One file per day Hourly messages Hourly messages Area covered by the GRIB must be wider than the one chosen for verification Area covered by the GRIB must be wider than the one chosen for verification FORECAST FORECAST Run 00 Run 00 Cumulated precipitation since the Cumulated precipitation since the beginning of the run beginning of the run Naming: Naming: mmm00_YYYYMMDD.grb (mmm=3 letter code for model name) (mmm=3 letter code for model name) OBSERVATION OBSERVATION Hourly precipitation in each message Hourly precipitation in each message Naming: Naming: YYYYMMDD.grb September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
STRUCTURE STRUCTURE VAST 3.0 is composed by 3 main sections, each one driven by a Bash script VAST 3.0 is composed by 3 main sections, each one driven by a Bash script Reads the input GRIB files and produces partial Reads the input GRIB files and produces partial results used by the other two scripts (TXT) results used by the other two scripts (TXT) VAST_3_0_main.sh Reads the partial files produced by the main script Reads the partial files produced by the main script and produces FSS plots for the required period and produces FSS plots for the required period make_fss_period_VAST_3_0.sh contingency_statistics_VAST_3_0.sh Reads the partial files produced by the main script Reads the partial files produced by the main script and produces plots of the chosen categorical score and produces plots of the chosen categorical score A Fortran script is associated to each Bash section to perform the elaboration. A Fortran script is associated to each Bash section to perform the elaboration. The configuration files needed by the Fortran scripts are produced and replaced The configuration files needed by the Fortran scripts are produced and replaced by the Bash scripts according to the user s requests. by the Bash scripts according to the user s requests. All Bash scripts require the user to input a list of option to produce the verification. All Bash scripts require the user to input a list of option to produce the verification. All options are mandatory (no All options are mandatory (no dafault dafault). ). September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
VAST_3_0_MAIN.SH / 1 VAST_3_0_MAIN.SH / 1 OPTIONS: OPTIONS: -s [YYYYMMDD] starting date -e [YYYYMMDD] ending date -m [aaa] -t [1t/3t] -d [d0/d1] -l [LT/24] -D [a] -a [a2/a3/new] verification area. a2-> Common area 2, 2.8 Km a3-> Common area 3, 2.8 Km new-> different area (insert required parameters in the script) -r [Y/N] remove or keep partial files (heavy on disk, not needed at the end of the verification) MODELS: MODELS: COSMO 2I ( COSMO 2I (- -m c22), ICON IT ( m c22), ICON IT (ico COSMO CE PL ( COSMO CE PL (ccp ccp), ICON PL ( ), ICON PL (ipl code of the model to be verified 2D verification/ 3D verification first 24h verification/second 24h verification Lead time (only compatible with 1t) /24h with 3 hourly cumulation single character to identify the verification ico), ICON 2I (i22), COSMO 1E (c1e), ), ICON 2I (i22), COSMO 1E (c1e), ipl), ICON D2 (id2), ICON IL (il4), ICON GR ( ), ICON D2 (id2), ICON IL (il4), ICON GR (igr igr) ) September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
VAST_3_0_MAIN.SH / 2 VAST_3_0_MAIN.SH / 2 OUTPUT: OUTPUT: The script produces TXT files with partial results to be used to create FSS The script produces TXT files with partial results to be used to create FSS and/or categorical scores aggregation and/or categorical scores aggregation For Fraction Skill Score calculation: For Fraction Skill Score calculation: fss_YYYYMMDD_mmmtt_dd_3hd_ts.txt (mmm=model, tt=1t/3t, dd=d0/d1, d=suffix, ts=03/06/09/ ) sample: sample: fss_20240304_c221t_d0_3h2_18.txt For categorical score calculation: For categorical score calculation: contingency_YYYYMMDD_mmmtt_dd_3hd_ts.txt sample: sample: contingency_20240304_c221t_d0_3h2_18.txt 2D or 3D verification: 2D or 3D verification: 1t 1t - -> 8 files > 8 files/day (LT or 24 option available) /day (LT or 24 option available) 3t 3t - -> 1 file/day (only 24 option available) > 1 file/day (only 24 option available) September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
MAKE_FSS_PERIOD_VAST_3_0.SH / 1 MAKE_FSS_PERIOD_VAST_3_0.SH / 1 OPTIONS: OPTIONS: -s [YYYYMMDD] starting date -e [YYYYMMDD] ending date -m [aaa] code of the model to be verified -t [1t/3t] 2D verification/ 3D verification -d [d0/d1] first 24h verification/second 24h verification -D [a] single character to identify the verification -a [a2/a3/new] verification area. a2-> Common area 2 a3-> Common area 3 new-> different area (the procedure will ask to insert the long name) -g [num] side of the smallest box (Km) (same of the one used in the main script) -l [LT/24] Lead time (only compatible with 1t) /24h with 3 hourly cumulation MODELS: MODELS: COSMO 2I ( COSMO 2I (- -m c22), ICON IT ( m c22), ICON IT (ico COSMO CE PL ( COSMO CE PL (ccp ccp), ICON PL ( ico), ICON 2I (i22), COSMO 1E (c1e), ), ICON 2I (i22), COSMO 1E (c1e), ), ICON PL (ipl ipl), ICON D2 (id2), ICON IL (il4), ICON GR ( ), ICON D2 (id2), ICON IL (il4), ICON GR (igr igr) ) September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
MAKE_FSS_PERIOD_VAST_3_0.SH / 3 MAKE_FSS_PERIOD_VAST_3_0.SH / 3 OUTPUT: OUTPUT: The script produces the Fractions Skill Score heatmap plots according to the The script produces the Fractions Skill Score heatmap plots according to the user s selection user s selection - -l 24 l 24 September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
MAKE_FSS_PERIOD_VAST_3_0.SH / 3 MAKE_FSS_PERIOD_VAST_3_0.SH / 3 OUTPUT: OUTPUT: The script produces the Fractions Skill Score heatmap plots according to the The script produces the Fractions Skill Score heatmap plots according to the user s selection user s selection - -l LT l LT September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
MAKE_FSS_PERIOD_VAST_3_0.SH / 3 MAKE_FSS_PERIOD_VAST_3_0.SH / 3 OUTPUT: OUTPUT: The script produces the Fractions Skill Score heatmap plots according to the The script produces the Fractions Skill Score heatmap plots according to the user s selection user s selection - -l LT l LT September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
CONTINGENCY_STATISTICS_VAST_3_0.SH / 1 CONTINGENCY_STATISTICS_VAST_3_0.SH / 1 OPTIONS: OPTIONS: -s [YYYYMMDD] starting date -e [YYYYMMDD] ending date -m [aaa] code of the model to be verified -t [1t/3t] 2D verification/ 3D verification -d [d0/d1] first 24h verification/second 24h verification -D [a] single character to identify the verification -a [a2/a3/new] verification area. a2-> Common area 2 a3-> Common area 3 new-> different area (the procedure will ask to insert the long name) -g [num] side of the smallest box (Km) (same of the one used in the main script) -l [LT/24] Lead time (only compatible with 1t) /24h with 3 hourly cumulation -i [FBI/ETS/FAR/POD] Categorical score to be plotted MODELS: MODELS: COSMO 2I ( COSMO 2I (- -m c22), ICON IT ( m c22), ICON IT (ico COSMO CE PL ( COSMO CE PL (ccp ccp), ICON PL ( ico), ICON 2I (i22), COSMO 1E (c1e), ), ICON 2I (i22), COSMO 1E (c1e), ), ICON PL (ipl ipl), ICON D2 (id2), ICON IL (il4), ICON GR ( ), ICON D2 (id2), ICON IL (il4), ICON GR (igr igr) ) September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
CONTINGENCY_STATISTICS_VAST_3_0.SH / 2 CONTINGENCY_STATISTICS_VAST_3_0.SH / 2 OUTPUT: OUTPUT: The script produces the heatmap plots for the chosen score according to the The script produces the heatmap plots for the chosen score according to the user s selection user s selection - -l 24 l 24 September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
CONTINGENCY_STATISTICS_VAST_3_0.SH / 3 CONTINGENCY_STATISTICS_VAST_3_0.SH / 3 OUTPUT: OUTPUT: The script produces the heatmap plots for the chosen score according to the The script produces the heatmap plots for the chosen score according to the user s selection user s selection - -l LT l LT September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach
THANK YOU FOR YOUR ATTENTION! Please Please send send all naima.vela@arpa.piemonte.it naima.vela@arpa.piemonte.it all your your comments comments and and questions questions to: to: September the 2 September the 2nd nd 2024 2024 - - 26 26th th COSMO General Meeting COSMO General Meeting - - Offenbach Offenbach