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

Diff of /MITgcm/pkg/gmredi/gmredi_calc_tensor.F

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

revision 1.12 by heimbach, Thu Nov 14 22:43:49 2002 UTC revision 1.13 by heimbach, Thu Nov 28 17:30:34 2002 UTC
# Line 155  C       Mask Iso-neutral slopes Line 155  C       Mask Iso-neutral slopes
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
# Line 186  C       Now we convert deltaH to a non-d Line 192  C       Now we convert deltaH to a non-d
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
# Line 232  C-- 2nd loop on k : compute Tensor Coeff Line 238  C-- 2nd loop on k : compute Tensor Coeff
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
# Line 255  C-    express the Tensor in term of Diff Line 261  C-    express the Tensor in term of Diff
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    
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22