/[lmdze]/trunk/phylmd/fisrtilp.f90
ViewVC logotype

Diff of /trunk/phylmd/fisrtilp.f90

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

revision 338 by guez, Thu Jun 13 14:40:06 2019 UTC revision 339 by guez, Thu Sep 26 17:08:42 2019 UTC
# Line 5  module fisrtilp_m Line 5  module fisrtilp_m
5  contains  contains
6    
7    SUBROUTINE fisrtilp(paprs, pplay, t, q, ptconv, ratqs, d_t, d_q, d_ql, rneb, &    SUBROUTINE fisrtilp(paprs, pplay, t, q, ptconv, ratqs, d_t, d_q, d_ql, rneb, &
8         radliq, rain, snow, pfrac_impa, pfrac_nucl, pfrac_1nucl, frac_impa, &         cldliq, rain, snow, pfrac_impa, pfrac_nucl, pfrac_1nucl, frac_impa, &
9         frac_nucl, prfl, psfl, rhcl)         frac_nucl, prfl, psfl, rhcl)
10    
11      ! From phylmd/fisrtilp.F, version 1.2, 2004/11/09 16:55:40      ! From phylmd/fisrtilp.F, version 1.2, 2004/11/09 16:55:40
# Line 38  contains Line 38  contains
38      REAL, INTENT (out):: d_ql(klon, klev) ! incrementation de l'eau liquide      REAL, INTENT (out):: d_ql(klon, klev) ! incrementation de l'eau liquide
39      REAL, INTENT (out):: rneb(klon, klev) ! fraction nuageuse      REAL, INTENT (out):: rneb(klon, klev) ! fraction nuageuse
40    
41      REAL, INTENT (out):: radliq(klon, klev)      REAL, INTENT (out):: cldliq(klon, klev)
42      ! eau liquide utilisee dans rayonnement      ! eau liquide utilisee dans rayonnement
43    
44      REAL, INTENT (out):: rain(klon) ! pluies (mm/s)      REAL, INTENT (out):: rain(klon) ! pluies (mm/s)
# Line 166  contains Line 166  contains
166            d_q(i, k) = 0.0            d_q(i, k) = 0.0
167            d_ql(i, k) = 0.0            d_ql(i, k) = 0.0
168            rneb(i, k) = 0.0            rneb(i, k) = 0.0
169            radliq(i, k) = 0.0            cldliq(i, k) = 0.0
170            frac_nucl(i, k) = 1.            frac_nucl(i, k) = 1.
171            frac_impa(i, k) = 1.            frac_impa(i, k) = 1.
172         END DO         END DO
# Line 344  contains Line 344  contains
344               zfice(i) = min(max(zfice(i), 0.0), 1.0)               zfice(i) = min(max(zfice(i), 0.0), 1.0)
345               zfice(i) = zfice(i)**nexpo               zfice(i) = zfice(i)**nexpo
346               zneb(i) = max(rneb(i, k), seuil_neb)               zneb(i) = max(rneb(i, k), seuil_neb)
347               radliq(i, k) = zoliq(i)/real(ninter+1)               cldliq(i, k) = zoliq(i)/real(ninter+1)
348            END IF            END IF
349         END DO         END DO
350    
# Line 375  contains Line 375  contains
375                  IF (zneb(i)==seuil_neb) ztot(i) = 0.0                  IF (zneb(i)==seuil_neb) ztot(i) = 0.0
376                  ztot(i) = min(max(ztot(i), 0.0), zoliq(i))                  ztot(i) = min(max(ztot(i), 0.0), zoliq(i))
377                  zoliq(i) = max(zoliq(i)-ztot(i), 0.0)                  zoliq(i) = max(zoliq(i)-ztot(i), 0.0)
378                  radliq(i, k) = radliq(i, k) + zoliq(i)/real(ninter+1)                  cldliq(i, k) = cldliq(i, k) + zoliq(i)/real(ninter+1)
379               END IF               END IF
380            END DO            END DO
381         END DO         END DO
Legend:
Removed from v.338  
changed lines
  Added in v.339
  ViewVC Help
Powered by ViewVC 1.1.21