[MITgcm-support] obcs with irregular boundary

[De Rydt, Jan]

Hi all,

I spent some time with the OBCS package trying to implement an irregular eastern boundary, i.e., OB_Ieast(j) is varying with j (see data.obcs below). I was unsuccessful in getting this to run, and the kind of errors I get look like this:

(PID.TID 0108.0001) *** ERROR *** S/R OBCS_CHECK: Inside Mask and OB locations disagree :
(PID.TID 0108.0001) *** ERROR ***  from insideMask (bi,bj=  1,  1) expects OBE=   16 @ j=    2 

What puzzles me is that the insideMaskFile seems to be defined from the OB list of indices, whereas the errors suggest that the interior mask is inconsistent with the OB list of indices...

Perhaps somebody has more insight in what is going wrong?

Thanks,
Jan


# ***************
# data.obcs
# ***************

&OBCS_PARM01
 OB_Iwest = 90*1.,
 OB_Ieast = -15.0,-15.0,-11.0,-10.0,-8.0,-7.0,-6.0,-5.0,-5.0,-4.0,-4.0,
            -3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,-3.0,
            -3.0,-4.0,-4.0,-5.0,-5.0,-6.0,-7.0,-8.0,-9.0,-10.0,-11.0,-12.0,
            -13.0,-14.0,-15.0,-15.0,-16.0,-17.0,-18.0,-18.0,-18.0,-19.0,-19.0,
            -21.0,-22.0,-23.0,-21.0,-21.0,-21.0,-21.0,-21.0,-21.0,-22.0,-22.0,
            -23.0,-24.0,-19.0,-18.0,-17.0,-16.0,-16.0,-15.0,-15.0,-14.0,-14.0,
            -14.0,-14.0,-14.0,-14.0,-14.0,-14.0,-13.0,-13.0,-13.0,-14.0,-14.0,
            -14.0,-14.0,-14.0,-15.0,-15.0,-16.0,-16.0,-17.0,-18.0,-19.0,-19.0,
#
 useOBCSprescribe = .true.,
#
 OBWtFile='OBWt.bin',
 OBWsFile='OBWs.bin',
#
 OBEtFile='OBEt.bin',
 OBEsFile='OBEs.bin',
 OBEuFile='OBEu.bin',
 OBEvFile='OBEv.bin',
#
 OBCS_balanceFacW = 1.,
 OBCS_balanceFacE = 0.,
#
&

--
British Antarctic Survey
High Cross, Madingley Road
Cambridge, CB3 0ET, UK
+44 (0)1223 221656

-- 
This message (and any attachments) is for the recipient only. NERC
is subject to the Freedom of Information Act 2000 and the contents
of this email and any reply you make may be disclosed by NERC unless
it is exempt from release under the Act. Any material supplied to
NERC may be stored in an electronic records management system.