155 |
ENDDO |
ENDDO |
156 |
|
|
157 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
158 |
CADJ STORE SlopeSqr(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE SlopeX(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
159 |
|
CADJ STORE SlopeY(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
160 |
|
CADJ STORE SlopeSqr(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
161 |
|
CADJ STORE taperFct(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
162 |
|
#ifdef GM_VISBECK_VARIABLE_K |
163 |
|
CADJ STORE dSigmaDrReal(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte |
164 |
|
#endif |
165 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
166 |
|
|
167 |
DO j=1-Oly+1,sNy+Oly-1 |
DO j=1-Oly+1,sNy+Oly-1 |
192 |
deltaH=deltaH/GM_Visbeck_depth |
deltaH=deltaH/GM_Visbeck_depth |
193 |
|
|
194 |
IF (K.eq.2) VisbeckK(i,j,bi,bj)=0. |
IF (K.eq.2) VisbeckK(i,j,bi,bj)=0. |
195 |
IF (Ssq(i,j).NE.0.) THEN |
IF ( Ssq(i,j).NE.0. .AND. dSigmaDrReal(i,j).NE.0. ) THEN |
196 |
N2= -Gravity*recip_RhoConst*dSigmaDrReal(i,j) |
N2= -Gravity*recip_RhoConst*dSigmaDrReal(i,j) |
197 |
SN=sqrt(Ssq(i,j)*N2) |
SN=sqrt(Ssq(i,j)*N2) |
198 |
VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH |
VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH |
238 |
|
|
239 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
240 |
kkey = (igmkey-1)*Nr + k |
kkey = (igmkey-1)*Nr + k |
241 |
#ifdef GM_NON_UNITY_DIAGONAL |
#ifdef GM_VISBECK_VARIABLE_K |
242 |
CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
243 |
CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
244 |
CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj, key=kkey, byte=isbyte |
261 |
& )*Kwz(i,j,k,bi,bj) |
& )*Kwz(i,j,k,bi,bj) |
262 |
ENDDO |
ENDDO |
263 |
ENDDO |
ENDDO |
264 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
265 |
|
#ifdef GM_VISBECK_VARIABLE_K |
266 |
|
CADJ STORE VisbeckK(:,:,bi,bj) = |
267 |
|
CADJ & comlev1_bibj, key=kkey, byte=isbyte |
268 |
|
#endif |
269 |
|
#endif |
270 |
|
|
271 |
#if ( defined (GM_NON_UNITY_DIAGONAL) || defined (GM_EXTRA_DIAGONAL) ) |
#if ( defined (GM_NON_UNITY_DIAGONAL) || defined (GM_EXTRA_DIAGONAL) ) |
272 |
|
|