pkg/kpp/kpp_transport_s.F

C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_transport_s.F,v 1.7 2001/05/29 14:01:38 adcroft Exp $
C $Name:  $

#include "KPP_OPTIONS.h"

      subroutine KPP_TRANSPORT_S (
     I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
     I     KappaRS,
     U     df )
C     /==========================================================\
C     | o SUBROUTINE KPP_TRANSPORT_S                             |
C     |   Add non local KPP transport term (ghat) to diffusive   |
C     |   salinity flux.                                         |
C     |==========================================================|
C     | The nonlocal transport term is nonzero only for scalars  |
C     | in unstable (convective) forcing conditions.             |
C     | Note: KappaRS(i,j,k) is defined at the top of grid cell  |
C     |       k while KPPghat(i,j,k,bi,bj) is defined at the     |
C     |       bottom of grid cell k.                             |
C     |       For K=1, KappaRS(i,j,k) = 0 which insures no flux  |
C     |       through the surface.                               |
C     \==========================================================/
      IMPLICIT NONE

C     == GLobal variables ==
#include "SIZE.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "FFIELDS.h"
#include "KPP.h"
#include "KPP_PARAMS.h"

C     == Routine arguments ==
C     iMin,iMax,jMin,  - Range of points for which calculation
C     jMax,bi,bj,k,km1   results will be set.
C     KappaRS          - Total diffusion in vertical for S.
C     df               - Diffusive flux component work array.
      INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
      _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
      _RL df     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

#if defined(ALLOW_KPP) && defined(KPP_GHAT)

C     == Local variables ==
C     I, J - Loop counters
      INTEGER I, J

      DO j=jMin,jMax
       DO i=iMin,iMax
        df(i,j) = df(i,j) - _rA(i,j,bi,bj) * 
     &         KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * drF(1) *
     &         SurfaceTendencyS(i,j,bi,bj)
       ENDDO
      ENDDO

#endif /* ALLOW_KPP and KPP_GHAT */

      RETURN
      END
1	jmc	1.8	C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_transport_s.F,v 1.7 2001/05/29 14:01:38 adcroft Exp $
2	adcroft	1.7	C $Name: $
3	adcroft	1.1
4			#include "KPP_OPTIONS.h"
5
6			subroutine KPP_TRANSPORT_S (
7			I iMin,iMax,jMin,jMax,bi,bj,k,km1,
8	adcroft	1.7	I KappaRS,
9	adcroft	1.1	U df )
10			C /==========================================================\
11			C \| o SUBROUTINE KPP_TRANSPORT_S \|
12			C \| Add non local KPP transport term (ghat) to diffusive \|
13			C \| salinity flux. \|
14			C \|==========================================================\|
15			C \| The nonlocal transport term is nonzero only for scalars \|
16			C \| in unstable (convective) forcing conditions. \|
17			C \| Note: KappaRS(i,j,k) is defined at the top of grid cell \|
18			C \| k while KPPghat(i,j,k,bi,bj) is defined at the \|
19			C \| bottom of grid cell k. \|
20			C \| For K=1, KappaRS(i,j,k) = 0 which insures no flux \|
21			C \| through the surface. \|
22			C \==========================================================/
23			IMPLICIT NONE
24
25			C == GLobal variables ==
26			#include "SIZE.h"
27			#include "DYNVARS.h"
28			#include "EEPARAMS.h"
29			#include "PARAMS.h"
30			#include "GRID.h"
31			#include "FFIELDS.h"
32			#include "KPP.h"
33			#include "KPP_PARAMS.h"
34
35			C == Routine arguments ==
36			C iMin,iMax,jMin, - Range of points for which calculation
37			C jMax,bi,bj,k,km1 results will be set.
38			C KappaRS - Total diffusion in vertical for S.
39			C df - Diffusive flux component work array.
40			INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
41			_RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
42			_RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43
44	heimbach	1.2	#if defined(ALLOW_KPP) && defined(KPP_GHAT)
45	adcroft	1.1
46			C == Local variables ==
47			C I, J - Loop counters
48			INTEGER I, J
49
50	heimbach	1.2	DO j=jMin,jMax
51			DO i=iMin,iMax
52	adcroft	1.7	df(i,j) = df(i,j) - _rA(i,j,bi,bj) *
53	jmc	1.8	& KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * drF(1) *
54	heimbach	1.2	& SurfaceTendencyS(i,j,bi,bj)
55	adcroft	1.1	ENDDO
56	heimbach	1.2	ENDDO
57	adcroft	1.1
58	heimbach	1.2	#endif /* ALLOW_KPP and KPP_GHAT */
59	adcroft	1.1
60			RETURN
61			END