--- trunk/phylmd/Orography/orolift.f 2014/03/05 14:57:53 82 +++ trunk/Sources/phylmd/Orography/orolift.f 2015/07/21 14:44:45 158 @@ -156,15 +156,8 @@ END IF 2110 CONTINUE - -200 CONTINUE - -!*********************************************************** - !* 3. COMPUTE MOUNTAIN LIFT -300 CONTINUE - DO 301 jl = 1, klon IF (ktest(jl)==1) THEN ztau(jl,klev+1) = -gklift*zrho(jl,klev+1)*2.*romega*2*pvaror(jl)*sin & @@ -182,10 +175,8 @@ !* -------------------- -400 CONTINUE - - DO 401 jk = 1, klev - DO 401 jl = 1, klon + DO jk = 1, klev + DO jl = 1, klon IF (ktest(jl)==1) THEN ztau(jl,jk) = ztau(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,14 +184,14 @@ ztau(jl,jk) = 0.0 ztav(jl,jk) = 0.0 END IF -401 CONTINUE + end DO + end DO !* 5. COMPUTE TENDENCIES. !* ------------------- IF (lifthigh) THEN -500 CONTINUE ! PRINT *,' DANS OROLIFT: 500' ! EXPLICIT SOLUTION AT ALL LEVELS @@ -217,8 +208,8 @@ ! PROJECT PERPENDICULARLY TO U NOT TO DESTROY ENERGY - DO 530 jk = 1, klev - DO 530 jl = 1, klon + DO jk = 1, klev + DO jl = 1, klon IF (ktest(jl)==1) THEN zslow = sqrt(pulow(jl)**2+pvlow(jl)**2) @@ -238,7 +229,8 @@ END IF END IF -530 CONTINUE + end DO + end DO ! 6. LOW LEVEL LIFT, SEMI IMPLICIT: ! ----------------------------------