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

Diff of /trunk/Sources/phylmd/nuage.f

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC 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 122  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 147  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.134  
changed lines
  Added in v.169
  ViewVC Help
Powered by ViewVC 1.1.21