42 |
integer ikppkey |
integer ikppkey |
43 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
44 |
|
|
|
EXTERNAL DIFFERENT_MULTIPLE |
|
|
LOGICAL DIFFERENT_MULTIPLE |
|
|
|
|
45 |
C !INPUT PARAMETERS: =================================================== |
C !INPUT PARAMETERS: =================================================== |
46 |
c Routine arguments |
c Routine arguments |
47 |
c bi, bj - array indices on which to apply calculations |
c bi, bj - array indices on which to apply calculations |
80 |
C Initialize local fields |
C Initialize local fields |
81 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
82 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
83 |
GH(I,J) = 0. |
GH(I,J) = 0. _d 0 |
84 |
GM(I,J) = 0. |
GM(I,J) = 0. _d 0 |
85 |
ENDDO |
ENDDO |
86 |
ENDDO |
ENDDO |
87 |
DO K = 1, Nr |
DO K = 1, Nr |
88 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
89 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
90 |
SH(I,J,K) = 0. |
SH(I,J,K) = 0. _d 0 |
91 |
SM(I,J,K) = 0. |
SM(I,J,K) = 0. _d 0 |
92 |
tke(I,J,K) = 0. |
tke(I,J,K) = 0. _d 0 |
93 |
ENDDO |
ENDDO |
94 |
ENDDO |
ENDDO |
95 |
ENDDO |
ENDDO |
129 |
C the meaning of shear and negative buoyancy frequency |
C the meaning of shear and negative buoyancy frequency |
130 |
DO J=jMin,jMax |
DO J=jMin,jMax |
131 |
DO I=iMin,iMax |
DO I=iMin,iMax |
132 |
GH(I,J) = 0. |
GH(I,J) = 0. _d 0 |
133 |
GM(I,J) = 0. |
GM(I,J) = 0. _d 0 |
134 |
ENDDO |
ENDDO |
135 |
ENDDO |
ENDDO |
136 |
C Find boundary length scale from energy weighted mean. |
C Find boundary length scale from energy weighted mean. |
149 |
C compute boundary length scale MYhbl |
C compute boundary length scale MYhbl |
150 |
DO J=jMin,jMax |
DO J=jMin,jMax |
151 |
DO I=iMin,iMax |
DO I=iMin,iMax |
152 |
IF ( GH(I,J) .EQ. 0. ) THEN |
IF ( GH(I,J) .EQ. 0. _d 0 ) THEN |
153 |
MYhbl(I,J,bi,bj) = 0. |
MYhbl(I,J,bi,bj) = 0. _d 0 |
154 |
ELSE |
ELSE |
155 |
MYhbl(I,J,bi,bj) = -GM(I,J)/GH(I,J)*MYhblScale |
MYhbl(I,J,bi,bj) = -GM(I,J)/GH(I,J)*MYhblScale |
156 |
ENDIF |
ENDIF |