[MITgcm-support] Stevens BC

[赤羽蓝]

Hi all,

I am running a regional ocean model with the OBCS package. When I use the
Stevens boundary conditions, the model always crashes with "stops due to
EXTREME Pot.Temp". The model is OK if I turn off the Stevens BC. I have
tried different inflows(outflow) at the Stevens boundary, which encounter
the same crashing. The sponge boundary layer didn't help either. I also
tried different relaxation time, no help.

Attaching my data.obcs file.

Any ideas will be grateful.

Lambert

# ***************
# Open boundaries
# ***************
 &OBCS_PARM01
# This flag turns off checking and fixing problematic topography across
# open boundaries.
 OB_Jnorth = 360*-1,
 OB_Jsouth = 360*1,
 OB_Iwest = 528*1,
 OB_Ieast = 528*-1,
 useOBCSprescribe=.TRUE.,
 useStevensWest = .TRUE.,
 useStevensNorth = .TRUE.,
# useOBCSsponge=.TRUE.,
 useOBCSbalance=.FALSE.,
# OBCS_balanceFacE = 1.,
# OBCS_balanceFacW = 1.,
# OBCS_balanceFacN = 1.,
# OBCS_balanceFacS = 1.,

 OBNsFile='sn.scs',
 OBNtFile='tn.scs',
 OBNvFile='vn.scs',
 OBNuFile='un.scs',

 OBSsFile='ss.scs',
 OBStFile='ts.scs',
 OBSvFile='vs.scs',
 OBSuFile='us.scs',

 OBEsFile='se.scs',
 OBEtFile='te.scs',
 OBEvFile='ve.scs',
 OBEuFile='ue.scs',

 OBWsFile='sw.scs',
 OBWtFile='tw.scs',
 OBWvFile='vw.scs',
 OBWuFile='uw.scs',
 &

# &OBCS_PARM03
# Urelaxobcsinner=864000,
# Urelaxobcsbound=864000,
# Vrelaxobcsinner=864000,
# Vrelaxobcsbound=864000,
# spongeThickness=5,
# &

 &OBCS_PARM04
 TrelaxStevens = 86400.,
 SrelaxStevens = 86400.,
 &
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210513/02b9d2b6/attachment.html>