32 aero_data, env_state, k)
43 real(kind=
dp),
intent(out) :: k
45 real(kind=
dp) :: v1, v2, d1, d2
68 real(kind=
dp),
intent(in) :: v1
70 real(kind=
dp),
intent(in) :: v2
76 real(kind=
dp),
intent(out) :: k_min
78 real(kind=
dp),
intent(out) :: k_max
81 integer,
parameter :: n_sample = 3
83 real(kind=
dp) :: d1, d2, d_min, d_max, k
87 d_min = minval(aero_data%density)
88 d_max = maxval(aero_data%density)
102 k_min = min(k_min, k)
103 k_max = max(k_max, k)
121 real(kind=
dp),
intent(in) :: v1
123 real(kind=
dp),
intent(in) :: d1
125 real(kind=
dp),
intent(in) :: v2
127 real(kind=
dp),
intent(in) :: d2
131 real(kind=
dp),
intent(in) :: temp
133 real(kind=
dp),
intent(out) :: bckernel
135 real(kind=
dp) :: rho, n_i, n_j, df, r0, n_i_inv_df, n_j_inv_df
137 rho = (v1 / (v1 + v2)) * d1 + (v2 / (v1 + v2)) * d2
140 df = aero_data%fractal%frac_dim
141 r0 = aero_data%fractal%prime_radius
142 n_i_inv_df = n_i**(1d0 / df)
143 n_j_inv_df = n_j**(1d0 / df)
144 bckernel = sqrt(6d0 *
const%boltzmann * temp * r0 / rho &
145 * (1d0 / n_i + 1d0 / n_j)) * (n_i_inv_df + n_j_inv_df)**2d0