[MITgcm-support] How to keep Prandtl number as one when 3D vertical eddy diffusivity is empolyed?

[Martin Losch]

Hi,

wouldn’t it be easier (and safer) to just add "diffKr(i,j,k,bi,bj)” to kappaRU/RV somewhere in calc_viscosity.F (or dynamics.F) in anology to what has been done in calc_3d_diffusitivity?

Martin

> On 19. Aug 2019, at 03:50, Ruan Xiaozhou <saberruan at hotmail.com> wrote:
> 
> Hello,
> 
> I have tried to do this and you could call CALC_3D_DIFFUSIVITY in the dynamics.F script and assign the 3D diffusivity field to viscosity. You’ll also need to define KappaRk in there too.
> 
> Cheers,
> Xiaozhou
> 
> 
> 
> #ifdef INCLUDE_CALC_DIFFUSIVITY_CALL
>       CALL CALC_3D_DIFFUSIVITY(
>      I         bi, bj, iMin, iMax, jMin, jMax,
>      I         GAD_TEMPERATURE, useGMredi, useKPP,
>      O         kappaRk,
>      I         myThid )
> #endif /* INCLUDE_CALC_DIFFUSIVITY_CALL */
> 
>       DO k = 1,Nr+1
>        ki = MIN(k,Nr)
> 
>        DO j = 1-OLy, sNy+OLy
>         DO i = 1-OLx, sNx+OLx
>          kappaRU(i,j,k) = kappaRk(i,j,ki)
>          kappaRV(i,j,k) = kappaRk(i,j,ki)
> 
> 
>         ENDDO
>        ENDDO
> 
>       ENDDO
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support