/[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.5 by heimbach, Mon Jan 29 20:09:23 2001 UTC revision 1.16 by jmc, Wed Jan 20 01:32:50 2010 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    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     maskC,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
16  C     |   salinity flux.                                         |  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     maskC            - Land mask for salt cells.  
 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
       _RS maskC  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _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    #ifdef ALLOW_GMREDI
62          _RL tmpFac
63    #endif
64    #ifdef ALLOW_SALT_PLUME
65          _RL tmpFac1
66    #endif
67    
68    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
69    
70    #ifdef ALLOW_SALT_PLUME
71          IF ( useSALT_PLUME ) THEN
72             tmpFac1 = 1. _d 0
73          ELSE
74             tmpFac1 = 0. _d 0
75          ENDIF
76    #endif /* ALLOW_SALT_PLUME */
77    
78    #ifdef ALLOW_GMREDI
79    
80          IF ( useGMRedi .AND. KPP_ghatUseTotalDiffus ) THEN
81            tmpFac = 1. _d 0
82          ELSE
83            tmpFac = 0. _d 0
84          ENDIF
85          DO j=jMin,jMax
86           DO i=iMin,iMax
87            df(i,j) = - _rA(i,j,bi,bj)
88         &        *( KPPdiffKzS(i,j,k,bi,bj) + tmpFac*Kwz(i,j,k,bi,bj) )
89         &        *KPPghat(i,j,km1,bi,bj)
90         &        *( surfaceForcingS(i,j,bi,bj)
91    #ifdef ALLOW_SALT_PLUME
92         &          + tmpFac1*saltPlumeFlux(i,j,bi,bj)*recip_rhoConst
93         &          *(1.-KPPplumefrac(i,j,bi,bj))
94    #endif /* ALLOW_SALT_PLUME */
95         &         )
96           ENDDO
97          ENDDO
98    
99    #else /* ALLOW_GMREDI */
100    
101        DO j=jMin,jMax        DO j=jMin,jMax
102         DO i=iMin,iMax         DO i=iMin,iMax
103          df(i,j) = df(i,j) - _rA(i,j,bi,bj) * maskC(i,j) *          df(i,j) = - _rA(i,j,bi,bj)
104       &         KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * delZ(1) *       &        *KPPdiffKzS(i,j,k,bi,bj)
105       &         SurfaceTendencyS(i,j,bi,bj)       &        *KPPghat(i,j,km1,bi,bj)
106         &        *( surfaceForcingS(i,j,bi,bj)
107    #ifdef ALLOW_SALT_PLUME
108         &          + tmpFac1*saltPlumeFlux(i,j,bi,bj)*recip_rhoConst
109         &          *(1.-KPPplumefrac(i,j,bi,bj))
110    #endif /* ALLOW_SALT_PLUME */
111         &         )
112         ENDDO         ENDDO
113        ENDDO        ENDDO
114    
115    #endif /* ALLOW_GMREDI */
116    
117  #endif /* ALLOW_KPP and KPP_GHAT */  #endif /* ALLOW_KPP and KPP_GHAT */
118    
119        RETURN        RETURN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22