[MITgcm-devel] obcs & ptracers

[Jean-Michel Campin]

Hi Martin,

I will check-in some modifications in obcs_init_variables.F
and when you have time, let me know what you think of it,
and then we can change it again/back if needed.

I noticed also a little problem:
all ptracers stuff should be within IF (usePTRACERS ) THEN / ENDIF.
I know that there are packages that don't work that way
(and if unfortunate), but this package does.
Right now, if I compile pkg ptracers but don't use it,
for instance PTRACERS_numInUse is not initialised
(I checked this, and managed to produce an error in
 obcs_init_variables.F).

Cheers,

Jean-Michel

On Wed, Dec 14, 2005 at 04:26:50PM +0100, Martin Losch wrote:
> Hi Jean-Michel,
> 
> I have to admit, that I never tested the restarts with obcs and 
> ptracers. I think you do need the to apply the open boundary conditions 
> at the very beginning, but not really after restarts (I was probably to 
> lazy to think about that). But the error is a different one: there are 
> exchanges missing for the ptracers (they are there for uVel, vVel, 
> theta, and salt). This is still the old problem, that in 
> obcs_prescribe_read, the "exchanges" for vertical slab are missing, so 
> that the obcs for ptracers are not correct on the halos. I should 
> probably look into that very soon, but I just don't have the time at 
> the moment :-(
> 
> Martin
>