/[MITgcm]/MITgcm/pkg/shap_filt/shap_filt_computvort.F
ViewVC logotype

Diff of /MITgcm/pkg/shap_filt/shap_filt_computvort.F

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

revision 1.1 by jmc, Fri Oct 7 00:21:07 2005 UTC revision 1.2 by jmc, Wed Nov 28 14:58:07 2007 UTC
# Line 58  C     == Local variables == Line 58  C     == Local variables ==
58  #endif /* ALLOW_EXCH2 */  #endif /* ALLOW_EXCH2 */
59  CEOP  CEOP
60    
61    #ifdef ALLOW_AUTODIFF_TAMC
62  C-    Initialisation :  C-    Initialisation :
63        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
64          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
65            vort(i,j)= 0.            vort(i,j)= 0.
66          ENDDO          ENDDO
67        ENDDO        ENDDO
68    #endif
69    
70  C-    replace Physical calc Div & Vort by computational one :  C-    replace Physical calc Div & Vort by computational one :
71        DO j=3-Oly,sNy+Oly-1        DO j=2-Oly,sNy+Oly
72          DO i=3-Olx,sNx+Olx-1          DO i=2-Olx,sNx+Olx
73            vort(i,j) = ( vFld(i,j)-vFld(i-1,j) )            vort(i,j) = ( vFld(i,j)-vFld(i-1,j) )
74       &              - ( uFld(i,j)-uFld(i,j-1) )       &              - ( uFld(i,j)-uFld(i,j-1) )
75            maskZ = (maskW(i,j,k,bi,bj)+maskW(i,j-1,k,bi,bj))            maskZ = (maskW(i,j,k,bi,bj)+maskW(i,j-1,k,bi,bj))
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22