[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