[MITgcm-support] Modeling of the surface waves

王刚 wangg at fio.org.cn
Sun Oct 28 21:33:11 EDT 2012


Dear all,

     I try to simulate the surface gravity waves in a small  water tank using non-hydrostatic model. The eastern boundary is sponge layers, and western boundary is prescribed periodical changing U and W. They are sinusoidal or cosinusoidal waves exponentially decayed with depth  [U=a*c*k*exp(k*z)*cos(k*x-sigma*t   and  w=-a*c*k*exp(k*z)*sin(k*x-sigma*t),   where  a=0.2m is the amplitude of sea lever,c=2m/s  is the phase velocity,k=2*pi/10m wavenumber,z=12.6m depth,sigma frequency]. However, the eta (sea surface lever) is very small (almost zero) in my result, and I can not find any clue of surface waves. Did I set wrong perameters?

. The following are my data 

# Continuous equation parameters
 &PARM01
 tRef =100*20,  
 sRef = 100*10,
 f0=0,           
 viscAz=1.0E-6,
 viscAh=1.0E-6,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 viscA4=0.E12,
 diffKhT=1.44E-7,
 diffKzT=1.44E-7,
#diffK4T=5.E3,
 diffKhS=1.44E-7,
 diffKzS=1.44E-7,
 beta=0.E-11,
 tAlpha=2.E-4,    
 sBeta =7.4E-4,
 gravity=9.81,
 gBaro=9.81,
 rigidLid=.FALSE.,
 nonHydrostatic=.TRUE.,
 implicitFreeSurface=.TRUE.,
 eosType='LINEAR',
 readBinaryPrec=64,
#- not safe to use globalFiles in multi-processors runs
#globalFiles=.TRUE.,
 writeBinaryPrec=64,
 writeStatePrec=64,
#saltStepping=.FALSE.,
 hFacMin=0.1,
 implicSurfPress=0.5,
 implicDiv2DFlow=0.5,
#nonlinFreeSurf=3,
#hFacInf=0.2,
#hFacSup=2.0,
 exactConserv=.TRUE.
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-10,
 cg3dMaxIters=100,
 cg3dTargetResidual=1.E-10,
 &

# Time stepping parameters
 &PARM03
 niter0=0,
 nTimeSteps=250000,
 deltaT=0.02,
 deltaTtracer=0.02,
 abEps=0.0001,
 pChkptFreq=500,
 chkptFreq=500,
 dumpFreq=1,
 monitorFreq=0.1,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=0.02,
 externForcingCycle=5,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
 delX=1000*0.5,
 delY=5,
 delZ=10*0.05,10*0.06,10*0.07,10*0.08,10*0.1,10*0.12,10*0.15,10*0.18,10*0.2,10*0.25,
&

# Input datasets
 &PARM05
 hydrogThetaFile='T.init',
#hydrogSaltFile='S.init',
 bathyFile='topog',
 &


  , and data.obcs  

        
# Open-boundaries
 &OBCS_PARM01
 OB_Ieast=-1,
 OB_Iwest=1,
 useOrlanskiEast =.TRUE.,
 useOrlanskiWest =.FALSE.,
 useOBCSbalance=.TRUE.,
 useOBCSprescribe=.TRUE.,
#useOBCSsponge=.TRUE.,
#OBEuFile='uEast.bin',
 OBWuFile='uWest.bin',
#OBWwFile='wWest.bin',
#OBWetaFile='etaWest.bin',
&

&OBCS_PARM02
 cmax=0.45D0,
 cvelTimeScale=5.D0,
&
 
# *****************************************
# Sponge Layer Parameters.
# *****************************************
 &OBCS_PARM03
# Urelaxobcsinner  = 0,
#Urelaxobcsbound  = 5,
#spongeThickness  = 20,
 &
# end of the parameters.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20121029/bca16371/attachment.htm>


More information about the MITgcm-support mailing list