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

Diff of /trunk/phylmd/nuage.f

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

trunk/phylmd/nuage.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/Sources/phylmd/nuage.f revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC
# Line 10  SUBROUTINE nuage (paprs, pplay, & Line 10  SUBROUTINE nuage (paprs, pplay, &
10    !    !
11    use dimens_m    use dimens_m
12    use dimphy    use dimphy
13      use nr_util, only: pi
14    use SUPHEC_M    use SUPHEC_M
15    IMPLICIT none    IMPLICIT none
16    !======================================================================    !======================================================================
# Line 55  SUBROUTINE nuage (paprs, pplay, & Line 56  SUBROUTINE nuage (paprs, pplay, &
56    PARAMETER (cetahb = 0.45, cetamb = 0.80)    PARAMETER (cetahb = 0.45, cetamb = 0.80)
57    !    !
58    INTEGER i, k    INTEGER i, k
59    REAL zflwp, zradef, zfice, zmsac    REAL zflwp, zfice
60    !    !
61    REAL radius, rad_froid, rad_chaud, rad_chau1, rad_chau2    REAL radius, rad_froid, rad_chaud, rad_chau1, rad_chau2
62    PARAMETER (rad_chau1=13.0, rad_chau2=9.0, rad_froid=35.0)    PARAMETER (rad_chau1=13.0, rad_chau2=9.0, rad_froid=35.0)
# Line 78  SUBROUTINE nuage (paprs, pplay, & Line 79  SUBROUTINE nuage (paprs, pplay, &
79    REAL re(klon, klev)       ! cloud droplet effective radius [um]    REAL re(klon, klev)       ! cloud droplet effective radius [um]
80    REAL sulfate_pi(klon, klev)  ! sulfate aerosol mass concentration [ug m-3] (pre-industrial value)    REAL sulfate_pi(klon, klev)  ! sulfate aerosol mass concentration [ug m-3] (pre-industrial value)
81    REAL cdnc_pi(klon, klev)     ! cloud droplet number concentration [m-3] (pi value)    REAL cdnc_pi(klon, klev)     ! cloud droplet number concentration [m-3] (pi value)
   REAL re_pi(klon, klev)       ! cloud droplet effective radius [um] (pi value)  
82    
83    REAL fl(klon, klev)  ! xliq * rneb (denominator to re ; fraction of liquid water clouds within the grid cell)    REAL fl(klon, klev)  ! xliq * rneb (denominator to re ; fraction of liquid water clouds within the grid cell)
84    
# Line 123  SUBROUTINE nuage (paprs, pplay, & Line 123  SUBROUTINE nuage (paprs, pplay, &
123             !             !
124             rad_chaud =  &             rad_chaud =  &
125                  1.1 * ( (pqlwp(i,k) * pplay(i,k) / (RD * T(i,k)) )   &                  1.1 * ( (pqlwp(i,k) * pplay(i,k) / (RD * T(i,k)) )   &
126                  / (4./3. * RPI * 1000. * cdnc(i,k)) )**(1./3.)                  / (4./3. * PI * 1000. * cdnc(i,k)) )**(1./3.)
127             !             !
128             ! Convert to um. CDR shall be at least 3 um.             ! Convert to um. CDR shall be at least 3 um.
129             !             !
# Line 148  SUBROUTINE nuage (paprs, pplay, & Line 148  SUBROUTINE nuage (paprs, pplay, &
148             ! ice cloud contribution) but using cdnc_pi instead of             ! ice cloud contribution) but using cdnc_pi instead of
149             ! cdnc.             ! cdnc.
150             radius = MAX(1.1e6 * ( (pqlwp(i,k)*pplay(i,k)/(RD*T(i,k)))   &             radius = MAX(1.1e6 * ( (pqlwp(i,k)*pplay(i,k)/(RD*T(i,k)))   &
151                  / (4./3.*RPI*1000.*cdnc_pi(i,k)) )**(1./3.),  &                  / (4./3.*PI*1000.*cdnc_pi(i,k)) )**(1./3.),  &
152                  3.) * (1.-zfice) + rad_froid * zfice                  3.) * (1.-zfice) + rad_froid * zfice
153             cldtaupi(i,k) = 3.0/2.0 * zflwp / radius             cldtaupi(i,k) = 3.0/2.0 * zflwp / radius
154    
Legend:
Removed from v.76  
changed lines
  Added in v.169
  ViewVC Help
Powered by ViewVC 1.1.21