LES experiments in ClimaAtmos

.buildkite/pipeline.yml

@@ -184,7 +184,6 @@ steps:
           slurm_mem: 16G
           slurm_gpus: 1
         soft_fail: true
-
   - group: "Plane Examples"
     steps:
       - label: ":computer: Agnesi linear hydrostatic mountain experiment (uniform)"

.gitignore

@@ -32,6 +32,8 @@ docs/src/config.md
 test/Manifest*.toml
 calibration/test/Manifest*.toml
 
+Artifacts.toml
+
 # ignore vscode artifacts
 *.vscode
 *.DS

config/default_configs/default_config.yml

@@ -145,6 +145,14 @@ vert_diff:
 hyperdiff:
   help: "Hyperdiffusion [`ClimaHyperdiffusion` (or `true`; default), `none` (or `false`)]"
   value: "CAM_SE"
+smagorinsky_lilly:
+  value: false
+c_smag: 
+  help: "Smagorinsky coefficient"
+  value: 0.2
+prandtl_turbulent_neutral:
+  help: "Turbulent Prandtl number for neutral stratification"
+  value: 0.333
 bubble:
   help: "Enable bubble correction for more accurate surface areas"
   value: true

config/model_configs/box_density_current_test.yml

@@ -1,16 +1,22 @@
+job_id: "box_density_current_test"
 dt_save_state_to_disk: "10secs"
 initial_condition: "DryDensityCurrentProfile"
 x_max: 51200.0
 z_elem: 45
-dt: "0.1secs"
-t_end: "10.0secs"
+dt: "0.3secs"
+t_end: "1000secs"
 discrete_hydrostatic_balance: true
-y_max: 6400.0
-y_elem: 15
+hyperdiff: "false"
+y_max: 51200.0
+y_elem: 45
 z_stretch: false
 x_elem: 45
 config: "box"
 z_max: 6400.0
+smagorinsky_lilly: true
+c_smag: 0.25
+prandtl_turbulent_neutral: 0.333
+netcdf_interpolation_num_points: [40, 40, 80]
 diagnostics:
-  - short_name: thetaa
-    period: 10secs
+  - short_name: [thetaa, wa, ua, ta, va, hfes, rhoa, rv, ha]
+    period: 30secs

config/model_configs/les_isdac_box.yml

@@ -15,6 +15,7 @@ hyperdiff: "false"
 apply_limiter: false
 smagorinsky_lilly: true
 c_smag: 0.20
+prandtl_turbulent_neutral: 0.333
 # time- and spatial discretization
 x_elem: 10
 x_max: 3.2e3
@@ -24,11 +25,13 @@ z_elem: 15
 z_max: 2.5e3
 z_stretch: false
 rayleigh_sponge: true
-toml: [toml/isdac_box.toml]  # sponge height
+toml: [toml/les_isdac.toml]  # sponge height
 ode_algo: "SSPKnoth"
 dt: "0.05secs" 
 t_end: "2mins" 
 dt_cloud_fraction: "10mins"
+restart_file: "/groups/esm/hervik/climaatmos-les-artifacts/les_isdac_day0.0.hdf5"
+# diagnostics
 dt_save_state_to_disk: "1mins"
 netcdf_interpolation_num_points: [30, 30, 150]
 diagnostics: