[MITgcm-support] SSH blow-up at Orlanski BC with small deltaT
Yilang Xu
yxu at whoi.edu
Thu Jun 6 18:52:33 EDT 2019
Hi everyone,
The subject title sounds crazy but I did encounter this problem. I set up a high-resolution experiment (1~2 km for grid size in x-y, and 10m in z) and initialized it with linear stratification. The boundaries are Orlanski with net flow balanced.
Initially I used deltaT=60s but the code blew up with NaNs due to CFL condition. Then I decrease deltaT to 30s, the run could come to a normal end, but SSH increased/decreased to unreasonable values at Orlanski boundaries. I then used deltaT=5s, and SSH seemed to have smaller anomalies, but still blew up at boundaries.
It seems that decreasing deltaT is the right way to ensure proper SSH behavior, and this did work for some cases with other kinds of initial stratifications. So I further decreased deltaT to 2s, but then I got much higher SSH blowups at boundaries than the case with deltaT=5s.
According to this old post http://mailman.mitgcm.org/pipermail/mitgcm-support/2013-January/008122.html , this might be related to viscAhGrid, which is already used in my model and is calculated as below:
viscAh = 0.25*L**2*viscAhGrid/deltaT
The formula above shows that decreasing deltaT increases viscAh. What is confusing to me is that decreasing deltaT can give me both better or worse results. I assume this is related to viscAhGrid, but there might be other SSH instability issues that matter (Orlanski, etc).
I have included the data file below. Really appreciate any help on this problem.
Thanks,
Yilang
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
Tref = 50*-1.9,
Sref = 34.0000, 34.0612, 34.1224, 34.1837, 34.2449,
34.3061, 34.3673, 34.4286, 34.4898, 34.5510,
34.6122, 34.6735, 34.7347, 34.7959, 34.8571,
34.9184, 34.9796, 35.0408, 35.1020, 35.1633,
35.2245, 35.2857, 35.3469, 35.4082, 35.4694,
35.5306, 35.5918, 35.6531, 35.7143, 35.7755,
35.8367, 35.8980, 35.9592, 36.0204, 36.0816,
36.1429, 36.2041, 36.2653, 36.3265, 36.3878,
36.4490, 36.5102, 36.5714, 36.6327, 36.6939,
36.7551, 36.8163, 36.8776, 36.9388, 37.0000,
viscAz=1.E-3,
# viscAh=10.0,
viscAhGrid=1.E-4,
no_slip_sides=.FALSE.,
no_slip_bottom=.FALSE.,
diffKhT=0.0,
diffKzT=1.E-6,
diffKhS=0.0,
diffKzS=1.E-6,
bottomDragQuadratic=2.5E-3,
eosType='JMD95Z',
HeatCapacity_cp = 3974.0,
rhoConst=1030.,
rhoNil=1030.,
gravity=9.81,
convertFW2Salt = 33.4,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
exactConserv=.TRUE.,
hFacMin=0.05,
nonHydrostatic=.FALSE.,
readBinaryPrec=64,
implicitDiffusion = .TRUE.,
implicitViscosity = .TRUE.,
saltAdvScheme=33,
saltVertAdvScheme =33,
tempAdvScheme=33,
tempVertAdvScheme=33,
staggerTimeStep=.TRUE.,
# useCDScheme = .TRUE.,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-13,
cg3dMaxIters=400,
cg3dTargetResidual=1.E-13,
&
# Time stepping parameters
&PARM03
nIter0=0,
nTimeSteps=1728000,
deltaT=2.0,
abEps=0.1,
cAdjFreq = 0.,
# tauCD = 400000.,
pChkptFreq=1728000.0,
chkptFreq=0.0,
dumpFreq=432000.0,
taveFreq=0.0,
monitorFreq=36000.,
monitorSelect=2,
&
# Gridding parameters
&PARM04
usingSphericalPolarGrid=.TRUE.,
ygOrigin = -80.0,
delXfile='dxfile.bin',
delYfile='dyfile.bin',
delZ=50*10.0,
&
# Input datasets
&PARM05
bathyFile='bathy.box',
# hydrogThetaFile='temp_3d_end_state.bin',
# hydrogSaltFile='salt_3d_end_state.bin',
# uVelInitFile='u_3d_end_state.bin',
# vVelInitFile='v_3d_end_state.bin',
&
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190606/76f9073c/attachment.html>
More information about the MITgcm-support
mailing list