DCMIP test cases

applications/gungho_model/example/configuration.nml

@@ -71,6 +71,7 @@ rotating=.true.,
 shallow=.false.,
 si_momentum_equation=.false.,
 theta_moist_source=.false.
+traditional_coriolis=.true.
 use_multires_coupling=.false.,
 use_physics=.false.,
 use_wavedynamics=.true.,

applications/lfric_atm/metadata/field_def_diags.xml

@@ -20,6 +20,7 @@
   <field id="rho_in_wth" name="rho_in_wth" long_name="dry_air_density_at_cell_interfaces" unit="kg m-3" grid_ref="full_level_face_grid" />
   <field id="wetrho_in_w3" name="wetrho_in_w3" long_name="wet_air_density_at_cell_centres" unit="kg m-3" grid_ref="half_level_face_grid" />
   <field id="wetrho_in_wth" name="wetrho_in_wth" long_name="wet_air_density_at_cell_interfaces" unit="kg m-3" grid_ref="full_level_face_grid" />
+  <field id="dbz" name="dbz" long_name="rainfall_radar_reflectivity" unit="1" grid_ref="full_level_face_grid" />
   <field id="exner_in_wth" name="exner_in_wth" long_name="exner_pressure_at_cell_interfaces" unit="1" grid_ref="full_level_face_grid" />
   <field id="exner_wth_n" name="exner_wth_n" long_name="start_of_timestep_exner_pressure_at_cell_interfaces" unit="1" grid_ref="full_level_face_grid" />
   <field id="u_in_w2h" name="u_in_w2h" long_name="eastward_wind_at_cell_faces" unit="m s-1" grid_ref="half_level_edge_grid" />
@@ -43,6 +44,7 @@
   <field id="w_in_w3_star" name="w_in_w3_star" long_name="physics_predictor_upward_air_velocity_at_cell_centres" unit="ms-1" grid_ref="half_level_face_grid" />
   <field id="height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_face_grid" />
   <field id="height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="full_level_face_grid" />
+  <field id="height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" standard_name="height_above_mean_sea_level" unit="m" grid_ref="half_level_edge_grid" />
   <!-- First line extracts level 0, second line removes vertical coordinate, cannot do as one -->
   <field id="zoomed__orography" field_ref="height_wth" grid_ref="full_level_0_face_zoom" />
   <field id="orography" name="orography" long_name="surface_altitude" standard_name="surface_altitude" unit="m" domain_ref="face" > zoomed__orography </field>

applications/lfric_atm/metadata/field_def_initial_diags.xml

@@ -9,6 +9,7 @@
     <field id="init_m_ci" name="m_ci" long_name="cloud_ice_mixing_ratio" unit="kg kg-1" grid_ref="full_level_face_grid" />
     <field id="init_m_s" name="m_s" long_name="snow_mixing_ratio"  unit="kg kg-1" grid_ref="full_level_face_grid" />
     <field id="init_m_g" name="m_g" long_name="graupel_mixing_ratio"  unit="kg kg-1" grid_ref="full_level_face_grid" />
+    <field id="init_dbz" name="dbz" long_name="rainfall_radar_reflectivity"  unit="1" grid_ref="full_level_face_grid" />
     <field id="init_rho" name="rho" long_name="air_density" standard_name="air_density" unit="kg m-3" grid_ref="half_level_face_grid" />
     <field id="init_xi1" name="xi1" long_name="eastward_vorticity" unit="s-1" grid_ref="half_level_face_grid" />
     <field id="init_xi2" name="xi2" long_name="northward_vorticity" unit="s-1" grid_ref="half_level_face_grid" />
@@ -19,8 +20,10 @@
     <field id="init_potential_vorticity" name="potential_vorticity" long_name="potential_vorticity" standard_name="potential_vorticity" unit="m2 s-1 K kg-1" grid_ref="half_level_face_grid" />
     <field id="init_divergence" name="divergence" long_name="divergence_of_wind" standard_name="divergence_of_wind" unit="s-1" grid_ref="half_level_face_grid" />
     <field id="init_exner" name="exner" long_name="exner_pressure" unit="1" grid_ref="half_level_face_grid" />
+    <field id="init_exner_in_wth" name="exner_in_wth" long_name="exner_pressure_at_cell_interfaces" unit="1" grid_ref="full_level_face_grid" />
     <field id="init_height_w3" name="height_w3" long_name="height_of_cell_centres_above_earth_radius" unit="m" grid_ref="half_level_face_grid" />
     <field id="init_height_wth" name="height_wth" long_name="height_of_cell_interfaces_above_earth_radius" unit="m" grid_ref="full_level_face_grid" />
+    <field id="init_height_w2h" name="height_w2h" long_name="height_of_cell_faces_above_earth_radius" unit="m" grid_ref="half_level_edge_grid" />
     <field id="init_ageofair" name="ageofair" long_name="ageofair" unit="s" grid_ref="half_level_face_grid" />
     <field id="init_area_at_msl" name="cell_area" long_name="grid_cell_area" standard_name="grid_cell_area" unit="m2" domain_ref="face"/>
   </field_group>

applications/lfric_atm/metadata/lfric_dictionary.xml

@@ -12,6 +12,7 @@
   <field id="m_ci" name="m_ci" long_name="m_ci" standard_name="cloud_ice_mixing_ratio" unit="kg kg-1" grid_ref="full_level_face_grid" />
   <field id="m_s" name="m_s" long_name="m_s" unit="kg kg-1" grid_ref="full_level_face_grid" />
   <field id="m_g" name="m_g" long_name="m_g" unit="kg kg-1" grid_ref="full_level_face_grid" />
+  <field id="dbz" name="dbz" long_name="rainfall_radar_reflectivity" unit="1" grid_ref="full_level_face_grid" />
   <field id="ageofair" name="ageofair" long_name="ageofair" unit="s" grid_ref="half_level_face_grid" />
   <field id="theta_ref" name="theta_ref" long_name="air_potential_temperature_ref" standard_name="air_potential_temperature_ref" unit="K" grid_ref="full_level_face_grid" />
   <field id="rho_ref" name="rho_ref" long_name="air_density_ref" standard_name="air_density_ref" unit="kg m-3" grid_ref="half_level_face_grid" />

dependencies.yaml

@@ -30,8 +30,8 @@ lfric_apps:
     ref:
 
 lfric_core:
-    source: git@github.com:MetOffice/lfric_core.git
-    ref: 2025.12.1
+    source: git@github.com:tommbendall/lfric_core.git
+    ref: TBendall/vn3.0_dcmip
 
 moci:
     source: git@github.com:MetOffice/moci.git

rose-stem/app/gungho_model/file/file_def_diags_nwp.xml

@@ -5,12 +5,16 @@
       <field field_ref="theta"/>
       <field field_ref="m_v"/>
       <field field_ref="m_cl"/>
+      <field field_ref="m_r"/>
       <field field_ref="rho"/>
       <field field_ref="u_in_w2h"/>
       <field field_ref="v_in_w2h"/>
       <field field_ref="w_in_wth"/>
       <field field_ref="divergence"/>
       <field field_ref="exner"/>
+      <field field_ref="exner_in_wth"/>
+      <field field_ref="dbz"/>
+      <field field_ref="xi3"/>
     </field_group>
   </file>
 

rose-stem/app/gungho_model/file/file_def_initial.xml

@@ -11,8 +11,10 @@
       <field field_ref="init_w_in_wth"/>
       <field field_ref="init_divergence"/>
       <field field_ref="init_exner"/>
+      <field field_ref="init_exner_in_wth"/>
       <field field_ref="init_height_w3"/>
       <field field_ref="init_height_wth"/>
+      <field field_ref="init_height_w2h"/>
     </field_group>
   </file>
 

rose-stem/app/gungho_model/opt/rose-app-breaking_gw.conf

@@ -0,0 +1,46 @@
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=40000.0
+method='uniform'
+number_of_layers=88
+
+[namelist:formulation]
+shallow=.false.
+use_physics=.true.
+
+[namelist:idealised]
+test='breaking_gw'
+
+[namelist:initial_wind]
+profile='breaking_gw'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=48
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='bell'
+
+[namelist:orography_bell_spherical]
+lambda_centre1=72.0
+lambda_centre2=140.0
+mountain_height=3000.0
+phi_centre1=45.0
+phi_centre2=45.0
+radius_lat=20.0
+radius_lon=3.5
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:section_choice]
+external_forcing=.true.

rose-stem/app/gungho_model/opt/rose-app-mountain_gap.conf

