Ricky Wong (mo-rickywong) · Apr 1, 2026 · Apr 3, 2026 · Apr 3, 2026 · Apr 21, 2026 · Apr 23, 2026
@@ -19,6 +19,8 @@ module adjt_compute_vorticity_alg_mod
   use fs_continuity_mod,                   only : W1, W2
   use function_space_collection_mod,       only : function_space_collection
 
+  use config_mod, only: config_type
+
   implicit none
 
   public
@@ -31,15 +33,16 @@ module adjt_compute_vorticity_alg_mod
   !>          Determined by testing the equality of inner products <Mx, Mx> and <AMx, x>,
   !>          where M is the tangent linear and A is the adjoint.
   !> @param[in]  mesh         The model mesh
-  subroutine adjt_compute_vorticity_alg( mesh )
+  subroutine adjt_compute_vorticity_alg(config, mesh)
 
     use compute_vorticity_alg_mod,           only : compute_vorticity_alg
     use adj_compute_vorticity_alg_mod,       only : adj_compute_vorticity_alg
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer, intent(in) :: mesh
+    type(config_type), intent(in) :: config
+    type(mesh_type),   intent(in) :: mesh
 
     ! Arguments for tl and adj calls
     type(field_type)            :: xi
@@ -101,7 +104,7 @@ module adjt_compute_vorticity_alg_mod
                  setval_x( u_input, u ) )
 
     ! Tangent linear
