model/src/convectively_mixtracer.F

C $Header: /u/gcmpack/MITgcm/model/src/convectively_mixtracer.F,v 1.2 2002/02/26 20:34:14 adcroft Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C !ROUTINE: CONVECTIVELY_MIXTRACER
C !INTERFACE:
      SUBROUTINE CONVECTIVELY_MIXTRACER(
     I                              bi,bj,k,weightA,weightB,
     U                              Tracer,
     I                              myThid)
C !DESCRIPTION:
C Mixes a tracer over two layers according to the weights pre-calculated
C as a function of stability.
C 
C Mixing is represented by:
C                       T(k-1) = T(k-1) + B * ( T(k) - T(k-1) )
C                       T(k)   = T(k)   + A * ( T(k-1) - T(k) )

C !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"

C !INPUT/OUTPUT PARAMETERS:
C     bi,bj,k - indices
C     weightA - weight for level K
C     weightB - weight for level K+1
C     myThid  - thread number
      INTEGER bi,bj,k
      _RL weightA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL weightB(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL Tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myThid

#ifdef INCLUDE_CONVECT_CALL

C !LOCAL VARIABLES:
C     i,j      :: Loop counter
C     delTrac  :: Difference between tracer in each layer
      INTEGER i,j
      _RL delTrac
CEOP

      DO j=1-Oly,sNy+Oly
       DO i=1-Olx,sNx+Olx

        delTrac=Tracer(i,j,k,bi,bj)-Tracer(i,j,k-1,bi,bj)
        Tracer(i,j,k-1,bi,bj)=Tracer(i,j,k-1,bi,bj)
     &                       +weightA(i,j)*delTrac
        Tracer(i,j,k,bi,bj)=Tracer(i,j,k,bi,bj)
     &                       -weightB(i,j)*delTrac

       ENDDO
      ENDDO

#endif /* INCLUDE_CONVECT_CALL */

      RETURN
      END
1	heimbach	1.2.4.1	C $Header: /u/gcmpack/MITgcm/model/src/convectively_mixtracer.F,v 1.2 2002/02/26 20:34:14 adcroft Exp $
2	adcroft	1.2	C $Name: $
3	adcroft	1.1
4			#include "CPP_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: CONVECTIVELY_MIXTRACER
8			C !INTERFACE:
9			SUBROUTINE CONVECTIVELY_MIXTRACER(
10			I bi,bj,k,weightA,weightB,
11			U Tracer,
12			I myThid)
13			C !DESCRIPTION:
14			C Mixes a tracer over two layers according to the weights pre-calculated
15			C as a function of stability.
16			C
17			C Mixing is represented by:
18			C T(k-1) = T(k-1) + B * ( T(k) - T(k-1) )
19			C T(k) = T(k) + A * ( T(k-1) - T(k) )
20
21	adcroft	1.2	C !USES:
22	adcroft	1.1	IMPLICIT NONE
23			#include "SIZE.h"
24			#include "EEPARAMS.h"
25
26	adcroft	1.2	C !INPUT/OUTPUT PARAMETERS:
27	adcroft	1.1	C bi,bj,k - indices
28			C weightA - weight for level K
29			C weightB - weight for level K+1
30			C myThid - thread number
31			INTEGER bi,bj,k
32			_RL weightA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
33			_RL weightB(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
34			_RL Tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
35			INTEGER myThid
36
37			#ifdef INCLUDE_CONVECT_CALL
38
39	adcroft	1.2	C !LOCAL VARIABLES:
40	adcroft	1.1	C i,j :: Loop counter
41			C delTrac :: Difference between tracer in each layer
42			INTEGER i,j
43			_RL delTrac
44			CEOP
45
46			DO j=1-Oly,sNy+Oly
47			DO i=1-Olx,sNx+Olx
48
49			delTrac=Tracer(i,j,k,bi,bj)-Tracer(i,j,k-1,bi,bj)
50			Tracer(i,j,k-1,bi,bj)=Tracer(i,j,k-1,bi,bj)
51			& +weightA(i,j)*delTrac
52			Tracer(i,j,k,bi,bj)=Tracer(i,j,k,bi,bj)
53			& -weightB(i,j)*delTrac
54
55			ENDDO
56			ENDDO
57
58			#endif /* INCLUDE_CONVECT_CALL */
59
60			RETURN
61			END