[MITgcm-support] ptracer surface boundary condition

[Samar Khatiwala]

Hi Jean-Michel

I have of course tried this but wasn't too happy with the results. If 
the restoring
is too strong I get over/undershoots. If its weak, the BC is not what I 
really
want to apply.

Thanks, Samar

On May 23, 2005, at 1:51 PM, Jean-Michel Campin wrote:

> Hi Samar,
>
> Generally, I put a strong restoring (short time-scale) and
> this works for me.
>
> Jean-Michel
>
> On Mon, May 23, 2005 at 01:46:59PM -0400, samar khatiwala wrote:
>> Hello
>>
>> I would like to perform some ptracer simulations in which the
>> ptracer concentration is prescribed at the surface. This could
>> be a time-dependent boundary condition. Unfortunately, I don't
>> see any simple way to use existing code to prescribe such a BC.
>> Simply setting pTracer(...,k=1,..)=BC does not work all that well.
>> To make this work robustly (and consistent with the time stepping) I 
>> have
>> to intervene at multiple locations in the code to (re)set pTracer and 
>> gPtr
>> to the BC.
>>
>> I was wondering if someone has figured out how to do this in a cleaner
>> and more consistent manner. Perhaps this could be incorporated into 
>> the
>> gchem package since it seems like a useful capability to have.
>>
>> Thanks
>> Samar
>>
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support