[MITgcm-support] open boundary prescribe
Matthew Mazloff
mmazloff at ucsd.edu
Thu Oct 20 12:32:52 EDT 2016
(1) What do you mean by “are not read”? What is the error?
(2) if you are just trying to add a tide you can do that with in data.obcs:
useOBCStides = .TRUE.,
and then give tidal period and amplitude and phase files.
MITgcm/verification/seaice_obcs
gives an example
Matt
> On Oct 20, 2016, at 12:33 AM, 伊藤海彦 <itohumihiko at eis.hokudai.ac.jp> wrote:
>
> 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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list