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 |
!====================================================================== |
!====================================================================== |
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 |
! |
! |
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 |
|
|