156 |
END IF |
END IF |
157 |
2110 CONTINUE |
2110 CONTINUE |
158 |
|
|
|
|
|
|
200 CONTINUE |
|
|
|
|
|
!*********************************************************** |
|
|
|
|
159 |
!* 3. COMPUTE MOUNTAIN LIFT |
!* 3. COMPUTE MOUNTAIN LIFT |
160 |
|
|
|
300 CONTINUE |
|
|
|
|
161 |
DO 301 jl = 1, klon |
DO 301 jl = 1, klon |
162 |
IF (ktest(jl)==1) THEN |
IF (ktest(jl)==1) THEN |
163 |
ztau(jl,klev+1) = -gklift*zrho(jl,klev+1)*2.*romega*2*pvaror(jl)*sin & |
ztau(jl,klev+1) = -gklift*zrho(jl,klev+1)*2.*romega*2*pvaror(jl)*sin & |
175 |
!* -------------------- |
!* -------------------- |
176 |
|
|
177 |
|
|
178 |
400 CONTINUE |
DO jk = 1, klev |
179 |
|
DO jl = 1, klon |
|
DO 401 jk = 1, klev |
|
|
DO 401 jl = 1, klon |
|
180 |
IF (ktest(jl)==1) THEN |
IF (ktest(jl)==1) THEN |
181 |
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) |
182 |
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) |
184 |
ztau(jl,jk) = 0.0 |
ztau(jl,jk) = 0.0 |
185 |
ztav(jl,jk) = 0.0 |
ztav(jl,jk) = 0.0 |
186 |
END IF |
END IF |
187 |
401 CONTINUE |
end DO |
188 |
|
end DO |
189 |
|
|
190 |
|
|
191 |
!* 5. COMPUTE TENDENCIES. |
!* 5. COMPUTE TENDENCIES. |
192 |
!* ------------------- |
!* ------------------- |
193 |
IF (lifthigh) THEN |
IF (lifthigh) THEN |
194 |
|
|
|
500 CONTINUE |
|
195 |
! PRINT *,' DANS OROLIFT: 500' |
! PRINT *,' DANS OROLIFT: 500' |
196 |
|
|
197 |
! EXPLICIT SOLUTION AT ALL LEVELS |
! EXPLICIT SOLUTION AT ALL LEVELS |
208 |
|
|
209 |
! PROJECT PERPENDICULARLY TO U NOT TO DESTROY ENERGY |
! PROJECT PERPENDICULARLY TO U NOT TO DESTROY ENERGY |
210 |
|
|
211 |
DO 530 jk = 1, klev |
DO jk = 1, klev |
212 |
DO 530 jl = 1, klon |
DO jl = 1, klon |
213 |
IF (ktest(jl)==1) THEN |
IF (ktest(jl)==1) THEN |
214 |
|
|
215 |
zslow = sqrt(pulow(jl)**2+pvlow(jl)**2) |
zslow = sqrt(pulow(jl)**2+pvlow(jl)**2) |
229 |
END IF |
END IF |
230 |
|
|
231 |
END IF |
END IF |
232 |
530 CONTINUE |
end DO |
233 |
|
end DO |
234 |
|
|
235 |
! 6. LOW LEVEL LIFT, SEMI IMPLICIT: |
! 6. LOW LEVEL LIFT, SEMI IMPLICIT: |
236 |
! ---------------------------------- |
! ---------------------------------- |