[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