physical_mixing() (Source File: physical_mixing.F90)

INTERFACE:

    subroutine physical_mixing(F,AH,diffusivity,pm3d,pm2d)

DESCRIPTION:

Here, the physical tracer variance decay for the tracer $F$, ${ D}^{\mbox{\scriptsize phys}}\left(\langle F \rangle^2 \right)$, due to horizontal and vertical mixing is calculated as proposed in Burchard and Rennau (2008):

$${ D}^{\mbox{\scriptsize phys}}\left(F^2 \right)= 2 K_h \left(\par... ...ght)^2 +2 K_h \left(\partial_y F \right)^2 +2 K_v \left(\partial_z F \right)^2.$$ (133)

USES:

    use domain,       only: imin,imax,jmin,jmax,kmax,H,au,av
 #if defined(SPHERICAL) || defined(CURVILINEAR)
    use domain, only: dxu,dyv
 #else
    use domain, only: dx,dy
 #endif
    use variables_3d, only: dt,nuh,hn,ssen
 
    IMPLICIT NONE

INPUT PARAMETERS:

    REALTYPE, intent(in)  :: F(I3DFIELD)
    REALTYPE, intent(in)  :: AH
    REALTYPE, intent(in)  :: diffusivity
   !INPUT PARAMETERS
    REALTYPE, intent(out) :: pm3d(I3DFIELD)
    REALTYPE, intent(out) :: pm2d(I2DFIELD)

REVISION HISTORY:

    Original author(s): Hannes Rennau

LOCAL VARIABLES:

    REALTYPE                  :: dupper,dlower
    integer                   :: i,j,k
    REALTYPE                  :: aux(I3DFIELD)