/[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 108 by guez, Tue Sep 16 14:00:41 2014 UTC trunk/Sources/phylmd/Orography/orolift.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 46  Line 46 
46  !*       0.2   LOCAL ARRAYS  !*       0.2   LOCAL ARRAYS
47  !              ------------  !              ------------
48        LOGICAL lifthigh        LOGICAL lifthigh
49        INTEGER klevm1, jl, ilevh, jk        INTEGER jl, jk
50        REAL zcons1, ztmst, zrtmst, zpi, zhgeo        REAL zcons1, ztmst, zpi, zhgeo
51        REAL zdelp, zslow, zsqua, zscav, zbet        REAL zdelp, zslow, zsqua, zscav, zbet
52        INTEGER iknub(klon), iknul(klon)        INTEGER iknub(klon), iknul(klon)
53        LOGICAL ll1(klon,klev+1)        LOGICAL ll1(klon,klev+1)
# Line 64  Line 64 
64    
65        IF (nlon/=klon .OR. nlev/=klev) STOP        IF (nlon/=klon .OR. nlev/=klev) STOP
66        zcons1 = 1./rd        zcons1 = 1./rd
       klevm1 = klev - 1  
67        ztmst = ptsphy        ztmst = ptsphy
       zrtmst = 1./ztmst  
68        zpi = acos(-1.)        zpi = acos(-1.)
69    
70        DO 1001 jl = 1, klon        DO 1001 jl = 1, klon
# Line 75  Line 73 
73          pvlow(jl) = 0.0          pvlow(jl) = 0.0
74          iknub(jl) = klev          iknub(jl) = klev
75          iknul(jl) = klev          iknul(jl) = klev
         ilevh = klev/3  
76          ll1(jl,klev+1) = .FALSE.          ll1(jl,klev+1) = .FALSE.
77          DO 1000 jk = 1, klev          DO 1000 jk = 1, klev
78            pvom(jl,jk) = 0.0            pvom(jl,jk) = 0.0
# Line 156  Line 153 
153          END IF          END IF
154  2110  CONTINUE  2110  CONTINUE
155    
   
 200   CONTINUE  
   
 !***********************************************************  
   
156  !*         3.      COMPUTE MOUNTAIN LIFT  !*         3.      COMPUTE MOUNTAIN LIFT
157    
 300   CONTINUE  
   
158        DO 301 jl = 1, klon        DO 301 jl = 1, klon
159          IF (ktest(jl)==1) THEN          IF (ktest(jl)==1) THEN
160            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 172 
172  !*                 --------------------  !*                 --------------------
173    
174    
 400   CONTINUE  
   
175        DO jk = 1, klev        DO jk = 1, klev
176          DO jl = 1, klon          DO jl = 1, klon
177            IF (ktest(jl)==1) THEN            IF (ktest(jl)==1) THEN
# Line 201  Line 189 
189  !*                 -------------------  !*                 -------------------
190        IF (lifthigh) THEN        IF (lifthigh) THEN
191    
 500     CONTINUE  
192  !     PRINT *,'  DANS OROLIFT: 500'  !     PRINT *,'  DANS OROLIFT: 500'
193    
194  !  EXPLICIT SOLUTION AT ALL LEVELS  !  EXPLICIT SOLUTION AT ALL LEVELS
Legend:
Removed from v.108  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21