184 |
|
|
185 |
400 CONTINUE |
400 CONTINUE |
186 |
|
|
187 |
DO 401 jk = 1, klev |
DO jk = 1, klev |
188 |
DO 401 jl = 1, klon |
DO jl = 1, klon |
189 |
IF (ktest(jl)==1) THEN |
IF (ktest(jl)==1) THEN |
190 |
ztau(jl,jk) = ztau(jl,klev+1)*paphm1(jl,jk)/paphm1(jl,klev+1) |
ztau(jl,jk) = ztau(jl,klev+1)*paphm1(jl,jk)/paphm1(jl,klev+1) |
191 |
ztav(jl,jk) = ztav(jl,klev+1)*paphm1(jl,jk)/paphm1(jl,klev+1) |
ztav(jl,jk) = ztav(jl,klev+1)*paphm1(jl,jk)/paphm1(jl,klev+1) |
193 |
ztau(jl,jk) = 0.0 |
ztau(jl,jk) = 0.0 |
194 |
ztav(jl,jk) = 0.0 |
ztav(jl,jk) = 0.0 |
195 |
END IF |
END IF |
196 |
401 CONTINUE |
end DO |
197 |
|
end DO |
198 |
|
|
199 |
|
|
200 |
!* 5. COMPUTE TENDENCIES. |
!* 5. COMPUTE TENDENCIES. |
218 |
|
|
219 |
! PROJECT PERPENDICULARLY TO U NOT TO DESTROY ENERGY |
! PROJECT PERPENDICULARLY TO U NOT TO DESTROY ENERGY |
220 |
|
|
221 |
DO 530 jk = 1, klev |
DO jk = 1, klev |
222 |
DO 530 jl = 1, klon |
DO jl = 1, klon |
223 |
IF (ktest(jl)==1) THEN |
IF (ktest(jl)==1) THEN |
224 |
|
|
225 |
zslow = sqrt(pulow(jl)**2+pvlow(jl)**2) |
zslow = sqrt(pulow(jl)**2+pvlow(jl)**2) |
239 |
END IF |
END IF |
240 |
|
|
241 |
END IF |
END IF |
242 |
530 CONTINUE |
end DO |
243 |
|
end DO |
244 |
|
|
245 |
! 6. LOW LEVEL LIFT, SEMI IMPLICIT: |
! 6. LOW LEVEL LIFT, SEMI IMPLICIT: |
246 |
! ---------------------------------- |
! ---------------------------------- |