71 |
|
|
72 |
C delta_Z( V ) @ interface k |
C delta_Z( V ) @ interface k |
73 |
c vmask_Km1=mask_Km1*maskS(i,j,Km1,bi,bj) |
c vmask_Km1=mask_Km1*maskS(i,j,Km1,bi,bj) |
74 |
vZm=(mask_Km1*vFld(I,J,Km1,bi,bj)-vFld(I,J,K,bi,bj)) |
vZm=(vFld(I,J,K,bi,bj)-mask_Km1*vFld(I,J,Km1,bi,bj))*rkSign |
|
c & *rkFac |
|
75 |
c2 & *recip_dRC(K) |
c2 & *recip_dRC(K) |
76 |
c IF (freeslip1) vZm=vZm*vmask_Km1 |
c IF (freeslip1) vZm=vZm*vmask_Km1 |
77 |
c IF (noslip1.AND.vmask_Km1.EQ.0.) vZm=vZm*2. |
c IF (noslip1.AND.vmask_Km1.EQ.0.) vZm=vZm*2. |
78 |
|
|
79 |
C delta_Z( V ) @ interface k+1 |
C delta_Z( V ) @ interface k+1 |
80 |
c vmask_Kp1=mask_Kp1*maskS(i,j,Kp1,bi,bj) |
c vmask_Kp1=mask_Kp1*maskS(i,j,Kp1,bi,bj) |
81 |
vZp=(vFld(I,J,K,bi,bj)-mask_Kp1*vFld(I,J,Kp1,bi,bj)) |
vZp=(mask_Kp1*vFld(I,J,Kp1,bi,bj)-vFld(I,J,K,bi,bj))*rkSign |
|
c & *rkFac |
|
82 |
c2 & *recip_dRC(Kp1) |
c2 & *recip_dRC(Kp1) |
83 |
c IF (freeslipK) vZp=vZp*vmask_Kp1 |
c IF (freeslipK) vZp=vZp*vmask_Kp1 |
84 |
c IF (noslipK.AND.vmask_Kp1.EQ.0.) vZp=vZp*2. |
c IF (noslipK.AND.vmask_Kp1.EQ.0.) vZp=vZp*2. |