/[lmdze]/trunk/Sources/phylmd/Orography/orolift.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Orography/orolift.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/Orography/orolift.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/phylmd/Orography/orolift.f revision 158 by guez, Tue Jul 21 14:44:45 2015 UTC
# Line 156  Line 156 
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 &
# Line 182  Line 175 
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)
# Line 193  Line 184 
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
# Line 217  Line 208 
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)
# Line 238  Line 229 
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  !  ----------------------------------  !  ----------------------------------
Legend:
Removed from v.82  
changed lines
  Added in v.158
  ViewVC Help
Powered by ViewVC 1.1.21