[MITgcm-support] 2D setup

[David Hebert]

Hello all,

I would like to run a very simple, 2D simulation with tides on E/W 
boundary and a specific T/S structure, and flat bottom. It seems that 
T/S are crashing after the first timestep (the U  continues to 
timestep). T/S are initialized good, but go to zero after 1st timestep. 
I'm wondering if I'm missing something regarding the obcs setup.

In OBCS_OPTIONS.h I have specified #undef ALLOW_OBCS_NORTH and #undef 
ALLOW_OBCS_SOUTH. Then in data.obcs I have specified boundary files for 
U/T/S on East/West only. Relevant settings from data.obcs:

&OBCS_PARM01
OB_Ieast=1*-1,
OB_Iwest=1*1,
useOBCSprescribe=.TRUE.,
OBWuFile='OBWu.bin',
OBWtFile='OBWt.bin',
OBWsFile='OBWs.bin',
OBEuFile='OBEu.bin',
OBEtFile='OBEt.bin',
OBEsFile='OBEs.bin',


I thought it would be periodic in N/S in this case, but perhaps I am 
mistaken. Is there a conflict with useOBCSprescribe and 2D, such as 
applying tref/sref on N/S boundary?

Thanks for any suggestions you might have regarding the setup.

David