[MITgcm-support] open boundary prescribe

[伊藤海彦]

Hi,Matt

I didn't realize there is a such useful option (useOBCStides = .TRUE.,) in
obsc!

I did it.
Thank you very much!

Umihiko ito

>(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_obcsgives 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20161022/8b65ea55/attachment-0001.htm>