85 |
_RL vDiss(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL vDiss(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
86 |
C I,J,K - Loop counters |
C I,J,K - Loop counters |
87 |
INTEGER i,j,k |
INTEGER i,j,k |
|
C rVelMaskOverride - Factor for imposing special surface boundary conditions |
|
|
C ( set according to free-surface condition ). |
|
|
C hFacROpen - Lopped cell factos used tohold fraction of open |
|
|
C hFacRClosed and closed cell wall. |
|
|
_RL rVelMaskOverride |
|
88 |
C xxxFac - On-off tracer parameters used for switching terms off. |
C xxxFac - On-off tracer parameters used for switching terms off. |
89 |
_RL ArDudrFac |
_RL ArDudrFac |
90 |
_RL phxFac |
_RL phxFac |
92 |
_RL ArDvdrFac |
_RL ArDvdrFac |
93 |
_RL phyFac |
_RL phyFac |
94 |
c _RL mtFacV |
c _RL mtFacV |
|
_RL wVelBottomOverride |
|
95 |
LOGICAL bottomDragTerms |
LOGICAL bottomDragTerms |
96 |
LOGICAL writeDiag |
LOGICAL writeDiag |
97 |
_RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
112 |
fVerV(1,1,kUp) = fVerV(1,1,kUp) |
fVerV(1,1,kUp) = fVerV(1,1,kUp) |
113 |
#endif |
#endif |
114 |
|
|
|
rVelMaskOverride=1. |
|
|
IF ( k .EQ. 1 ) rVelMaskOverride=freeSurfFac |
|
|
wVelBottomOverride=1. |
|
|
IF (k.EQ.Nr) wVelBottomOverride=0. |
|
115 |
writeDiag = DIFFERENT_MULTIPLE(diagFreq, myTime, |
writeDiag = DIFFERENT_MULTIPLE(diagFreq, myTime, |
116 |
& myTime-deltaTClock) |
& myTime-deltaTClock) |
117 |
|
|
204 |
|
|
205 |
IF (momViscosity) THEN |
IF (momViscosity) THEN |
206 |
C Calculate del^2 u and del^2 v for bi-harmonic term |
C Calculate del^2 u and del^2 v for bi-harmonic term |
207 |
IF (viscA4.NE.0. |
IF ( (viscA4.NE.0. .AND. no_slip_sides) |
208 |
|
& .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0. |
209 |
& .OR. viscA4Grid.NE.0. |
& .OR. viscA4Grid.NE.0. |
210 |
& .OR. viscC4leith.NE.0. |
& .OR. viscC4leith.NE.0. |
211 |
& ) THEN |
& ) THEN |