Skip to content

Add option to horizontally homogenize the Stokes drift when used via … #967

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
merged 2 commits into from
Sep 9, 2025
Merged

Add option to horizontally homogenize the Stokes drift when used via … #967

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
07be9de
Add option to horizontally homogenize the Stokes drift when used via …
Sep 2, 2025
53e2942
Add variable description in new method for horizontally averaging Sto…
Sep 9, 2025
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
16 changes: 16 additions & 0 deletions src/user/MOM_wave_interface.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ module MOM_wave_interface
use MOM_io, only : file_exists, get_var_sizes, read_variable
use MOM_io, only : vardesc, var_desc
use MOM_safe_alloc, only : safe_alloc_ptr
use MOM_spatial_means, only : global_area_mean
use MOM_time_manager, only : time_type, operator(+), operator(/)
use MOM_unit_scaling, only : unit_scale_type
use MOM_variables, only : thermo_var_ptrs, surface
Expand Down Expand Up @@ -66,6 +67,7 @@ module MOM_wave_interface
logical, public :: Stokes_DDT = .false. !< Developmental:
!! True if Stokes d/dt is used
logical, public :: Passive_Stokes_DDT = .false. !< Keeps Stokes_DDT on, but doesn't affect dynamics
logical :: Homogenize_Surfbands !< True to homogenize surface band Stokes drift in the horizontal

real, allocatable, dimension(:,:,:), public :: &
Us_x !< 3d zonal Stokes drift profile [L T-1 ~> m s-1]
Expand Down Expand Up @@ -441,6 +443,11 @@ subroutine MOM_wave_interface_init(time, G, GV, US, param_file, CS, diag)
"A layer thickness below which the cell-center Stokes drift is used instead of "//&
"the cell average. This is only used if WAVE_INTERFACE_ANSWER_DATE < 20230101.", &
units="m", default=0.1, scale=US%m_to_Z, do_not_log=(CS%answer_date>=20230101))
call get_param(param_file, mdl, "HOMOGENIZE_SURFBANDS", CS%Homogenize_Surfbands, &
"A logical which causes the code to horizontally homogenize the surface band "//&
"Stokes drift, which is needed in column mode to avoid round-off differences. "//&
"At present it only works with DATAOVERRIDE, and is not coded for COUPLER.",&
default=.false.)
call get_param(param_file, mdl, "SURFBAND_SOURCE", TMPSTRING2, &
"Choice of SURFACE_BANDS data mode, valid options include: \n"//&
" DATAOVERRIDE - Read from NetCDF using FMS DataOverride. \n"//&
Expand Down Expand Up @@ -1071,6 +1078,7 @@ subroutine Surface_Bands_by_data_override(Time, G, GV, US, CS)
character(len=48) :: dim_name(4) ! The names of the dimensions of the variable.
character(len=20) :: varname ! The name of an input variable for data override.
real :: PI ! 3.1415926535... [nondim]
real :: avgx, avgy ! The global averages of temp_x and temp_y [L T-1 ~> m s-1]
logical :: wavenumber_exists
integer :: ndims, b, i, j

Expand Down Expand Up @@ -1151,6 +1159,14 @@ subroutine Surface_Bands_by_data_override(Time, G, GV, US, CS)
endif
enddo
enddo
if (CS%Homogenize_Surfbands) then
avgx = global_area_mean(temp_x, G)
avgy = global_area_mean(temp_y, G)
do j = G%jsd,G%jed ; do i = G%Isd,G%Ied ; if (G%mask2dT(i,j) > 0.0) then
temp_y(i,j) = avgy
temp_x(i,j) = avgx
endif ; enddo ; enddo
endif

! Interpolate to u/v grids
do j = G%jsc,G%jec
Expand Down