/[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.6 by jmc, Thu Jun 9 15:59:20 2005 UTC revision 1.7 by jmc, Wed Jun 22 00:33:14 2005 UTC
# Line 71  C       Transport at interface k+1 (here Line 71  C       Transport at interface k+1 (here
71    
72  C delta_Z( U )  @ interface k  C delta_Z( U )  @ interface k
73  c       umask_Km1=mask_Km1*maskW(i,j,Km1,bi,bj)  c       umask_Km1=mask_Km1*maskW(i,j,Km1,bi,bj)
74          uZm=(mask_Km1*uFld(I,J,Km1,bi,bj)-uFld(I,J,K,bi,bj))          uZm=(uFld(I,J,K,bi,bj)-mask_Km1*uFld(I,J,Km1,bi,bj))*rkSign
 c    &      *rkFac  
75  c2   &      *recip_dRC(K)  c2   &      *recip_dRC(K)
76  c       IF (freeslip1) uZm=uZm*umask_Km1  c       IF (freeslip1) uZm=uZm*umask_Km1
77  c       IF (noslip1.AND.umask_Km1.EQ.0.) uZm=uZm*2.  c       IF (noslip1.AND.umask_Km1.EQ.0.) uZm=uZm*2.
78    
79  C delta_Z( U )  @ interface k+1  C delta_Z( U )  @ interface k+1
80  c       umask_Kp1=mask_Kp1*maskW(i,j,Kp1,bi,bj)  c       umask_Kp1=mask_Kp1*maskW(i,j,Kp1,bi,bj)
81          uZp=(uFld(I,J,K,bi,bj)-mask_Kp1*uFld(I,J,Kp1,bi,bj))          uZp=(mask_Kp1*uFld(I,J,Kp1,bi,bj)-uFld(I,J,K,bi,bj))*rkSign
 c    &      *rkFac  
82  c2   &      *recip_dRC(Kp1)  c2   &      *recip_dRC(Kp1)
83  c       IF (freeslipK) uZp=uZp*umask_Kp1  c       IF (freeslipK) uZp=uZp*umask_Kp1
84  c       IF (noslipK.AND.umask_Kp1.EQ.0.) uZp=uZp*2.  c       IF (noslipK.AND.umask_Kp1.EQ.0.) uZp=uZp*2.
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22