118 |
if (K.eq.2) VisbeckK(i,j,bi,bj)=0. |
if (K.eq.2) VisbeckK(i,j,bi,bj)=0. |
119 |
Calt? if (dSigmaDrReal(i,j).NE.0.) then |
Calt? if (dSigmaDrReal(i,j).NE.0.) then |
120 |
Calt? N2=(-Gravity*recip_Rhonil)*dSigmaDrReal(i,j) |
Calt? N2=(-Gravity*recip_Rhonil)*dSigmaDrReal(i,j) |
121 |
if (dRdSigmaLtd(i,j).NE.0.) then |
if ( dRdSigmaLtd(i,j).NE.0. .AND. Ssq.NE.0. ) then |
122 |
N2=(-Gravity*recip_Rhonil)/dRdSigmaLtd(i,j) |
N2=(-Gravity*recip_Rhonil)/dRdSigmaLtd(i,j) |
123 |
SN=sqrt(Ssq*N2) |
SN=sqrt(Ssq*N2) |
124 |
VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH |
VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH |