Skip to content

Import of high order turbulence closure scheme #734

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
Shusuke Nishimoto (mo-snishimoto) wants to merge 12 commits into
base: main
Choose a base branch
from
Open

Import of high order turbulence closure scheme #734

Show file tree
Hide file tree
Changes from 8 commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
41fff2c
Import 3DTE scheme (UM source codes are just copies of um@vn14.2).
mo-snishimoto Jul 29, 2026
aed80b9
Show minimum necessary modifications to UM's source code.
mo-snishimoto Jul 29, 2026
6a5b436
Convert real_umphys into r_bl & Apply LFRic code styling script (lfri…
mo-snishimoto Jul 29, 2026
dfe0b36
Additional modification so that 3DTE scheme can work in single precis…
mo-snishimoto Jul 29, 2026
1606343
Merge adv and con tracer into pos tracer and create gen tracer.
mo-snishimoto Jul 29, 2026
672999a
Implement full 3DTE Mk1, including TKE based leonard term.
mo-snishimoto Jul 29, 2026
0e029dc
Add new tasks in rose-stem suite
mo-snishimoto Aug 11, 2026
5e3914b
Update KGO.
mo-snishimoto Aug 25, 2026
b249835
Review comment
mo-snishimoto Aug 25, 2026
93c9d93
Sci/Tech review comment.
mo-snishimoto Sep 8, 2026
58bda9b
Merge upstream/main into the branch and resolve conflicts.
mo-snishimoto Sep 9, 2026
5ef66cc
Update KGO.
mo-snishimoto Sep 9, 2026
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 applications/jedi_lfric_tests/example/configuration.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -389,7 +389,7 @@ ffsl_splitting=5*1,
ffsl_unity_3d=.false.,
ffsl_vertical_order=2,2,1,2,2,
field_names='density','potential_temperature','wind','moisture',
'con_tracer',
'pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=5*1,
Expand Down
2 changes: 1 addition & 1 deletion applications/jedi_lfric_tests/example_forecast/configuration.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -337,7 +337,7 @@ ffsl_splitting=5*1,
ffsl_unity_3d=.false.,
ffsl_vertical_order=2,2,1,2,2,
field_names='density','potential_temperature','wind','moisture',
'con_tracer',
'pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=5*1,
Expand Down
2 changes: 1 addition & 1 deletion applications/jules/example/configuration.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -804,7 +804,7 @@ ffsl_splitting=6*1,
ffsl_unity_3d=.true.,
ffsl_vertical_order=2,2,1,2,2,2,
field_names='density','potential_temperature','wind','moisture',
'adv_tracer','con_tracer',
'gen_tracer','pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=6*2,
Expand Down
3 changes: 2 additions & 1 deletion applications/lfric_atm/example/configuration.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,7 @@ prime_mesh_name='dynamics',
topology='fully_periodic',
/
&blayer
bl_scheme='9c',
Comment thread
Adrian-Lock marked this conversation as resolved.
Outdated
a_ent_2=0.056,
a_ent_shr=1.6,
bl_levels=50,
Expand Down Expand Up @@ -1000,7 +1001,7 @@ ffsl_splitting=6*1,
ffsl_unity_3d=.true.,
ffsl_vertical_order=2,2,1,2,2,2,
field_names='density','potential_temperature','wind','moisture',
'adv_tracer','con_tracer',
'gen_tracer','pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=6*2,
Expand Down
19 changes: 18 additions & 1 deletion applications/lfric_atm/metadata/field_def_diags.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -664,7 +664,6 @@
<field id="turbulence__zhsc" name="zhsc" long_name="height_at_top_of_decoupled_stratocumulus_layer" unit="m" domain_ref="face" />
<field id="turbulence__oblen" name="oblen" long_name="obukhov_length" unit="m" domain_ref="face" />
<field id="turbulence__bl_weight_1dbl" name="bl_weight_1dbl" long_name="Blending_weight_in_boundary_layer" unit="1" domain_ref="face" />
<field id="turbulence__tke" name="tke_bl" long_name="turbulent_kinetic_energy" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="turbulence__wvar" name="wvar" long_name="vertical_velocity_variance" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="turbulence__dsldzm" name="dsldzm" long_name="liquid_temperature_gradient" unit="kg m-1" grid_ref="full_level_face_grid" />
<field id="turbulence__mix_len_bm" name="mix_len_bm" long_name="Turb length-scale for bimodal" unit="m" grid_ref="full_level_face_grid" />
Expand All @@ -673,6 +672,7 @@
<field id="zoomed__rib" field_ref="turbulence__gradrinr" grid_ref="full_level_0_face_zoom" />
<field id="turbulence__rib" name="rib" long_name="lowest_layer_bulk_richardson_number" unit="1" domain_ref="face"> zoomed__rib </field>
<field id="turbulence__rhokh" name="rhokh_bl" long_name="scalar_eddy_diffusivity_times_density" unit="kg m-1 s-1" grid_ref="half_level_face_grid" />
<field id="turbulence__rhokm" name="rhokm_bl" long_name="momentum_eddy_diffusivity_times_density" unit="kg m-1 s-1" grid_ref="full_level_face_grid" />
<field id="turbulence__heat_flux_bl" name="heat_flux_bl" long_name="upward_heat_flux_in_air" unit="W m-2" grid_ref="half_level_face_grid" />
<field id="turbulence__moist_flux_bl" name="moist_flux_bl" long_name="upward_moisture_flux_in_air" unit="kg m-2 s-1" grid_ref="half_level_face_grid" />
<field id="turbulence__taux" name="taux" long_name="eastward_wind_stress" unit="N m-2" grid_ref="half_level_face_grid" />
Expand All @@ -685,6 +685,23 @@
<field id="turbulence__ent_t_frac_dsc" name="ent_t_frac_dsc" long_name="fraction_of_time_decoupled_stratocumulus_inversion_is_above_level" unit="1" domain_ref="face" axis_ref="entrainment_levels" />
<field id="turbulence__ent_zrzi" name="ent_zrzi" long_name="level_height_as_fraction_of_sml_inversion_height_above_mixed_layer_base" unit="1" domain_ref="face" axis_ref="entrainment_levels" />
<field id="turbulence__ent_zrzi_dsc" name="ent_zrzi_dsc" long_name="level_height_as_fraction_of_dsc_inversion_height_above_mixed_layer_base" unit="1" domain_ref="face" axis_ref="entrainment_levels" />
<field id="turbulence__tke" name="tke_bl" long_name="turbulent_kinetic_energy" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="turbulence__tsq" name="tsq_bl" long_name="second_moment_of_liquid_water_potential_temperature" unit="K2" grid_ref="full_level_face_grid" />
<field id="turbulence__qsq" name="qsq_bl" long_name="second_moment_of_total_water_mass_ratio" unit="kg2 kg-2" grid_ref="full_level_face_grid" />
<field id="turbulence__cov" name="cov_bl" long_name="covariance_of_liquid_water_potential_temperature_and_total_water_mass_ratio" unit="K kg kg-1" grid_ref="full_level_face_grid" />
<field id="turbulence__zhpar_shcu" name="zhpar_shcu" long_name="height_for_non-gradient_buoyancy_flux_in_TKE_scheme" unit="m" domain_ref="face" />
Comment thread
Adrian-Lock marked this conversation as resolved.
Outdated
<field id="turbulence__master_length" name="master_length" long_name="turbulent_length_scale_in_TKE_scheme" unit="m" grid_ref="full_level_face_grid" />
<field id="turbulence__rhogamu" name="rhogamu_bl" long_name="non-local_eastward_wind_stress_in_TKE_scheme" unit="N m-2" grid_ref="full_level_face_grid" />
<field id="turbulence__rhogamv" name="rhogamv_bl" long_name="non-local_northward_wind_stress_in_TKE_scheme" unit="N m-2" grid_ref="full_level_face_grid" />
<field id="turbulence__rhogamt" name="rhogamt_bl" long_name="non-local_sensible_heat_flux_in_TKE_scheme" unit="W m-2" grid_ref="half_level_face_grid" />
<field id="turbulence__rhogamq" name="rhogamq_bl" long_name="non-local_latent_heat_flux_in_TKE_scheme" unit="W m-2" grid_ref="half_level_face_grid" />
<field id="turbulence__tke_shr_prod" name="tke_shr_prod" long_name="production_rate_of_TKE_by_shear" unit="m2 s-3" grid_ref="full_level_face_grid" />
<field id="turbulence__tke_boy_prod" name="tke_boy_prod" long_name="production_rate_of_TKE_by_buoyancy" unit="m2 s-3" grid_ref="full_level_face_grid" />
<field id="turbulence__tke_dissp" name="tke_dissp" long_name="dissipation_rate_of_TKE" unit="m2 s-3" grid_ref="full_level_face_grid" />
<field id="turbulence__sm25" name="sm25" long_name="stability_function_for_momentum_in_TKE_scheme" unit="1" grid_ref="full_level_face_grid" />
<field id="turbulence__sh25" name="sh25" long_name="stability_function_for_scalar_in_TKE_scheme" unit="1" grid_ref="full_level_face_grid" />
<field id="turbulence__dbdz" name="dbdz" long_name="vertical_gradient_of_buoyancy" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="turbulence__dvdzm" name="dvdzm" long_name="modulus_of_wind_shear" unit="s-1" grid_ref="full_level_face_grid" />
<!-- diffusion diagnostics -->
<field id="diffusion__dtheta_smag" name="dtheta_smag" long_name="theta_increment_from_smagorinsky_mixing" unit="K" grid_ref="full_level_face_grid" />
<field id="diffusion__dmv_smag" name="dmv_smag" long_name="vapour_increment_from_smagorinsky_mixing" unit="kg kg-1" grid_ref="full_level_face_grid" />
Expand Down
5 changes: 5 additions & 0 deletions applications/lfric_atm/metadata/lfric_dictionary.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -269,7 +269,12 @@
<field id="ng_mphys" name="ng_mphys" long_name="number_concentration_of_graupel_after_microphysics" unit="kg-1" grid_ref="full_level_face_grid" />
<!-- turbulence fields -->
<field id="zh" name="zh" long_name="boundary_layer_depth" unit="m" domain_ref="face" />
<field id="zhpar_shcu" name="zhpar_shcu" long_name="height_for_non-gradient_buoyancy_flux" unit="m" domain_ref="face" />
<field id="wvar" name="wvar" long_name="vertical_velocity_variance" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="tke_bl" name="tke_bl" long_name="turbulent_kinetic_energy" unit="m2 s-2" grid_ref="full_level_face_grid" />
<field id="tsq_bl" name="tsq_bl" long_name="second_moment_of_liquid_water_potential_temperature" unit="K2" grid_ref="full_level_face_grid" />
<field id="qsq_bl" name="qsq_bl" long_name="second_moment_of_total_water_mass_ratio" unit="kg2 kg-2" grid_ref="full_level_face_grid" />
<field id="cov_bl" name="cov_bl" long_name="covariance_of_liquid_water_potential_temperature_and_total_water_mass_ratio" unit="K kg kg-1" grid_ref="full_level_face_grid" />
<!-- surface fields -->
<field id="z0msea" name="z0msea" long_name="sea_surface_roughness_length" unit="m" domain_ref="face" />
<field id="surface_conductance" name="surface_conductance" long_name="surface_conductance" unit="m s-1" domain_ref="face" />
Expand Down
2 changes: 1 addition & 1 deletion applications/ngarch/example/configuration_bl.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -508,7 +508,7 @@ ffsl_splitting=6*1,
ffsl_unity_3d=.true.,
ffsl_vertical_order=2,2,1,2,2,2,
field_names='density','potential_temperature','wind','moisture',
'adv_tracer','con_tracer',
'gen_tracer','pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=6*2,
Expand Down
2 changes: 1 addition & 1 deletion applications/ngarch/example/configuration_casim.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -569,7 +569,7 @@ ffsl_splitting=6*1,
ffsl_unity_3d=.true.,
ffsl_vertical_order=2,2,1,2,2,2,
field_names='density','potential_temperature','wind','moisture',
'adv_tracer','con_tracer',
'gen_tracer','pos_tracer',
fv_horizontal_order=2,
fv_vertical_order=2,
horizontal_method=6*2,
Expand Down
Loading
Loading