[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