/[MITgcm]/MITgcm/pkg/gmredi/gmredi_calc_psi_b.F
ViewVC logotype

Contents of /MITgcm/pkg/gmredi/gmredi_calc_psi_b.F

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

Revision 1.10 - (show annotations) (download)
Fri May 30 02:50:16 2008 UTC (16 years ago) by gforget
Branch: MAIN
CVS Tags: checkpoint62c, checkpoint62q, checkpoint62p, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint60, checkpoint61, checkpoint62, checkpoint62b, checkpoint61f, checkpoint61n, checkpoint61q, checkpoint61e, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.9: +8 -1 lines 
o bridging the gap between eddy stress and GM.
  -> eddyTau is replaced with eddyPsi (eddyTau = f x rho0 x eddyPsi)
      along with a change in CPP option (now ALLOW_EDDYPSI).
  -> when using GM w/ GM_AdvForm:
      The total eddy streamfunction (Psi = eddyPsi + K x Slope)
      is applied either in the tracer Eq. or in momentum Eq.
      depending on data.gmredi (intro. GM_InMomAsStress).
  -> ALLOW_EDDYPSI_CONTROL for estimation purpose.
  The key modifications are in model/src/taueddy_external_forcing.F
  pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F
1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_calc_psi_b.F,v 1.9 2008/03/28 18:48:05 heimbach Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 CStartOfInterface
7 SUBROUTINE GMREDI_CALC_PSI_B(
8 I bi, bj, iMin, iMax, jMin, jMax,
9 I sigmaX, sigmaY, sigmaR,
10 I ldd97_LrhoW, ldd97_LrhoS,
11 I myThid )
12 C /==========================================================\
13 C | SUBROUTINE GMREDI_CALC_PSI_B |
14 C | o Calculate stream-functions for GM bolus velocity |
15 C |==========================================================|
16 C \==========================================================/
17 IMPLICIT NONE
18
19 C == Global variables ==
20 #include "SIZE.h"
21 #include "GRID.h"
22 #include "DYNVARS.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "GMREDI.h"
26 #include "FFIELDS.h"
27
28 #ifdef ALLOW_AUTODIFF_TAMC
29 #include "tamc.h"
30 #include "tamc_keys.h"
31 #endif /* ALLOW_AUTODIFF_TAMC */
32
33 C == Routine arguments ==
34 C
35 _RL sigmaX(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
36 _RL sigmaY(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
37 _RL sigmaR(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
38 _RL ldd97_LrhoW(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
39 _RL ldd97_LrhoS(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
40 INTEGER bi,bj,iMin,iMax,jMin,jMax
41 INTEGER myThid
42 CEndOfInterface
43
44 #ifdef ALLOW_GMREDI
45 #ifdef GM_BOLUS_ADVEC
46
47 C == Local variables ==
48 INTEGER i,j,k, km1
49 _RL SlopeX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
50 _RL SlopeY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
51 _RL dSigmaDrW(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
52 _RL dSigmaDrS(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
53 _RL taperX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
54 _RL taperY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
55
56 C- Initialization : <= done in S/R gmredi_init
57
58 #ifdef ALLOW_AUTODIFF_TAMC
59 act1 = bi - myBxLo(myThid)
60 max1 = myBxHi(myThid) - myBxLo(myThid) + 1
61 act2 = bj - myByLo(myThid)
62 max2 = myByHi(myThid) - myByLo(myThid) + 1
63 act3 = myThid - 1
64 max3 = nTx*nTy
65 act4 = ikey_dynamics - 1
66 igmkey = (act1 + 1) + act2*max1
67 & + act3*max1*max2
68 & + act4*max1*max2*max3
69 #endif /* ALLOW_AUTODIFF_TAMC */
70
71 #ifdef ALLOW_AUTODIFF_TAMC
72 # ifdef GM_VISBECK_VARIABLE_K
73 CADJ STORE VisbeckK(:,:,bi,bj) = comlev1_bibj, key=igmkey, byte=isbyte
74 # endif
75 #endif
76 IF (GM_AdvForm) THEN
77 DO k=2,Nr
78 km1 = k-1
79
80 #ifdef ALLOW_AUTODIFF_TAMC
81 kkey = (igmkey-1)*Nr + k
82 DO j=1-Oly,sNy+Oly
83 DO i=1-Olx,sNx+Olx
84 SlopeX(i,j) = 0. _d 0
85 SlopeY(i,j) = 0. _d 0
86 dSigmaDrW(i,j) = 0. _d 0
87 dSigmaDrS(i,j) = 0. _d 0
88 ENDDO
89 ENDDO
90 #endif
91
92 C Gradient of Sigma below U and V points
93 DO j=1-Oly,sNy+Oly
94 DO i=1-Olx+1,sNx+Olx
95 SlopeX(i,j)=op5*( sigmaX(i,j,km1)+sigmaX(i,j,k) )
96 & *maskW(i,j,k,bi,bj)
97 dSigmaDrW(i,j)=op5*( sigmaR(i-1,j,k)+sigmaR(i,j,k) )
98 & *maskW(i,j,k,bi,bj)
99 ENDDO
100 ENDDO
101 DO j=1-Oly+1,sNy+Oly
102 DO i=1-Olx,sNx+Olx
103 SlopeY(i,j)=op5*( sigmaY(i,j,km1)+sigmaY(i,j,k) )
104 & *maskS(i,j,k,bi,bj)
105 dSigmaDrS(i,j)=op5*( sigmaR(i,j-1,k)+sigmaR(i,j,k) )
106 & *maskS(i,j,k,bi,bj)
107 ENDDO
108 ENDDO
109
110 C Calculate slopes , taper and/or clip
111 CALL GMREDI_SLOPE_PSI(
112 O taperX, taperY,
113 U SlopeX, SlopeY,
114 U dSigmaDrW, dSigmaDrS,
115 I ldd97_LrhoW, ldd97_LrhoS, rF(k), k,
116 I bi, bj, myThid )
117
118 #ifdef ALLOW_AUTODIFF_TAMC
119 CADJ STORE SlopeX(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
120 CADJ STORE SlopeY(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
121 CADJ STORE taperX(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
122 CADJ STORE taperY(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
123 #endif /* ALLOW_AUTODIFF_TAMC */
124
125 C- Compute the 2 stream-function Components ( GM bolus vel.)
126 DO j=1-Oly,sNy+Oly
127 DO i=1-Olx+1,sNx+Olx
128 GM_PsiX(i,j,k,bi,bj) = SlopeX(i,j)*taperX(i,j)
129 #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
130 & *( kapgm(i,j,k,bi,bj)
131 #else
132 & *( GM_background_K
133 #endif
134 #ifdef GM_VISBECK_VARIABLE_K
135 & +op5*(VisbeckK(i-1,j,bi,bj)+VisbeckK(i,j,bi,bj))
136 #endif
137 & )*maskW(i,j,k,bi,bj)
138 #ifdef ALLOW_EDDYPSI
139 & +eddyPsiX(i,j,k,bi,bj)*maskW(i,j,k,bi,bj)
140 #endif
141 ENDDO
142 ENDDO
143 DO j=1-Oly+1,sNy+Oly
144 DO i=1-Olx,sNx+Olx
145 GM_PsiY(i,j,k,bi,bj) = SlopeY(i,j)*taperY(i,j)
146 #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
147 & *( kapgm(i,j,k,bi,bj)
148 #else
149 & *( GM_background_K
150 #endif
151 #ifdef GM_VISBECK_VARIABLE_K
152 & +op5*(VisbeckK(i,j-1,bi,bj)+VisbeckK(i,j,bi,bj))
153 #endif
154 & )*maskS(i,j,k,bi,bj)
155 #ifdef ALLOW_EDDYPSI
156 & +eddyPsiY(i,j,k,bi,bj)*maskS(i,j,k,bi,bj)
157 #endif
158 ENDDO
159 ENDDO
160
161 C----- end of loop on level k
162 ENDDO
163
164 ENDIF
165 #endif /* GM_BOLUS_ADVEC */
166 #endif /* ALLOW_GMREDI */
167
168 RETURN
169 END
  ViewVC Help
Powered by ViewVC 1.1.22