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

Diff of /trunk/phylmd/nuage.f90

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

revision 336 by guez, Thu Jun 13 14:40:06 2019 UTC revision 337 by guez, Mon Sep 16 16:54:50 2019 UTC
# Line 4  module nuage_m Line 4  module nuage_m
4    
5  contains  contains
6    
7    SUBROUTINE nuage (paprs, pplay, t, pqlwp, pclc, pcltau, pclemi, pch, pcl, &    SUBROUTINE nuage (paprs, pplay, t, pqlwp, pclc, cldtau, pclemi, pch, pcl, &
8         pcm, pct, pctlwp)         pcm, pct, pctlwp)
9    
10      ! From LMDZ4/libf/phylmd/nuage.F, version 1.1.1.1, 2004/05/19 12:53:07      ! From LMDZ4/libf/phylmd/nuage.F, version 1.1.1.1, 2004/05/19 12:53:07
# Line 28  contains Line 28  contains
28      REAL, intent(inout):: pclc(klon, klev)      REAL, intent(inout):: pclc(klon, klev)
29      ! couverture nuageuse pour le rayonnement (0 \`a 1)      ! couverture nuageuse pour le rayonnement (0 \`a 1)
30    
31      REAL, intent(out):: pcltau(klon, klev) ! \'epaisseur optique des nuages      REAL, intent(out):: cldtau(klon, klev) ! \'epaisseur optique des nuages
32      real, intent(out):: pclemi(klon, klev) ! \'emissivit\'e des nuages (0 \`a 1)      real, intent(out):: pclemi(klon, klev) ! \'emissivit\'e des nuages (0 \`a 1)
33      REAL pch(klon), pcl(klon), pcm(klon), pct(klon), pctlwp(klon)      REAL pch(klon), pcl(klon), pcm(klon), pct(klon), pctlwp(klon)
34    
# Line 69  contains Line 69  contains
69    
70            radius = rad_chaud * (1.-zfice) + rad_froid * zfice            radius = rad_chaud * (1.-zfice) + rad_froid * zfice
71            coef = coef_chau * (1.-zfice) + coef_froi * zfice            coef = coef_chau * (1.-zfice) + coef_froi * zfice
72            pcltau(i, k) = 3.0/2.0 * zflwp / radius            cldtau(i, k) = 3.0/2.0 * zflwp / radius
73            pclemi(i, k) = 1.0 - EXP(- coef * zflwp)            pclemi(i, k) = 1.0 - EXP(- coef * zflwp)
74            lo = (pclc(i, k) <= seuil_neb)            lo = (pclc(i, k) <= seuil_neb)
75            IF (lo) pclc(i, k) = 0.0            IF (lo) pclc(i, k) = 0.0
76            IF (lo) pcltau(i, k) = 0.0            IF (lo) cldtau(i, k) = 0.0
77            IF (lo) pclemi(i, k) = 0.0            IF (lo) pclemi(i, k) = 0.0
78         END DO         END DO
79      END DO      END DO
Legend:
Removed from v.336  
changed lines
  Added in v.337
  ViewVC Help
Powered by ViewVC 1.1.21