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

Contents 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 - (show annotations) (download)
Thu Aug 25 16:22:17 2005 UTC (18 years, 6 months ago) by dimitri
Branch: release1_50yr
Changes since 1.1: +62 -0 lines
adding ecco1x1 verification/global_with_CFC11 experiment

1 C $Header:
2 C $Name: release1_50yr $
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