Update surface albedo in TUV-x prior to calculating rate constants

schemes/musica/micm/musica_ccpp_micm.F90

@@ -16,7 +16,7 @@ module musica_ccpp_micm
 
 contains
 
-  !> Register MICM constituent properties with the CCPP
+  !> Registers MICM constituent properties with the CCPP
   subroutine micm_register(solver_type, num_grid_cells, constituent_props, errmsg, errcode)
     use ccpp_constituent_prop_mod, only: ccpp_constituent_properties_t
     use musica_micm,               only: Rosenbrock, RosenbrockStandardOrder
@@ -76,7 +76,7 @@ subroutine micm_register(solver_type, num_grid_cells, constituent_props, errmsg,
 
   end subroutine micm_register
 
-  !> Intitialize MICM
+  !> Initializes MICM
   subroutine micm_init(errmsg, errcode)
     character(len=512), intent(out) :: errmsg
     integer,            intent(out) :: errcode
@@ -86,7 +86,7 @@ subroutine micm_init(errmsg, errcode)
 
   end subroutine micm_init
 
-  !> Solve chemistry at the current time step
+  !> Solves chemistry at the current time step
   subroutine micm_run(time_step, temperature, pressure, dry_air_density, &
                       user_defined_rate_parameters, constituents, errmsg, errcode)
     use musica_micm,   only: solver_stats_t
@@ -124,7 +124,7 @@ subroutine micm_run(time_step, temperature, pressure, dry_air_density, &
 
   end subroutine micm_run
 
-  !> Finalize MICM
+  !> Finalizes MICM
   subroutine micm_final(errmsg, errcode)
     character(len=512), intent(out) :: errmsg
     integer,            intent(out) :: errcode

schemes/musica/musica_ccpp.F90

@@ -26,14 +26,17 @@ end subroutine musica_ccpp_register
   !> \section arg_table_musica_ccpp_init Argument Table
   !! \htmlinclude musica_ccpp_init.html
   subroutine musica_ccpp_init(vertical_layer_dimension, vertical_interface_dimension, &
-                              errmsg, errcode)
-    integer,            intent(in)  :: vertical_layer_dimension     ! (count)
-    integer,            intent(in)  :: vertical_interface_dimension ! (count)
+                              photolysis_wavelength_grid_interfaces, errmsg, errcode)
+    use ccpp_kinds, only : kind_phys
+
+    integer,            intent(in)  :: vertical_layer_dimension                 ! (count)
+    integer,            intent(in)  :: vertical_interface_dimension             ! (count)
+    real(kind_phys),    intent(in)  :: photolysis_wavelength_grid_interfaces(:) ! m
     character(len=512), intent(out) :: errmsg
     integer,            intent(out) :: errcode
 
     call tuvx_init(vertical_layer_dimension, vertical_interface_dimension, &
-                   errmsg, errcode)
+                   photolysis_wavelength_grid_interfaces, errmsg, errcode)
     if (errcode /= 0) return
     call micm_init(errmsg, errcode)
     if (errcode /= 0) return
@@ -43,13 +46,14 @@ end subroutine musica_ccpp_init
   !> \section arg_table_musica_ccpp_run Argument Table
   !! \htmlinclude musica_ccpp_run.html
   !!
-  !! The standard name for the variable 'surface_tempearture' is
-  !! `blackbody_temperature_at_surface` because this is what we have as
-  !! the standard name for `cam_in%ts`, whcih represents the same quantity.
+  !! The standard name for the variable 'surface_temperature' is
+  !! 'blackbody_temperature_at_surface' because this is what we have as
+  !! the standard name for 'cam_in%ts', whcih represents the same quantity.
   subroutine musica_ccpp_run(time_step, temperature, pressure, dry_air_density, constituent_props, &
                              constituents, geopotential_height_wrt_surface_at_midpoint,            &
                              geopotential_height_wrt_surface_at_interface, surface_temperature,    &
-                             surface_geopotential, standard_gravitational_acceleration, errmsg, errcode)
+                             surface_geopotential, surface_albedo,                                 &
+                             standard_gravitational_acceleration, errmsg, errcode)
     use musica_ccpp_micm_util,     only: reshape_into_micm_arr, reshape_into_ccpp_arr
     use musica_ccpp_micm_util,     only: convert_to_mol_per_cubic_meter, convert_to_mass_mixing_ratio
     use ccpp_constituent_prop_mod, only: ccpp_constituent_prop_ptr_t
@@ -67,6 +71,7 @@ subroutine musica_ccpp_run(time_step, temperature, pressure, dry_air_density, co
     real(kind_phys),    intent(in)    :: geopotential_height_wrt_surface_at_interface(:,:) ! m (column, interface)
     real(kind_phys),    intent(in)    :: surface_temperature(:)                            ! K
     real(kind_phys),    intent(in)    :: surface_geopotential(:)                           ! m2 s-2
+    real(kind_phys),    intent(in)    :: surface_albedo                                    ! unitless
     real(kind_phys),    intent(in)    :: standard_gravitational_acceleration               ! m s-2
     character(len=512), intent(out)   :: errmsg
     integer,            intent(out)   :: errcode
@@ -92,8 +97,8 @@ subroutine musica_ccpp_run(time_step, temperature, pressure, dry_air_density, co
     call tuvx_run(temperature, dry_air_density,                 &
                   geopotential_height_wrt_surface_at_midpoint,  &
                   geopotential_height_wrt_surface_at_interface, &
-                  surface_temperature,                          &
-                  surface_geopotential,                         &
+                  surface_temperature, surface_geopotential,    &
+                  surface_albedo,                               &
                   standard_gravitational_acceleration,          &
                   photolysis_rate_constants,                    &
                   errmsg, errcode)

schemes/musica/musica_ccpp.meta

@@ -19,6 +19,12 @@
   type = integer
   dimensions = ()
   intent = in
+[ photolysis_wavelength_grid_interfaces ]
+  standard_name = photolysis_wavelength_grid_interfaces
+  units = m
+  type = real | kind = kind_phys
+  dimensions = (photolysis_wavelength_grid_interface_dimension)
+  intent = in
 [ errmsg ]
   standard_name = ccpp_error_message
   units = none
@@ -95,6 +101,12 @@
   units = m2 s-2
   dimensions = (horizontal_loop_extent)
   intent = in
+[ surface_albedo ]
+  standard_name = surface_albedo_due_to_UV_and_VIS_direct
+  type = real | kind = kind_phys
+  units = None
+  dimensions = ()
+  intent = in
 [ standard_gravitational_acceleration ]
   standard_name = standard_gravitational_acceleration
   units = m s-2