115 |
! for liquid water clouds: |
! for liquid water clouds: |
116 |
IF (ok_aie) THEN |
IF (ok_aie) THEN |
117 |
cdnc(i, k) = 10.**(bl95_b0 + bl95_b1 & |
cdnc(i, k) = 10.**(bl95_b0 + bl95_b1 & |
118 |
* log10(MAX(sulfate(i, k), 1e-4))) * 1.e6 |
* log10(MAX(sulfate(i, k), 1e-4)) + 6.) |
119 |
cdnc_pi(i, k) = 10.**(bl95_b0 + bl95_b1 & |
cdnc_pi(i, k) = 10.**(bl95_b0 + bl95_b1 & |
120 |
* log10(MAX(sulfate_pi(i, k), 1e-4))) * 1e6 |
* log10(MAX(sulfate_pi(i, k), 1e-4)) + 6.) |
121 |
|
|
122 |
! Restrict to interval [20, 1000] cm^3: |
! Restrict to interval [20, 1000] cm-3: |
123 |
cdnc(i, k) = MIN(1000e6, MAX(20e6, cdnc(i, k))) |
cdnc(i, k) = MIN(1000e6, MAX(20e6, cdnc(i, k))) |
124 |
cdnc_pi(i, k) = MIN(1000e6, MAX(20e6, cdnc_pi(i, k))) |
cdnc_pi(i, k) = MIN(1000e6, MAX(20e6, cdnc_pi(i, k))) |
125 |
|
|