[MITgcm-devel] explicit vertical mixing

[Samar Khatiwala]

Hi Martin,

Yes the check in ini_parms.F is one of the things that has to be  
commented out. There are several other
statements/calls in gad_calc_rhs.F and thermodynamics.F that I had to  
comment out/block. But it seems
to be working now. I had just hoped that there would be some magical  
combination of parameters (this
being the mitgcmUV) that would save me the effort.

Actually, and somewhat paradoxically, I'm doing this because I want  
to take veryyyyy long time steps
(like 1 year).

Samar

On Nov 9, 2007, at 9:39 AM, Martin Losch wrote:

> Hi Samar,
> didn't look at the code, but you can probably just comment out the  
> stop statement and deal with the large diffusivities (of the ivdc  
> parameter) by having a very small timestep, right?
>
> Martin
>
> On 7 Nov 2007, at 21:51, Samar Khatiwala wrote:
>
>> Hi
>>
>> Is it possible to treat vertical mixing explicitly in the MIT  
>> model? I understand that for most (99.999999%) of
>> problems this would be a crazy thing to want to do. But it so  
>> happens that I have a compelling need to
>> treat advection-diffusion completely explicitly, but WITHOUT using  
>> convective adjustment. Furthermore, I'd
>> like to do this only for ptracers and leave T and S alone (but  
>> ignore that complication for now).
>>
>> So I'd like to turn off implicitDiffusion, but still make use of  
>> the ivdc parameter in computing the vertical diffusivity.
>> From what I can make out, there seems to be no combination of  
>> parameters that would allow me to do this.
>> I think I know how to hack the code to make this happen, but  
>> perhaps someone here has a better idea (i.e.,
>> requiring less work).
>>
>> Thanks, Samar
>>
>>
>>
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel