[MITgcm-support] Obcs problems
Rutkiewicz M.B.
mbr1x07 at soton.ac.uk
Thu Jun 19 11:41:14 EDT 2008
Hi!
I tried that and I still get the same error.
Maria
On 19/06/2008 14:23, "Matthew Mazloff" <mmazloff at MIT.EDU> wrote:
Hi Maria,
I think its possible that the OBCS_OPTIONS.h in your build folder is
not linked to the one in your code folder. Try compiling again after
either cleaning out your build folder or at least deleting
OBCS_OPTIONS from it and relinking to the one in your code folder.
-matt
On Jun 19, 2008, at 5:10 AM, Rutkiewicz M.B. wrote:
> Hello!
> I am trying to add obcs to my configuration but I am having some
> problems.
> I am running a zonal channel configuration with solid boundaries
> North and
> South. Now I am trying to open the Northern boundary by prescribing a
> temperature profile, which I have created with gendata.m.
>
> I added obcs to packages.conf (OBCS_OPTIONS.h being included in my
> directory
> code/) and the compilation was fine, but when I try running the
> executable,
> I get the following error message:
>
> OBCS_CHECK: ERROR: useOBCSprescribe = .TRUE. for
> OBCS_CHECK: undef ALLOW_OBCS_PRESCRIBE
> ABNORMAL END: S/R OBCS_CHECK
>
> I am not sure what I have done wrong, I switched on define
> ALLOW_OBCS_NORTH
> and define ALLOW_OBCS_PRESCRIBE in OBCS_OPTIONS.h and see my data and
> data.obcs files below!
>
> What have I missed?
>
> Cheers,
> Maria
>
> data:
> # ====================
> # | Model parameters |
> # ====================
> #
> # Continuous equation parameters
> &PARM01
> #
> # gives the reference temperature profile.
> #
> tRef = 10*2.0,
> sRef = 10*34.5,
> # horizontal and vertical eddy viscosity
> viscAh=1000.,
> viscAz=1.E-4,
> no_slip_sides=.FALSE.,
> no_slip_bottom=.TRUE.,
> # viscA4=0.E12,
> # horizontal and vertical eddy diffusivity
> diffKhT=1000.,
> diffKzT=1.E-4,
> f0=-1e-4,
> beta=1.E-11,
> tAlpha=2.4e-4,
> sBeta =7.6e-4,
> gbaro=9.81,
> gravity=9.81,
> rigidLid=.FALSE.,
> implicitFreeSurface=.TRUE.,
> eosType='LINEAR',
> nonHydrostatic=.FALSE.,
> useSingleCpuIO = .TRUE.,
> readBinaryPrec=64,
> &
>
> # Elliptic solver parameters
> &PARM02
> cg2dMaxIters=1000,
> cg2dTargetResidual=1.E-7,
> cg3dMaxIters=40,
> cg3dTargetResidual=1.E-9,
> &
>
> # Time stepping parameters
> &PARM03
> nIter0=0,
> nTimeSteps=8766,
> deltaT=3600.0,
> abEps=0.1,
> tauThetaClimRelax = 2629800.,
> cAdjFreq = 3600.0,
>
> # create a pickup file each day
> # pChkpt: *permanent pickup file
> pChkptFreq=31557600.0,
> #alternating pickup file: each 1/6th day:
> chkptFreq=0.0,
> dumpFreq=2629800.0,
> monitorFreq=0,
> # outputTypesInclusive=.FALSE.,
> &
>
> # Gridding parameters
> &PARM04
> usingCartesianGrid=.TRUE.,
> usingSphericalPolarGrid=.FALSE.,
> delX=30*50000.,
> delY=15*50000.,
> delZ=10*300.,
> &
>
> # Input datasets
> &PARM05
> # to specify the initial temperature gradient
> thetaClimFile='mytemp',
> bathyFile='topog.channel',
> zonalWindFile='mywind',
> &
>
>
> data.obcs:
>
> # Open-boundaries
> &OBCS_PARM01
>
> #specify northern open boundary
> OB_Jnorth=30*-1,
> # OB_Jsouth=80*1,
> # OB_Ieast=42*-1,
> # OB_Iwest=42*1,
>
> #useOrlanskiNorth=.TRUE.,
> #useOrlanskiSouth=.TRUE.,
> #useOrlanskiEast=.TRUE.,
> #useOrlanskiWest=.FALSE.,
>
> useOBCSprescribe = .TRUE.,
> OBNtFile='myobn_t',
>
> # OBNuFile = 'OBmeridU.bin',
> # OBSuFile = 'OBmeridU.bin',
> # OBWuFile = 'OBzonalU.bin',
> # OBEuFile = 'OBzonalU.bin',
> # OBWsFile = 'OBzonalS.bin',
> # OBWptrFile(1) = 'OBzonalS.bin',
> &
>
> # Orlanski parameters
> &OBCS_PARM02
> #Cmax=0.45,
> #cVelTimeScale=1000.,
> &
>
> --
> Miss Maria Rutkiewicz
>
> PhD Student Physical Oceanography
> Room 186/07
> National Oceanography Centre, Southampton
> University of Southampton Waterfront Campus
> European Way
> Southampton SO14 3ZH
> United Kingdom
>
> email: mbr1x07 at soton.ac.uk
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list