/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_v_del2v.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_v_del2v.F

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

revision 1.6 by mlosch, Wed Apr 5 15:43:15 2006 UTC revision 1.7 by jmc, Tue Dec 5 05:30:38 2006 UTC
# Line 53  C     Zonal flux d/dx V Line 53  C     Zonal flux d/dx V
53       &   *_recip_dxV(i,j,bi,bj)       &   *_recip_dxV(i,j,bi,bj)
54       &   *(vFld(i,j)-vFld(i-1,j))       &   *(vFld(i,j)-vFld(i-1,j))
55  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
56       &   *sqcosFacV(J,bi,bj)       &   *sqCosFacV(J,bi,bj)
57  #endif  #endif
58    c    &   *deepFacC(k)        ! dyU scaling factor
59    c    &   *recip_deepFacC(k)  ! recip_dxV scaling factor
60         ENDDO         ENDDO
61        ENDDO        ENDDO
62    
# Line 66  C     Meridional flux d/dy V Line 68  C     Meridional flux d/dy V
68       &   *_recip_dyF(i,j,bi,bj)       &   *_recip_dyF(i,j,bi,bj)
69       &   *(vFld(i,j+1)-vFld(i,j))       &   *(vFld(i,j+1)-vFld(i,j))
70  #if (defined (ISOTROPIC_COS_SCALING) && defined (COSINEMETH_III))  #if (defined (ISOTROPIC_COS_SCALING) && defined (COSINEMETH_III))
71       &   *sqcosFacU(J,bi,bj)       &   *sqCosFacU(J,bi,bj)
72  #endif  #endif
73    c    &   *deepFacC(k)        ! dxF scaling factor
74    c    &   *recip_deepFacC(k)  ! recip_dyF scaling factor
75         ENDDO         ENDDO
76        ENDDO        ENDDO
77    
78  C     del^2 V  C     del^2 V
79        DO j=1-Oly+1,sNy+Oly-1        DO j=1-Oly+1,sNy+Oly-1
80         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
81          del2v(i,j) = recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)          del2v(i,j) =
82       &  *recip_rAs(i,j,bi,bj)       &   recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)
83         &  *recip_rAs(i,j,bi,bj)*recip_deepFac2C(k)
84       &  *( fZon(i+1,j)    - fZon(i, j )       &  *( fZon(i+1,j)    - fZon(i, j )
85       &    +fMer( i ,j)    - fMer(i,j-1)       &    +fMer( i ,j)    - fMer(i,j-1)
86       &   )*_maskS(i,j,k,bi,bj)       &   )*_maskS(i,j,k,bi,bj)
# Line 90  C-- No-slip BCs impose a drag at walls.. Line 95  C-- No-slip BCs impose a drag at walls..
95          hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)          hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)
96          del2v(i,j) = del2v(i,j)          del2v(i,j) = del2v(i,j)
97       &  -_recip_hFacS(i,j,k,bi,bj)*recip_drF(k)       &  -_recip_hFacS(i,j,k,bi,bj)*recip_drF(k)
98       &  *recip_rAs(i,j,bi,bj)       &   *recip_rAs(i,j,bi,bj)*recip_deepFac2C(k)
99       &   *( hFacZClosedW*dyU( i ,j,bi,bj)       &   *( hFacZClosedW*dyU( i ,j,bi,bj)
100       &      *_recip_dxV( i ,j,bi,bj)       &      *_recip_dxV( i ,j,bi,bj)
101       &     +hFacZClosedE*dyU(i+1,j,bi,bj)       &     +hFacZClosedE*dyU(i+1,j,bi,bj)
102       &      *_recip_dxV(i+1,j,bi,bj)       &      *_recip_dxV(i+1,j,bi,bj)
103       &    )*drF(k)*2.*vFld(i,j)       &    )*drF(k)*2.*vFld(i,j)
104       &      *_maskS(i,j,k,bi,bj)       &     *_maskS(i,j,k,bi,bj)
105         ENDDO         ENDDO
106        ENDDO        ENDDO
107        ENDIF        ENDIF
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22