New clubb external Nov '23.

Author: bstephens82

Skx_module.F90

@@ -75,7 +75,7 @@ subroutine Skx_func( nz, ngrdcol, xp2, xp3, &
     !$acc parallel loop gang vector collapse(2) default(present)
     do k = 1, nz
       do i = 1, ngrdcol
-        Skx(i,k) = xp3(i,k) * sqrt(( xp2(i,k) + Skx_denom_tol )**(-3))
+        Skx(i,k) = xp3(i,k) * sqrt( xp2(i,k) + Skx_denom_tol )**(-3)
       end do
     end do
     !$acc end parallel loop

advance_clubb_core_module.F90

@@ -562,6 +562,7 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
                                         l_use_thvm_in_bv_freq, &
                                         l_use_shear_Richardson, &
                                         l_modify_limiters_for_cnvg_test, & 
+                                        stats_metadata, &
                                         stats_zm, & 
                                         Cx_fnc_Richardson )
 
@@ -596,8 +597,7 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
         ibv_efold
 
     use stats_variables, only: &
-        ishear_sqd, &           ! Variable(s)
-        l_stats_samp
+        stats_metadata_type
 
     use stats_type_utilities, only: &
         stat_update_var      ! Procedure
@@ -606,7 +606,7 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
 
     implicit none
 
-    !!------------------------------ Constant Parameters !------------------------------
+    !------------------------------ Constant Parameters ------------------------------
     real( kind = core_rknd ), parameter :: &
       Richardson_num_divisor_threshold = 1.0e-6_core_rknd, &
       Cx_fnc_Richardson_below_ground_value = one
@@ -621,7 +621,7 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
                                           ! https://github.com/larson-group/clubb/issues/965#issuecomment-1119816722
                                           ! for a plot on how output behaves with varying min_max_smth_mag
 
-    !------------------------------ Input ------------------------------
+    !------------------------------ Input Variables ------------------------------
     integer, intent(in) :: &
       nz, &
       ngrdcol
@@ -658,6 +658,9 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
     logical, intent(in) :: &
       l_modify_limiters_for_cnvg_test
 
+    type (stats_metadata_type), intent(in) :: &
+      stats_metadata
+
     !------------------------------ InOut Variable ------------------------------
     type (stats), target, dimension(ngrdcol), intent(inout) :: &
       stats_zm
@@ -739,10 +742,10 @@ subroutine compute_Cx_fnc_Richardson( nz, ngrdcol, gr, &
     end do
     !$acc end parallel loop
 
-    if ( l_stats_samp ) then
+    if ( stats_metadata%l_stats_samp ) then
       !$acc update host(shear_sqd)
       do i = 1, ngrdcol
-        call stat_update_var( ishear_sqd, shear_sqd(i,:), & ! intent(in)
+        call stat_update_var( stats_metadata%ishear_sqd, shear_sqd(i,:), & ! intent(in)
                               stats_zm(i) )               ! intent(inout)
       end do
     end if