/[MITgcm]/MITgcm/pkg/kpp/kpp_transport_s.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_transport_s.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.9 by heimbach, Sat Jul 13 03:12:30 2002 UTC revision 1.14 by dimitri, Tue Dec 18 09:05:59 2007 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "KPP_OPTIONS.h"  #include "KPP_OPTIONS.h"
5    #ifdef ALLOW_GMREDI
6    #include "GMREDI_OPTIONS.h"
7    #endif
8    
9        subroutine KPP_TRANSPORT_S (        subroutine KPP_TRANSPORT_S (
10       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
11       I     KappaRS,       O     df,
12       U     df )       I     myTime, myIter, myThid )
13  C     /==========================================================\  C     /==========================================================\
14  C     | o SUBROUTINE KPP_TRANSPORT_S                             |  C     | o SUBROUTINE KPP_TRANSPORT_S                             |
15  C     |   Add non local KPP transport term (ghat) to diffusive   |  C     |   Add non local KPP transport term (ghat) to diffusive   |
# Line 14  C     |   salinity flux. Line 17  C     |   salinity flux.
17  C     |==========================================================|  C     |==========================================================|
18  C     | The nonlocal transport term is nonzero only for scalars  |  C     | The nonlocal transport term is nonzero only for scalars  |
19  C     | in unstable (convective) forcing conditions.             |  C     | in unstable (convective) forcing conditions.             |
20  C     | Note: KappaRS(i,j,k) is defined at the top of grid cell  |  C     | Note: KPPdiffKzS(-,k) is defined at the top of grid cell |
21  C     |       k while KPPghat(i,j,k,bi,bj) is defined at the     |  C     |       k while KPPghat(i,j,k,bi,bj) is defined at the     |
22  C     |       bottom of grid cell k.                             |  C     |       bottom of grid cell k.                             |
23  C     |       For K=1, KappaRS(i,j,k) = 0 which insures no flux  |  C     |       For K=1, KPPdiffKzS(-,k) = 0 which insures no flux |
24  C     |       through the surface.                               |  C     |       through the surface.                               |
25  C     \==========================================================/  C     \==========================================================/
26        IMPLICIT NONE        IMPLICIT NONE
27    
28  C     == GLobal variables ==  C     == GLobal variables ==
29  #include "SIZE.h"  #include "SIZE.h"
 #include "DYNVARS.h"  
30  #include "EEPARAMS.h"  #include "EEPARAMS.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32  #include "GRID.h"  #include "GRID.h"
33  #include "FFIELDS.h"  #include "FFIELDS.h"
34  #include "KPP.h"  #include "KPP.h"
35  #include "KPP_PARAMS.h"  #include "KPP_PARAMS.h"
36    #ifdef ALLOW_GMREDI
37    # include "GMREDI.h"
38    #endif
39    #ifdef ALLOW_SALT_PLUME
40    # include "SALT_PLUME.h"
41    #endif
42    
43  C     == Routine arguments ==  C     == Routine arguments ==
44  C     iMin,iMax,jMin,  - Range of points for which calculation  C     iMin,iMax,jMin,  - Range of points for which calculation
45  C     jMax,bi,bj,k,km1   results will be set.  C     jMax,bi,bj,k,km1   results will be set.
 C     KappaRS          - Total diffusion in vertical for S.  
46  C     df               - Diffusive flux component work array.  C     df               - Diffusive flux component work array.
47    C     myTime :: Current time in simulation
48    C     myIter :: Current iteration number in simulation
49    C     myThid :: My Thread Id. number
50        INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1        INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
       _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)  
51        _RL df     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL df     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52          _RL     myTime
53          INTEGER myIter
54          INTEGER myThid
55    
56  #if defined(ALLOW_KPP) && defined(KPP_GHAT)  #if defined(ALLOW_KPP) && defined(KPP_GHAT)
57    
58  C     == Local variables ==  C     == Local variables ==
59  C     I, J - Loop counters  C     I, J - Loop counters
60        INTEGER I, J        INTEGER I, J
61          _RL tmpFac
62    
63    #ifdef ALLOW_GMREDI
64    
65          IF ( useGMRedi .AND. KPP_ghatUseTotalDiffus ) THEN
66            tmpFac = 1. _d 0
67          ELSE
68            tmpFac = 0. _d 0
69          ENDIF
70        DO j=jMin,jMax        DO j=jMin,jMax
71         DO i=iMin,iMax         DO i=iMin,iMax
72          df(i,j) = df(i,j) - _rA(i,j,bi,bj) *          df(i,j) = - _rA(i,j,bi,bj)
73       &         KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * drF(1) *       &        *( KPPdiffKzS(i,j,k,bi,bj) + tmpFac*Kwz(i,j,k,bi,bj) )
74       &         SurfaceTendencyS(i,j,bi,bj)       &        *KPPghat(i,j,km1,bi,bj)
75         &        *( surfaceForcingS(i,j,bi,bj)
76    #ifdef ALLOW_SALT_PLUME
77         &          + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst
78         &          *(1.-KPPplumefrac(i,j,bi,bj))
79    #endif /* ALLOW_SALT_PLUME */
80         &         )
81         ENDDO         ENDDO
82        ENDDO        ENDDO
83    
84    #else /* ALLOW_GMREDI */
85    
86          DO j=jMin,jMax
87           DO i=iMin,iMax
88            df(i,j) = - _rA(i,j,bi,bj)
89         &        *KPPdiffKzS(i,j,k,bi,bj)
90         &        *KPPghat(i,j,km1,bi,bj)
91         &        *( surfaceForcingS(i,j,bi,bj)
92    #ifdef ALLOW_SALT_PLUME
93         &          + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst
94         &          *(1.-KPPplumefrac(i,j,bi,bj))
95    #endif /* ALLOW_SALT_PLUME */
96         &         )
97           ENDDO
98          ENDDO
99    
100    #endif /* ALLOW_GMREDI */
101    
102    
103  #endif /* ALLOW_KPP and KPP_GHAT */  #endif /* ALLOW_KPP and KPP_GHAT */
104    
105        RETURN        RETURN
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22