[MITgcm-devel] obcs balance

[Martin Losch]

Hi all,
do you agree with this (obcs_calc.F)?
> [...], it seems strange that  the loop ranges of the loop for the  
> calculation of the average flux (the first loop in each section),  
> to be substracted from the velocity field in the second loop, range  
> from I=1-Olx,sNx+Olx (for north/south) and J=1-Oly,sNy+Oly (for  
> east/west).
> Naively I would think that they should range only over the  
> computational domain *without* the overlaps, so I=1,sNx and  
> J=1,sNy, but I may have overlooked something.
It seems to reduce 名 徐 's problem (how on earth do you transcribe  
this name?) dramatically, but since none of the verification  
experiments test this code my changing this in the repository may  
have unexpected effects for some users of thie code. Should I go  
ahead and do it anyway (unless you think it's correct)?

Martin