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):

$\displaystyle { 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)