[MITgcm-support] OBCS north/south and periodic east/west?

S. Bachman sb965 at cam.ac.uk
Fri Apr 11 15:05:48 EDT 2014


Hi all,

I seem to be having trouble getting MITgcm to run with OBCS enabled on 
the northern and southern boundaries and just leaving the east and west 
periodic.  I've undefined ALLOW_OBCS_EAST and ALLOW_OBCS_WEST in 
OBCS_OPTIONS.h, and am reading in the N/S boundaries from file.  My 
domain is 240x25x60 divided into 10 tiles (24 points each) in the zonal 
direction.  My data.obcs looks like this:

# ***************
# Open boundaries
# ***************
  &OBCS_PARM01
  OB_Jnorth=   24*-1,
  OB_Jsouth=    24*1,

#
  useOrlanskiNorth=.FALSE.,
  useOBCSsponge=.FALSE.,
  useOBCSbalance=.FALSE.,
  useOBCSprescribe=.TRUE.,
#
  OBNuFile='UVEL_north.bin',
  OBNvFile='VVEL_north.bin',
  OBNtFile='THETA_north.bin',
  OBNsFile='SALT_north.bin',

  OBSuFile='UVEL_south.bin',
  OBSvFile='VVEL_south.bin',
  OBStFile='THETA_south.bin',
  OBSsFile='SALT_south.bin',

but I'm getting an error at startup

(PID.TID 0000.0001) *** ERROR *** S/R OBCS_CHECK: Inside Mask and OB 
locations disagree :
(PID.TID 0000.0001) *** ERROR ***  from insideMask (bi,bj=  1,  1) 
expects OBE=    1 @ j=    1
(PID.TID 0000.0001) *** ERROR ***  from insideMask (bi,bj=  1,  1) 
expects OBW=   24 @ j=    1
(PID.TID 0000.0001) *** ERROR ***  from insideMask (bi,bj=  1,  1) 
expects OBE=    1 @ j=   25
(PID.TID 0000.0001) *** ERROR ***  from insideMask (bi,bj=  1,  1) 
expects OBW=   24 @ j=   25
(PID.TID 0000.0001) *** ERROR *** OBCS_CHECK:     8 errors in OB 
location vs Mask
(PID.TID 0000.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run

This is really confusing because OBCS isn't even enabled on the E/W 
boundaries, so I'm not sure whether I should be adding something to my 
data.obcs or if something else is going on here.  I also noticed that it 
is calling for OBE and OBW at the corner points of the tile only, if 
that means anything.

Any help would be greatly appreciated, thanks!
Scott




More information about the MITgcm-support mailing list