1 | ! |
---|
2 | ! $Id: nucleation_tstep_mod.F90 3690 2020-05-30 14:14:30Z nfevrier $ |
---|
3 | ! |
---|
4 | MODULE nucleation_tstep_mod |
---|
5 | |
---|
6 | CONTAINS |
---|
7 | |
---|
8 | SUBROUTINE nucleation_rate(rhoa,t_seri,pplay,rh,a_xm,b_xm,c_xm,nucl_rate,ntot_n,x_n) |
---|
9 | |
---|
10 | USE aerophys |
---|
11 | USE infotrac_phy |
---|
12 | USE YOMCST, ONLY : RPI, RD, RMD, RKBOL, RNAVO |
---|
13 | |
---|
14 | IMPLICIT NONE |
---|
15 | |
---|
16 | ! input variables |
---|
17 | LOGICAL, PARAMETER :: flag_new_nucl=.TRUE. |
---|
18 | REAL rhoa ! H2SO4 number density [molecules/cm3] |
---|
19 | REAL t_seri ! temperature (K) |
---|
20 | REAL pplay ! pressure (Pa) |
---|
21 | REAL rh ! relative humidity (between 0 and 1) |
---|
22 | REAL a_xm, b_xm, c_xm |
---|
23 | |
---|
24 | ! output variables |
---|
25 | REAL nucl_rate |
---|
26 | REAL ntot_n ! total number of molecules in the critical cluster for neutral nucleation |
---|
27 | REAL x_n ! mole fraction of H2SO4 in the critical cluster for neutral nucleation |
---|
28 | |
---|
29 | ! local variables |
---|
30 | REAL jnuc_n ! nucleation rate in 1/cm3s (10^-7-10^10 1/cm3s) for neutral nucleation |
---|
31 | REAL rc_n ! radius of the critical cluster in nm for neutral nucleation |
---|
32 | REAL na_n ! sulfuric acid molecules in the neutral critical cluster (NOT IN USE) |
---|
33 | LOGICAL kinetic_n ! true if kinetic neutral nucleation (NOT IN USE) |
---|
34 | LOGICAL kinetic_i ! true if kinetic ion-induced nucleation (NOT IN USE) |
---|
35 | REAL rhoatres ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation (NOT IN USE) |
---|
36 | REAL VH2SO4mol |
---|
37 | REAL ntot_i, x_i, jnuc_i, rc_i, na_i, n_i ! quantities for charged nucleation (NOT IN USE) |
---|
38 | REAL csi ! Ion condensation sink (s-1) NOT IN USE |
---|
39 | REAL airn ! Air molecule concentration in (cm-3) NOT IN USE |
---|
40 | REAL ipr ! Ion pair production rate (cm-3 s-1) NOT IN USE |
---|
41 | |
---|
42 | ! call nucleation routine |
---|
43 | IF (.NOT.flag_new_nucl) THEN |
---|
44 | ! Use older routine from Hanna Vehkamäki (FMI) |
---|
45 | CALL binapara(t_seri,rh,rhoa,jnuc_n,x_n,ntot_n,rc_n) |
---|
46 | ! when total number of molecules is too small |
---|
47 | ! then set jnuc_n to collision rate of two H2SO4 molecules (following personal communication of Ulrike Niemeier and Hanna Vehkamäki) |
---|
48 | IF (ntot_n < 4.0) THEN |
---|
49 | VH2SO4mol=mH2SO4mol/(1.E-3*(a_xm+t_seri*(b_xm+t_seri*c_xm))) !cm3 |
---|
50 | jnuc_n = rhoa**2. *(3./4.*RPI)**(1./6.) *(12.*RKBOL*t_seri/mH2SO4mol)**0.5 & |
---|
51 | & *100.*(2.*VH2SO4mol**(1./3.))**2. !1/(cm3s) |
---|
52 | ntot_n=2.0 |
---|
53 | x_n=1.0 |
---|
54 | ENDIF |
---|
55 | ELSE |
---|
56 | ! Use new routine from Anni Maattanen (LATMOS) |
---|
57 | csi=0.0 ! no charged nucleation for now |
---|
58 | ipr=-1.0 ! dummy value to make sure charged nucleation does not occur |
---|
59 | airn=0.0 ! NOT IN USE |
---|
60 | ! airn=pplay/t_seri/RD/1.E3*RNAVO/RMD ! molec cm-3 (for future use, to be confirmed) |
---|
61 | CALL newbinapara(t_seri,rh,rhoa,csi,airn,ipr,jnuc_n,ntot_n,jnuc_i,ntot_i, & |
---|
62 | & x_n,x_i,na_n,na_i,rc_n,rc_i,n_i,kinetic_n,kinetic_i,rhoatres) |
---|
63 | ENDIF |
---|
64 | |
---|
65 | ! convert jnuc_n from particles/cm3/s to kg(H2SO4)/kgA/s |
---|
66 | nucl_rate=jnuc_n*ntot_n*x_n*mH2SO4mol/(pplay/t_seri/RD/1.E6) |
---|
67 | |
---|
68 | END SUBROUTINE nucleation_rate |
---|
69 | |
---|
70 | !-------------------------------------------------------------------------------------------------- |
---|
71 | |
---|
72 | SUBROUTINE nucleation_part(nucl_rate,ntot,x,dt,Vbin,tr_seri) |
---|
73 | |
---|
74 | USE aerophys |
---|
75 | USE infotrac_phy |
---|
76 | |
---|
77 | IMPLICIT NONE |
---|
78 | |
---|
79 | ! input variable |
---|
80 | REAL nucl_rate |
---|
81 | REAL ntot ! total number of molecules in the critical cluster |
---|
82 | REAL x ! mole raction of H2SO4 in the critical cluster |
---|
83 | REAL dt |
---|
84 | REAL Vbin(nbtr_bin) |
---|
85 | |
---|
86 | ! output variable |
---|
87 | REAL tr_seri(nbtr) |
---|
88 | |
---|
89 | ! local variable |
---|
90 | INTEGER k |
---|
91 | REAL Vnew |
---|
92 | REAL ff(nbtr_bin) |
---|
93 | |
---|
94 | ! dry volume of nucleated particle (at reference temperature) |
---|
95 | Vnew=mH2SO4mol*ntot*x/dens_aer_dry |
---|
96 | |
---|
97 | ! compute distribution factor for particles of intermediate size (from Jacobson 1994, equation 13) |
---|
98 | ff(:)=0.0 |
---|
99 | DO k=1, nbtr_bin |
---|
100 | ! CK 20160531: bug fix for first bin |
---|
101 | IF (k.LE.nbtr_bin-1) THEN |
---|
102 | IF (Vbin(k).LE.Vnew.AND.Vnew.LT.Vbin(k+1)) THEN |
---|
103 | ff(k)= Vbin(k)/Vnew*(Vbin(k+1)-Vnew)/(Vbin(k+1)-Vbin(k)) |
---|
104 | ENDIF |
---|
105 | ENDIF |
---|
106 | IF (k.EQ.1.AND.Vnew.LE.Vbin(k)) THEN |
---|
107 | ff(k)= 1. |
---|
108 | ENDIF |
---|
109 | IF (k.GT.1) THEN |
---|
110 | IF (Vbin(k-1).LT.Vnew.AND.Vnew.LT.Vbin(k)) THEN |
---|
111 | ff(k)= 1.-ff(k-1) |
---|
112 | ENDIF |
---|
113 | ENDIF |
---|
114 | IF (k.EQ.nbtr_bin.AND.Vnew.GE.Vbin(k)) THEN |
---|
115 | ff(k)= 1. |
---|
116 | ENDIF |
---|
117 | ENDDO |
---|
118 | ! correction of ff for volume conservation |
---|
119 | DO k=1, nbtr_bin |
---|
120 | ff(k)=ff(k)*Vnew/Vbin(k) |
---|
121 | ENDDO |
---|
122 | |
---|
123 | ! add nucleated H2SO4 to corresponding aerosol bin |
---|
124 | DO k=1, nbtr_bin |
---|
125 | tr_seri(k+nbtr_sulgas)=tr_seri(k+nbtr_sulgas)+nucl_rate*dt/(mH2SO4mol*ntot*x)*ff(k) |
---|
126 | ENDDO |
---|
127 | |
---|
128 | END SUBROUTINE nucleation_part |
---|
129 | |
---|
130 | !--------------------------------------------------------------------------------------------------- |
---|
131 | |
---|
132 | SUBROUTINE binapara(pt,prh,rhoa_in,jnuc,x,ntot,rc) |
---|
133 | |
---|
134 | |
---|
135 | ! Fortran 90 subroutine binapara |
---|
136 | ! |
---|
137 | ! Calculates parametrized values of nucleation rate, |
---|
138 | ! mole fraction of sulphuric acid |
---|
139 | ! total number of particles, and the radius of the critical cluster |
---|
140 | ! in H2O-H2SO4 system IF temperature, saturatio ratio of water and |
---|
141 | ! sulfuric acid concentration are given. |
---|
142 | ! |
---|
143 | ! Copyright (C) 2002 Hanna Vehkamäki |
---|
144 | ! |
---|
145 | ! Division of Atmospheric Sciences |
---|
146 | ! Department of Physical Sciences |
---|
147 | ! P.O. Box 64 |
---|
148 | ! FIN-00014 University of Helsinki |
---|
149 | ! Finland |
---|
150 | ! |
---|
151 | ! hanna.vehkamaki@helsinki.fi |
---|
152 | |
---|
153 | IMPLICIT NONE |
---|
154 | |
---|
155 | REAL :: pt,t !temperature in K (190.15-300.15K) |
---|
156 | REAL :: prh,rh !saturation ratio of water (0.0001-1) |
---|
157 | REAL,INTENT(IN) :: rhoa_in !sulfuric acid concentration in 1/cm3 (10^4-10^11 1/cm3) |
---|
158 | REAL,INTENT(OUT) :: jnuc !nucleation rate in 1/cm3s (10^-7-10^10 1/cm3s) |
---|
159 | REAL,INTENT(OUT) :: ntot !total number of molecules in the critical cluster (ntot>4) |
---|
160 | REAL,INTENT(OUT) :: x ! mole fraction of H2SO4 in the critical cluster |
---|
161 | REAL,INTENT(OUT) :: rc !radius of the critical cluster in nm |
---|
162 | REAL :: rhotres ! threshold concentration of h2so4 (1/cm^3) |
---|
163 | ! which produces nucleation rate 1/(cm^3 s) as a function of rh and t |
---|
164 | REAL rhoa |
---|
165 | |
---|
166 | ! CK: use intermediate variable to avoid overwriting |
---|
167 | t=pt |
---|
168 | rh=prh |
---|
169 | rhoa=rhoa_in |
---|
170 | |
---|
171 | IF (t < 190.15) THEN |
---|
172 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): temperature < 190.15 K, using 190.15 K (T=',t,'K)' |
---|
173 | t=190.15 |
---|
174 | ENDIF |
---|
175 | |
---|
176 | IF (t > 300.15) THEN |
---|
177 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): temperature > 300.15 K, using 300.15 K' |
---|
178 | t=300.15 |
---|
179 | ENDIF |
---|
180 | |
---|
181 | IF (rh < 0.0001) THEN |
---|
182 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): saturation ratio of water < 0.0001, using 0.0001' |
---|
183 | rh=0.0001 |
---|
184 | ENDIF |
---|
185 | |
---|
186 | IF (rh > 1.0) THEN |
---|
187 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): saturation ratio of water > 1 using 1' |
---|
188 | ! print *, 'rh=',rh |
---|
189 | rh=1.0 |
---|
190 | ENDIF |
---|
191 | |
---|
192 | IF (rhoa < 1.e4) THEN |
---|
193 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): sulfuric acid concentration < 1e4 1/cm3, using 1e4 1/cm3' |
---|
194 | rhoa=1.e4 |
---|
195 | ENDIF |
---|
196 | |
---|
197 | IF (rhoa > 1.e11) THEN |
---|
198 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev,'): sulfuric acid concentration > 1e11 1/cm3, using 1e11 1/cm3' |
---|
199 | rhoa=1.e11 |
---|
200 | ENDIF |
---|
201 | |
---|
202 | x= 0.7409967177282139 - 0.002663785665140117*t + 0.002010478847383187*LOG(rh) & |
---|
203 | & - 0.0001832894131464668*t*LOG(rh) + 0.001574072538464286*LOG(rh)**2 & |
---|
204 | & - 0.00001790589121766952*t*LOG(rh)**2 + 0.0001844027436573778*LOG(rh)**3 & |
---|
205 | & - 1.503452308794887E-6*t*LOG(rh)**3 - 0.003499978417957668*LOG(rhoa) & |
---|
206 | & + 0.0000504021689382576*t*LOG(rhoa) |
---|
207 | |
---|
208 | jnuc= 0.1430901615568665 + 2.219563673425199*t - 0.02739106114964264*t**2 + & |
---|
209 | & 0.00007228107239317088*t**3 + 5.91822263375044/x + & |
---|
210 | & 0.1174886643003278*LOG(rh) + 0.4625315047693772*t*LOG(rh) - & |
---|
211 | & 0.01180591129059253*t**2*LOG(rh) + & |
---|
212 | & 0.0000404196487152575*t**3*LOG(rh) + (15.79628615047088*LOG(rh))/x - & |
---|
213 | & 0.215553951893509*LOG(rh)**2 - 0.0810269192332194*t*LOG(rh)**2 + & |
---|
214 | & 0.001435808434184642*t**2*LOG(rh)**2 - & |
---|
215 | & 4.775796947178588E-6*t**3*LOG(rh)**2 - & |
---|
216 | & (2.912974063702185*LOG(rh)**2)/x - 3.588557942822751*LOG(rh)**3 + & |
---|
217 | & 0.04950795302831703*t*LOG(rh)**3 - & |
---|
218 | & 0.0002138195118737068*t**2*LOG(rh)**3 + & |
---|
219 | & 3.108005107949533E-7*t**3*LOG(rh)**3 - & |
---|
220 | & (0.02933332747098296*LOG(rh)**3)/x + & |
---|
221 | & 1.145983818561277*LOG(rhoa) - & |
---|
222 | & 0.6007956227856778*t*LOG(rhoa) + & |
---|
223 | & 0.00864244733283759*t**2*LOG(rhoa) - & |
---|
224 | & 0.00002289467254710888*t**3*LOG(rhoa) - & |
---|
225 | & (8.44984513869014*LOG(rhoa))/x + & |
---|
226 | & 2.158548369286559*LOG(rh)*LOG(rhoa) + & |
---|
227 | & 0.0808121412840917*t*LOG(rh)*LOG(rhoa) - & |
---|
228 | & 0.0004073815255395214*t**2*LOG(rh)*LOG(rhoa) - & |
---|
229 | & 4.019572560156515E-7*t**3*LOG(rh)*LOG(rhoa) + & |
---|
230 | & (0.7213255852557236*LOG(rh)*LOG(rhoa))/x + & |
---|
231 | & 1.62409850488771*LOG(rh)**2*LOG(rhoa) - & |
---|
232 | & 0.01601062035325362*t*LOG(rh)**2*LOG(rhoa) + & |
---|
233 | & 0.00003771238979714162*t**2*LOG(rh)**2*LOG(rhoa) + & |
---|
234 | & 3.217942606371182E-8*t**3*LOG(rh)**2*LOG(rhoa) - & |
---|
235 | & (0.01132550810022116*LOG(rh)**2*LOG(rhoa))/x + & |
---|
236 | & 9.71681713056504*LOG(rhoa)**2 - & |
---|
237 | & 0.1150478558347306*t*LOG(rhoa)**2 + & |
---|
238 | & 0.0001570982486038294*t**2*LOG(rhoa)**2 + & |
---|
239 | & 4.009144680125015E-7*t**3*LOG(rhoa)**2 + & |
---|
240 | & (0.7118597859976135*LOG(rhoa)**2)/x - & |
---|
241 | & 1.056105824379897*LOG(rh)*LOG(rhoa)**2 + & |
---|
242 | & 0.00903377584628419*t*LOG(rh)*LOG(rhoa)**2 - & |
---|
243 | & 0.00001984167387090606*t**2*LOG(rh)*LOG(rhoa)**2 + & |
---|
244 | & 2.460478196482179E-8*t**3*LOG(rh)*LOG(rhoa)**2 - & |
---|
245 | & (0.05790872906645181*LOG(rh)*LOG(rhoa)**2)/x - & |
---|
246 | & 0.1487119673397459*LOG(rhoa)**3 + & |
---|
247 | & 0.002835082097822667*t*LOG(rhoa)**3 - & |
---|
248 | & 9.24618825471694E-6*t**2*LOG(rhoa)**3 + & |
---|
249 | & 5.004267665960894E-9*t**3*LOG(rhoa)**3 - & |
---|
250 | & (0.01270805101481648*LOG(rhoa)**3)/x |
---|
251 | jnuc=EXP(jnuc) !1/(cm3s) |
---|
252 | |
---|
253 | ntot =-0.002954125078716302 - 0.0976834264241286*t + 0.001024847927067835*t**2 - 2.186459697726116E-6*t**3 - & |
---|
254 | & 0.1017165718716887/x - 0.002050640345231486*LOG(rh) - 0.007585041382707174*t*LOG(rh) + & |
---|
255 | & 0.0001926539658089536*t**2*LOG(rh) - 6.70429719683894E-7*t**3*LOG(rh) - & |
---|
256 | & (0.2557744774673163*LOG(rh))/x + 0.003223076552477191*LOG(rh)**2 + 0.000852636632240633*t*LOG(rh)**2 - & |
---|
257 | & 0.00001547571354871789*t**2*LOG(rh)**2 + 5.666608424980593E-8*t**3*LOG(rh)**2 + & |
---|
258 | & (0.03384437400744206*LOG(rh)**2)/x + 0.04743226764572505*LOG(rh)**3 - & |
---|
259 | & 0.0006251042204583412*t*LOG(rh)**3 + 2.650663328519478E-6*t**2*LOG(rh)**3 - & |
---|
260 | & 3.674710848763778E-9*t**3*LOG(rh)**3 - (0.0002672510825259393*LOG(rh)**3)/x - & |
---|
261 | & 0.01252108546759328*LOG(rhoa) + 0.005806550506277202*t*LOG(rhoa) - & |
---|
262 | & 0.0001016735312443444*t**2*LOG(rhoa) + 2.881946187214505E-7*t**3*LOG(rhoa) + & |
---|
263 | & (0.0942243379396279*LOG(rhoa))/x - 0.0385459592773097*LOG(rh)*LOG(rhoa) - & |
---|
264 | & 0.0006723156277391984*t*LOG(rh)*LOG(rhoa) + 2.602884877659698E-6*t**2*LOG(rh)*LOG(rhoa) + & |
---|
265 | & 1.194163699688297E-8*t**3*LOG(rh)*LOG(rhoa) - (0.00851515345806281*LOG(rh)*LOG(rhoa))/x - & |
---|
266 | & 0.01837488495738111*LOG(rh)**2*LOG(rhoa) + 0.0001720723574407498*t*LOG(rh)**2*LOG(rhoa) - & |
---|
267 | & 3.717657974086814E-7*t**2*LOG(rh)**2*LOG(rhoa) - & |
---|
268 | & 5.148746022615196E-10*t**3*LOG(rh)**2*LOG(rhoa) + & |
---|
269 | & (0.0002686602132926594*LOG(rh)**2*LOG(rhoa))/x - 0.06199739728812199*LOG(rhoa)**2 + & |
---|
270 | & 0.000906958053583576*t*LOG(rhoa)**2 - 9.11727926129757E-7*t**2*LOG(rhoa)**2 - & |
---|
271 | & 5.367963396508457E-9*t**3*LOG(rhoa)**2 - (0.007742343393937707*LOG(rhoa)**2)/x + & |
---|
272 | & 0.0121827103101659*LOG(rh)*LOG(rhoa)**2 - 0.0001066499571188091*t*LOG(rh)*LOG(rhoa)**2 + & |
---|
273 | & 2.534598655067518E-7*t**2*LOG(rh)*LOG(rhoa)**2 - & |
---|
274 | & 3.635186504599571E-10*t**3*LOG(rh)*LOG(rhoa)**2 + & |
---|
275 | & (0.0006100650851863252*LOG(rh)*LOG(rhoa)**2)/x + 0.0003201836700403512*LOG(rhoa)**3 - & |
---|
276 | & 0.0000174761713262546*t*LOG(rhoa)**3 + 6.065037668052182E-8*t**2*LOG(rhoa)**3 - & |
---|
277 | & 1.421771723004557E-11*t**3*LOG(rhoa)**3 + (0.0001357509859501723*LOG(rhoa)**3)/x |
---|
278 | ntot=EXP(ntot) |
---|
279 | |
---|
280 | rc=EXP(-1.6524245+0.42316402*x+0.33466487*LOG(ntot)) !nm |
---|
281 | |
---|
282 | IF (jnuc < 1.E-7) THEN |
---|
283 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev'): nucleation rate < 1E-7/cm3s, using 0.0/cm3s,' |
---|
284 | jnuc=0.0 |
---|
285 | ENDIF |
---|
286 | |
---|
287 | IF (jnuc > 1.e10) THEN |
---|
288 | ! print *,'Warning (ilon=',ilon,'ilev=',ilev, & |
---|
289 | ! & '): nucleation rate > 1e10/cm3s, using 1e10/cm3s' |
---|
290 | jnuc=1.e10 |
---|
291 | ENDIF |
---|
292 | |
---|
293 | rhotres=EXP( -279.2430007512709 + 11.73439886096903*rh + 22700.92970508331/t & |
---|
294 | & - (1088.644983466801*rh)/t + 1.144362942094912*t & |
---|
295 | & - 0.03023314602163684*rh*t - 0.001302541390154324*t**2 & |
---|
296 | & - 6.386965238433532*LOG(rh) + (854.980361026715*LOG(rh))/t & |
---|
297 | & + 0.00879662256826497*t*LOG(rh)) !1/cm3 |
---|
298 | |
---|
299 | RETURN |
---|
300 | |
---|
301 | END SUBROUTINE binapara |
---|
302 | |
---|
303 | !--------------------------------------------------------------------------------------------------- |
---|
304 | |
---|
305 | SUBROUTINE newbinapara(t,satrat,rhoa,csi,airn,ipr,jnuc_n_real,ntot_n_real,jnuc_i_real,ntot_i_real, & |
---|
306 | & x_n_real,x_i_real,na_n_real,na_i_real,rc_n_real,rc_i_real,n_i_real, & |
---|
307 | & kinetic_n,kinetic_i,rhoatres_real) |
---|
308 | |
---|
309 | ! Fortran 90 subroutine newbinapara |
---|
310 | ! |
---|
311 | ! Calculates parametrized values for neutral and ion-induced sulfuric acid-water particle formation rate |
---|
312 | ! of critical clusters, |
---|
313 | ! number of particle in the critical clusters, the radii of the critical clusters |
---|
314 | ! in H2O-H2SO4-ion system if temperature, saturation ratio of water, sulfuric acid concentration, |
---|
315 | ! and, optionally, either condensation sink due to pre-existing particle and ion pair production rate, |
---|
316 | ! or atmospheric concentration of negative ions are given. |
---|
317 | ! |
---|
318 | ! The code calculates also the kinetic limit and the particle formation rate |
---|
319 | ! above this limit (in which case we set ntot=1 and na=1) |
---|
320 | ! |
---|
321 | ! Copyright (C)2018 Määttänen et al. 2018 |
---|
322 | ! |
---|
323 | ! anni.maattanen@latmos.ipsl.fr |
---|
324 | ! joonas.merikanto@fmi.fi |
---|
325 | ! hanna.vehkamaki@helsinki.fi |
---|
326 | ! |
---|
327 | ! References |
---|
328 | ! A. Määttänen, J. Merikanto, H. Henschel, J. Duplissy, R. Makkonen, |
---|
329 | ! I. K. Ortega and H. Vehkamäki (2018), New parameterizations for |
---|
330 | ! neutral and ion-induced sulfuric acid-water particle formation in |
---|
331 | ! nucleation and kinetic regimes, J. Geophys. Res. Atmos., 122, doi:10.1002/2017JD027429. |
---|
332 | ! |
---|
333 | ! Brasseur, G., and A. Chatel (1983), paper presented at the 9th Annual Meeting of the |
---|
334 | ! European Geophysical Society, Leeds, Great Britain, August 1982. |
---|
335 | ! |
---|
336 | ! Dunne, Eimear M., et al.(2016), Global atmospheric particle formation from CERN CLOUD measurements, |
---|
337 | ! Science 354.6316, 1119-1124. |
---|
338 | ! |
---|
339 | |
---|
340 | USE aerophys |
---|
341 | USE YOMCST, ONLY : RPI, RKBOL |
---|
342 | |
---|
343 | IMPLICIT NONE |
---|
344 | |
---|
345 | !---------------------------------------------------- |
---|
346 | |
---|
347 | !Global intent in |
---|
348 | REAL,INTENT(IN) :: t ! temperature in K |
---|
349 | REAL,INTENT(IN) :: satrat ! saturatio ratio of water (between zero and 1) |
---|
350 | REAL,INTENT(IN) :: rhoa ! sulfuric acid concentration in 1/cm3 |
---|
351 | REAL,INTENT(IN) :: csi ! Ion condensation sink (s-1) |
---|
352 | REAL,INTENT(IN) :: airn ! Air molecule concentration in (cm-3) |
---|
353 | REAL,INTENT(IN) :: ipr ! Ion pair production rate (cm-3 s-1) |
---|
354 | !Global intent out |
---|
355 | REAL,INTENT(OUT) :: jnuc_n_real ! Neutral nucleation rate in 1/cm3s (J>10^-7 1/cm3s) |
---|
356 | REAL,INTENT(OUT) :: ntot_n_real ! total number of molecules in the neutral critical cluster |
---|
357 | REAL,INTENT(OUT) :: jnuc_i_real ! Charged nucleation rate in 1/cm3s (J>10^-7 1/cm3s) |
---|
358 | REAL,INTENT(OUT) :: ntot_i_real ! total number of molecules in the charged critical cluster |
---|
359 | REAL,INTENT(OUT) :: x_n_real ! mole fraction of H2SO4 in the neutral critical cluster |
---|
360 | REAL,INTENT(OUT) :: x_i_real ! mole fraction of H2SO4 in the charged critical cluster |
---|
361 | ! (note that x_n=x_i in nucleation regime) |
---|
362 | REAL,INTENT(OUT) :: na_n_real ! sulfuric acid molecules in the neutral critical cluster |
---|
363 | REAL,INTENT(OUT) :: na_i_real ! sulfuric molecules in the charged critical cluster |
---|
364 | REAL,INTENT(OUT) :: rc_n_real ! radius of the charged critical cluster in nm |
---|
365 | REAL,INTENT(OUT) :: rc_i_real ! radius of the charged critical cluster in nm |
---|
366 | REAL,INTENT(OUT) :: n_i_real ! number of ion pairs in air (cm-3) |
---|
367 | REAL,INTENT(OUT) :: rhoatres_real ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation |
---|
368 | LOGICAL,INTENT(OUT) :: kinetic_n ! true if kinetic neutral nucleation |
---|
369 | LOGICAL,INTENT(OUT) :: kinetic_i ! true if kinetic ion-induced nucleation |
---|
370 | |
---|
371 | ! Local |
---|
372 | DOUBLE PRECISION :: jnuc_n ! Neutral nucleation rate in 1/cm3s (J>10^-7 1/cm3s) |
---|
373 | DOUBLE PRECISION :: ntot_n ! total number of molecules in the neutral critical cluster |
---|
374 | DOUBLE PRECISION :: jnuc_i ! Charged nucleation rate in 1/cm3s (J>10^-7 1/cm3s) |
---|
375 | DOUBLE PRECISION :: ntot_i ! total number of molecules in the charged critical cluster |
---|
376 | DOUBLE PRECISION :: x_n ! mole fraction of H2SO4 in the neutral critical cluster |
---|
377 | DOUBLE PRECISION :: x_i ! mole fraction of H2SO4 in the charged critical cluster |
---|
378 | ! (note that x_n=x_i in nucleation regime) |
---|
379 | DOUBLE PRECISION :: na_n ! sulfuric acid molecules in the neutral critical cluster |
---|
380 | DOUBLE PRECISION :: na_i ! sulfuric molecules in the charged critical cluster |
---|
381 | DOUBLE PRECISION :: rc_n ! radius of the charged critical cluster in nm |
---|
382 | DOUBLE PRECISION :: rc_i ! radius of the charged critical cluster in nm |
---|
383 | DOUBLE PRECISION :: n_i ! number of ion pairs in air (cm-3) |
---|
384 | DOUBLE PRECISION :: rhoatres ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation |
---|
385 | DOUBLE PRECISION :: x ! mole fraction of H2SO4 in the critical cluster |
---|
386 | DOUBLE PRECISION :: satratln ! bounded water saturation ratio for neutral case (between 5.E-6 - 1.0) |
---|
387 | DOUBLE PRECISION :: satratli ! bounded water saturation ratio for ion-induced case (between 1.E-7 - 0.95) |
---|
388 | DOUBLE PRECISION :: rhoaln ! bounded concentration of h2so4 for neutral case (between 10^10 - 10^19 m-3) |
---|
389 | DOUBLE PRECISION :: rhoali ! bounded concentration of h2so4 for ion-induced case (between 10^10 - 10^22 m-3) |
---|
390 | DOUBLE PRECISION :: tln ! bounded temperature for neutral case (between 165-400 K) |
---|
391 | DOUBLE PRECISION :: tli ! bounded temperature for ion-induced case (195-400 K) |
---|
392 | DOUBLE PRECISION :: kinrhotresn ! threshold sulfuric acid for neutral kinetic nucleation |
---|
393 | DOUBLE PRECISION :: kinrhotresi ! threshold sulfuric acid for ion-induced kinetic nucleation |
---|
394 | DOUBLE PRECISION :: jnuc_i1 ! Ion-induced rate for n_i=1 cm-3 |
---|
395 | DOUBLE PRECISION :: xloss ! Ion loss rate |
---|
396 | DOUBLE PRECISION :: recomb ! Ion-ion recombination rate |
---|
397 | |
---|
398 | !--- 0) Initializations: |
---|
399 | |
---|
400 | kinetic_n=.FALSE. |
---|
401 | kinetic_i=.FALSE. |
---|
402 | jnuc_n=0.0 |
---|
403 | jnuc_i=0.0 |
---|
404 | ntot_n=0.0 |
---|
405 | ntot_i=0.0 |
---|
406 | na_n=0.0 |
---|
407 | na_i=0.0 |
---|
408 | rc_n=0.0 |
---|
409 | rc_i=0.0 |
---|
410 | x=0.0 |
---|
411 | x_n=0.0 |
---|
412 | x_i=0.0 |
---|
413 | satratln=satrat |
---|
414 | satratli=satrat |
---|
415 | rhoaln=rhoa |
---|
416 | rhoali=rhoa |
---|
417 | tln=t |
---|
418 | tli=t |
---|
419 | n_i=0.0 |
---|
420 | |
---|
421 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
422 | |
---|
423 | !Boundary values according to parameterization limits |
---|
424 | |
---|
425 | !Temperature bounds |
---|
426 | IF (t.LE.165.) THEN |
---|
427 | ! print *,'Warning: temperature < 165.0 K, using 165.0 K in neutral nucleation calculation' |
---|
428 | tln=165.0 |
---|
429 | ENDIF |
---|
430 | IF (t.LE.195.) THEN |
---|
431 | ! print *,'Warning: temperature < 195.0 K, using 195.0 K in ion-induced nucleation calculation' |
---|
432 | tli=195.0 |
---|
433 | ENDIF |
---|
434 | IF (t.GE.400.) THEN |
---|
435 | ! print *,'Warning: temperature > 400. K, using 400. K in nucleation calculations' |
---|
436 | tln=400. |
---|
437 | tli=400. |
---|
438 | ENDIF |
---|
439 | |
---|
440 | ! Saturation ratio bounds |
---|
441 | IF (satrat.LT.1.E-7) THEN |
---|
442 | ! print *,'Warning: saturation ratio of water < 1.E-7, using 1.E-7 in ion-induced nucleation calculation' |
---|
443 | satratli=1.E-7 |
---|
444 | ENDIF |
---|
445 | IF (satrat.LT.1.E-5) THEN |
---|
446 | ! print *,'Warning: saturation ratio of water < 1.E-5, using 1.E-5 in neutral nucleation calculation' |
---|
447 | satratln=1.E-5 |
---|
448 | ENDIF |
---|
449 | IF (satrat.GT.0.95) THEN |
---|
450 | ! print *,'Warning: saturation ratio of water > 0.95, using 0.95 in ion-induced nucleation calculation' |
---|
451 | satratli=0.95 |
---|
452 | ENDIF |
---|
453 | IF (satrat.GT.1.0) THEN |
---|
454 | ! print *,'Warning: saturation ratio of water > 1 using 1 in neutral nucleation calculation' |
---|
455 | satratln=1.0 |
---|
456 | ENDIF |
---|
457 | |
---|
458 | ! Sulfuric acid concentration bounds |
---|
459 | IF (rhoa.LE.1.E4) THEN |
---|
460 | ! print *,'Warning: sulfuric acid < 1e4 1/cm3, using 1e4 1/cm3 in nucleation calculation' |
---|
461 | rhoaln=1.E4 |
---|
462 | rhoali=1.E4 |
---|
463 | ENDIF |
---|
464 | IF (rhoa.GT.1.E13) THEN |
---|
465 | ! print *,'Warning: sulfuric acid > 1e13 1/cm3, using 1e13 1/cm3 in neutral nucleation calculation' |
---|
466 | rhoaln=1.E13 |
---|
467 | ENDIF |
---|
468 | IF (rhoa.GT.1.E16) THEN |
---|
469 | ! print *,'Warning: sulfuric acid concentration > 1e16 1/cm3, using 1e16 1/cm3 in ion-induced nucleation calculation' |
---|
470 | rhoali=1.E16 |
---|
471 | ENDIF |
---|
472 | |
---|
473 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
474 | |
---|
475 | !Critical cluster composition (valid for both cases, bounds not used here) |
---|
476 | x_n= 7.9036365428891719E-1 - 2.8414059650092153E-3*tln + 1.4976802556584141E-2*LOG(satratln) & |
---|
477 | & - 2.4511581740839115E-4*tln*LOG(satratln) + 3.4319869471066424E-3 *LOG(satratln)**2 & |
---|
478 | & - 2.8799393617748428E-5*tln*LOG(satratln)**2 + 3.0174314126331765E-4*LOG(satratln)**3 & |
---|
479 | & - 2.2673492408841294E-6*tln*LOG(satratln)**3 - 4.3948464567032377E-3*LOG(rhoaln) & |
---|
480 | & + 5.3305314722492146E-5*tln*LOG(rhoaln) |
---|
481 | x_i= 7.9036365428891719E-1 - 2.8414059650092153E-3*tli + 1.4976802556584141E-2*LOG(satratli) & |
---|
482 | & - 2.4511581740839115E-4*tli*LOG(satratli) + 3.4319869471066424E-3 *LOG(satratli)**2 & |
---|
483 | & - 2.8799393617748428E-5*tli*LOG(satratli)**2 + 3.0174314126331765E-4*LOG(satratli)**3 & |
---|
484 | & - 2.2673492408841294E-6*tli*LOG(satratli)**3 - 4.3948464567032377E-3*LOG(rhoali) & |
---|
485 | & + 5.3305314722492146E-5*tli*LOG(rhoali) |
---|
486 | |
---|
487 | x_n=MIN(MAX(x_n,1.E-30),1.) |
---|
488 | x_i=MIN(MAX(x_i,1.E-30),1.) |
---|
489 | |
---|
490 | !Neutral nucleation |
---|
491 | |
---|
492 | !Kinetic limit check |
---|
493 | IF (satratln .GE. 1.E-2 .AND. satratln .LE. 1.) THEN |
---|
494 | kinrhotresn=EXP(7.8920778706888086E+1 + 7.3665492897447082*satratln - 1.2420166571163805E+4/tln & |
---|
495 | & + (-6.1831234251470971E+2*satratln)/tln - 2.4501159970109945E-2*tln & |
---|
496 | & -1.3463066443605762E-2*satratln*tln + 8.3736373989909194E-06*tln**2 & |
---|
497 | & -1.4673887785408892*LOG(satratln) + (-3.2141890006517094E+1*LOG(satratln))/tln & |
---|
498 | & + 2.7137429081917556E-3*tln*LOG(satratln)) !1/cm3 |
---|
499 | IF (kinrhotresn.LT.rhoaln) kinetic_n=.TRUE. |
---|
500 | ENDIF |
---|
501 | |
---|
502 | IF (satratln .GE. 1.E-4 .AND. satratln .LT. 1.E-2) THEN |
---|
503 | kinrhotresn=EXP(7.9074383049843647E+1 - 2.8746005462158347E+1*satratln - 1.2070272068458380E+4/tln & |
---|
504 | & + (-5.9205040320056632E+3*satratln)/tln - 2.4800372593452726E-2*tln & |
---|
505 | & -4.3983007681295948E-2*satratln*tln + 2.5943854791342071E-5*tln**2 & |
---|
506 | & -2.3141363245211317*LOG(satratln) + (9.9186787997857735E+1*LOG(satratln))/tln & |
---|
507 | & + 5.6819382556144681E-3*tln*LOG(satratln)) !1/cm3 |
---|
508 | IF (kinrhotresn.LT.rhoaln) kinetic_n=.TRUE. |
---|
509 | ENDIF |
---|
510 | |
---|
511 | IF (satratln .GE. 5.E-6 .AND. satratln .LT. 1.E-4) THEN |
---|
512 | kinrhotresn=EXP(8.5599712000361677E+1 + 2.7335119660796581E+3*satratln - 1.1842350246291651E+4/tln & |
---|
513 | & + (-1.2439843468881438E+6*satratln)/tln - 5.4536964974944230E-2*tln & |
---|
514 | & + 5.0886987425326087*satratln*tln + 7.1964722655507067E-5*tln**2 & |
---|
515 | & -2.4472627526306372*LOG(satratln) + (1.7561478001423779E+2*LOG(satratln))/tln & |
---|
516 | & + 6.2640132818141811E-3*tln*LOG(satratln)) !1/cm3 |
---|
517 | IF (kinrhotresn.LT.rhoaln) kinetic_n=.TRUE. |
---|
518 | ENDIF |
---|
519 | |
---|
520 | IF (kinetic_n) THEN |
---|
521 | ! Dimer formation rate |
---|
522 | jnuc_n=1.E6*(2.*0.3E-9)**2.*SQRT(8.*RPI*RKBOL*(1./mH2SO4mol+1./mH2SO4mol))/2.*SQRT(t)*rhoa**2. |
---|
523 | ! jnuc_n=1.E6*(2.*0.3E-9)**2.*SQRT(8.*3.141593*1.38E-23*(1./(1.661E-27*98.07)+1./(1.661E-27*98.07)))/2.*SQRT(t)*rhoa**2. |
---|
524 | ntot_n=1. !set to 1 |
---|
525 | na_n=1. ! The critical cluster contains one molecules but the produced cluster contains 2 molecules |
---|
526 | x_n=na_n/ntot_n ! so also set this to 1 |
---|
527 | rc_n=0.3E-9 |
---|
528 | ELSE |
---|
529 | jnuc_n= 2.1361182605986115E-1 + 3.3827029855551838*tln -3.2423555796175563E-2*tln**2 + & |
---|
530 | & 7.0120069477221989E-5*tln**3 +8.0286874752695141/x_n + & |
---|
531 | & (-2.6939840579762231E-1)*LOG(satratln) +1.6079879299099518*tln*LOG(satratln) + & |
---|
532 | & (-1.9667486968141933E-2)*tln**2*LOG(satratln) + & |
---|
533 | & 5.5244755979770844E-5*tln**3*LOG(satratln) + (7.8884704837892468*LOG(satratln))/x_n + & |
---|
534 | & 4.6374659198909596*LOG(satratln)**2 - 8.2002809894792153E-2*tln*LOG(satratln)**2 + & |
---|
535 | & 8.5077424451172196E-4*tln**2*LOG(satratln)**2 + & |
---|
536 | & (-2.6518510168987462E-6)*tln**3*LOG(satratln)**2 + & |
---|
537 | & (-1.4625482500575278*LOG(satratln)**2)/x_n - 5.2413002989192037E-1*LOG(satratln)**3 + & |
---|
538 | & 5.2755117653715865E-3*tln*LOG(satratln)**3 + & |
---|
539 | & (-2.9491061332113830E-6)*tln**2*LOG(satratln)**3 + & |
---|
540 | & (-2.4815454194486752E-8)*tln**3*LOG(satratln)**3 + & |
---|
541 | & (-5.2663760117394626E-2*LOG(satratln)**3)/x_n + & |
---|
542 | & 1.6496664658266762*LOG(rhoaln) + & |
---|
543 | & (-8.0809397859218401E-1)*tln*LOG(rhoaln) + & |
---|
544 | & 8.9302927091946642E-3*tln**2*LOG(rhoaln) + & |
---|
545 | & (-1.9583649496497497E-5)*tln**3*LOG(rhoaln) + & |
---|
546 | & (-8.9505572676891685*LOG(rhoaln))/x_n + & |
---|
547 | & (-3.0025283601622881E+1)*LOG(satratln)*LOG(rhoaln) + & |
---|
548 | & 3.0783365644763633E-1*tln*LOG(satratln)*LOG(rhoaln) + & |
---|
549 | & (-7.4521756337984706E-4)*tln**2*LOG(satratln)*LOG(rhoaln) + & |
---|
550 | & (-5.7651433870681853E-7)*tln**3*LOG(satratln)*LOG(rhoaln) + & |
---|
551 | & (1.2872868529673207*LOG(satratln)*LOG(rhoaln))/x_n + & |
---|
552 | & (-6.1739867501526535E-1)*LOG(satratln)**2*LOG(rhoaln) + & |
---|
553 | & 7.2347385705333975E-3*tln*LOG(satratln)**2*LOG(rhoaln) + & |
---|
554 | & (-3.0640494530822439E-5)*tln**2*LOG(satratln)**2*LOG(rhoaln) + & |
---|
555 | & 6.5944609194346214E-8*tln**3*LOG(satratln)**2*LOG(rhoaln) + & |
---|
556 | & (-2.8681650332461055E-2*LOG(satratln)**2*LOG(rhoaln))/x_n + & |
---|
557 | & 6.5213802375160306*LOG(rhoaln)**2 + & |
---|
558 | & (-4.7907162004793016E-2)*tln*LOG(rhoaln)**2 + & |
---|
559 | & (-1.0727890114215117E-4)*tln**2*LOG(rhoaln)**2 + & |
---|
560 | & 5.6401818280534507E-7*tln**3*LOG(rhoaln)**2 + & |
---|
561 | & (5.4113070888923009E-1*LOG(rhoaln)**2)/x_n + & |
---|
562 | & 5.2062808476476330E-1*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
563 | & (-6.0696882500824584E-3)*tln*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
564 | & 2.3851383302608477E-5*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
565 | & (-1.5243837103067096E-8)*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
566 | & (-5.6543192378015687E-2*LOG(satratln)*LOG(rhoaln)**2)/x_n + & |
---|
567 | & (-1.1630806410696815E-1)*LOG(rhoaln)**3 + & |
---|
568 | & 1.3806404273119610E-3*tln*LOG(rhoaln)**3 + & |
---|
569 | & (-2.0199865087650833E-6)*tln**2*LOG(rhoaln)**3 + & |
---|
570 | & (-3.0200284885763192E-9)*tln**3*LOG(rhoaln)**3 + & |
---|
571 | & (-6.9425267104126316E-3*LOG(rhoaln)**3)/x_n |
---|
572 | jnuc_n=EXP(jnuc_n) |
---|
573 | |
---|
574 | ntot_n =-3.5863435141979573E-3 - 1.0098670235841110E-1*tln + 8.9741268319259721E-4*tln**2 - 1.4855098605195757E-6*tln**3 & |
---|
575 | & - 1.2080330016937095E-1/x_n + 1.1902674923928015E-3*LOG(satratln) - 1.9211358507172177E-2*tln*LOG(satratln) + & |
---|
576 | & 2.4648094311204255E-4*tln**2*LOG(satratln) - 7.5641448594711666E-7*tln**3*LOG(satratln) + & |
---|
577 | & (-2.0668639384228818E-02*LOG(satratln))/x_n - 3.7593072011595188E-2*LOG(satratln)**2 + & |
---|
578 | & 9.0993182774415718E-4 *tln*LOG(satratln)**2 + & |
---|
579 | & (-9.5698412164297149E-6)*tln**2*LOG(satratln)**2 + 3.7163166416110421E-8*tln**3*LOG(satratln)**2 + & |
---|
580 | & (1.1026579525210847E-2*LOG(satratln)**2)/x_n + 1.1530844115561925E-2 *LOG(satratln)**3 + & |
---|
581 | & (-1.8083253906466668E-4)*tln*LOG(satratln)**3 + 8.0213604053330654E-7*tln**2*LOG(satratln)**3 + & |
---|
582 | & (-8.5797885383051337E-10)*tln**3*LOG(satratln)**3 + (1.0243693899717402E-3*LOG(satratln)**3)/x_n + & |
---|
583 | & (-1.7248695296299649E-2)*LOG(rhoaln) + 1.1294004162437157E-2*tln*LOG(rhoaln) + & |
---|
584 | & (-1.2283640163189278E-4)*tln**2*LOG(rhoaln) + 2.7391732258259009E-7*tln**3*LOG(rhoaln) + & |
---|
585 | & (6.8505583974029602E-2*LOG(rhoaln))/x_n +2.9750968179523635E-1*LOG(satratln)*LOG(rhoaln) + & |
---|
586 | & (-3.6681154503992296E-3)*tln*LOG(satratln)*LOG(rhoaln) + 1.0636473034653114E-5*tln**2*LOG(satratln)*LOG(rhoaln)+ & |
---|
587 | & 5.8687098466515866E-9*tln**3*LOG(satratln)*LOG(rhoaln) + (-5.2028866094191509E-3*LOG(satratln)*LOG(rhoaln))/x_n+ & |
---|
588 | & 7.6971988880587231E-4*LOG(satratln)**2*LOG(rhoaln) - 2.4605575820433763E-5*tln*LOG(satratln)**2*LOG(rhoaln) + & |
---|
589 | & 2.3818484400893008E-7*tln**2*LOG(satratln)**2*LOG(rhoaln) + & |
---|
590 | & (-8.8474102392445200E-10)*tln**3*LOG(satratln)**2*LOG(rhoaln) + & |
---|
591 | & (-1.6640566678168968E-4*LOG(satratln)**2*LOG(rhoaln))/x_n - 7.7390093776705471E-2*LOG(rhoaln)**2 + & |
---|
592 | & 5.8220163188828482E-4*tln*LOG(rhoaln)**2 + 1.2291679321523287E-6*tln**2*LOG(rhoaln)**2 + & |
---|
593 | & (-7.4690997508075749E-9)*tln**3*LOG(rhoaln)**2 + (-5.6357941220497648E-3*LOG(rhoaln)**2)/x_n + & |
---|
594 | & (-4.7170109625089768E-3)*LOG(satratln)*LOG(rhoaln)**2 + 6.9828868534370193E-5*tln*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
595 | & (-3.1738912157036403E-7)*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
596 | & 2.3975538706787416E-10*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & |
---|
597 | & (4.2304213386288567E-4*LOG(satratln)*LOG(rhoaln)**2)/x_n + 1.3696520973423231E-3*LOG(rhoaln)**3 + & |
---|
598 | & (-1.6863387574788199E-5)*tln*LOG(rhoaln)**3 + 2.7959499278844516E-8*tln**2*LOG(rhoaln)**3 + & |
---|
599 | & 3.9423927013227455E-11*tln**3*LOG(rhoaln)**3 + (8.6136359966337272E-5*LOG(rhoaln)**3)/x_n |
---|
600 | ntot_n=EXP(ntot_n) |
---|
601 | |
---|
602 | rc_n=EXP(-22.378268374023630+0.44462953606125100*x_n+0.33499495707849131*LOG(ntot_n)) !in meters |
---|
603 | |
---|
604 | na_n=x_n*ntot_n |
---|
605 | IF (na_n .LT. 1.) THEN |
---|
606 | print *, 'Warning: number of acid molecules < 1 in nucleation regime, setting na_n=1' |
---|
607 | na_n=1.0 |
---|
608 | ENDIF |
---|
609 | ENDIF |
---|
610 | |
---|
611 | ! Set the neutral nucleation rate to 0.0 if less than 1.0E-7 |
---|
612 | IF (jnuc_n.LT.1.E-7) THEN |
---|
613 | jnuc_n=0.0 |
---|
614 | ENDIF |
---|
615 | |
---|
616 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
617 | |
---|
618 | ! Threshold neutral nucleation rate (j > 1/cm3s) parameterization (can be commented out if not needed) |
---|
619 | IF (tln .GE. 310.) THEN |
---|
620 | rhoatres=EXP(-2.8220714121794250 + 1.1492362322651116E+1*satratln -3.3034839106184218E+3/tln & |
---|
621 | & + (-7.1828571490168133E+2*satratln)/tln + 1.4649510835204091E-1*tln & |
---|
622 | & -3.0442736551916524E-2*satratln*tln -9.3258567137451497E-5*tln**2 & |
---|
623 | & -1.1583992506895649E+1*LOG(satratln) + (1.5184848765906165E+3*LOG(satratln))/tln & |
---|
624 | & + 1.8144983916747057E-2*tln*LOG(satratln)) !1/cm3 |
---|
625 | ENDIF |
---|
626 | |
---|
627 | IF (tln .GT. 190. .AND. tln .LT. 310.) THEN |
---|
628 | rhoatres=EXP(-3.1820396091231999E+2 + 7.2451289153199676*satratln + 2.6729355170089486E+4/tln & |
---|
629 | & + (-7.1492506076423069E+2*satratln)/tln + 1.2617291148391978*tln & |
---|
630 | & - 1.6438112080468487E-2*satratln*tln -1.4185518234553220E-3*tln**2 & |
---|
631 | & -9.2864597847386694*LOG(satratln) + (1.2607421852455602E+3*LOG(satratln))/tln & |
---|
632 | & + 1.3324434472218746E-2*tln*LOG(satratln)) !1/cm3 |
---|
633 | ENDIF |
---|
634 | |
---|
635 | IF (tln .LT. 185. .AND. tln .GT. 155.) THEN |
---|
636 | rhoatres=1.1788859232398459E+5 - 1.0244255702550814E+4*satratln + & |
---|
637 | & 4.6815029684321962E+3*satratln**2 -1.6755952338499657E+2*tln |
---|
638 | ENDIF |
---|
639 | |
---|
640 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
641 | |
---|
642 | ! Ion-induced nucleation parameterization |
---|
643 | |
---|
644 | IF (ipr.GT.0.0) THEN ! if the ion production rate is above zero |
---|
645 | |
---|
646 | ! Calculate the ion induced nucleation rate wrt. concentration of 1 ion/cm3 |
---|
647 | |
---|
648 | kinrhotresi = 5.3742280876674478E1 - 6.6837931590012266E-3 *LOG(satratli)**(-2) & |
---|
649 | & - 1.0142598385422842E-01 * LOG(satratli)**(-1) - 6.4170597272606873E+00 * LOG(satratli) & |
---|
650 | & - 6.4315798914824518E-01 * LOG(satratli)**2 - 2.4428391714772721E-02 * LOG(satratli)**3 & |
---|
651 | & - 3.5356658734539019E-04 * LOG(satratli)**4 + 2.5400015099140506E-05 * tli * LOG(satratli)**(-2) & |
---|
652 | & - 2.7928900816637790E-04 * tli * LOG(satratli)**(-1) + 4.4108573484923690E-02 * tli * LOG(satratli) & |
---|
653 | & + 6.3943789012475532E-03 * tli * LOG(satratli)**(2) + 2.3164296174966580E-04 * tli * LOG(satratli)**(3) & |
---|
654 | & + 3.0372070669934950E-06 * tli * LOG(satratli)**4 + 3.8255873977423475E-06 * tli**2 * LOG(satratli)**(-1) & |
---|
655 | & - 1.2344793083561629E-04 * tli**2 * LOG(satratli) - 1.7959048869810192E-05 * tli**2 * LOG(satratli)**(2) & |
---|
656 | & - 3.2165622558722767E-07 * tli**2 * LOG(satratli)**3 - 4.7136923780988659E-09 * tli**3 * LOG(satratli)**(-1) & |
---|
657 | & + 1.1873317184482216E-07 * tli**3 * LOG(satratli) + 1.5685860354866621E-08 * tli**3 * LOG(satratli)**2 & |
---|
658 | & - 1.4329645891059557E+04 * tli**(-1) + 1.3842599842575321E-01 * tli & |
---|
659 | & - 4.1376265912842938E-04 * tli**(2) + 3.9147639775826004E-07 * tli**3 |
---|
660 | |
---|
661 | kinrhotresi=EXP(kinrhotresi) !1/cm3 |
---|
662 | |
---|
663 | IF (kinrhotresi.LT.rhoali) kinetic_i=.true. |
---|
664 | |
---|
665 | IF (kinetic_i) THEN |
---|
666 | ! jnuc_i1=1.0E6*(0.3E-9 + 0.487E-9)**2.*SQRT(8.*3.141593*1.38E-23*(1./(1.661E-27*98.07)+1./(1.661E-27*98.07)))* & |
---|
667 | jnuc_i1=1.0E6*(0.3E-9 + 0.487E-9)**2.*SQRT(8.*RPI*RKBOL*(1./mH2SO4mol+1./mH2SO4mol))* & |
---|
668 | & SQRT(tli)*rhoali !1/cm3s |
---|
669 | ntot_i=1. !set to 1 |
---|
670 | na_i=1. |
---|
671 | x_i=na_i/ntot_i ! so also set this to 1 |
---|
672 | rc_i=0.487E-9 |
---|
673 | ELSE |
---|
674 | jnuc_i1 = 3.0108954259038608E+01+tli*6.1176722090512577E+01+(tli**2)*8.7240333618891663E-01+(tli**3)* & |
---|
675 | & (-4.6191788649375719E-03)+(tli**(-1))*8.3537059107024481E-01 + & |
---|
676 | & (1.5028549216690628E+01+tli*(-1.9310989753720623E-01)+(tli**2)*8.0155514634860480E-04+(tli**3)* & |
---|
677 | & (-1.0832730707799128E-06)+(tli**(-1))*1.7577660457989019)*(LOG(satratli)**(-2)) + & |
---|
678 | & (-2.0487870170216488E-01 + tli * 1.3263949252910405E-03 + (tli**2) *(-8.4195688402450274E-06) + & |
---|
679 | & (tli**3)*1.6154895940993287E-08 + (tli**(-1))*3.8734212545203874E+01) * (LOG(satratli)**(-2)*LOG(rhoali)) + & |
---|
680 | & (1.4955918863858371 + tli * 9.2290004245522454E+01 + (tli**2) *(-8.9006965195392618E-01) + & |
---|
681 | & (tli**3) * 2.2319123411013099E-03 + (tli**(-1)) * 4.0180079996840852E-03) * & |
---|
682 | & (LOG(satratli)**(-1) * LOG(rhoali)**(-1)) + & |
---|
683 | & (7.9018031228561085 + tli *(-1.1649433968658949E+01) + (tli**2) * 1.1400827854910951E-01 + & |
---|
684 | & (tli**3) *(-3.1941526492127755E-04) + (tli**(-1)) *(-3.7662115740271446E-01)) * (LOG(satratli)**(-1)) + & |
---|
685 | & (1.5725237111225979E+02 + tli *(-1.0051649979836277) + (tli**2) * 1.1866484014507624E-03 + & |
---|
686 | & (tli**3) * 7.3557614998540389E-06 + (tli**(-1)) * 2.6270197023115189) * (LOG(satratli)**(-1) * LOG(rhoali)) + & |
---|
687 | & (-1.6973840122470968E+01 + tli * 1.1258423691432135E-01 + (tli**2) *(-2.9850139351463793E-04) + (tli**3) * & |
---|
688 | & 1.4301286324827064E-07 + (tli**(-1)) * 1.3163389235253725E+01) * (LOG(satratli)**(-1) * LOG(rhoali)**2) + & |
---|
689 | & (-1.0399591631839757 + tli * 2.7022055588257691E-03 + (tli**2) *(-2.1507467231330936E-06) + (tli**3) * & |
---|
690 | & 3.8059489037584171E-10 + (tli**(-1)) * 1.5000492788553410E+02) * (LOG(satratli)**(-1) * LOG(rhoali)**3) + & |
---|
691 | & (1.2250990965305315 + tli * 3.0495946490079444E+01 + (tli**2) * 2.1051563135187106E+01 + (tli**3) * & |
---|
692 | & (-8.2200682916580878E-02) + (tli**(-1)) * 2.9965871386685029E-02) * (LOG(rhoali)**(-2)) + & |
---|
693 | & (4.8281605955680433 + tli * 1.7346551710836445E+02 + (tli**2) *(-1.0113602140796010E+01) + (tli**3) * & |
---|
694 | & 3.7482518458685089E-02 + (tli**(-1)) *(-1.4449998158558205E-01)) * (LOG(rhoali)**(-1)) + & |
---|
695 | & (2.3399230964451237E+02 + tli *(-2.3099267235261948E+01) + (tli**2) * 8.0122962140916354E-02 + & |
---|
696 | & (tli**3) * 6.1542576994557088E-05 + (tli**(-1)) * 5.3718413254843007) * (LOG(rhoali)) + & |
---|
697 | & (1.0299715519499360E+02 + tli *(-6.4663357203364136E-02) + (tli**2) *(-2.0487150565050316E-03) + & |
---|
698 | & (tli**3) * 8.7935289055530897E-07 + (tli**(-1)) * 3.6013204601215229E+01) * (LOG(rhoali)**2) + & |
---|
699 | & (-3.5452115439584042 + tli * 1.7083445731159330E-02 + (tli**2) *(-1.2552625290862626E-05) + (tli**3) * & |
---|
700 | & 1.2968447449182847E-09 + (tli**(-1)) * 1.5748687512056560E+02) * (LOG(rhoali)**3) + & |
---|
701 | & (2.2338490119517975 + tli * 1.0229410216045540E+02 + (tli**2) *(-3.2103611955174052) + (tli**3) * & |
---|
702 | & 1.3397152304977591E-02 + (tli**(-1)) *(-2.4155187776460030E-02)) * (LOG(satratli)* LOG(rhoali)**(-2)) + & |
---|
703 | & (3.7592282990713963 + tli *(-1.5257988769009816E+02) + (tli**2) * 2.6113805420558802 + (tli**3) * & |
---|
704 | & (-9.0380721653694363E-03) + (tli**(-1)) *(-1.3974197138171082E-01)) * (LOG(satratli)* LOG(rhoali)**(-1)) + & |
---|
705 | & (1.8293600730573988E+01 + tli * 1.8344728606002992E+01 + (tli**2) *(-4.0063363221106751E-01) + (tli**3) & |
---|
706 | & * 1.4842749371258522E-03 + (tli**(-1)) * 1.1848846003282287) * (LOG(satratli)) + & |
---|
707 | & (-1.7634531623032314E+02 + tli * 4.9011762441271278 + (tli**2) *(-1.3195821562746339E-02) + (tli**3) * & |
---|
708 | & (-2.8668619526430859E-05) + (tli**(-1)) *(-2.9823396976393551E-01)) * (LOG(satratli)* LOG(rhoali)) + & |
---|
709 | & (-3.2944043694275727E+01 + tli * 1.2517571921051887E-01 + (tli**2) * 8.3239769771186714E-05 + (tli**3) * & |
---|
710 | & 2.8191859341519507E-07 + (tli**(-1)) *(-2.7352880736682319E+01)) * (LOG(satratli)* LOG(rhoali)**2) + & |
---|
711 | & (-1.1451811137553243 + tli * 2.0625997485732494E-03 + (tli**2) *(-3.4225389469233624E-06) + (tli**3) * & |
---|
712 | & 4.4437613496984567E-10 + (tli**(-1)) * 1.8666644332606754E+02) * (LOG(satratli)* LOG(rhoali)**3) + & |
---|
713 | & (3.2270897099493567E+01 + tli * 7.7898447327513687E-01 + (tli**2) *(-6.5662738484679626E-03) + (tli**3) * & |
---|
714 | & 3.7899330796456790E-06 + (tli**(-1)) * 7.1106427501756542E-01) * (LOG(satratli)**2 * LOG(rhoali)**(-1)) + & |
---|
715 | & (-2.8901906781697811E+01 + tli *(-1.5356398793054860) + (tli**2) * 1.9267271774384788E-02 + (tli**3) * & |
---|
716 | & (-5.3886270475516162E-05) + (tli**(-1)) * 5.0490415975693426E-01) * (LOG(satratli)**2) + & |
---|
717 | & (3.3365683645733924E+01 + tli *(-3.6114561564894537E-01) + (tli**2) * 9.2977354471929262E-04 + (tli**3) * & |
---|
718 | & 1.9549769069511355E-07 + (tli**(-1)) *(-8.8865930095112855)) * (LOG(satratli)**2 * LOG(rhoali)) + & |
---|
719 | & (2.4592563042806375 + tli *(-8.3227071743101084E-03) + (tli**2) * 8.2563338043447783E-06 + (tli**3) * & |
---|
720 | & (-8.4374976698593496E-09) + (tli**(-1)) *(-2.0938173949893473E+02)) * (LOG(satratli)**2 * LOG(rhoali)**2) + & |
---|
721 | & (4.4099823444352317E+01 + tli * 2.5915665826835252 + (tli**2) *(-1.6449091819482634E-02) + (tli**3) * & |
---|
722 | & 2.6797249816144721E-05 + (tli**(-1)) * 5.5045672663909995E-01)* satratli |
---|
723 | jnuc_i1=EXP(jnuc_i1) |
---|
724 | |
---|
725 | ntot_i = ABS((-4.8324296064013375E+04 + tli * 5.0469120697428906E+02 + (tli**2) *(-1.1528940488496042E+00) + & |
---|
726 | & (tli**(-1)) *(-8.6892744676239192E+02) + (tli**(3)) * 4.0030302028120469E-04) + & |
---|
727 | & (-6.7259105232039847E+03 + tli * 1.9197488157452008E+02 + (tli**2) *(-1.3602976930126354E+00) + & |
---|
728 | & (tli**(-1)) *(-1.1212637938360332E+02) + (tli**(3)) * 2.8515597265933207E-03) * & |
---|
729 | & LOG(satratli)**(-2) * LOG(rhoali)**(-2) + & |
---|
730 | & (2.6216455217763342E+02 + tli *(-2.3687553252750821E+00) + (tli**2) * 7.4074554767517521E-03 + & |
---|
731 | & (tli**(-1)) *(-1.9213956820114927E+03) + (tli**(3)) *(-9.3839114856129453E-06)) * LOG(satratli)**(-2) + & |
---|
732 | & (3.9652478944137344E+00 + tli * 1.2469375098256536E-02 + (tli**2) *(-9.9837754694045633E-05) + (tli**(-1)) * & |
---|
733 | & (-5.1919499210175138E+02) + (tli**(3)) * 1.6489001324583862E-07) * LOG(satratli)**(-2) * LOG(rhoali) + & |
---|
734 | & (2.4975714429096206E+02 + tli * 1.7107594562445172E+02 + (tli**2) *(-7.8988711365135289E-01) + (tli**(-1)) * & |
---|
735 | & (-2.2243599782483177E+01) + (tli**(3)) *(-1.6291523004095427E-04)) * LOG(satratli)**(-1) * LOG(rhoali)**(-2) +& |
---|
736 | & (-8.9270715592533611E+02 + tli * 1.2053538883338946E+02 + (tli**2) *(-1.5490408828541018E+00) + (tli**(-1)) * & |
---|
737 | & (-1.1243275579419826E+01) + (tli**(3)) * 4.8053105606904655E-03) * LOG(satratli)**(-1) * LOG(rhoali)**(-1) + & |
---|
738 | & (7.6426441642091631E+03 + tli *(-7.1785462414656578E+01) + (tli**2) * 2.3851864923199523E-01 + (tli**(-1)) * & |
---|
739 | & 8.5591775688708395E+01 + (tli**(3)) *(-3.7000473243342858E-04)) * LOG(satratli)**(-1) + & |
---|
740 | & (-5.1516826398607911E+01 + tli * 9.1385720811460558E-01 + (tli**2) *(-3.5477100262158974E-03) + & |
---|
741 | & (tli**(-1)) * 2.7545544507625586E+03 + (tli**(3)) * 5.4708262093640928E-06) * LOG(satratli)**(-1) * LOG(rhoali) + & |
---|
742 | & (-3.0386767129196176E+02 + tli *(-1.1033438883583569E+04) + (tli**2) * 8.1296859732896067E+01 + (tli**(-1)) * & |
---|
743 | & 1.2625883141097162E+01 + (tli**(3)) *(-1.2728497822219101E-01)) * LOG(rhoali)**(-2) + & |
---|
744 | & (-3.3763494256461472E+03 + tli * 3.1916579136391006E+03 + (tli**2) *(-2.7234339474441143E+01) + (tli**(-1)) * & |
---|
745 | & (-2.1897653262707397E+01) + (tli**(3)) * 5.1788505812259071E-02) * LOG(rhoali)**(-1) + & |
---|
746 | & (-1.8817843873687068E+03 + tli * 4.3038072285882070E+00 + (tli**2) * 6.6244087689671860E-03 + (tli**(-1)) * & |
---|
747 | & (-2.7133073605696295E+03) + (tli**(3)) *(-1.7951557394285043E-05)) * LOG(rhoali) + & |
---|
748 | & (-1.7668827539244447E+02 + tli * 4.8160932330629913E-01 + (tli**2) *(-6.3133007671100293E-04) + (tli**(-1)) * & |
---|
749 | & 2.5631774669873157E+04 + (tli**(3)) * 4.1534484127873519E-07) * LOG(rhoali)**(2) + & |
---|
750 | & (-1.6661835889222382E+03 + tli * 1.3708900504682877E+03 + (tli**2) *(-1.7919060052198969E+01) + (tli**(-1)) * & |
---|
751 | & (-3.5145029804436405E+01) + (tli**(3)) * 5.1047240947371224E-02) * LOG(satratli)* LOG(rhoali)**(-2) + & |
---|
752 | & (1.0843549363030939E+04 + tli *(-7.3557073636139577E+01) + (tli**2) * 1.2054625131778862E+00 + (tli**(-1)) * & |
---|
753 | & 1.9358737917864391E+02 + (tli**(3)) *(-4.2871620775911338E-03)) * LOG(satratli)* LOG(rhoali)**(-1) + & |
---|
754 | & (-2.4269802549752835E+03 + tli * 1.1348265061941714E+01 + (tli**2) *(-5.0430423939495157E-02) + (tli**(-1)) * & |
---|
755 | & 2.3709874548950634E+03 + (tli**(3)) * 1.4091851828620244E-04) * LOG(satratli) + & |
---|
756 | & (5.2745372575251588E+02 + tli *(-2.6080675912627314E+00) + (tli**2) * 5.6902218056670145E-03 + (tli**(-1)) * & |
---|
757 | & (-3.2149319482897838E+04) + (tli**(3)) *(-5.4121996056745853E-06)) * LOG(satratli)* LOG(rhoali) + & |
---|
758 | & (-1.6401959518360403E+01 + tli * 2.4322962162439640E-01 + (tli**2) * 1.1744366627725344E-03 + (tli**(-1)) * & |
---|
759 | & (-8.2694427518413195E+03) + (tli**(3)) *(-5.0028379203873102E-06))* LOG(satratli)**(2) + & |
---|
760 | & (-2.7556572017167782E+03 + tli * 4.9293344495058264E+01 + (tli**2) *(-2.6503456520676050E-01) + (tli**(-1)) * & |
---|
761 | & 1.2130698030982167E+03 + (tli**(3)) * 4.3530610668042957E-04)* LOG(satratli)**2 * LOG(rhoali)**(-1) + & |
---|
762 | & (-6.3419182228959192E+00 + tli * 4.0636212834605827E-02 + (tli**2) *(-1.0450112687842742E-04) + (tli**(-1)) * & |
---|
763 | & 3.1035882189759656E+02 + (tli**(3)) * 9.4328418657873500E-08)* LOG(satratli)**(-3) + & |
---|
764 | & (3.0189213304689042E+03 + tli *(-2.3804654203861684E+01) + (tli**2) * 6.8113013411972942E-02 + (tli**(-1)) * & |
---|
765 | & 6.3112071081188913E+02 + (tli**(3)) *(-9.4460854261685723E-05))* (satratli) * LOG(rhoali) + & |
---|
766 | & (1.1924791930673702E+04 + tli *(-1.1973824959206000E+02) + (tli**2) * 1.6888713097971020E-01 + (tli**(-1)) * & |
---|
767 | & 1.8735938211539585E+02 + (tli**(3)) * 5.0974564680442852E-04)* (satratli) + & |
---|
768 | & (3.6409071302482083E+01 + tli * 1.7919859306449623E-01 + (tli**2) *(-1.0020116255895206E-03) + (tli**(-1)) * & |
---|
769 | & (-8.3521083354432303E+03) + (tli**(3)) * 1.5879900546795635E-06)* satratli * LOG(rhoali)**(2)) |
---|
770 | |
---|
771 | rc_i = (-3.6318550637865524E-08 + tli * 2.1740704135789128E-09 + (tli**2) * & |
---|
772 | & (-8.5521429066506161E-12) + (tli**3) *(-9.3538647454573390E-15)) + & |
---|
773 | & (2.1366936839394922E-08 + tli *(-2.4087168827395623E-10) + (tli**2) * 8.7969869277074319E-13 + & |
---|
774 | & (tli**3) *(-1.0294466881303291E-15))* LOG(satratli)**(-2) * LOG(rhoali)**(-1) + & |
---|
775 | & (-7.7804007761164303E-10 + tli * 1.0327058173517932E-11 + (tli**2) *(-4.2557697639692428E-14) + & |
---|
776 | & (tli**3) * 5.4082507061618662E-17)* LOG(satratli)**(-2) + & |
---|
777 | & (3.2628927397420860E-12 + tli *(-7.6475692919751066E-14) + (tli**2) * 4.1985816845259788E-16 + & |
---|
778 | & (tli**3) *(-6.2281395889592719E-19))* LOG(satratli)**(-2) * LOG(rhoali) + & |
---|
779 | & (2.0442205540818555E-09 + tli * 4.0441858911249830E-08 + (tli**2) *(-3.3423487629482825E-10) + & |
---|
780 | & (tli**3) * 6.8000404742985678E-13)* LOG(satratli)**(-1) * LOG(rhoali)**(-2) + & |
---|
781 | & (1.8381489183824627E-08 + tli *(-8.9853322951518919E-09) + (tli**2) * 7.5888799566036185E-11 + & |
---|
782 | & (tli**3) *(-1.5823457864755549E-13))* LOG(satratli)**(-1) * LOG(rhoali)**(-1) + & |
---|
783 | & (1.1795760639695057E-07 + tli *(-8.1046722896375875E-10) + (tli**2) * 9.1868604369041857E-14 + & |
---|
784 | & (tli**3) * 4.7882428237444610E-15)* LOG(satratli)**(-1) + & |
---|
785 | & (-4.4028846582545952E-09 + tli * 4.6541269232626618E-11 + (tli**2) *(-1.1939929984285194E-13) + & |
---|
786 | & (tli**3) * 2.3602037016614437E-17)* LOG(satratli)**(-1) * LOG(rhoali) + & |
---|
787 | & (2.7885056884209128E-11 + tli *(-4.5167129624119121E-13) + (tli**2) * 1.6558404997394422E-15 + & |
---|
788 | & (tli**3) *(-1.2037336621218054E-18))* LOG(satratli)**(-1) * LOG(rhoali)**2 + & |
---|
789 | & (-2.3719627171699983E-09 + tli *(-1.5260127909292053E-07) + (tli**2) * 1.7177017944754134E-09 + & |
---|
790 | & (tli**3) *(-4.7031737537526395E-12))* LOG(rhoali)**(-2) + & |
---|
791 | & (-5.6946433724699646E-09 + tli * 8.4629788237081735E-09 + (tli**2) *(-1.7674135187061521E-10) + & |
---|
792 | & (tli**3) * 6.6236547903091862E-13)* LOG(rhoali)**(-1) + & |
---|
793 | & (-2.2808617930606012E-08 + tli * 1.4773376696847775E-10 + (tli**2) *(-1.3076953119957355E-13) + & |
---|
794 | & (tli**3) * 2.3625301497914000E-16)* LOG(rhoali) + & |
---|
795 | & (1.4014269939947841E-10 + tli *(-2.3675117757377632E-12) + (tli**2) * 5.1514033966707879E-15 + & |
---|
796 | & (tli**3) *(-4.8864233454747856E-18))* LOG(rhoali)**2 + & |
---|
797 | & (6.5464943868885886E-11 + tli * 1.6494354816942769E-08 + (tli**2) *(-1.7480097393483653E-10) + & |
---|
798 | & (tli**3) * 4.7460075628523984E-13)* LOG(satratli)* LOG(rhoali)**(-2) + & |
---|
799 | & (8.4737893183927871E-09 + tli *(-6.0243327445597118E-09) + (tli**2) * 5.8766070529814883E-11 + & |
---|
800 | & (tli**3) *(-1.4926748560042018E-13))* LOG(satratli)* LOG(rhoali)**(-1) + & |
---|
801 | & (1.0761964135701397E-07 + tli *(-1.0142496009071148E-09) + (tli**2) * 2.1337312466519190E-12 + & |
---|
802 | & (tli**3) * 1.6376014957685404E-15)* LOG(satratli) + & |
---|
803 | & (-3.5621571395968670E-09 + tli * 4.1175339587760905E-11 + (tli**2) *(-1.3535372357998504E-13) + & |
---|
804 | & (tli**3) * 8.9334219536920720E-17)* LOG(satratli)* LOG(rhoali) + & |
---|
805 | & (2.0700482083136289E-11 + tli *(-3.9238944562717421E-13) + (tli**2) * 1.5850961422040196E-15 + & |
---|
806 | & (tli**3) *(-1.5336775610911665E-18))* LOG(satratli)* LOG(rhoali)**2 + & |
---|
807 | & (1.8524255464416206E-09 + tli *(-2.1959816152743264E-11) + (tli**2) *(-6.4478119501677012E-14) + & |
---|
808 | & (tli**3) * 5.5135243833766056E-16)* LOG(satratli)**2 * LOG(rhoali)**(-1) + & |
---|
809 | & (1.9349488650922679E-09 + tli *(-2.2647295919976428E-11) + (tli**2) * 9.2917479748268751E-14 + & |
---|
810 | & (tli**3) *(-1.2741959892173170E-16))* LOG(satratli)**2 + & |
---|
811 | & (2.1484978031650972E-11 + tli *(-9.3976642475838013E-14) + (tli**2) *(-4.8892738002751923E-16) + & |
---|
812 | & (tli**3) * 1.4676120441783832E-18)* LOG(satratli)**2 * LOG(rhoali) + & |
---|
813 | & (6.7565715216420310E-13 + tli *(-3.5421162549480807E-15) + (tli**2) *(-3.4201196868693569E-18) + & |
---|
814 | & (tli**3) * 2.2260187650412392E-20)* LOG(satratli)**3 * LOG(rhoali) |
---|
815 | |
---|
816 | na_i=x_i*ntot_i |
---|
817 | IF (na_i .LT. 1.) THEN |
---|
818 | ! print *, 'Warning: number of acid molecules < 1 in nucleation regime, setting na_n=1' |
---|
819 | na_n=1.0 |
---|
820 | ENDIF |
---|
821 | ENDIF |
---|
822 | |
---|
823 | jnuc_i=jnuc_i1 |
---|
824 | ! Ion loss rate (1/s) |
---|
825 | xloss=csi+jnuc_i |
---|
826 | |
---|
827 | ! Recombination (here following Brasseur and Chatel, 1983) |
---|
828 | recomb=6.0E-8*SQRT(300./tli)+6.0E-26*airn*(300./tli)**4 |
---|
829 | |
---|
830 | ! Small ion concentration in air (1/cm3) (following Dunne et al., 2016) |
---|
831 | ! max function is to avoid n_i to go practically zero at very high J_ion |
---|
832 | n_i=MAX(0.01,(SQRT(xloss**2.0+4.0*recomb*ipr)-xloss)/(2.0*recomb)) |
---|
833 | |
---|
834 | ! Ion-induced nucleation rate |
---|
835 | ! Min function is to ensure that max function above does not cause J_ion to overshoot |
---|
836 | jnuc_i=MIN(ipr,n_i*jnuc_i1) |
---|
837 | ! Set the ion-induced nucleation rate to 0.0 if less than 1.0E-7 |
---|
838 | IF (jnuc_i.LT.1.E-7) THEN |
---|
839 | jnuc_i=0.0 |
---|
840 | ENDIF |
---|
841 | |
---|
842 | ENDIF |
---|
843 | |
---|
844 | !--conversion from double precision to float in case the model is run in single precision |
---|
845 | jnuc_n_real = jnuc_n |
---|
846 | ntot_n_real = ntot_n |
---|
847 | jnuc_i_real = jnuc_i |
---|
848 | ntot_i_real = ntot_i |
---|
849 | x_n_real = x_n |
---|
850 | x_i_real = x_i |
---|
851 | na_n_real = na_n |
---|
852 | na_i_real = na_i |
---|
853 | rc_n_real = rc_n |
---|
854 | rc_i_real = rc_i |
---|
855 | n_i_real = n_i |
---|
856 | rhoatres_real = rhoatres |
---|
857 | |
---|
858 | END SUBROUTINE newbinapara |
---|
859 | |
---|
860 | END MODULE nucleation_tstep_mod |
---|