Skip to content

Cricital latitude modification of Henyey #599

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
marshallward merged 3 commits into from
Apr 22, 2024
Merged

Cricital latitude modification of Henyey #599

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
851e17b
Proto-type cricital latitude modification of Henyey
adcroft Jul 14, 2023
ad33544
Merge branch 'dev/gfdl' into no-arctic-henyey
marshallward Apr 19, 2024
dbab09f
Merge branch 'dev/gfdl' into no-arctic-henyey
marshallward Apr 22, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 7 additions & 0 deletions src/parameterizations/vertical/MOM_bkgnd_mixing.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -59,6 +59,8 @@ module MOM_bkgnd_mixing
real :: N0_2Omega !< ratio of the typical Buoyancy frequency to
!! twice the Earth's rotation period, used with the
!! Henyey scaling from the mixing [nondim]
real :: Henyey_max_lat !< A latitude poleward of which the Henyey profile
!! is returned to the minimum diffusivity [degN]
real :: prandtl_bkgnd !< Turbulent Prandtl number used to convert
!! vertical background diffusivity into viscosity [nondim]
real :: Kd_tanh_lat_scale !< A nondimensional scaling for the range of
Expand Down Expand Up @@ -282,6 +284,10 @@ subroutine bkgnd_mixing_init(Time, G, GV, US, param_file, diag, CS, physical_OBL
call get_param(param_file, mdl, "OMEGA", CS%omega, &
"The rotation rate of the earth.", &
units="s-1", default=7.2921e-5, scale=US%T_to_s)
call get_param(param_file, mdl, "HENYEY_MAX_LAT", CS%Henyey_max_lat, &
"A latitude poleward of which the Henyey profile "//&
"is returned to the minimum diffusivity", &
units="degN", default=95.0)
endif

call get_param(param_file, mdl, "KD_TANH_LAT_FN", CS%Kd_tanh_lat_fn, &
Expand Down Expand Up @@ -447,6 +453,7 @@ subroutine calculate_bkgnd_mixing(h, tv, N2_lay, Kd_lay, Kd_int, Kv_bkgnd, j, G,
I_x30 = 2.0 / invcosh(CS%N0_2Omega*2.0) ! This is evaluated at 30 deg.
do i=is,ie
abs_sinlat = abs(sin(G%geoLatT(i,j)*deg_to_rad))
if (abs(G%geoLatT(i,j))>CS%Henyey_max_lat) abs_sinlat = min_sinlat
Kd_sfc(i) = max(CS%Kd_min, CS%Kd * &
((abs_sinlat * invcosh(CS%N0_2Omega / max(min_sinlat, abs_sinlat))) * I_x30) )
enddo
Expand Down