From e04ddb951bae7dfc17759bece02ad492a7474316 Mon Sep 17 00:00:00 2001 From: Kate Hedstrom Date: Wed, 10 Jan 2024 14:37:24 -0900 Subject: [PATCH] little_h can go to zero in MLE diagnostic --- config_src/external/GFDL_ocean_BGC/generic_tracer.F90 | 2 +- src/parameterizations/lateral/MOM_mixed_layer_restrat.F90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/config_src/external/GFDL_ocean_BGC/generic_tracer.F90 b/config_src/external/GFDL_ocean_BGC/generic_tracer.F90 index 37f3499901..5e26384a40 100644 --- a/config_src/external/GFDL_ocean_BGC/generic_tracer.F90 +++ b/config_src/external/GFDL_ocean_BGC/generic_tracer.F90 @@ -78,7 +78,7 @@ subroutine generic_tracer_source(Temp,Salt,tv,rho_dzt,dzt,hblt_depth,ilb,jlb,tau integer, intent(in) :: jlb !< Lower bounds of y extent of input arrays on data domain real, dimension(ilb:,jlb:,:), intent(in) :: Temp !< Potential temperature [deg C] real, dimension(ilb:,jlb:,:), intent(in) :: Salt !< Salinity [psu] - type(thermo_var_ptrs), intent(in) :: tv + type(thermo_var_ptrs), intent(in) :: tv !< structure containing pointers to available thermodynamic fields real, dimension(ilb:,jlb:,:), intent(in) :: rho_dzt !< Mass per unit area of each layer [kg m-2] real, dimension(ilb:,jlb:,:), intent(in) :: dzt !< Ocean layer thickness [m] real, dimension(ilb:,jlb:), intent(in) :: hblt_depth !< Boundary layer depth [m] diff --git a/src/parameterizations/lateral/MOM_mixed_layer_restrat.F90 b/src/parameterizations/lateral/MOM_mixed_layer_restrat.F90 index 645ff0fa7a..8ccdf4488d 100644 --- a/src/parameterizations/lateral/MOM_mixed_layer_restrat.F90 +++ b/src/parameterizations/lateral/MOM_mixed_layer_restrat.F90 @@ -944,7 +944,7 @@ subroutine mixedlayer_restrat_Bodner(CS, G, GV, US, h, uhtr, vhtr, tv, forces, d (G%CoriolisBu(I-1,J) + G%CoriolisBu(I,J-1))))**2, absurdly_small_freq2) lf_bodner_diag(i,j) = 0.25 * (( CS%mstar * u_star3 + CS%nstar * w_star3 )**two_thirds & ! (this line m2 s-2) * ( US%m_to_L * GV%m_to_H * US%T_to_s**2 )) & ! [L H s2 m-2 T-2 ~> 1 or kg m-3] - / (f2_h * little_h(i,j))! [T-2 H ~> m s-2] + / (f2_h * (max(little_h(i,j), GV%Angstrom_H)))! [T-2 H ~> m s-2] enddo ; enddo