From 76baf3e62c611e793b626c733671c456f70b17ec Mon Sep 17 00:00:00 2001
From: "Takumi.Matsunobu" <Takumi.Matsunobu@physik.uni-muenchen.de>
Date: Wed, 16 Aug 2023 15:33:45 +0200
Subject: [PATCH] perturbations can be configured by conf.yml
---
conf/real-from-d2-ana/simulation.yml | 4 ++++
.../real-from-d2-ana/prepare_namelist.py | 22 +++++++++++++++++--
2 files changed, 24 insertions(+), 2 deletions(-)
diff --git a/conf/real-from-d2-ana/simulation.yml b/conf/real-from-d2-ana/simulation.yml
index ed8083f..04f9537 100644
--- a/conf/real-from-d2-ana/simulation.yml
+++ b/conf/real-from-d2-ana/simulation.yml
@@ -33,3 +33,7 @@ simulation:
parent_folder: /dss/dsskcsfs01/pn34ca/pn34ca-dss-0007/DATA-OP/
+
+ uncertanty:
+ psp: False
+ mpp: False # must be capitalised
\ No newline at end of file
diff --git a/templates/real-from-d2-ana/prepare_namelist.py b/templates/real-from-d2-ana/prepare_namelist.py
index 059c5d7..2eb9e2d 100644
--- a/templates/real-from-d2-ana/prepare_namelist.py
+++ b/templates/real-from-d2-ana/prepare_namelist.py
@@ -56,6 +56,23 @@ master_namelist_string = """
# Compute difference in seconds
checkpoint_time = int((Chunk_END_DATE - Chunk_START_DATE).total_seconds())
+# Turn on or off PSP perturbations
+if %simulation.uncertanty.psp% :
+ inwp_psp = 1
+else:
+ inwp_psp = 0
+
+# Microphysics Parameter Perturbations
+ccn_type = 7
+cloudnue = 1
+
+MPPCOMB = ((7,1), (6, 0), (6, 8), (9, 0), (9,8))
+IBC_MPP_COMBLIST = [0, 1, 4, 4, 0, 3, 4, 1, 2, 0, 2, 0, 1, 3, 1, 3, 2, 4, 2, 3]
+if %simulation.uncertanty.mpp% :
+ ccn_type = MPPCOMB[IBC_MPP_COMBLIST[int(MEMBER[1:])-1]][0]
+ cloudnue = MPPCOMB[IBC_MPP_COMBLIST[int(MEMBER[1:])-1]][1]
+
+
# TODO: Is that really necessary?
# Add 10 minutes to allow the model to write the restarts
# Chunk_END_DATE = Chunk_END_DATE + timedelta(minutes=10)
@@ -80,14 +97,15 @@ atmosphere_namelist_replacements = {
"latbc_boundary_grid": "%simulation.lateral_boundary_grid_filename%",
},
"nwp_phy_nml": {
+ "inwp_psp" : inwp_psp,
"psp_rnd_seed": int(MEMBER[1:]),
#"psp_apply_ddt_core": ".false.", # True if vertical wind perturbations should be applied at the explicit dynamical core timestep
"inwp_gscp": 4 , # 1: default, 2: graupel scheme for convection-permitting scales
},
"twomom_mcrph_nml": {
- "ccn_type" : 7,
- "cloudnue" : 0,
+ "ccn_type" : ccn_type,
+ "cloudnue" : cloudnue,
"rainnue" : 0,
},
--
GitLab