[MITgcm-support] problem of open boundary
Tongya Liu
liutongya13 at mails.ucas.ac.cn
Mon Mar 30 22:43:07 EDT 2015
Hi all,
I met a very strange prolem in my model recently. My model terrain is 109-135E&0-30N(the grid is 360*300). The northern,eastern and southern boundary are open boundaries. HYCOM flow feild is used to supply the open boundary velocity. But there is a strong northward flow just close to the eastern boundary, whose velocity and direction is quite different from the eastern boundary. It seems that the boundary input didn't work in this place. This strange flow has influenced the interior circulation. I adjusted the boundary velocity and the boundary conditions, but it still exist.
I show the result in this web: http://photo.blog.sina.com.cn/photo/1955901475/0028mLjZgy6R7I55jUIf0#pic
Please give me some advices. Thanks very much.
Tongya Liu
data:
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
viscAr=1.77E-3,
viscAh=100,
diffKhT=0,
diffKrT=3E-5,
diffKhS=0,
diffKrS=3E-5,
useFullLeith=.TRUE.,
viscC4Leith=1.5,
viscC4Leithd=1.5,
viscA4GridMax=1.0,
useAreaViscLength=.TRUE.,
useAnisotropicViscAGridMax=.TRUE.,
tempAdvScheme=33,
saltAdvScheme=33,
eosType = 'JMD95Z',
multiDimAdvection=.TRUE.,
vectorInvariantMomentum=.TRUE.,
no_slip_sides=.FALSE.,
no_slip_bottom=.TRUE.,
bottomDragQuadratic=0.0025,
implicitDiffusion=.TRUE.,
implicitFreeSurface=.TRUE.,
implicitViscosity=.TRUE.,
staggerTimeStep=.TRUE.,
implicitFreeSurface=.TRUE.,
exactConserv=.TRUE.,
hFacMin=0.3,
hFacMindr=50.,
hFacInf=0.1,
hFacSup=5.,
readBinaryPrec=64,
#- not safe to use globalFiles in multi-processors runs
globalFiles=.TRUE.,
useSingleCpuIO=.TRUE.,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=500,
cg2dTargetResidual=1.E-5,
&
# Time stepping parameters
&PARM03
nIter0= 0,
nTimeSteps = 25920,
#nTimeSteps = 200,
#startTime=0.,
#endTime=311040000.,
basetime=0,
deltaT=1200.0,
abEps=0.1,
momDissip_In_AB=.FALSE.,
dumpFreq=864000,
pChkptFreq=31104000.,
chkptFreq=2592000.,
monitorSelect=2,
monitorFreq=864000.,
&
# Gridding parameters
&PARM04
usingSphericalPolarGrid=.TRUE.,
ygOrigin=0.0,
xgOrigin=99.,
dxSpacing=0.1,
dySpacing=0.1,
delR=30.,30.,30.,30.,30.,30.,
30.,30.,30.,30.,30.,30.,
30.,30.,30.,30.,30.,30.,
30.,30.,50.,50.,100.,200.,
300.,500.,600.,800.,800.,1000.,
&
# Input datasets
&PARM05
bathyFile='topog.bin',
hydrogThetaFile='Temp.bin',
hydrogSaltFile='Salt.bin',
&
data.obcs
# Open-boundaries
&OBCS_PARM01
# open boundaries.
OBCSfixTopo=.FALSE.,
OB_Jnorth= 360*-1,
OB_Jsouth= 360*1,
OB_Ieast = 300*-1,
useOBCSprescribe = .TRUE.,
#useOBCSsponge=.TRUE.,
useOBCSbalance=.TRUE.,
OBCS_balanceFacN= 1.,
OBCS_balanceFacS= 1,
OBCS_balanceFacE= 1.,
OBNuFile = 'HY_NorthU.bin',
OBSuFile = 'HY_SouthU.bin',
OBEuFile = 'HY_EastU.bin',
OBNvFile = 'HY_NorthV.bin',
OBSvFile = 'HY_SouthV.bin',
OBEvFile = 'HY_EastV.bin',
OBNsFile = 'snorth.bin',
OBNtFile = 'tnorth.bin',
OBEsFile = 'seast.bin',
OBEtFile = 'teast.bin',
OBSsFile = 'ssouth.bin',
OBStFile = 'tsouth.bin',
OBCS_monitorFreq= 2592000.,
OBCS_monSelect = 1,
&
# *****************************************
# Sponge Layer Parameters.
# *****************************************
&OBCS_PARM03
Urelaxobcsinner=432000.E0,
Urelaxobcsbound=43200.E0,
Vrelaxobcsinner=432000.E0,
Vrelaxobcsbound=43200.E0,
spongeThickness=30,
&
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20150331/747c5cdb/attachment.htm>
More information about the MITgcm-support
mailing list