1 | PROGRAM ammonium_sulfate |
---|
2 | !--based on Tang and Munkelwitz JGR 1994 and Tang JGR 1997 |
---|
3 | IMPLICIT NONE |
---|
4 | REAL xms !--solute mass fraction |
---|
5 | REAL aw !-water activity (RH here) |
---|
6 | REAL molality !-molality |
---|
7 | REAL rhop !-dry droplet density |
---|
8 | REAL nr !--real refractive index |
---|
9 | REAL RR, Rwater, Ramm !--molar refraction |
---|
10 | PARAMETER (Rwater=3.717, Ramm=23.50) |
---|
11 | REAL VV !--molal volume |
---|
12 | REAL Mwater, Mamm |
---|
13 | PARAMETER (Mwater=18.02, Mamm=132.14) |
---|
14 | REAL y !--mole fraction |
---|
15 | PARAMETER (rhop=1.769) !--g cm-3 |
---|
16 | REAL B0, B1, B2, B3, B4 |
---|
17 | PARAMETER (B0=110.65495, B1=-367.59197,B2=504.62934) |
---|
18 | PARAMETER (B3=-315.43839,B4=67.70824) |
---|
19 | REAL C1, C2, C3, C4 |
---|
20 | PARAMETER (C1=-2.715E-3,C2=3.113E-5,C3=-2.336E-6,C4=1.412E-8) |
---|
21 | REAL A0, A1, A2, A3 |
---|
22 | PARAMETER (A0=0.9971, A1=5.92E-3, A2=-5.036E-6, A3=1.024E-8) |
---|
23 | c------------------------- |
---|
24 | INTEGER nbre_rh,IRH,wve |
---|
25 | INTEGER nbre_dry,ii |
---|
26 | PARAMETER(nbre_rh=12,nbre_dry=4) |
---|
27 | REAL DELTA |
---|
28 | REAL RH_tab(nbre_rh) |
---|
29 | DATA RH_tab/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./ |
---|
30 | INTEGER Nwvmax |
---|
31 | PARAMETER(Nwvmax=100) |
---|
32 | c |
---|
33 | c--wavelength in m |
---|
34 | REAL lambda_min, lambda_max, lambda(Nwvmax) |
---|
35 | PARAMETER (lambda_min=0.2E-6, lambda_max=5.0E-6) |
---|
36 | |
---|
37 | REAL n_rr(nbre_rh,Nwvmax) |
---|
38 | REAL n_ii(nbre_rh,Nwvmax) |
---|
39 | REAL ratio(nbre_rh) !--diameter growth factor |
---|
40 | REAL rhod(nbre_rh) !--density |
---|
41 | REAL n_r_exact(nbre_rh) !--refractive index |
---|
42 | c |
---|
43 | DO wve=1,Nwvmax |
---|
44 | lambda(wve)= |
---|
45 | . lambda_min+FLOAT(wve-1)*(lambda_max-lambda_min)/FLOAT(Nwvmax-1) |
---|
46 | ENDDO |
---|
47 | c |
---|
48 | c--xms is the varying variable as per Tang and Munkelwitz 1994 |
---|
49 | OPEN (unit=4,file='ri_sul_v2') |
---|
50 | OPEN (unit=10,file='growth_sul_v2') |
---|
51 | c |
---|
52 | DO IRH=1,nbre_rh |
---|
53 | c |
---|
54 | IF (IRH.LE.nbre_dry) THEN |
---|
55 | c |
---|
56 | ratio(IRH)=1.0 |
---|
57 | rhod(IRH)=rhop |
---|
58 | n_r_exact(IRH)=1.521 |
---|
59 | c |
---|
60 | ELSE |
---|
61 | c |
---|
62 | aw=rh_tab(IRH)/100. !--water activity |
---|
63 | molality=B0+aw*(B1+aw*(B2+aw*(B3+aw*B4))) !--molality |
---|
64 | xms=(molality*Mamm/1000.)/(molality*Mamm/1000.+1.)*100. !--solute mass fraction |
---|
65 | c |
---|
66 | rhod(IRH)=A0+xms*(A1+xms*(A2+xms*A3)) !--density |
---|
67 | ratio(IRH)=(rhop/rhod(IRH))**(1./3.) * (xms/100.)**(-1./3.) !--growth factor |
---|
68 | y=Mwater*xms/100./(Mamm*(1.-xms/100.)+xms/100.*Mwater) |
---|
69 | VV=1./rhod(IRH)*(Mwater+(Mamm-Mwater)*y) |
---|
70 | RR=Rwater+(Ramm-Rwater)*y |
---|
71 | nr=(2.*RR/VV+1.)/(1.-RR/VV) |
---|
72 | n_r_exact(IRH)=sqrt(nr) |
---|
73 | c |
---|
74 | ENDIF |
---|
75 | c |
---|
76 | ENDDO |
---|
77 | c |
---|
78 | write(10,20)'RH', (rh_tab(IRH),IRH=1,nbre_rh) |
---|
79 | write(10,10)"growth factor",(ratio(IRH),IRH=1,nbre_rh) |
---|
80 | write(10,10)"density factor",(rhod(IRH)/rhod(1),IRH=1,nbre_rh) |
---|
81 | write(10,10)'rho @ RH', (rhod(IRH),IRH=1,nbre_rh) |
---|
82 | write(10,10)'nr @ RH', (n_r_exact(IRH),IRH=1,nbre_rh) |
---|
83 | 10 FORMAT(A14, F6.3,11(',',F6.3)) |
---|
84 | 20 FORMAT(A14, F6.1,11(',',F6.1)) |
---|
85 | |
---|
86 | c |
---|
87 | c-----Computation of n_r at different wavelengths using |
---|
88 | c n_r(lambda)=n_r(0.589)-0.03(lambda-0.589) from |
---|
89 | c Palmer and Williams (1975) and Kent et al. (1983). |
---|
90 | c Above computed n_r_exact is assumed to be at 0.589 nm |
---|
91 | c |
---|
92 | DO IRH=1,nbre_rh |
---|
93 | DO wve=1,Nwvmax |
---|
94 | c |
---|
95 | n_rr(IRH,wve)=n_r_exact(irh)-0.03*(lambda(wve)*1.E6-0.589) |
---|
96 | n_ii(IRH,wve)=0.0 |
---|
97 | c |
---|
98 | WRITE(4,*) RH_tab(IRH), lambda(wve), n_rr(IRH,wve), n_ii(IRH,wve) |
---|
99 | c |
---|
100 | ENDDO |
---|
101 | ENDDO |
---|
102 | c |
---|
103 | 63 FORMAT(F10.6,E13.6) |
---|
104 | c |
---|
105 | END |
---|