MOM6: Fix interior obcs (#814)

- NOAA-GFDL/MOM6@5e4f97bf4 Fix interior obcs (#814)
NOAA-GFDL · Feb 26, 2025 · 19369da · 19369da
1 parent 74bc126
commit 19369da
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 1 deletion.
diff --git a/ocean_only/DOME/MOM_parameter_doc.all b/ocean_only/DOME/MOM_parameter_doc.all
@@ -298,6 +298,9 @@ OBC_RAMP_TIMESCALE = 1.0        !   [days] default = 1.0
                                 ! If RAMP_OBCS is true, this sets the ramping timescale.
 OBC_TIDE_N_CONSTITUENTS = 0     ! default = 0
                                 ! Number of tidal constituents being added to the open boundary.
+EXTERIOR_OBC_BUG = True         !   [Boolean] default = True
+                                ! If true, recover a bug in barotropic solver and other routines when boundary
+                                ! contitions interior to the domain are used.
 OBC_SEGMENT_001 = "J=N,I=110:100,SIMPLE" !
                                 ! Documentation needs to be dynamic?????
 OBC_REMAPPING_SCHEME = "PLM"    ! default = "PLM"

diff --git a/ocean_only/circle_obcs/MOM_parameter_doc.all b/ocean_only/circle_obcs/MOM_parameter_doc.all
@@ -283,6 +283,9 @@ OBC_RAMP_TIMESCALE = 1.0        !   [days] default = 1.0
                                 ! If RAMP_OBCS is true, this sets the ramping timescale.
 OBC_TIDE_N_CONSTITUENTS = 0     ! default = 0
                                 ! Number of tidal constituents being added to the open boundary.
+EXTERIOR_OBC_BUG = True         !   [Boolean] default = True
+                                ! If true, recover a bug in barotropic solver and other routines when boundary
+                                ! contitions interior to the domain are used.
 OBC_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI" !
                                 ! Documentation needs to be dynamic?????
 OBC_SEGMENT_002 = "J=0,I=0:N,FLATHER,ORLANSKI" !

diff --git a/src/MOM6 b/src/MOM6
+58 −8		src/core/MOM_barotropic.F90
+0 −3		src/core/MOM_boundary_update.F90
+2 −1		src/core/MOM_dynamics_split_RK2.F90
+1 −1		src/core/MOM_dynamics_split_RK2b.F90
+19 −6		src/core/MOM_open_boundary.F90
+35 −0		src/tracer/MOM_tracer_advect.F90