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

Contents of /MITgcm/verification/global_with_CFC11/code1x1/calc_gtr1.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, 7 months ago) by dimitri
Branch: release1_50yr
Changes since 1.1: +114 -0 lines
adding ecco1x1 verification/global_with_CFC11 experiment

1 C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/global_with_CFC11/code1x1/Attic/calc_gtr1.F,v 1.1.2.1 2005/08/25 16:22:17 dimitri Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: CALC_GTR1
8 C !INTERFACE:
9 SUBROUTINE CALC_GTR1(
10 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
11 I xA,yA,uTrans,vTrans,rTrans,maskUp,
12 I KappaRT,
13 U fVerT,
14 I myTime,myIter,myThid )
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | SUBROUTINE CALC_GTR1
18 C | o Calculate the passive tracer tendency terms.
19 C *==========================================================*
20
21 C !USES:
22 IMPLICIT NONE
23 C == GLobal variables ==
24 #include "SIZE.h"
25 #include "DYNVARS.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28 #include "GAD.h"
29 #ifdef ALLOW_PASSIVE_TRACER
30 #include "TR1.h"
31 #endif
32
33 C !INPUT/OUTPUT PARAMETERS:
34 C == Routine arguments ==
35 C fVerT :: Flux of temperature (T) in the vertical
36 C direction at the upper(U) and lower(D) faces of a cell.
37 C maskUp :: Land mask used to denote base of the domain.
38 C xA :: Tracer cell face area normal to X
39 C yA :: Tracer cell face area normal to X
40 C uTrans :: Zonal volume transport through cell face
41 C vTrans :: Meridional volume transport through cell face
42 C rTrans :: Vertical volume transport through cell face
43 C bi, bj, iMin, iMax, jMin, jMax :: Range of points for which calculation
44 C results will be set.
45 C myThid - Instance number for this innvocation of CALC_GT
46 _RL fVerT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
47 _RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
48 _RS yA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49 _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
50 _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51 _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52 _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53 _RL KappaRT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
54 INTEGER k,kUp,kDown,kM1
55 INTEGER bi,bj,iMin,iMax,jMin,jMax
56 _RL myTime
57 INTEGER myIter
58 INTEGER myThid
59
60 CEOP
61
62 #ifdef ALLOW_PASSIVE_TRACER
63 INTEGER i,j
64
65 #ifdef ALLOW_AUTODIFF_TAMC
66 C-- only the kUp part of fverT is set in this subroutine
67 C-- the kDown is still required
68 fVerT(1,1,kDown) = fVerT(1,1,kDown)
69 #endif
70
71 #ifdef INCLUDE_TR_FORCING_CODE
72 C-- External thermal forcing term(s)
73 CALL EXTERNAL_FORCING_TR(
74 I iMin,iMax,jMin,jMax,bi,bj,k,
75 I myTime,myThid)
76 #endif /* INCLUDE_TR_FORCING_CODE */
77
78 CALL GAD_CALC_RHS(
79 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
80 I xA,yA,uTrans,vTrans,rTrans,maskUp,
81 I diffKhT, diffK4T, KappaRT, tr1,
82 I GAD_TR1, tracerAdvScheme,
83 U fVerT, gTr1,
84 I myThid )
85
86 DO j=jMin,jMax
87 DO i=iMin,iMax
88 gTr1(i,j,1,bi,bj) = gTr1(i,j,1,bi,bj) +
89 & surfaceTendencyTr1(i,j,bi,bj)
90 ENDDO
91 ENDDO
92
93 IF ( tracerAdvScheme.EQ.ENUM_CENTERED_2ND
94 & .OR.tracerAdvScheme.EQ.ENUM_UPWIND_3RD
95 & .OR.tracerAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN
96 CALL ADAMS_BASHFORTH2(
97 I bi, bj, K,
98 U gTr1, gTr1nm1,
99 I myIter, myThid )
100 ENDIF
101
102 #ifdef NONLIN_FRSURF
103 IF (nonlinFreeSurf.GT.0) THEN
104 CALL FREESURF_RESCALE_G(
105 I bi, bj, K,
106 U gTr1,
107 I myThid )
108 ENDIF
109 #endif /* NONLIN_FRSURF */
110
111 #endif /* ALLOW_PASSIVE_TRACER */
112
113 RETURN
114 END

  ViewVC Help
Powered by ViewVC 1.1.22