[MITgcm-support] OBCS north/south and periodic east/west?
Scott Bachman
sb965 at cam.ac.uk
Mon Apr 14 06:07:25 EDT 2014
Yun's solution fixed my problem. Thank you for your replies Dimitris and Yun!
Scott
On Apr 11, 2014, at 10:08 PM, Dimitris Menemenlis <dmenemenlis at gmail.com> wrote:
> north and southern boundaries are 25, right, not 250?
>
> On Apr 11, 2014, at 12:19 PM, Yun Xu <yunx at uci.edu> wrote:
>
>> 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
More information about the MITgcm-support
mailing list