pkg/ptracers/ptracers_convect.F

C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_convect.F,v 1.4 2007/01/09 22:27:12 jmc Exp $
C $Name:  $

#include "PTRACERS_OPTIONS.h"

CBOP
C !ROUTINE: PTRACERS_CONVECT

C !INTERFACE: ==========================================================
      SUBROUTINE PTRACERS_CONVECT( bi,bj,k,weightA,weightB,myThid )

C !DESCRIPTION:
C     do passive tracers convection

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PTRACERS_SIZE.h"
#include "PTRACERS_PARAMS.h"
#include "PTRACERS_FIELDS.h"

C !INPUT PARAMETERS: ===================================================
C  bi,bj,k :: tile indices
C  weightA :: weight for level K-1
C  weightB :: weight for level K
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)
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  none

#ifdef ALLOW_PTRACERS

C !LOCAL VARIABLES: ====================================================
C  iTracer :: tracer index
      INTEGER iTracer
CEOP

C Loop over tracers
      DO iTracer=1,PTRACERS_numInUse

       CALL CONVECTIVELY_MIXTRACER(
     I         bi,bj,k,weightA,weightB,
     U         pTracer(1-Olx,1-Oly,1,1,1,iTracer),
     I         myThid)

C End of tracer loop
      ENDDO

#endif /* ALLOW_PTRACERS */

      RETURN
      END
1	jmc	1.5	C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_convect.F,v 1.4 2007/01/09 22:27:12 jmc Exp $
2	jmc	1.2	C $Name: $
3	adcroft	1.1
4			#include "PTRACERS_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: PTRACERS_CONVECT
8
9			C !INTERFACE: ==========================================================
10			SUBROUTINE PTRACERS_CONVECT( bi,bj,k,weightA,weightB,myThid )
11
12			C !DESCRIPTION:
13	jmc	1.4	C do passive tracers convection
14	adcroft	1.1
15			C !USES: ===============================================================
16			IMPLICIT NONE
17			#include "SIZE.h"
18			#include "EEPARAMS.h"
19	jmc	1.3	#include "PTRACERS_SIZE.h"
20	jmc	1.5	#include "PTRACERS_PARAMS.h"
21			#include "PTRACERS_FIELDS.h"
22	adcroft	1.1
23			C !INPUT PARAMETERS: ===================================================
24			C bi,bj,k :: tile indices
25			C weightA :: weight for level K-1
26			C weightB :: weight for level K
27			C myThid :: thread number
28			INTEGER bi,bj,k
29			_RL weightA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30			_RL weightB(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
31			INTEGER myThid
32
33			C !OUTPUT PARAMETERS: ==================================================
34			C none
35
36			#ifdef ALLOW_PTRACERS
37
38			C !LOCAL VARIABLES: ====================================================
39			C iTracer :: tracer index
40			INTEGER iTracer
41			CEOP
42
43			C Loop over tracers
44			DO iTracer=1,PTRACERS_numInUse
45
46			CALL CONVECTIVELY_MIXTRACER(
47			I bi,bj,k,weightA,weightB,
48	jmc	1.2	U pTracer(1-Olx,1-Oly,1,1,1,iTracer),
49	adcroft	1.1	I myThid)
50
51			C End of tracer loop
52			ENDDO
53
54			#endif /* ALLOW_PTRACERS */
55
56			RETURN
57			END