[3385] | 1 | PROGRAM sea_salt |
---|
| 2 | !--based on Tang and Munkelwitz JGR 1994 and Tang JGR 1997 |
---|
| 3 | !--then correct for observed growth factor |
---|
| 4 | IMPLICIT NONE |
---|
| 5 | REAL rhop, rhow, nrw, niw, nrp, nip, xx |
---|
| 6 | PARAMETER (rhop=2.16,rhow=1.0) !--g cm-3 |
---|
| 7 | INTEGER nbre_RH, Nmode |
---|
| 8 | PARAMETER(nbre_rh=12, Nmode=1) |
---|
| 9 | REAL RH_tab(nbre_RH) |
---|
| 10 | DATA RH_tab/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./ |
---|
| 11 | |
---|
| 12 | c--observed growth factor AS |
---|
| 13 | c--Carrico et al. Hygroscopic growth behaviour of a carbon-dominated |
---|
| 14 | C--aerosol in Yosemite Park, Atm. Env. 39, 1393-1404, 2005 |
---|
| 15 | REAL growth(nbre_RH,Nmode) |
---|
| 16 | DATA growth/ |
---|
| 17 | . 1.0,1.0,1.0,1.0,1.0,1.05,1.10,1.20,1.29,1.30,1.31,1.40/ !--AS |
---|
| 18 | c |
---|
| 19 | c--wavelength in m |
---|
| 20 | INTEGER Nwvmax, Nwvhale, Nwvdry, wve, wv |
---|
| 21 | PARAMETER (Nwvmax=200,Nwvhale=60,Nwvdry=61) |
---|
| 22 | REAL lambda(Nwvmax) |
---|
| 23 | REAL lambda_hale(Nwvhale), n_r_hale(Nwvhale), n_i_hale(Nwvhale) |
---|
| 24 | REAL lambda_dry(Nwvdry), n_r_dry(Nwvdry), n_i_dry(Nwvdry) |
---|
| 25 | c |
---|
| 26 | INTEGER IRH, mode |
---|
| 27 | REAL lambda_min, lambda_max |
---|
| 28 | PARAMETER (lambda_min=0.2E-6, lambda_max=5.0E-6) |
---|
| 29 | c |
---|
| 30 | REAL n_r(nbre_rh,Nmode,Nwvmax) |
---|
| 31 | REAL n_i(nbre_rh,Nmode,Nwvmax) |
---|
| 32 | c |
---|
| 33 | REAL ratio(nbre_rh,Nmode) !--diameter growth factor |
---|
| 34 | REAL rhod(nbre_rh,Nmode) !--density |
---|
| 35 | c |
---|
| 36 | DO wve=1,Nwvmax |
---|
| 37 | lambda(wve)= |
---|
| 38 | . lambda_min+FLOAT(wve-1)*(lambda_max-lambda_min)/FLOAT(Nwvmax-1) |
---|
| 39 | ENDDO |
---|
| 40 | c |
---|
| 41 | c--reading water refractive index |
---|
| 42 | OPEN(unit=11,file='r_water_hale.dat') |
---|
| 43 | DO wve=1,Nwvhale |
---|
| 44 | READ(11,*) lambda_hale(wve), n_r_hale(wve) |
---|
| 45 | ENDDO |
---|
| 46 | CLOSE(11) |
---|
| 47 | c print *,'n_r_hale=', n_r_hale |
---|
| 48 | c |
---|
| 49 | OPEN(unit=11,file='i_water_hale.dat') |
---|
| 50 | DO wve=1,Nwvhale |
---|
| 51 | READ(11,*) lambda_hale(wve), n_i_hale(wve) |
---|
| 52 | ENDDO |
---|
| 53 | CLOSE(11) |
---|
| 54 | c print *,'n_i_hale=', n_i_hale |
---|
| 55 | c |
---|
| 56 | c--reading dry bc refractive index |
---|
| 57 | OPEN(unit=11,file='r_bcsoot_she.dat') |
---|
| 58 | DO wve=1,Nwvdry |
---|
| 59 | READ(11,*) lambda_dry(wve), n_r_dry(wve) |
---|
| 60 | ENDDO |
---|
| 61 | CLOSE(11) |
---|
| 62 | c print *,'n_r_dry=', n_r_dry |
---|
| 63 | c |
---|
| 64 | OPEN(unit=11,file='i_bcsoot_she.dat') |
---|
| 65 | DO wve=1,Nwvdry |
---|
| 66 | READ(11,*) lambda_dry(wve), n_i_dry(wve) |
---|
| 67 | ENDDO |
---|
| 68 | CLOSE(11) |
---|
| 69 | c print *,'n_i_dry=', n_i_dry |
---|
| 70 | c |
---|
| 71 | OPEN (unit=11,file='growth_bc_AS_v2') |
---|
| 72 | c |
---|
| 73 | OPEN (unit=21,file='ri_bc_AS_v2') |
---|
| 74 | c |
---|
| 75 | DO mode=1, Nmode |
---|
| 76 | c |
---|
| 77 | DO IRH=1, nbre_RH |
---|
| 78 | c |
---|
| 79 | xx=1./growth(IRH,mode)**3. |
---|
| 80 | rhod(IRH,mode)=xx*rhop+(1.-xx)*rhow |
---|
| 81 | c |
---|
| 82 | DO wve=1,Nwvmax |
---|
| 83 | c |
---|
| 84 | c-plus proche voisin |
---|
| 85 | nrw=n_r_hale(1) |
---|
| 86 | niw=n_i_hale(1) |
---|
| 87 | DO wv=1, Nwvhale |
---|
| 88 | IF (lambda_hale(wv)/1.e9.LT.lambda(wve)) THEN |
---|
| 89 | nrw=n_r_hale(wv) |
---|
| 90 | niw=n_i_hale(wv) |
---|
| 91 | ENDIF |
---|
| 92 | ENDDO |
---|
| 93 | |
---|
| 94 | nrp=n_r_dry(1) |
---|
| 95 | nip=n_i_dry(1) |
---|
| 96 | DO wv=1, Nwvdry |
---|
| 97 | IF (lambda_dry(wv)/1.e9.LT.lambda(wve)) THEN |
---|
| 98 | nrp=n_r_dry(wv) |
---|
| 99 | nip=n_i_dry(wv) |
---|
| 100 | ENDIF |
---|
| 101 | ENDDO |
---|
| 102 | c |
---|
| 103 | n_r(IRH,mode,wve) = xx*nrp +(1.-xx)*nrw |
---|
| 104 | n_i(IRH,mode,wve) = xx*nip +(1.-xx)*niw |
---|
| 105 | c |
---|
| 106 | WRITE(20+mode,*) RH_tab(IRH), lambda(wve), |
---|
| 107 | . n_r(IRH,mode,wve), n_i(IRH,mode,wve) |
---|
| 108 | c |
---|
| 109 | ENDDO !-wve |
---|
| 110 | ENDDO !-IRH |
---|
| 111 | c |
---|
| 112 | WRITE(10+mode,20)'RH', (RH_tab(IRH),IRH=1,nbre_rh) |
---|
| 113 | WRITE(10+mode,10)"growth factor",(growth(IRH,mode),IRH=1,nbre_rh) |
---|
| 114 | WRITE(10+mode,10)"density factor", |
---|
| 115 | . (rhod(IRH,mode)/rhod(1,mode),IRH=1,nbre_rh) |
---|
| 116 | WRITE(10+mode,10)'rho @ RH', (rhod(IRH,mode),IRH=1,nbre_rh) |
---|
| 117 | c |
---|
| 118 | ENDDO |
---|
| 119 | |
---|
| 120 | 10 FORMAT(A14, F6.3,11(',',F6.3)) |
---|
| 121 | 20 FORMAT(A14, F6.1,11(',',F6.1)) |
---|
| 122 | 63 FORMAT(F10.6,E13.6) |
---|
| 123 | |
---|
| 124 | CLOSE(11) |
---|
| 125 | CLOSE(12) |
---|
| 126 | CLOSE(13) |
---|
| 127 | |
---|
| 128 | CLOSE(21) |
---|
| 129 | CLOSE(22) |
---|
| 130 | CLOSE(23) |
---|
| 131 | c |
---|
| 132 | END |
---|