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