/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vi_u_vertshear.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vi_u_vertshear.F

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

revision 1.3 by jmc, Thu Sep 6 17:59:36 2001 UTC revision 1.4 by jmc, Mon Sep 10 14:14:36 2001 UTC
# Line 59  c       wBarXp=0.5*(wFld(I,J,Kp1,bi,bj)+ Line 59  c       wBarXp=0.5*(wFld(I,J,Kp1,bi,bj)+
59  c    &         *mask_Kp1  c    &         *mask_Kp1
60    
61  C       Transport at interface k  C       Transport at interface k
62          wBarXm=0.5*(wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)          wBarXm=0.5*(
63       &             +wFld(I-1,J,K,bi,bj)*rA(i-1,j,bi,bj))       &    wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(I,J,Km1,bi,bj)
64         &   +wFld(I-1,J,K,bi,bj)*rA(i-1,j,bi,bj)*maskC(I-1,J,Km1,bi,bj)
65         &             )*mask_Km1
66    
67  C       Transport at interface k+1  C       Transport at interface k+1 (here wFld is already masked)
68          wBarXp=0.5*(wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj)          wBarXp=0.5*(
69       &             +wFld(I-1,J,Kp1,bi,bj)*rA(i-1,j,bi,bj))*mask_Kp1       &    wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj)
70         &   +wFld(I-1,J,Kp1,bi,bj)*rA(i-1,j,bi,bj)
71         &             )*mask_Kp1
72    
73  C delta_Z( U )  @ interface k  C delta_Z( U )  @ interface k
74          umask_Km1=mask_Km1*maskW(i,j,Km1,bi,bj)  c       umask_Km1=mask_Km1*maskW(i,j,Km1,bi,bj)
75          uZm=(umask_Km1*uFld(I,J,Km1,bi,bj)-uFld(I,J,K,bi,bj))          uZm=(mask_Km1*uFld(I,J,Km1,bi,bj)-uFld(I,J,K,bi,bj))
76  c2   &      *recip_dRC(K)  c2   &      *recip_dRC(K)
77          IF (freeslip1) uZm=uZm*umask_Km1  c       IF (freeslip1) uZm=uZm*umask_Km1
78          IF (noslip1.AND.umask_Km1.EQ.0.) uZm=uZm*2.  c       IF (noslip1.AND.umask_Km1.EQ.0.) uZm=uZm*2.
79    
80  C delta_Z( U )  @ interface k+1  C delta_Z( U )  @ interface k+1
81          umask_Kp1=mask_Kp1*maskW(i,j,Kp1,bi,bj)  c       umask_Kp1=mask_Kp1*maskW(i,j,Kp1,bi,bj)
82          uZp=(uFld(I,J,K,bi,bj)-umask_Kp1*uFld(I,J,Kp1,bi,bj))          uZp=(uFld(I,J,K,bi,bj)-mask_Kp1*uFld(I,J,Kp1,bi,bj))
83  c2   &      *recip_dRC(Kp1)  c2   &      *recip_dRC(Kp1)
84          IF (freeslipK) uZp=uZp*umask_Kp1  c       IF (freeslipK) uZp=uZp*umask_Kp1
85          IF (noslipK.AND.umask_Kp1.EQ.0.) uZp=uZp*2.  c       IF (noslipK.AND.umask_Kp1.EQ.0.) uZp=uZp*2.
86    
87  c1      IF (upwindShear) THEN  c1      IF (upwindShear) THEN
88  c1       wBarXZ=0.5*( wBarXm + wBarXp )  c1       wBarXZ=0.5*( wBarXm + wBarXp )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22