Skip to content

Ensemble products per species #78

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
pirmink wants to merge 5 commits into
base: main
Choose a base branch
from
Open

Ensemble products per species #78

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
4592a21
Plot some ensemble statistics per species.
pirmink Nov 14, 2025
bac0e6a
Correct error in preset for shapefiles.
pirmink Nov 25, 2025
722880d
Add species_id to filename to discern plot for different species.
pirmink Nov 26, 2025
c684a4e
Bump patch level of version.
pirmink Dec 1, 2025
35ac3db
Correct naming of percentile shapefiles, consistent naming of concent…
pirmink Dec 1, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion pyproject.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
[tool.poetry]
name = "pyflexplot"
version = "2.2.4"
version = "2.2.5"
description = "PyFlexPlot - Visualize and post-process FLEXPART dispersion simulation results stored in NetCDF format"
authors = ["Stefan Ruedisuehli <stefan.ruedisuehli@env.ethz.ch>", "Pirmin Kaufmann <pirmin.kaufmann@meteoswiss.ch>"]
packages = [
Expand Down
12 changes: 11 additions & 1 deletion src/pyflexplot/data/presets/opr/icon-ch2-eps/all_pdf.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@ ens_member_id = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 1
model = "ICON-CH2-EPS"
lang = "de"
outfile_time_format = "%Y%m%dT%H%M"
combine_species = true
combine_species = false

[_base._concentration]
plot_variable = "concentration"
Expand All @@ -29,21 +29,25 @@ plot_variable = "affected_area"
integrate = true
level = 0
time = -1
combine_species = true

[_base._concentration."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_probability_{domain}.pdf"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._concentration._integr."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_probability_{domain}.pdf"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._tot_deposition."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_probability_{domain}.pdf"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._affected_area."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_probability_{domain}.pdf"
Expand All @@ -55,6 +59,7 @@ outfile = "{release_start}_{release_site}_FLEXPART-{model}_probability_{domain}.
plot_type = "multipanel"
multipanel_param = "time"
time = [2, 4, 8, -1]
combine_species = true

[_base._concentration._integr."_pctl95+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_percentile_{domain}.pdf"
Expand Down Expand Up @@ -93,9 +98,11 @@ time = 0

[_base._concentration._cloud."_arrival+"]
ens_variable = "cloud_arrival_time"
combine_species = true

[_base._concentration._cloud."_departure+"]
ens_variable = "cloud_departure_time"
combine_species = true

[_base._concentration._integr."_min+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_cloud_statistics_{domain}.pdf"
Expand Down Expand Up @@ -136,14 +143,17 @@ ens_variable = "mean"
[_base._concentration._integr."_median+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_cloud_statistics_{domain}.pdf"
ens_variable = "median"

[_base._tot_deposition."_median+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_cloud_statistics_{domain}.pdf"
ens_variable = "median"

[_base._affected_area."_median+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_cloud_statistics_{domain}.pdf"
ens_variable = "median"

["**".full]
domain = "full"

["**".ch]
domain = "ch"
10 changes: 8 additions & 2 deletions src/pyflexplot/data/presets/opr/icon-ch2-eps/all_png.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,8 +7,8 @@ ens_member_id = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 1
model = "ICON-CH2-EPS"
lang = "de"
outfile_time_format = "%Y%m%dT%H%M"
combine_species = true
outfile = "forecast-iconch2eps-dispersion~{base_time}~{plot_variable}-{ens_variable}~{domain}~{release_site}~1~{time_step}.png"
combine_species = false
outfile = "forecast-iconch2eps-dispersion~{base_time}~{plot_variable}-{ens_variable}~{domain}~{release_site}~{species_id}~{time_step}.png"

[_base._concentration]
plot_variable = "concentration"
Expand All @@ -30,18 +30,22 @@ plot_variable = "affected_area"
integrate = true
level = 0
time = -1
combine_species = true

[_base._concentration."_probability+"]
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._concentration._integr."_probability+"]
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._tot_deposition."_probability+"]
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._affected_area."_probability+"]
ens_variable = "probability"
Expand All @@ -66,9 +70,11 @@ time = 0

[_base._concentration._cloud."_arrival+"]
ens_variable = "cloud_arrival_time"
combine_species = true

[_base._concentration._cloud."_departure+"]
ens_variable = "cloud_departure_time"
combine_species = true

[_base._concentration._integr."_min+"]
ens_variable = "minimum"
Expand Down
38 changes: 22 additions & 16 deletions src/pyflexplot/data/presets/opr/icon-ch2-eps/all_shp.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@ ens_member_id = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 1
model = "ICON-CH2-EPS"
lang = "de"
outfile_time_format = "%Y%m%dT%H%M"
combine_species = true
combine_species = false

[_base._concentration]
plot_variable = "concentration"
Expand All @@ -29,39 +29,43 @@ plot_variable = "affected_area"
integrate = true
level = 0
time = -1
combine_species = true

[_base._concentration."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit-Konzentration_{species_id}_{time_step}.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit_Konzentration_{species_id}_{time_step}.shp"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._concentration._integr."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit-integrierte-Luftaktivitaet_{species_id}.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._tot_deposition."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit-Deposition_{species_id}.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit_Deposition_{species_id}.shp"
ens_variable = "probability"
ens_param_thr = 0
combine_species = true

[_base._affected_area."_probability+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit-beaufschl_Gebiet.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wahrscheinlichkeit_beaufschl_Gebiet.shp"
ens_variable = "probability"
ens_param_thr = 0

[_base._concentration._integr."_pctl95+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "percentile"
ens_param_pctl = 95

[_base._tot_deposition."_pctl95+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil_Deposition_{species_id}.shp"
ens_variable = "percentile"
ens_param_pctl = 95

[_base._affected_area."_pctl95+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_95-Perzentil_beaufschl_Gebiet.shp"
ens_variable = "percentile"
ens_param_pctl = 95

Expand All @@ -73,53 +77,55 @@ time = 0
[_base._concentration._cloud."_arrival+"]
ens_variable = "cloud_arrival_time"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wolkenankunftszeit.shp"
combine_species = true

[_base._concentration._cloud."_departure+"]
ens_variable = "cloud_departure_time"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Wolkenabzugszeit.shp"
combine_species = true

[_base._concentration._integr."_min+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Minimum_integrierte-Luftaktivitaet.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Minimum_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "minimum"

[_base._tot_deposition."_min+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Minimum_Deposition.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Minimum_Deposition_{species_id}.shp"
ens_variable = "minimum"

[_base._affected_area."_min+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Minimum_beaufschl_Gebiet.shp"
ens_variable = "minimum"

[_base._concentration._integr."_max+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Maximum_integrierte-Luftaktivitaet.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Maximum_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "maximum"

[_base._tot_deposition."_max+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Maximum_Deposition.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Maximum_Deposition_{species_id}.shp"
ens_variable = "maximum"

[_base._affected_area."_max+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Maximum_beaufschl_Gebiet.shp"
ens_variable = "maximum"

[_base._concentration._integr."_mean+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Mittel_integrierte-Luftaktivitaet.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Mittel_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "mean"

[_base._tot_deposition."_mean+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Mittel_Deposition.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Mittel_Deposition_{species_id}.shp"
ens_variable = "mean"

# [_base._affected_area."_mean+"] # for affected area, mean is same area as max -> skip mean
# outfile = "{release_start}_{release_site}_FLEXPART-{model}_Mittel_beaufschl_Gebiet.shp"
# ens_variable = "mean"

[_base._concentration._integr."_median+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Median_integrierte-Luftaktivitaet.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Median_integrierte-Luftaktivitaet_{species_id}.shp"
ens_variable = "median"

[_base._tot_deposition."_median+"]
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Median_Deposition.shp"
outfile = "{release_start}_{release_site}_FLEXPART-{model}_Median_Deposition_{species_id}.shp"
ens_variable = "median"

[_base._affected_area."_median+"]
Expand Down