@@ -0,0 +1,55 @@
+[namelist:damping_layer]
+dl_base=10000.0
+dl_str=0.15
+dl_type='standard'
+
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=20000.0
+method='dcmip_mountain'
+number_of_layers=57
+
+[namelist:formulation]
+dlayer_on=.false.
+shallow=.true.
+traditional_coriolis=.true.
+use_physics=.true.
+
+[namelist:idealised]
+test='horizontal_mountain'
+
+[namelist:initial_pressure]
+method='balanced'
+
+[namelist:initial_wind]
+profile='horizontal_mountain'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=120
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='dcmip_gap'
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:planet]
+cp=1004.64
+gravity=9.80616
+omega=7.29211E-5
+p_zero=100000.0
+rd=287.04
+scaling_factor=20.0
+
+[namelist:section_choice]
+external_forcing=.true.

rose-stem/app/gungho_model/opt/rose-app-mountain_gap_fail.conf

@@ -0,0 +1,50 @@
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=20000.0
+method='dcmip_mountain'
+number_of_layers=57
+
+[namelist:formulation]
+dlayer_on=.false.
+shallow=.true.
+traditional_coriolis=.true.
+use_physics=.true.
+
+[namelist:idealised]
+test='horizontal_mountain'
+
+[namelist:initial_pressure]
+method='balanced'
+
+[namelist:initial_wind]
+profile='horizontal_mountain'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=160
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='dcmip_gap'
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:planet]
+cp=1004.64
+gravity=9.80616
+omega=7.29211E-5
+p_zero=100000.0
+rd=287.04
+scaling_factor=20.0
+
+[namelist:section_choice]
+external_forcing=.true.

rose-stem/app/gungho_model/opt/rose-app-mountain_gap_no_rotation.conf

@@ -0,0 +1,50 @@
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=20000.0
+method='dcmip_mountain'
+number_of_layers=57
+
+[namelist:formulation]
+dlayer_on=.false.
+rotating=.false.
+shallow=.true.
+use_physics=.true.
+
+[namelist:idealised]
+test='horizontal_mountain'
+
+[namelist:initial_pressure]
+method='balanced'
+
+[namelist:initial_wind]
+profile='horizontal_mountain'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=160
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='dcmip_gap'
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:planet]
+cp=1004.64
+gravity=9.80616
+omega=0.0
+p_zero=100000.0
+rd=287.04
+scaling_factor=20.0
+
+[namelist:section_choice]
+external_forcing=.true.

rose-stem/app/gungho_model/opt/rose-app-mountain_vortex.conf

@@ -0,0 +1,50 @@
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=20000.0
+method='dcmip_mountain'
+number_of_layers=57
+
+[namelist:formulation]
+dlayer_on=.false.
+shallow=.true.
+traditional_coriolis=.true.
+use_physics=.true.
+
+[namelist:idealised]
+test='horizontal_mountain'
+
+[namelist:initial_pressure]
+method='balanced'
+
+[namelist:initial_wind]
+profile='horizontal_mountain'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=160
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='dcmip_vortex'
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:planet]
+cp=1004.64
+gravity=9.80616
+omega=7.29211E-5
+p_zero=100000.0
+rd=287.04
+scaling_factor=20.0
+
+[namelist:section_choice]
+external_forcing=.true.

rose-stem/app/gungho_model/opt/rose-app-mountain_vortex_no_rotation.conf

@@ -0,0 +1,50 @@
+[namelist:external_forcing]
+hs_random=.false.
+wind_forcing='rayleigh_friction'
+
+[namelist:extrusion]
+domain_height=20000.0
+method='dcmip_mountain'
+number_of_layers=57
+
+[namelist:formulation]
+dlayer_on=.false.
+rotating=.false.
+shallow=.true.
+use_physics=.true.
+
+[namelist:idealised]
+test='horizontal_mountain'
+
+[namelist:initial_pressure]
+method='balanced'
+
+[namelist:initial_wind]
+profile='horizontal_mountain'
+
+[namelist:initialization]
+zero_w2v_wind=.true.
+
+[namelist:io]
+diagnostic_frequency=160
+
+[namelist:orography]
+orog_init_option='analytic'
+profile='dcmip_vortex'
+
+[namelist:physics]
+limit_drag_incs=.false.
+sample_physics_scalars=.true.
+sample_physics_winds=.true.
+sample_physics_winds_correction=.false.
+
+[namelist:planet]
+cp=1004.64
+gravity=9.80616
+omega=0.0
+p_zero=100000.0
+rd=287.04
+scaling_factor=20.0
+
+[namelist:section_choice]
+external_forcing=.true.