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

[Ruan Xiaozhou]

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190819/252327ba/attachment.html>