[MITgcm-devel] obcs & ptracers

[Martin Losch]

Hi Jean-Michel,

I noticed that, too, and I have now included the usePTRACERS in 
obcs_init_variables.F (but there still may be problems elsewhere).

Martin

PS. I won't be around too much in the next days, but I'll try to check 
out what you (will) have done.

On Dec 14, 2005, at 6:03 PM, Jean-Michel Campin wrote:

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