[MITgcm-support] open boundary prescribe
伊藤海彦
itohumihiko at eis.hokudai.ac.jp
Thu Oct 20 03:33:38 EDT 2016
Hi all,
I'm trying to prescribe boundary condition to East boundary.
But it's yet to work well.
I set useOBCSprescribe=.TRUE.,,but boundary condition files are not read.
So I think the format of the files is something wrong.
Could you tell me how to make boundary condition files in gendata.m
The script of my gendata.m to make these files for uVel is this.
nx=200;
ny=1;
nz=50;
t=10;
ieee='b';
prec='real*8';
for k=1,t
uEast=sin((2.0*3.14/t)*k)*ones(ny,nz,k);
end
fid=fopen('OBEastU.bin','w',ieee); fwrite(fid,uEast,prec); fclose(fid);
and this is my data, data.obcs,data.exf,
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef=50*0.5,
sRef=50*34.06,
bottomDragLinear=0.E-4,
#bottomDragQuadratic=2.6E-1,
no_slip_sides=.FALSE.,
no_slip_bottom=.FALSE.,
viscA4=0.0E4,
viscAh=25.E-3,
viscAz=1.E-3,
diffK4T=0.E4,
diffKhT=25.E-3,
diffKzT=1.E-3,
diffK4S=0.E4,
diffKhS=1.E-2,
diffKzS=1.E-5,
f0=0.937E-4,
fPrime=1.12E-4,
use3dCoriolis=.TRUE.,
beta=0.E-11,
eosType='MDJWF',
tAlpha=2.E-4,
sBeta =0.E-4,
rigidLid=.TRUE.,
#implicitFreeSurface=.TRUE.,
#implicSurfPress=0,
#implicDiv2Dflow=0,
hFacMin=0.05,
#use3Dsolver=.TRUE.,
nonHydrostatic=.TRUE.,
readBinaryPrec=64,
writeBinaryPrec=64,
#- not safe to use globalFiles in multi-processors runs
#globalFiles=.TRUE.,
tempAdvScheme=33,
staggerTimeStep=.TURE.,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=5000,
# cg2dUseMinResSol=1,
cg2dTargetResidual=1.E-12,
cg3dMaxIters=500,
cg3dTargetResidual=1.E-8,
&
# Time stepping parameters
&PARM03
nIter0=0,
nTimeSteps=20000,
#nTimeSteps=8640,
deltaT=60.0,
abEps=0.01,
pChkptFreq=0.0,
chkptFreq=0.0,
dumpFreq=6000.0,
#taveFreq=864000.0,
monitorSelect=2,
monitorFreq=6000.,
# for time dependent open boundary conditions, activate the following 3
lines:
periodicExternalForcing=.TRUE.,
externForcingPeriod=36000.,
externForcingCycle =36000.,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
#delXfile='dx.bin',
#dYspacing=200.,
delX=200*50.0,
delY=1*50.0,
delZ=50*10.0,
&
# Input datasets
&PARM05
bathyFile='topog_manaita_xz_0004.slope',
hydrogThetaFile='T_manaita_xz_0005.init',
#surfQfile='Qnet.forcing',
#hydrogThetaFile='T.pickup',
#uVelInitFile='U.pickup',
#pSurfInitFile='Eta.pickup',
# Open-boundaries
&OBCS_PARM01
OB_Ieast=1*-1,
#OB_Iwest=1*1,
#OB_Jnorth=200*-1,
#OB_Jsouth=200*1,
useOBCSprescribe=.TRUE.,
OBEuFile = 'OBEastU.bin',
#OBEvFile = 'OBEastV.bin',
#OBEwFile = 'OBEastW.bin',
#OBEtFile = 'OBEastT.bin',
#OBEsFile = 'OBEastS.bin',
#
#useOrlanskiEast=.TRUE.,
#useOrlanskiWest=.TRUE.,
#useOrlanskiNorth=.TRUE.,
#useOrlanskiSouth=.TRUE.,
#
useOBCSbalance=.TRUE.,
OBCS_balanceFacE = 0,
OBCS_balanceFacW = 0,
OBCS_balanceFacN = 1,
OBCS_balanceFacS = 1.,
OBCS_monitorFreq=6000.,
&
#
# *********************
# External Forcing Data
# *********************
&EXF_NML_01
#
useExfCheckRange = .TRUE.,
repeatPeriod = 36000.0,
exf_iprec = 64,
#
&
# *********************
&EXF_NML_02
&
# *********************
&EXF_NML_03
&
# *********************
&EXF_NML_04
&
# *********************
&EXF_NML_OBCS
obcsEstartdate1 = 20000101,
obcsEstartdate2 = 00000,
obcsEperiod = 3600.0,
&
Many thanks to any assistance!
Umihiko ito
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20161020/925657bc/attachment-0001.htm>
More information about the MITgcm-support
mailing list