/[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.15 by jmc, Sun Jan 12 21:35:27 2003 UTC revision 1.16 by heimbach, Tue Jan 21 19:34:13 2003 UTC
# Line 152  C       Mask Iso-neutral slopes Line 152  C       Mask Iso-neutral slopes
152        ENDDO        ENDDO
153    
154  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
155    CADJ STORE SlopeX(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
156    CADJ STORE SlopeY(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
157  CADJ STORE SlopeSqr(:,:)     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE SlopeSqr(:,:)     = comlev1_bibj_k, key=kkey, byte=isbyte
158    CADJ STORE dsigmadrreal(:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
159    CADJ STORE taperfct(:,:)     = comlev1_bibj_k, key=kkey, byte=isbyte
160  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
161    
162        DO j=1-Oly+1,sNy+Oly-1        DO j=1-Oly+1,sNy+Oly-1
# Line 216  C-     Limit range that KapGM can take Line 220  C-     Limit range that KapGM can take
220          ENDDO          ENDDO
221         ENDDO         ENDDO
222        ENDIF        ENDIF
223    cph( NEW
224    #ifdef ALLOW_AUTODIFF_TAMC
225    CADJ STORE VisbeckK(:,:,bi,bj) = comlev1_bibj, key=igmkey, byte=isbyte
226    #endif
227    cph)
228  #endif /* GM_VISBECK_VARIABLE_K */  #endif /* GM_VISBECK_VARIABLE_K */
229    
230    
# Line 229  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  #if (defined (GM_NON_UNITY_DIAGONAL) || \
242         defined (GM_VISBECK_VARIABLE_K))
243  CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE Kwx(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
244  CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE Kwy(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
245  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
# Line 241  C-    express the Tensor in term of Diff Line 251  C-    express the Tensor in term of Diff
251         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
252          Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K          Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K
253  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
254       &          + VisbeckK(i,j,bi,bj)*(1.+GM_skewflx)           &          + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx)    
255  #endif  #endif
256          Kwx(i,j,k,bi,bj)= Kgm_tmp*Kwx(i,j,k,bi,bj)          Kwx(i,j,k,bi,bj)= Kgm_tmp*Kwx(i,j,k,bi,bj)
257          Kwy(i,j,k,bi,bj)= Kgm_tmp*Kwy(i,j,k,bi,bj)          Kwy(i,j,k,bi,bj)= Kgm_tmp*Kwy(i,j,k,bi,bj)
# Line 284  C     Calculate slopes for use in tensor Line 294  C     Calculate slopes for use in tensor
294       O             SlopeSqr, taperFct,       O             SlopeSqr, taperFct,
295       I             bi, bj, myThid )       I             bi, bj, myThid )
296    
297    cph( NEW
298    #ifdef ALLOW_AUTODIFF_TAMC
299    cph(
300    CADJ STORE taperfct(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
301    cph)
302    #endif /* ALLOW_AUTODIFF_TAMC */
303    cph)
304    
305  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
306          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
307           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
# Line 297  C     Calculate slopes for use in tensor Line 315  C     Calculate slopes for use in tensor
315           ENDDO           ENDDO
316          ENDDO          ENDDO
317  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
318  # ifndef GM_TAPER_ORIG_CLIPPING  # ifdef GM_EXCLUDE_CLIPPING
319  CADJ STORE Kux(:,:,k,bi,bj)  = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE Kux(:,:,k,bi,bj)  = comlev1_bibj_k, key=kkey, byte=isbyte
320  # endif  # endif
321  #endif  #endif
# Line 357  C     Calculate slopes for use in tensor Line 375  C     Calculate slopes for use in tensor
375       O             SlopeSqr, taperFct,       O             SlopeSqr, taperFct,
376       I             bi, bj, myThid )       I             bi, bj, myThid )
377    
378    cph(
379    #ifdef ALLOW_AUTODIFF_TAMC
380    cph(
381    CADJ STORE taperfct(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
382    cph)
383    #endif /* ALLOW_AUTODIFF_TAMC */
384    cph)
385    
386  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
387          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
388           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
# Line 370  C     Calculate slopes for use in tensor Line 396  C     Calculate slopes for use in tensor
396           ENDDO           ENDDO
397          ENDDO          ENDDO
398  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
399  # ifndef GM_TAPER_ORIG_CLIPPING  # ifdef GM_EXCLUDE_CLIPPING
400  CADJ STORE Kvy(:,:,k,bi,bj)  = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE Kvy(:,:,k,bi,bj)  = comlev1_bibj_k, key=kkey, byte=isbyte
401  # endif  # endif
402  #endif  #endif
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22