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

Yun Xu yunx at uci.edu
Fri Apr 11 15:19:34 EDT 2014


Hi Scott, 

try:
OB_Jnorth=   240*-1,
OB_Jsouth=    240*1,

Yun

On Apr 11, 2014, at 12:05 PM, S. Bachman <sb965 at cam.ac.uk> wrote:

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

--
Yun Xu <yunx at uci.edu>
Department of Earth System Science 
University of California, Irvine
240P Rowland Hall
Irvine, CA, 92697








-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1014 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140411/db93fcf7/attachment.p7s>


More information about the MITgcm-support mailing list