[MITgcm-support] OBCS_CHECK

[Martin Losch]

Hi Nadine,

in this case it makes sense to have a closer look what exactly causes the error. In obcs_check.F (and also in your STDERR.* files) you’ll find that this error is caused by a conflict between CPP-flags and the actual domain specification. You’ll have to look for something like

OBCS_CHECK: #undef ALLOW_OBCS_NORTH and
OBCS_CHECK: tile bi,bj=',bi,bj, ' has Northern OB

That is likely to help you find the problem. Without this information it’s difficult to help.

Martin


> On 20. Dec 2019, at 16:30, Nadine Steiger <Nadine.Steiger at uib.no> wrote:
> 
> Hei,
> 
> I am trying to change my bed topography of a so far working setup, but get this error message:
> 
> (PID.TID 0100.0001) *** ERROR *** OBCS_CHECK:    10 errors in tile OB set-up
> (PID.TID 0100.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run
> 
> My data.obcs looks as follows:
> # Open-boundaries
>  &OBCS_PARM01
> # This flag turns off checking and fixing problematic topography across
> # open boundaries.
>  OBCSfixTopo=.FALSE.,
>  OB_Jnorth= 1000*-1,
>  OB_Jsouth= 1000*1,
>  OB_Ieast = 480*-1,
>  OB_Iwest = 480*1,
>  useOBCSprescribe = .TRUE.,
>  OBNvFile = 'OB_NorthV.bin',
>  OBSvFile = 'OB_SouthV.bin',
>  OBCS_monitorFreq= 1200.,
>  OBCS_monSelect = 1,
>  &
> 
> I found an earlier email thread where it was suggested to describe the open boundaries from reading in hFacW/S
> (http://mailman.mitgcm.org/pipermail/mitgcm-support/2014-September/009468.html)
> 
> In this thread someone writes that:
> in fortran, the array sections are accessed (after reading in), as
> >       OBE = int(hFacW(NNX  ,1:NNY,1))
> >       OBW = int(hFacW(2    ,1:NNY,1))
> >       OBN = int(hFacS(1:NNX,NNY  ,1))
> >       OBS = int(hFacS(1:NNX,2    ,1))
> 
> 
> So, I retried it, this time prescribing the boundaries as
> OB_Jnorth = 1*0,719*-1,280*0,
> OB_Jsouth = 1*0,719*1,280*0,
> # OB_Ieast = 480*-1,
> # OB_Iwest = 480*1,
> where I got the values from hFacS.
> 
> In addition, I put a wall at the western and eastern boundary so i wrote in OBCS_OPTIONS.h:
> #undef ALLOW_OBCS_EAST
> #undef ALLOW_OBCS_WEST
> and commented OB_Ieast and OB_Iwest in data.obcs.
> 
> Even with this more detailed description of OB_Jnorth/south, I still get the same error message.
> 
> What strikes me is, that it worked before and I really only changed the bed topography (which is very idealized) slightly.
> The new topography works when I turn off the ice shelf package, the old topography worked with the ice shelf package.
> 
> What am I missing here?
> 
> Cheers,
> Nadine
> 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support