Add constituent tendency updater

schemes/dry_adiabatic_adjust/dadadj.meta

@@ -119,6 +119,7 @@
   type = real | kind = kind_phys
   dimensions = (horizontal_loop_extent, vertical_layer_dimension)
   intent = out
+  constituent = True
 [ dadpdf ]
   standard_name = binary_indicator_for_dry_adiabatic_adjusted_grid_cell
   units = fraction

schemes/utilities/physics_tendency_updaters.F90

@@ -9,6 +9,7 @@ module physics_tendency_updaters
   public :: apply_tendency_of_northward_wind_run
   public :: apply_heating_rate_run
   public :: apply_tendency_of_air_temperature_run
+  public :: apply_constituent_tendencies_run
 
 CONTAINS
 
@@ -117,4 +118,29 @@ subroutine apply_tendency_of_air_temperature_run(nz, t_tend, temp, dtdT_total,
 
    end subroutine apply_tendency_of_air_temperature_run
 
+   !> \section arg_table_apply_constituent_tendencies_run Argument Table
+   !!! \htmlinclude apply_constituent_tendencies_run.html
+   subroutine apply_constituent_tendencies_run(nz, const_tend, const, dt, errcode, errmsg)
+      ! Dummy arguments
+      integer,            intent(in)    :: nz                 ! Num vertical layers
+      real(kind_phys),    intent(inout) :: const_tend(:,:,:)  ! constituent tendency array
+      real(kind_phys),    intent(inout) :: const(:,:,:)       ! constituent state array
+      real(kind_phys),    intent(in)    :: dt                 ! physics time step
+      integer,            intent(out)   :: errcode
+      character(len=512), intent(out)   :: errmsg
+
+      ! Local variables
+      integer :: klev
+
+      errcode = 0
+      errmsg = ''
+
+      do klev = 1, nz
+         const(:, klev, :) = const(:, klev, :) + (const_tend(:, klev, :) * dt)
+      end do
+
+      const_tend = 0._kind_phys
+
+   end subroutine apply_constituent_tendencies_run
+
 end module physics_tendency_updaters

schemes/utilities/physics_tendency_updaters.meta

@@ -215,4 +215,53 @@
   type = character | kind = len=512
   dimensions = ()
   intent = out
+#####################################################################
+[ccpp-table-properties]
+  name = apply_constituent_tendencies
+  type = scheme
+[ccpp-arg-table]
+  name = apply_constituent_tendencies_run
+  type = scheme
+[ nz ]
+  standard_name = vertical_layer_dimension
+  long_name = Number of vertical layers
+  units = count
+  type = integer
+  dimensions = ()
+  intent = in
+[ const_tend ]
+  standard_name = ccpp_constituent_tendencies
+  long_name = ccpp constituent tendencies
+  units = none
+  type = real | kind = kind_phys
+  dimensions = (horizontal_loop_extent, vertical_layer_dimension, number_of_ccpp_constituents)
+  intent = inout
+[ const ]
+  standard_name = ccpp_constituents
+  long_name = ccpp constituents
+  units = none
+  type = real | kind = kind_phys
+  dimensions = (horizontal_loop_extent, vertical_layer_dimension, number_of_ccpp_constituents)
+  intent = inout
+[ dt ]
+  standard_name = timestep_for_physics
+  long_name = time step
+  units = s
+  dimensions = ()
+  type = real | kind = kind_phys
+  intent = in
+[ errcode ]
+  standard_name = ccpp_error_code
+  long_name = Error flag for error handling in CCPP
+  units = 1
+  type = integer
+  dimensions = ()
+  intent = out
+[ errmsg ]
+  standard_name = ccpp_error_message
+  long_name = Error message for error handling in CCPP
+  units = none
+  type = character | kind = len=512
+  dimensions = ()
+  intent = out
 #########################################################

test/test_suites/suite_dry_adiabatic_adjust.xml

@@ -4,7 +4,7 @@
   <group name="physics_before_coupler">
     <!-- Dry Adiabatic Adjustment -->
     <scheme>dadadj</scheme>
-    <scheme>dadadj_apply_qv_tendency</scheme>
+    <scheme>apply_constituent_tendencies</scheme>
     <scheme>apply_heating_rate</scheme>
     <scheme>qneg</scheme>
     <scheme>geopotential_temp</scheme>