[MITgcm-support] Obcs problems

Fri Jun 20 02:54:12 EDT 2008

Maria,

the other possibility is that you misspelled ALLOW_OBCS_PRESCRIBE in  
your OBCS_OPTIONS.h (this happens me all the time).

suggestions (this is how I normally try to track down a problem like  
yours):
1. check spelling of cpp-flags (because the cpp processor will not  
complain about that)
2. in your build directory, check where are the files in question  
linked from (is your OBCS_OPTIONS.h really the one in your code- 
directory)
3. delete obcs_check.f (lower case "f"!!!) and run "make  
obcs_check.f" to rebuild it. check (also by comparison to  
obcs_check.F) if this part of the code:
> #ifndef ALLOW_OBCS_PRESCRIBE
>       IF (useOBCSprescribe) THEN
>        WRITE(msgBuf,'(A)')
>      &  'OBCS_CHECK: ERROR: useOBCSprescribe = .TRUE. for'
>        CALL PRINT_ERROR( msgBuf , 1)
>        WRITE(msgBuf,'(A)')
>      &  'OBCS_CHECK: undef ALLOW_OBCS_PRESCRIBE'
>        CALL PRINT_ERROR( msgBuf , 1)
>        STOP 'ABNORMAL END: S/R OBCS_CHECK'
>       ENDIF
> #endif
is still in obcs_check.f and if *your* customized version of  
OBCS_OPTIONS.h has really been included into obcs_check.f (beginning  
of the file, you need to add some comments into the file to be able  
to see that).

Repeat 1-3 until the above code snipped does not appear in  
obcs_check.f any more. Then everything is OK and you should be able  
to get a working executable with "make makefile && make CLEAN && make  
depend && make".

Martin

On 19 Jun 2008, at 17:41, Rutkiewicz M.B. wrote:

> 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
>
>
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support