/[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.3 by jmc, Tue May 12 19:56:36 2009 UTC
# Line 29  C     == Global variables === Line 29  C     == Global variables ===
29  #include "PARAMS.h"  #include "PARAMS.h"
30  #include "GRID.h"  #include "GRID.h"
31  #ifdef ALLOW_EXCH2  #ifdef ALLOW_EXCH2
32    #include "W2_EXCH2_SIZE.h"
33  #include "W2_EXCH2_TOPOLOGY.h"  #include "W2_EXCH2_TOPOLOGY.h"
 #include "W2_EXCH2_PARAMS.h"  
34  #endif /* ALLOW_EXCH2 */  #endif /* ALLOW_EXCH2 */
35    
36  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# 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.3
  ViewVC Help
Powered by ViewVC 1.1.22