/[MITgcm]/MITgcm/verification/global_with_CFC11/code1x1/kpp_transport_tr1.F
ViewVC logotype

Annotation of /MITgcm/verification/global_with_CFC11/code1x1/kpp_transport_tr1.F

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


Revision 1.1.2.1 - (hide annotations) (download)
Thu Aug 25 16:22:17 2005 UTC (18 years, 7 months ago) by dimitri
Branch: release1_50yr
Changes since 1.1: +62 -0 lines
adding ecco1x1 verification/global_with_CFC11 experiment

1 dimitri 1.1.2.1 C
2     C $Name: $
3    
4     #include "KPP_OPTIONS.h"
5    
6     subroutine KPP_TRANSPORT_TR1 (
7     I iMin,iMax,jMin,jMax,bi,bj,k,km1,
8     I KappaRS,
9     U df )
10     C /==========================================================\
11     C | o SUBROUTINE KPP_TRANSPORT_TR1 |
12     C | Add non local KPP transport term (ghat) to diffusive |
13     C | passive tracer 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     #include "TR1.h"
35    
36     C == Routine arguments ==
37     C iMin,iMax,jMin, - Range of points for which calculation
38     C jMax,bi,bj,k,km1 results will be set.
39     C KappaRS - Total diffusion in vertical for S.
40     C df - Diffusive flux component work array.
41     INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
42     _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
43     _RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44    
45     #if defined(ALLOW_KPP) && defined(KPP_GHAT)
46    
47     C == Local variables ==
48     C I, J - Loop counters
49     INTEGER I, J
50    
51     DO j=jMin,jMax
52     DO i=iMin,iMax
53     df(i,j) = df(i,j) - _rA(i,j,bi,bj) *
54     & KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * drF(1) *
55     & SurfaceTendencyTr1(i,j,bi,bj)
56     ENDDO
57     ENDDO
58    
59     #endif /* ALLOW_KPP and KPP_GHAT */
60    
61     RETURN
62     END

  ViewVC Help
Powered by ViewVC 1.1.22