-    call compute_vorticity_alg( xi, u, weak_form_u )
+    call compute_vorticity_alg(config, xi, u, weak_form_u)
 
     ! < Mx, Mx >
     call invoke( x_innerproduct_x( xi_inner_prod, xi ), &
@@ -130,7 +133,7 @@ module adjt_compute_vorticity_alg_mod
     u_u_input_inner_prod = 0.0_r_def
 
     ! Adjoint
-    call adj_compute_vorticity_alg( xi, u, weak_form_u )
+    call adj_compute_vorticity_alg(config, xi, u, weak_form_u)
 
     ! < AMx, x >
     call invoke( x_innerproduct_y( xi_xi_input_inner_prod,         &

@@ -6,6 +6,7 @@
 !> @brief Module containing adjoint test for atl_derive_exner_from_eos_alg
 module atlt_derive_exner_from_eos_alg_mod
 
+  use config_mod,                          only : config_type
   use field_mod,                           only : field_type
   use function_space_mod,                  only : function_space_type
   use constants_mod,                       only : r_def, i_def
@@ -33,16 +34,18 @@ module atlt_derive_exner_from_eos_alg_mod
   !> @details Passes if adjoint is transpose of tangent linear.
   !>          Determined by testing the equality of inner products <Mx, Mx> and <AMx, x>,
   !>          where M is the tangent linear and A is the adjoint.
+  !> @param[in]  config       Application configuration object
   !> @param[in]  mesh         The model mesh
-  subroutine atlt_derive_exner_from_eos_alg( mesh )
+  subroutine atlt_derive_exner_from_eos_alg(config, mesh)
 
     use tl_derive_exner_from_eos_alg_mod,    only : tl_derive_exner_from_eos
     use atl_derive_exner_from_eos_alg_mod,   only : atl_derive_exner_from_eos
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer, intent(in) :: mesh
+    type(config_type), intent(in) :: config
+    type(mesh_type),   intent(in) :: mesh
 
     ! Arguments for tl and atl calls
     ! Form of state and rhs is [u,theta,rho,exner]
@@ -111,7 +114,7 @@ module atlt_derive_exner_from_eos_alg_mod
                  setval_random( ls_moist_dyn ) )
 
     ! Tangent linear
-    call tl_derive_exner_from_eos( state, moist_dyn, ls_state, ls_moist_dyn )
+    call tl_derive_exner_from_eos( config, state, moist_dyn, ls_state, ls_moist_dyn )
 
     ! < Mx, Mx >
     call invoke( x_innerproduct_x( state_inner_prod(igh_t), state(igh_t) ),           &
@@ -151,7 +154,7 @@ module atlt_derive_exner_from_eos_alg_mod
     moist_dyn_moist_dyn_input_inner_prod = 0.0_r_def
 
     ! Adjoint
-    call atl_derive_exner_from_eos( state, moist_dyn, ls_state, ls_moist_dyn )
+    call atl_derive_exner_from_eos( config, state, moist_dyn, ls_state, ls_moist_dyn )
 
     ! < AMx, x >
     call invoke( x_innerproduct_y( state_state_input_inner_prod(igh_t),           &

@@ -6,6 +6,7 @@
 !> @brief Module containing adjoint test for atl_rhs_alg
 module atlt_rhs_alg_mod
 
+  use config_mod,                          only : config_type
   use sci_assign_field_random_range_alg_mod, &
                                            only : assign_field_random_range
   use field_mod,                           only : field_type
@@ -41,18 +42,20 @@ module atlt_rhs_alg_mod
   !> @details Passes if adjoint is transpose of tangent linear.
   !>          Determined by testing the equality of inner products <Mx, Mx> and <AMx, x>,
   !>          where M is the tangent linear and A is the adjoint.
+  !> @param[in]  config       Application configuration object
   !> @param[in]  mesh         The model mesh
   !> @param[in]  model_clock  The model clock
-  subroutine atlt_rhs_alg( mesh, model_clock )
+  subroutine atlt_rhs_alg( config, mesh, model_clock )
 
-    use tl_rhs_alg_mod,                      only : tl_rhs_alg
-    use atl_rhs_alg_mod,                     only : atl_rhs_alg
+    use tl_rhs_alg_mod,  only: tl_rhs_alg
+    use atl_rhs_alg_mod, only: atl_rhs_alg
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer, intent(in) :: mesh
-    type(model_clock_type),   intent(in) :: model_clock
+    type(config_type),      intent(in) :: config
+    type(mesh_type),        intent(in) :: mesh
+    type(model_clock_type), intent(in) :: model_clock
 
     ! Arguments for tl and atl calls
     ! Form of state and rhs is [u,theta,rho,exner]
@@ -206,9 +209,9 @@ module atlt_rhs_alg_mod
 
 
     ! Tangent linear
-    call tl_rhs_alg( rhs, alpha_dt, base_state, state, moist_dyn, &
-                     ls_state, ls_moist_dyn, compute_eos,         &
-                     dlayer_rhs, model_clock )
+    call tl_rhs_alg( config, rhs, alpha_dt, base_state, state, moist_dyn, &
+                     ls_state, ls_moist_dyn, compute_eos, dlayer_rhs,     &
+                     model_clock )
 
     ! < Mx, Mx >
     call invoke( x_innerproduct_x( base_state_inner_prod(igh_u), base_state(igh_u) ), &
@@ -310,9 +313,9 @@ module atlt_rhs_alg_mod
     rhs_rhs_input_inner_prod(:) = 0.0_r_def
 
     ! Adjoint
-    call atl_rhs_alg( rhs, alpha_dt, base_state, state, moist_dyn, &
-                      ls_state, ls_moist_dyn, compute_eos,         &
-                      dlayer_rhs, model_clock )
+    call atl_rhs_alg( config, rhs, alpha_dt, base_state, state, moist_dyn, &
+                      ls_state, ls_moist_dyn, compute_eos, dlayer_rhs,     &
+                      model_clock )
 
     ! < AMx, x >
     call invoke( x_innerproduct_y( base_state_base_state_input_inner_prod(igh_u), &

@@ -7,6 +7,7 @@
 !> @brief Module containing standard adjoint test for adj_convert_hdiv_field_kernel
 module adjt_sci_convert_hdiv_field_alg_mod
 
+  use config_mod,                           only: config_type
   use field_mod,                            only: field_type
   use function_space_mod,                   only: function_space_type
   use mesh_mod,                             only: mesh_type
@@ -30,18 +31,20 @@ module adjt_sci_convert_hdiv_field_alg_mod
 
   !=============================================================================
   !> @brief Adjoint test for adj_convert_hdiv_field_kernel
+  !> @param[in] config    Application configuration object
   !> @param[in] mesh      A mesh object
   !> @param[in] chi       The coordinate field
   !> @param[in] panel_id  The panel id
-  subroutine adjt_sci_convert_hdiv_field_alg(mesh, chi, panel_id)
+  subroutine adjt_sci_convert_hdiv_field_alg(config, mesh, chi, panel_id)
 
     use sci_convert_hdiv_field_kernel_mod,    only: convert_hdiv_field_kernel_type
     use invoke_adj_cvt_hdiv_field_kernel_mod, only: invoke_adj_convert_hdiv_field_kernel
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer,       intent(in) :: mesh
+    type(config_type),              intent(in) :: config
+    type(mesh_type),                intent(in) :: mesh
     type(field_type), dimension(3), intent(in) :: chi
     type(field_type),               intent(in) :: panel_id
 
@@ -74,6 +77,12 @@ module adjt_sci_convert_hdiv_field_alg_mod
     real(kind=r_def), parameter :: overall_tolerance = 1500.0_r_def
     real(kind=r_def), parameter :: eps = 1e-30_r_def
 
+    integer(i_def) :: coord_system
+    real(r_def)    :: scaled_radius
+
+    coord_system  = config%finite_element%coord_system()
+    scaled_radius = config%planet%scaled_radius()
+
     ! Initialise arguments
     vector_space_w3_ptr => function_space_collection%get_fs(mesh, element_order_h, &
                                                             element_order_v, W3)
@@ -96,7 +105,10 @@ module adjt_sci_convert_hdiv_field_alg_mod
     ! Final initialisation and call TL kernel
     call invoke( setval_random(computational_field), &
                  setval_X(computational_field_input, computational_field), &
-                 convert_hdiv_field_kernel_type(physical_field3, computational_field, chi, panel_id) )
+                 convert_hdiv_field_kernel_type(physical_field3, computational_field, &
+                                                chi, panel_id, mesh%geometry(), &
+                                                mesh%topology(), &
+                                                coord_system, scaled_radius) )
 
     ! Get dot product <M dx, M dx>
     write(log_scratch_space, *) "adjt_sci_convert_hdiv_field inner products:"
@@ -118,7 +130,10 @@ module adjt_sci_convert_hdiv_field_alg_mod
     call invoke(inc_a_times_X( computational_field_sf, computational_field ))
 
     ! Call adjoint
-    call invoke_adj_convert_hdiv_field_kernel(physical_field3, computational_field, chi, panel_id)
+    call invoke_adj_convert_hdiv_field_kernel(physical_field3, computational_field, &
+                                              chi, panel_id, mesh%geometry(),       &
+                                              mesh%topology(), coord_system,        &
+                                              scaled_radius)
 
     ! Calculate dot product <AM dx, dx>
     inner2 = 0.0_r_def

@@ -7,7 +7,7 @@
 !> @brief Module containing standard adjoint tests for the w3/wth to/from w2 interpolations
 module adjt_interpolation_alg_mod
 
-  use base_mesh_config_mod,          only : geometry, topology
+  use config_mod,                    only : config_type
   use field_mod,                     only : field_type
   use function_space_mod,            only : function_space_type
   use mesh_mod,                      only : mesh_type
@@ -32,16 +32,18 @@ module adjt_interpolation_alg_mod
 
   !=============================================================================
   !> @brief Adjoint test for w2 to w3/wth interpolations
-  !> @param[in] mesh  A mesh object
-  subroutine adjt_interp_w2_to_w3wth_alg( mesh )
+  !> @param[in] config  Application configuration object
+  !> @param[in] mesh    A mesh object
+  subroutine adjt_interp_w2_to_w3wth_alg( config, mesh )
 
     use interpolation_alg_mod,         only : interp_w2_to_w3wth_alg
     use adj_interpolation_alg_mod,     only : adj_interp_w2_to_w3wth_alg
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer, intent(in) :: mesh
+    type(config_type), intent(in) :: config
+    type(mesh_type),   intent(in) :: mesh
 
     ! Arguments for forward and adjoint calls
     type(field_type)                   :: u_in_w2
@@ -108,7 +110,7 @@ module adjt_interpolation_alg_mod
                  setval_random(v_in_w3), setval_x(v_in_w3_inp, v_in_w3), &
                  setval_random(w_in_wth), setval_x(w_in_wth_inp, w_in_wth) )
 
-    call interp_w2_to_w3wth_alg( u_in_w2, u_in_w3, v_in_w3, w_in_wth )
+    call interp_w2_to_w3wth_alg( config, u_in_w2, u_in_w3, v_in_w3, w_in_wth )
 
     call invoke( x_innerproduct_x(u2_inner_prod, u_in_w2), &
                  x_innerproduct_x(u3_inner_prod, u_in_w3), &
@@ -148,7 +150,7 @@ module adjt_interpolation_alg_mod
     v3_v3_inp_inner_prod = 0.0_r_def
     wt_wt_inp_inner_prod = 0.0_r_def
 
-    call adj_interp_w2_to_w3wth_alg( u_in_w2, u_in_w3, v_in_w3, w_in_wth )
+    call adj_interp_w2_to_w3wth_alg( config, u_in_w2, u_in_w3, v_in_w3, w_in_wth )
 
     call invoke( x_innerproduct_y( u2_u2_inp_inner_prod, u_in_w2, u_in_w2_inp ), &
                  x_innerproduct_y( u3_u3_inp_inner_prod, u_in_w3, u_in_w3_inp ), &
@@ -178,16 +180,18 @@ module adjt_interpolation_alg_mod
 
   !=============================================================================
   !> @brief Adjoint test for w3/wth to w2 interpolations
-  !> @param[in] mesh  A mesh object
-  subroutine adjt_interp_w3wth_to_w2_alg( mesh )
+  !> @param[in] config  Application configuration object
+  !> @param[in] mesh    A mesh object
+  subroutine adjt_interp_w3wth_to_w2_alg( config, mesh )
 
     use interpolation_alg_mod,         only : interp_w3wth_to_w2_alg
     use adj_interpolation_alg_mod,     only : adj_interp_w3wth_to_w2_alg
 
     implicit none
 
     ! Arguments
-    type(mesh_type), pointer, intent(in) :: mesh
+    type(config_type), intent(in) :: config
+    type(mesh_type),   intent(in) :: mesh
 
     ! Arguments for forward and adjoint calls
     type(field_type)                   :: u_in_w2
@@ -254,8 +258,7 @@ module adjt_interpolation_alg_mod
                  setval_random(v_in_w3), setval_x(v_in_w3_inp, v_in_w3), &
                  setval_random(w_in_wth), setval_x(w_in_wth_inp, w_in_wth) )
 
-    call interp_w3wth_to_w2_alg( u_in_w2, u_in_w3, v_in_w3, w_in_wth, &
-                                 geometry, topology )
+    call interp_w3wth_to_w2_alg( config, u_in_w2, u_in_w3, v_in_w3, w_in_wth )
 
     call invoke( x_innerproduct_x(u2_inner_prod, u_in_w2), &
                  x_innerproduct_x(u3_inner_prod, u_in_w3), &
@@ -295,8 +298,7 @@ module adjt_interpolation_alg_mod
     v3_v3_inp_inner_prod = 0.0_r_def
     wt_wt_inp_inner_prod = 0.0_r_def
 
-    call adj_interp_w3wth_to_w2_alg( u_in_w2, u_in_w3, v_in_w3, w_in_wth, &
-                                     geometry, topology )
+    call adj_interp_w3wth_to_w2_alg( config, u_in_w2, u_in_w3, v_in_w3, w_in_wth )
 
     call invoke( x_innerproduct_y( u2_u2_inp_inner_prod, u_in_w2, u_in_w2_inp ), &
                  x_innerproduct_y( u3_u3_inp_inner_prod, u_in_w3, u_in_w3_inp ), &

@@ -6,7 +6,9 @@
 !> @brief Module containing adjoint test for adj_w3h_adv_upd_lookup_lookup_kernel
 module adjt_w3h_adv_upd_lookup_alg_mod
 
+
   use adj_trans_lookup_cache_mod,      only : adj_trans_lookup_cache_type
+  use config_mod,                      only : config_type
   use integer_field_mod,               only : integer_field_type
   use r_tran_field_mod,                only : r_tran_field_type
   use fs_continuity_mod,               only : W3, W2
@@ -36,17 +38,19 @@ module adjt_w3h_adv_upd_lookup_alg_mod
   !> @details Passes if adjoint is transpose of tangent linear.
   !>          Determined by testing the equality of inner products <Mx, Mx> and <AMx, x>,
   !>          where M is the tangent linear and A is the adjoint.
+  !> @param[in]  config                  Application configuration object
   !> @param[in]  mesh                    The model mesh
   !> @param[in]  adj_lookup_table_cache  The lookup table cache
-  subroutine adjt_w3h_adv_upd_lookup_alg( mesh, adj_lookup_table_cache )
+  subroutine adjt_w3h_adv_upd_lookup_alg( config, mesh, adj_lookup_table_cache )
 
     use w3h_advective_update_kernel_mod,     only : w3h_advective_update_kernel_type
     use invoke_adj_w3h_adv_upd_lookup_mod,   only : invoke_adj_w3h_adv_upd_lookup_kernel_type
 
     implicit none
 
     ! Arguments
-    type(mesh_type),          pointer, intent(in) :: mesh
+    type(config_type),                 intent(in) :: config
+    type(mesh_type),                   intent(in) :: mesh
     type(adj_trans_lookup_cache_type), intent(in) :: adj_lookup_table_cache
 
     ! Internal arguments for kernels
@@ -84,8 +88,9 @@ module adjt_w3h_adv_upd_lookup_alg_mod
     real(kind=r_tran)            :: machine_tol
     real(kind=r_tran)            :: relative_diff
 
-    nullify(m3_inv, lookup_w3h_adv_upd, lookup_w3h_adv_upd_field, set_count_w3h_adv_upd_field, &
-            vector_space_w3_md_ptr, vector_space_w3_ptr, vector_space_w2_ptr)
+    nullify(m3_inv, lookup_w3h_adv_upd, lookup_w3h_adv_upd_field, &
+            set_count_w3h_adv_upd_field, vector_space_w3_md_ptr,  &
+            vector_space_w3_ptr, vector_space_w2_ptr)
 
     vector_space_w3_md_ptr => function_space_collection%get_fs(                &
         mesh, element_order_h, element_order_v, W3, ndata, ndata_first=.false. &
@@ -101,7 +106,7 @@ module adjt_w3h_adv_upd_lookup_alg_mod
     call wind%initialise( vector_space=vector_space_w2_ptr, name='wind' )
     call advective_increment_input%initialise( vector_space=vector_space_w3_ptr, name='advective_increment_input' )
     call tracer_input%initialise( vector_space=vector_space_w3_md_ptr, name='tracer_input' )
-    m3_inv => get_inverse_mass_matrix_fe( W3, mesh%get_id() )
+    m3_inv => get_inverse_mass_matrix_fe( config, mesh, W3 )
     advective_increment_inner_prod = 0.0_r_def
     tracer_inner_prod = 0.0_r_def
 

@@ -148,7 +148,7 @@ contains
 
     call vector_mx%set_scalar(0.0_r_def)
     call vector_amx%set_scalar(0.0_r_def)
-    call mixed_op%apply( vector_x, vector_mx )
+    call mixed_op%apply( config, vector_x, vector_mx )
 
     ! --------------------------------------------------------------------------
     ! Compute <Mx,Mx>
@@ -186,7 +186,7 @@ contains
 
     adj_mixed_op = adj_mixed_operator_type()
 
-    call adj_mixed_op%apply( vector_mx, vector_amx )
+    call adj_mixed_op%apply( config, vector_mx, vector_amx )
 
     ! --------------------------------------------------------------------------
     ! Compute <AMx,x>