Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer
- Files:
-
- 11 edited
-
aer_sedimnt.F90 (modified) (1 diff)
-
coagulate.F90 (modified) (2 diffs)
-
cond_evap_tstep_mod.F90 (modified) (4 diffs)
-
micphy_tstep.F90 (modified) (6 diffs)
-
miecalc_aer.F90 (modified) (1 diff)
-
nucleation_tstep_mod.F90 (modified) (13 diffs)
-
so2_to_h2so4.F90 (modified) (1 diff)
-
stratdistrib.F90 (modified) (2 diffs)
-
stratemit.F90 (modified) (6 diffs)
-
sulfate_aer_mod.F90 (modified) (12 diffs)
-
traccoag_mod.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/aer_sedimnt.F90
r5082 r5087 122 122 !compute budg_sed_part as sum over bins in kg(S)/m2/s 123 123 budg_sed_part(JL)=budg_sed_part(JL)+ZRHO*ZAERONWM1(JL,nb)*ZVAER(JL,1,nb)*(mSatom/mH2SO4mol) & 124 &*dens_aer_dry*4./3.*RPI*(mdw(nb)/2.)**3124 *dens_aer_dry*4./3.*RPI*(mdw(nb)/2.)**3 125 125 ENDDO 126 126 ENDDO -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/coagulate.F90
r5082 r5087 197 197 DO i=1, nbtr_bin 198 198 delta(i)=((2.*radiuswet(i)+mfppar(i))**3.-(4.*radiuswet(i)**2.+mfppar(i)**2.)**1.5)/ & 199 &(6.*radiuswet(i)*mfppar(i))-2.*radiuswet(i)199 (6.*radiuswet(i)*mfppar(i))-2.*radiuswet(i) 200 200 ENDDO 201 201 … … 208 208 num=4.*RPI*(radiuswet(i)+radiuswet(j))*(Di(i)+Di(j)) 209 209 denom=(radiuswet(i)+radiuswet(j))/(radiuswet(i)+radiuswet(j)+sqrt(delta(i)**2.+delta(j)**2.))+ & 210 &4.*(Di(i)+Di(j))/(sqrt(thvelpar(i)**2.+thvelpar(j)**2.)*(radiuswet(i)+radiuswet(j)))210 4.*(Di(i)+Di(j))/(sqrt(thvelpar(i)**2.+thvelpar(j)**2.)*(radiuswet(i)+radiuswet(j))) 211 211 beta(i,j)=num/denom 212 212 ! -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
r5081 r5087 10 10 11 11 SUBROUTINE condens_evapor_rate_kelvin(R2SO4G,t_seri,pplay,R2SO4, & 12 &DENSO4,f_r_wet,R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR)12 DENSO4,f_r_wet,R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR) 13 13 ! 14 14 ! INPUT: … … 140 140 ! SENFELD 141 141 FL(IK)=FL1*RRSI_wet(IK)**2*( 1.0 +RKNUD ) & 142 &/( 1.0 +ALPHA/(2.0*RKNUD) +RKNUD )142 /( 1.0 +ALPHA/(2.0*RKNUD) +RKNUD ) 143 143 ! TURCO 144 144 ! RL= (4.0/3.0 +0.71/RKNUD)/(1.0+1.0/RKNUD) … … 166 166 !******************************************************************** 167 167 SUBROUTINE condens_evapor_rate(R2SO4G,t_seri,pplay,ACTSO4,R2SO4, & 168 &DENSO4,f_r_wet,FL,ASO4,DNDR)168 DENSO4,f_r_wet,FL,ASO4,DNDR) 169 169 ! 170 170 ! INPUT: … … 263 263 ! SENFELD 264 264 FL(IK)=FL1*RRSI_wet(IK)**2*( 1.0 +RKNUD ) & 265 &/( 1.0 +ALPHA/(2.0*RKNUD) +RKNUD )265 /( 1.0 +ALPHA/(2.0*RKNUD) +RKNUD ) 266 266 ! TURCO 267 267 ! RL= (4.0/3.0 +0.71/RKNUD)/(1.0+1.0/RKNUD) -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/micphy_tstep.F90
r5082 r5087 51 51 !coefficients for H2SO4 density parametrization used for nucleation if ntot<4 52 52 a_xm = 0.7681724 + 1.*(2.1847140 + 1.*(7.1630022 + 1.*(-44.31447 + & 53 &1.*(88.75606 + 1.*(-75.73729 + 1.*23.43228)))))53 1.*(88.75606 + 1.*(-75.73729 + 1.*23.43228))))) 54 54 b_xm = 1.808225e-3 + 1.*(-9.294656e-3 + 1.*(-0.03742148 + 1.*(0.2565321 + & 55 &1.*(-0.5362872 + 1.*(0.4857736 - 1.*0.1629592)))))55 1.*(-0.5362872 + 1.*(0.4857736 - 1.*0.1629592))))) 56 56 c_xm = -3.478524e-6 + 1.*(1.335867e-5 + 1.*(5.195706e-5 + 1.*(-3.717636e-4 + & 57 &1.*(7.990811e-4 + 1.*(-7.458060e-4 + 1.*2.58139e-4 )))))57 1.*(7.990811e-4 + 1.*(-7.458060e-4 + 1.*2.58139e-4 ))))) 58 58 59 59 IF(.not.flag_new_strat_compo) THEN … … 85 85 ! convert tr_seri(GASH2SO4) (in kg/kgA) to H2SO4 number density (in molecules/cm3) 86 86 rhoa=tr_seri(ilon,ilev,id_H2SO4_strat) & 87 &*pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mH2SO4mol87 *pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mH2SO4mol 88 88 ! compute nucleation rate in kg(H2SO4)/kgA/s 89 89 CALL nucleation_rate(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev),rh(ilon,ilev), & 90 &a_xm,b_xm,c_xm,nucl_rate,ntot,x)90 a_xm,b_xm,c_xm,nucl_rate,ntot,x) 91 91 !NL - add nucleation box (if flag on) 92 92 IF (flag_nuc_rate_box) THEN … … 102 102 f_r_wetik(:) = f_r_wetB(ilon,ilev,:) 103 103 CALL condens_evapor_rate_kelvin(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev), & 104 &R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), &105 &R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR)104 R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), & 105 R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR) 106 106 ELSE 107 107 CALL condens_evapor_rate(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev), & 108 &ACTSO4(ilon,ilev),R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), &109 &FL,ASO4,DNDR)108 ACTSO4(ilon,ilev),R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), & 109 FL,ASO4,DNDR) 110 110 ENDIF 111 111 ! Compute H2SO4 saturate vapor for big particules … … 134 134 ! compute fluxes as diagnostic in [kg(S)/m2/layer/s] (now - for evap and + for cond) 135 135 budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & 136 &*cond_evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys136 *cond_evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys 137 137 budg_3D_nucl(ilon,ilev)=budg_3D_nucl(ilon,ilev)+mSatom/mH2SO4mol & 138 &*nucl_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys138 *nucl_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys 139 139 ! update time step 140 140 PDT=PDT-dt … … 142 142 ! convert tr_seri(GASH2SO4) (in kg/kgA) to H2SO4 number density (in molecules/cm3) 143 143 rhoa=tr_seri(ilon,ilev,id_H2SO4_strat) & 144 &*pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mH2SO4mol144 *pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mH2SO4mol 145 145 ! compute cond/evap rate in kg(H2SO4)/kgA/s (now only evap for pdtphys) 146 146 IF(flag_new_strat_compo) THEN 147 147 CALL condens_evapor_rate_kelvin(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev), & 148 &R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), &149 &R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR)148 R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), & 149 R2SO4ik,DENSO4ik,f_r_wetik,FL,ASO4,DNDR) 150 150 ELSE 151 151 CALL condens_evapor_rate(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev), & 152 &ACTSO4(ilon,ilev),R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), &153 &FL,ASO4,DNDR)152 ACTSO4(ilon,ilev),R2SO4(ilon,ilev),DENSO4(ilon,ilev),f_r_wet(ilon,ilev), & 153 FL,ASO4,DNDR) 154 154 ENDIF 155 155 ! limit evaporation (negative FL) over one physics time step to H2SO4 content of the droplet … … 170 170 ! compute fluxes as diagnostic in [kg(S)/m2/layer/s] (now - for evap and + for cond) 171 171 budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & 172 &*evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG172 *evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG 173 173 ! compute vertically integrated flux due to the net effect of nucleation and condensation/evaporation 174 174 budg_h2so4_to_part(ilon)=budg_h2so4_to_part(ilon)+(H2SO4_init-tr_seri(ilon,ilev,id_H2SO4_strat)) & 175 &*mSatom/mH2SO4mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys175 *mSatom/mH2SO4mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys 176 176 ENDIF 177 177 ENDDO -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/miecalc_aer.F90
r5082 r5087 258 258 Nwv=1 259 259 dlambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)= & 260 &lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)- &261 &lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv+1)260 lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)- & 261 lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv+1) 262 262 DO Nwv=2, NwvmaxLW-1 263 263 dlambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)= & 264 &(lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv-1)- &265 &lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv+1))/2.264 (lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv-1)- & 265 lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv+1))/2. 266 266 ENDDO 267 267 Nwv=NwvmaxLW 268 268 dlambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)= & 269 &lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv-1)- &270 &lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv)269 lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv-1)- & 270 lambda_int(NwvmaxSW+nwave_sw+nwave_lw+Nwv) 271 271 272 272 IF (refr_ind_interpol) THEN -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/nucleation_tstep_mod.F90
r5082 r5087 49 49 VH2SO4mol=mH2SO4mol/(1.E-3*(a_xm+t_seri*(b_xm+t_seri*c_xm))) !cm3 50 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)51 *100.*(2.*VH2SO4mol**(1./3.))**2. !1/(cm3s) 52 52 ntot_n=2.0 53 53 x_n=1.0 … … 60 60 ! airn=pplay/t_seri/RD/1.E3*RNAVO/RMD ! molec cm-3 (for future use, to be confirmed) 61 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)62 x_n,x_i,na_n,na_i,rc_n,rc_i,n_i,kinetic_n,kinetic_i) 63 63 ENDIF 64 64 … … 197 197 198 198 x= 0.7409967177282139 - 0.002663785665140117*t + 0.002010478847383187*LOG(rh) & 199 &- 0.0001832894131464668*t*LOG(rh) + 0.001574072538464286*LOG(rh)**2 &200 &- 0.00001790589121766952*t*LOG(rh)**2 + 0.0001844027436573778*LOG(rh)**3 &201 &- 1.503452308794887E-6*t*LOG(rh)**3 - 0.003499978417957668*LOG(rhoa) &202 &+ 0.0000504021689382576*t*LOG(rhoa)199 - 0.0001832894131464668*t*LOG(rh) + 0.001574072538464286*LOG(rh)**2 & 200 - 0.00001790589121766952*t*LOG(rh)**2 + 0.0001844027436573778*LOG(rh)**3 & 201 - 1.503452308794887E-6*t*LOG(rh)**3 - 0.003499978417957668*LOG(rhoa) & 202 + 0.0000504021689382576*t*LOG(rhoa) 203 203 204 204 jnuc= 0.1430901615568665 + 2.219563673425199*t - 0.02739106114964264*t**2 + & 205 &0.00007228107239317088*t**3 + 5.91822263375044/x + &206 &0.1174886643003278*LOG(rh) + 0.4625315047693772*t*LOG(rh) - &207 &0.01180591129059253*t**2*LOG(rh) + &208 &0.0000404196487152575*t**3*LOG(rh) + (15.79628615047088*LOG(rh))/x - &209 &0.215553951893509*LOG(rh)**2 - 0.0810269192332194*t*LOG(rh)**2 + &210 &0.001435808434184642*t**2*LOG(rh)**2 - &211 &4.775796947178588E-6*t**3*LOG(rh)**2 - &212 &(2.912974063702185*LOG(rh)**2)/x - 3.588557942822751*LOG(rh)**3 + &213 &0.04950795302831703*t*LOG(rh)**3 - &214 &0.0002138195118737068*t**2*LOG(rh)**3 + &215 &3.108005107949533E-7*t**3*LOG(rh)**3 - &216 &(0.02933332747098296*LOG(rh)**3)/x + &217 &1.145983818561277*LOG(rhoa) - &218 &0.6007956227856778*t*LOG(rhoa) + &219 &0.00864244733283759*t**2*LOG(rhoa) - &220 &0.00002289467254710888*t**3*LOG(rhoa) - &221 &(8.44984513869014*LOG(rhoa))/x + &222 &2.158548369286559*LOG(rh)*LOG(rhoa) + &223 &0.0808121412840917*t*LOG(rh)*LOG(rhoa) - &224 &0.0004073815255395214*t**2*LOG(rh)*LOG(rhoa) - &225 &4.019572560156515E-7*t**3*LOG(rh)*LOG(rhoa) + &226 &(0.7213255852557236*LOG(rh)*LOG(rhoa))/x + &227 &1.62409850488771*LOG(rh)**2*LOG(rhoa) - &228 &0.01601062035325362*t*LOG(rh)**2*LOG(rhoa) + &229 &0.00003771238979714162*t**2*LOG(rh)**2*LOG(rhoa) + &230 &3.217942606371182E-8*t**3*LOG(rh)**2*LOG(rhoa) - &231 &(0.01132550810022116*LOG(rh)**2*LOG(rhoa))/x + &232 &9.71681713056504*LOG(rhoa)**2 - &233 &0.1150478558347306*t*LOG(rhoa)**2 + &234 &0.0001570982486038294*t**2*LOG(rhoa)**2 + &235 &4.009144680125015E-7*t**3*LOG(rhoa)**2 + &236 &(0.7118597859976135*LOG(rhoa)**2)/x - &237 &1.056105824379897*LOG(rh)*LOG(rhoa)**2 + &238 &0.00903377584628419*t*LOG(rh)*LOG(rhoa)**2 - &239 &0.00001984167387090606*t**2*LOG(rh)*LOG(rhoa)**2 + &240 &2.460478196482179E-8*t**3*LOG(rh)*LOG(rhoa)**2 - &241 &(0.05790872906645181*LOG(rh)*LOG(rhoa)**2)/x - &242 &0.1487119673397459*LOG(rhoa)**3 + &243 &0.002835082097822667*t*LOG(rhoa)**3 - &244 &9.24618825471694E-6*t**2*LOG(rhoa)**3 + &245 &5.004267665960894E-9*t**3*LOG(rhoa)**3 - &246 &(0.01270805101481648*LOG(rhoa)**3)/x205 0.00007228107239317088*t**3 + 5.91822263375044/x + & 206 0.1174886643003278*LOG(rh) + 0.4625315047693772*t*LOG(rh) - & 207 0.01180591129059253*t**2*LOG(rh) + & 208 0.0000404196487152575*t**3*LOG(rh) + (15.79628615047088*LOG(rh))/x - & 209 0.215553951893509*LOG(rh)**2 - 0.0810269192332194*t*LOG(rh)**2 + & 210 0.001435808434184642*t**2*LOG(rh)**2 - & 211 4.775796947178588E-6*t**3*LOG(rh)**2 - & 212 (2.912974063702185*LOG(rh)**2)/x - 3.588557942822751*LOG(rh)**3 + & 213 0.04950795302831703*t*LOG(rh)**3 - & 214 0.0002138195118737068*t**2*LOG(rh)**3 + & 215 3.108005107949533E-7*t**3*LOG(rh)**3 - & 216 (0.02933332747098296*LOG(rh)**3)/x + & 217 1.145983818561277*LOG(rhoa) - & 218 0.6007956227856778*t*LOG(rhoa) + & 219 0.00864244733283759*t**2*LOG(rhoa) - & 220 0.00002289467254710888*t**3*LOG(rhoa) - & 221 (8.44984513869014*LOG(rhoa))/x + & 222 2.158548369286559*LOG(rh)*LOG(rhoa) + & 223 0.0808121412840917*t*LOG(rh)*LOG(rhoa) - & 224 0.0004073815255395214*t**2*LOG(rh)*LOG(rhoa) - & 225 4.019572560156515E-7*t**3*LOG(rh)*LOG(rhoa) + & 226 (0.7213255852557236*LOG(rh)*LOG(rhoa))/x + & 227 1.62409850488771*LOG(rh)**2*LOG(rhoa) - & 228 0.01601062035325362*t*LOG(rh)**2*LOG(rhoa) + & 229 0.00003771238979714162*t**2*LOG(rh)**2*LOG(rhoa) + & 230 3.217942606371182E-8*t**3*LOG(rh)**2*LOG(rhoa) - & 231 (0.01132550810022116*LOG(rh)**2*LOG(rhoa))/x + & 232 9.71681713056504*LOG(rhoa)**2 - & 233 0.1150478558347306*t*LOG(rhoa)**2 + & 234 0.0001570982486038294*t**2*LOG(rhoa)**2 + & 235 4.009144680125015E-7*t**3*LOG(rhoa)**2 + & 236 (0.7118597859976135*LOG(rhoa)**2)/x - & 237 1.056105824379897*LOG(rh)*LOG(rhoa)**2 + & 238 0.00903377584628419*t*LOG(rh)*LOG(rhoa)**2 - & 239 0.00001984167387090606*t**2*LOG(rh)*LOG(rhoa)**2 + & 240 2.460478196482179E-8*t**3*LOG(rh)*LOG(rhoa)**2 - & 241 (0.05790872906645181*LOG(rh)*LOG(rhoa)**2)/x - & 242 0.1487119673397459*LOG(rhoa)**3 + & 243 0.002835082097822667*t*LOG(rhoa)**3 - & 244 9.24618825471694E-6*t**2*LOG(rhoa)**3 + & 245 5.004267665960894E-9*t**3*LOG(rhoa)**3 - & 246 (0.01270805101481648*LOG(rhoa)**3)/x 247 247 jnuc=EXP(jnuc) !1/(cm3s) 248 248 249 249 ntot =-0.002954125078716302 - 0.0976834264241286*t + 0.001024847927067835*t**2 - 2.186459697726116E-6*t**3 - & 250 &0.1017165718716887/x - 0.002050640345231486*LOG(rh) - 0.007585041382707174*t*LOG(rh) + &251 &0.0001926539658089536*t**2*LOG(rh) - 6.70429719683894E-7*t**3*LOG(rh) - &252 &(0.2557744774673163*LOG(rh))/x + 0.003223076552477191*LOG(rh)**2 + 0.000852636632240633*t*LOG(rh)**2 - &253 &0.00001547571354871789*t**2*LOG(rh)**2 + 5.666608424980593E-8*t**3*LOG(rh)**2 + &254 &(0.03384437400744206*LOG(rh)**2)/x + 0.04743226764572505*LOG(rh)**3 - &255 &0.0006251042204583412*t*LOG(rh)**3 + 2.650663328519478E-6*t**2*LOG(rh)**3 - &256 &3.674710848763778E-9*t**3*LOG(rh)**3 - (0.0002672510825259393*LOG(rh)**3)/x - &257 &0.01252108546759328*LOG(rhoa) + 0.005806550506277202*t*LOG(rhoa) - &258 &0.0001016735312443444*t**2*LOG(rhoa) + 2.881946187214505E-7*t**3*LOG(rhoa) + &259 &(0.0942243379396279*LOG(rhoa))/x - 0.0385459592773097*LOG(rh)*LOG(rhoa) - &260 &0.0006723156277391984*t*LOG(rh)*LOG(rhoa) + 2.602884877659698E-6*t**2*LOG(rh)*LOG(rhoa) + &261 &1.194163699688297E-8*t**3*LOG(rh)*LOG(rhoa) - (0.00851515345806281*LOG(rh)*LOG(rhoa))/x - &262 &0.01837488495738111*LOG(rh)**2*LOG(rhoa) + 0.0001720723574407498*t*LOG(rh)**2*LOG(rhoa) - &263 &3.717657974086814E-7*t**2*LOG(rh)**2*LOG(rhoa) - &264 &5.148746022615196E-10*t**3*LOG(rh)**2*LOG(rhoa) + &265 &(0.0002686602132926594*LOG(rh)**2*LOG(rhoa))/x - 0.06199739728812199*LOG(rhoa)**2 + &266 &0.000906958053583576*t*LOG(rhoa)**2 - 9.11727926129757E-7*t**2*LOG(rhoa)**2 - &267 &5.367963396508457E-9*t**3*LOG(rhoa)**2 - (0.007742343393937707*LOG(rhoa)**2)/x + &268 &0.0121827103101659*LOG(rh)*LOG(rhoa)**2 - 0.0001066499571188091*t*LOG(rh)*LOG(rhoa)**2 + &269 &2.534598655067518E-7*t**2*LOG(rh)*LOG(rhoa)**2 - &270 &3.635186504599571E-10*t**3*LOG(rh)*LOG(rhoa)**2 + &271 &(0.0006100650851863252*LOG(rh)*LOG(rhoa)**2)/x + 0.0003201836700403512*LOG(rhoa)**3 - &272 &0.0000174761713262546*t*LOG(rhoa)**3 + 6.065037668052182E-8*t**2*LOG(rhoa)**3 - &273 &1.421771723004557E-11*t**3*LOG(rhoa)**3 + (0.0001357509859501723*LOG(rhoa)**3)/x250 0.1017165718716887/x - 0.002050640345231486*LOG(rh) - 0.007585041382707174*t*LOG(rh) + & 251 0.0001926539658089536*t**2*LOG(rh) - 6.70429719683894E-7*t**3*LOG(rh) - & 252 (0.2557744774673163*LOG(rh))/x + 0.003223076552477191*LOG(rh)**2 + 0.000852636632240633*t*LOG(rh)**2 - & 253 0.00001547571354871789*t**2*LOG(rh)**2 + 5.666608424980593E-8*t**3*LOG(rh)**2 + & 254 (0.03384437400744206*LOG(rh)**2)/x + 0.04743226764572505*LOG(rh)**3 - & 255 0.0006251042204583412*t*LOG(rh)**3 + 2.650663328519478E-6*t**2*LOG(rh)**3 - & 256 3.674710848763778E-9*t**3*LOG(rh)**3 - (0.0002672510825259393*LOG(rh)**3)/x - & 257 0.01252108546759328*LOG(rhoa) + 0.005806550506277202*t*LOG(rhoa) - & 258 0.0001016735312443444*t**2*LOG(rhoa) + 2.881946187214505E-7*t**3*LOG(rhoa) + & 259 (0.0942243379396279*LOG(rhoa))/x - 0.0385459592773097*LOG(rh)*LOG(rhoa) - & 260 0.0006723156277391984*t*LOG(rh)*LOG(rhoa) + 2.602884877659698E-6*t**2*LOG(rh)*LOG(rhoa) + & 261 1.194163699688297E-8*t**3*LOG(rh)*LOG(rhoa) - (0.00851515345806281*LOG(rh)*LOG(rhoa))/x - & 262 0.01837488495738111*LOG(rh)**2*LOG(rhoa) + 0.0001720723574407498*t*LOG(rh)**2*LOG(rhoa) - & 263 3.717657974086814E-7*t**2*LOG(rh)**2*LOG(rhoa) - & 264 5.148746022615196E-10*t**3*LOG(rh)**2*LOG(rhoa) + & 265 (0.0002686602132926594*LOG(rh)**2*LOG(rhoa))/x - 0.06199739728812199*LOG(rhoa)**2 + & 266 0.000906958053583576*t*LOG(rhoa)**2 - 9.11727926129757E-7*t**2*LOG(rhoa)**2 - & 267 5.367963396508457E-9*t**3*LOG(rhoa)**2 - (0.007742343393937707*LOG(rhoa)**2)/x + & 268 0.0121827103101659*LOG(rh)*LOG(rhoa)**2 - 0.0001066499571188091*t*LOG(rh)*LOG(rhoa)**2 + & 269 2.534598655067518E-7*t**2*LOG(rh)*LOG(rhoa)**2 - & 270 3.635186504599571E-10*t**3*LOG(rh)*LOG(rhoa)**2 + & 271 (0.0006100650851863252*LOG(rh)*LOG(rhoa)**2)/x + 0.0003201836700403512*LOG(rhoa)**3 - & 272 0.0000174761713262546*t*LOG(rhoa)**3 + 6.065037668052182E-8*t**2*LOG(rhoa)**3 - & 273 1.421771723004557E-11*t**3*LOG(rhoa)**3 + (0.0001357509859501723*LOG(rhoa)**3)/x 274 274 ntot=EXP(ntot) 275 275 … … 285 285 286 286 rhotres=EXP( -279.2430007512709 + 11.73439886096903*rh + 22700.92970508331/t & 287 &- (1088.644983466801*rh)/t + 1.144362942094912*t &288 &- 0.03023314602163684*rh*t - 0.001302541390154324*t**2 &289 &- 6.386965238433532*LOG(rh) + (854.980361026715*LOG(rh))/t &290 &+ 0.00879662256826497*t*LOG(rh)) !1/cm3287 - (1088.644983466801*rh)/t + 1.144362942094912*t & 288 - 0.03023314602163684*rh*t - 0.001302541390154324*t**2 & 289 - 6.386965238433532*LOG(rh) + (854.980361026715*LOG(rh))/t & 290 + 0.00879662256826497*t*LOG(rh)) !1/cm3 291 291 292 292 RETURN … … 297 297 298 298 SUBROUTINE newbinapara(t,satrat,rhoa,csi,airn,ipr,jnuc_n_real,ntot_n_real,jnuc_i_real,ntot_i_real, & 299 & x_n_real,x_i_real,na_n_real,na_i_real,rc_n_real,rc_i_real,n_i_real, &300 &kinetic_n,kinetic_i)299 x_n_real,x_i_real,na_n_real,na_i_real,rc_n_real,rc_i_real,n_i_real, & 300 kinetic_n,kinetic_i) 301 301 302 302 ! Fortran 90 subroutine newbinapara … … 457 457 !Critical cluster composition (valid for both cases, bounds not used here) 458 458 x_n= 7.9036365428891719E-1 - 2.8414059650092153E-3*tln + 1.4976802556584141E-2*LOG(satratln) & 459 &- 2.4511581740839115E-4*tln*LOG(satratln) + 3.4319869471066424E-3 *LOG(satratln)**2 &460 &- 2.8799393617748428E-5*tln*LOG(satratln)**2 + 3.0174314126331765E-4*LOG(satratln)**3 &461 &- 2.2673492408841294E-6*tln*LOG(satratln)**3 - 4.3948464567032377E-3*LOG(rhoaln) &462 &+ 5.3305314722492146E-5*tln*LOG(rhoaln)459 - 2.4511581740839115E-4*tln*LOG(satratln) + 3.4319869471066424E-3 *LOG(satratln)**2 & 460 - 2.8799393617748428E-5*tln*LOG(satratln)**2 + 3.0174314126331765E-4*LOG(satratln)**3 & 461 - 2.2673492408841294E-6*tln*LOG(satratln)**3 - 4.3948464567032377E-3*LOG(rhoaln) & 462 + 5.3305314722492146E-5*tln*LOG(rhoaln) 463 463 x_i= 7.9036365428891719E-1 - 2.8414059650092153E-3*tli + 1.4976802556584141E-2*LOG(satratli) & 464 &- 2.4511581740839115E-4*tli*LOG(satratli) + 3.4319869471066424E-3 *LOG(satratli)**2 &465 &- 2.8799393617748428E-5*tli*LOG(satratli)**2 + 3.0174314126331765E-4*LOG(satratli)**3 &466 &- 2.2673492408841294E-6*tli*LOG(satratli)**3 - 4.3948464567032377E-3*LOG(rhoali) &467 &+ 5.3305314722492146E-5*tli*LOG(rhoali)464 - 2.4511581740839115E-4*tli*LOG(satratli) + 3.4319869471066424E-3 *LOG(satratli)**2 & 465 - 2.8799393617748428E-5*tli*LOG(satratli)**2 + 3.0174314126331765E-4*LOG(satratli)**3 & 466 - 2.2673492408841294E-6*tli*LOG(satratli)**3 - 4.3948464567032377E-3*LOG(rhoali) & 467 + 5.3305314722492146E-5*tli*LOG(rhoali) 468 468 469 469 x_n=MIN(MAX(x_n,1.E-30),1.) … … 475 475 IF (satratln >= 1.E-2 .AND. satratln <= 1.) THEN 476 476 kinrhotresn=EXP(7.8920778706888086E+1 + 7.3665492897447082*satratln - 1.2420166571163805E+4/tln & 477 &+ (-6.1831234251470971E+2*satratln)/tln - 2.4501159970109945E-2*tln &478 &-1.3463066443605762E-2*satratln*tln + 8.3736373989909194E-06*tln**2 &479 &-1.4673887785408892*LOG(satratln) + (-3.2141890006517094E+1*LOG(satratln))/tln &480 & + 2.7137429081917556E-3*tln*LOG(satratln)) !1/cm3477 + (-6.1831234251470971E+2*satratln)/tln - 2.4501159970109945E-2*tln & 478 -1.3463066443605762E-2*satratln*tln + 8.3736373989909194E-06*tln**2 & 479 -1.4673887785408892*LOG(satratln) + (-3.2141890006517094E+1*LOG(satratln))/tln & 480 + 2.7137429081917556E-3*tln*LOG(satratln)) !1/cm3 481 481 IF (kinrhotresn<rhoaln) kinetic_n=.TRUE. 482 482 ENDIF … … 484 484 IF (satratln >= 1.E-4 .AND. satratln < 1.E-2) THEN 485 485 kinrhotresn=EXP(7.9074383049843647E+1 - 2.8746005462158347E+1*satratln - 1.2070272068458380E+4/tln & 486 &+ (-5.9205040320056632E+3*satratln)/tln - 2.4800372593452726E-2*tln &487 &-4.3983007681295948E-2*satratln*tln + 2.5943854791342071E-5*tln**2 &488 &-2.3141363245211317*LOG(satratln) + (9.9186787997857735E+1*LOG(satratln))/tln &489 &+ 5.6819382556144681E-3*tln*LOG(satratln)) !1/cm3486 + (-5.9205040320056632E+3*satratln)/tln - 2.4800372593452726E-2*tln & 487 -4.3983007681295948E-2*satratln*tln + 2.5943854791342071E-5*tln**2 & 488 -2.3141363245211317*LOG(satratln) + (9.9186787997857735E+1*LOG(satratln))/tln & 489 + 5.6819382556144681E-3*tln*LOG(satratln)) !1/cm3 490 490 IF (kinrhotresn<rhoaln) kinetic_n=.TRUE. 491 491 ENDIF … … 493 493 IF (satratln >= 5.E-6 .AND. satratln < 1.E-4) THEN 494 494 kinrhotresn=EXP(8.5599712000361677E+1 + 2.7335119660796581E+3*satratln - 1.1842350246291651E+4/tln & 495 &+ (-1.2439843468881438E+6*satratln)/tln - 5.4536964974944230E-2*tln &496 &+ 5.0886987425326087*satratln*tln + 7.1964722655507067E-5*tln**2 &497 &-2.4472627526306372*LOG(satratln) + (1.7561478001423779E+2*LOG(satratln))/tln &498 &+ 6.2640132818141811E-3*tln*LOG(satratln)) !1/cm3495 + (-1.2439843468881438E+6*satratln)/tln - 5.4536964974944230E-2*tln & 496 + 5.0886987425326087*satratln*tln + 7.1964722655507067E-5*tln**2 & 497 -2.4472627526306372*LOG(satratln) + (1.7561478001423779E+2*LOG(satratln))/tln & 498 + 6.2640132818141811E-3*tln*LOG(satratln)) !1/cm3 499 499 IF (kinrhotresn<rhoaln) kinetic_n=.TRUE. 500 500 ENDIF … … 509 509 ELSE 510 510 jnuc_n= 2.1361182605986115E-1 + 3.3827029855551838*tln -3.2423555796175563E-2*tln**2 + & 511 &7.0120069477221989E-5*tln**3 +8.0286874752695141/x_n + &512 &(-2.6939840579762231E-1)*LOG(satratln) +1.6079879299099518*tln*LOG(satratln) + &513 &(-1.9667486968141933E-2)*tln**2*LOG(satratln) + &514 &5.5244755979770844E-5*tln**3*LOG(satratln) + (7.8884704837892468*LOG(satratln))/x_n + &515 &4.6374659198909596*LOG(satratln)**2 - 8.2002809894792153E-2*tln*LOG(satratln)**2 + &516 &8.5077424451172196E-4*tln**2*LOG(satratln)**2 + &517 &(-2.6518510168987462E-6)*tln**3*LOG(satratln)**2 + &518 &(-1.4625482500575278*LOG(satratln)**2)/x_n - 5.2413002989192037E-1*LOG(satratln)**3 + &519 &5.2755117653715865E-3*tln*LOG(satratln)**3 + &520 &(-2.9491061332113830E-6)*tln**2*LOG(satratln)**3 + &521 &(-2.4815454194486752E-8)*tln**3*LOG(satratln)**3 + &522 &(-5.2663760117394626E-2*LOG(satratln)**3)/x_n + &523 &1.6496664658266762*LOG(rhoaln) + &524 &(-8.0809397859218401E-1)*tln*LOG(rhoaln) + &525 &8.9302927091946642E-3*tln**2*LOG(rhoaln) + &526 &(-1.9583649496497497E-5)*tln**3*LOG(rhoaln) + &527 &(-8.9505572676891685*LOG(rhoaln))/x_n + &528 &(-3.0025283601622881E+1)*LOG(satratln)*LOG(rhoaln) + &529 &3.0783365644763633E-1*tln*LOG(satratln)*LOG(rhoaln) + &530 &(-7.4521756337984706E-4)*tln**2*LOG(satratln)*LOG(rhoaln) + &531 &(-5.7651433870681853E-7)*tln**3*LOG(satratln)*LOG(rhoaln) + &532 &(1.2872868529673207*LOG(satratln)*LOG(rhoaln))/x_n + &533 &(-6.1739867501526535E-1)*LOG(satratln)**2*LOG(rhoaln) + &534 &7.2347385705333975E-3*tln*LOG(satratln)**2*LOG(rhoaln) + &535 &(-3.0640494530822439E-5)*tln**2*LOG(satratln)**2*LOG(rhoaln) + &536 &6.5944609194346214E-8*tln**3*LOG(satratln)**2*LOG(rhoaln) + &537 &(-2.8681650332461055E-2*LOG(satratln)**2*LOG(rhoaln))/x_n + &538 &6.5213802375160306*LOG(rhoaln)**2 + &539 &(-4.7907162004793016E-2)*tln*LOG(rhoaln)**2 + &540 &(-1.0727890114215117E-4)*tln**2*LOG(rhoaln)**2 + &541 &5.6401818280534507E-7*tln**3*LOG(rhoaln)**2 + &542 &(5.4113070888923009E-1*LOG(rhoaln)**2)/x_n + &543 &5.2062808476476330E-1*LOG(satratln)*LOG(rhoaln)**2 + &544 &(-6.0696882500824584E-3)*tln*LOG(satratln)*LOG(rhoaln)**2 + &545 &2.3851383302608477E-5*tln**2*LOG(satratln)*LOG(rhoaln)**2 + &546 &(-1.5243837103067096E-8)*tln**3*LOG(satratln)*LOG(rhoaln)**2 + &547 &(-5.6543192378015687E-2*LOG(satratln)*LOG(rhoaln)**2)/x_n + &548 &(-1.1630806410696815E-1)*LOG(rhoaln)**3 + &549 &1.3806404273119610E-3*tln*LOG(rhoaln)**3 + &550 &(-2.0199865087650833E-6)*tln**2*LOG(rhoaln)**3 + &551 &(-3.0200284885763192E-9)*tln**3*LOG(rhoaln)**3 + &552 &(-6.9425267104126316E-3*LOG(rhoaln)**3)/x_n511 7.0120069477221989E-5*tln**3 +8.0286874752695141/x_n + & 512 (-2.6939840579762231E-1)*LOG(satratln) +1.6079879299099518*tln*LOG(satratln) + & 513 (-1.9667486968141933E-2)*tln**2*LOG(satratln) + & 514 5.5244755979770844E-5*tln**3*LOG(satratln) + (7.8884704837892468*LOG(satratln))/x_n + & 515 4.6374659198909596*LOG(satratln)**2 - 8.2002809894792153E-2*tln*LOG(satratln)**2 + & 516 8.5077424451172196E-4*tln**2*LOG(satratln)**2 + & 517 (-2.6518510168987462E-6)*tln**3*LOG(satratln)**2 + & 518 (-1.4625482500575278*LOG(satratln)**2)/x_n - 5.2413002989192037E-1*LOG(satratln)**3 + & 519 5.2755117653715865E-3*tln*LOG(satratln)**3 + & 520 (-2.9491061332113830E-6)*tln**2*LOG(satratln)**3 + & 521 (-2.4815454194486752E-8)*tln**3*LOG(satratln)**3 + & 522 (-5.2663760117394626E-2*LOG(satratln)**3)/x_n + & 523 1.6496664658266762*LOG(rhoaln) + & 524 (-8.0809397859218401E-1)*tln*LOG(rhoaln) + & 525 8.9302927091946642E-3*tln**2*LOG(rhoaln) + & 526 (-1.9583649496497497E-5)*tln**3*LOG(rhoaln) + & 527 (-8.9505572676891685*LOG(rhoaln))/x_n + & 528 (-3.0025283601622881E+1)*LOG(satratln)*LOG(rhoaln) + & 529 3.0783365644763633E-1*tln*LOG(satratln)*LOG(rhoaln) + & 530 (-7.4521756337984706E-4)*tln**2*LOG(satratln)*LOG(rhoaln) + & 531 (-5.7651433870681853E-7)*tln**3*LOG(satratln)*LOG(rhoaln) + & 532 (1.2872868529673207*LOG(satratln)*LOG(rhoaln))/x_n + & 533 (-6.1739867501526535E-1)*LOG(satratln)**2*LOG(rhoaln) + & 534 7.2347385705333975E-3*tln*LOG(satratln)**2*LOG(rhoaln) + & 535 (-3.0640494530822439E-5)*tln**2*LOG(satratln)**2*LOG(rhoaln) + & 536 6.5944609194346214E-8*tln**3*LOG(satratln)**2*LOG(rhoaln) + & 537 (-2.8681650332461055E-2*LOG(satratln)**2*LOG(rhoaln))/x_n + & 538 6.5213802375160306*LOG(rhoaln)**2 + & 539 (-4.7907162004793016E-2)*tln*LOG(rhoaln)**2 + & 540 (-1.0727890114215117E-4)*tln**2*LOG(rhoaln)**2 + & 541 5.6401818280534507E-7*tln**3*LOG(rhoaln)**2 + & 542 (5.4113070888923009E-1*LOG(rhoaln)**2)/x_n + & 543 5.2062808476476330E-1*LOG(satratln)*LOG(rhoaln)**2 + & 544 (-6.0696882500824584E-3)*tln*LOG(satratln)*LOG(rhoaln)**2 + & 545 2.3851383302608477E-5*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & 546 (-1.5243837103067096E-8)*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & 547 (-5.6543192378015687E-2*LOG(satratln)*LOG(rhoaln)**2)/x_n + & 548 (-1.1630806410696815E-1)*LOG(rhoaln)**3 + & 549 1.3806404273119610E-3*tln*LOG(rhoaln)**3 + & 550 (-2.0199865087650833E-6)*tln**2*LOG(rhoaln)**3 + & 551 (-3.0200284885763192E-9)*tln**3*LOG(rhoaln)**3 + & 552 (-6.9425267104126316E-3*LOG(rhoaln)**3)/x_n 553 553 jnuc_n=EXP(jnuc_n) 554 554 555 555 ntot_n =-3.5863435141979573E-3 - 1.0098670235841110E-1*tln + 8.9741268319259721E-4*tln**2 - 1.4855098605195757E-6*tln**3 & 556 &- 1.2080330016937095E-1/x_n + 1.1902674923928015E-3*LOG(satratln) - 1.9211358507172177E-2*tln*LOG(satratln) + &557 &2.4648094311204255E-4*tln**2*LOG(satratln) - 7.5641448594711666E-7*tln**3*LOG(satratln) + &558 &(-2.0668639384228818E-02*LOG(satratln))/x_n - 3.7593072011595188E-2*LOG(satratln)**2 + &559 &9.0993182774415718E-4 *tln*LOG(satratln)**2 + &560 &(-9.5698412164297149E-6)*tln**2*LOG(satratln)**2 + 3.7163166416110421E-8*tln**3*LOG(satratln)**2 + &561 &(1.1026579525210847E-2*LOG(satratln)**2)/x_n + 1.1530844115561925E-2 *LOG(satratln)**3 + &562 &(-1.8083253906466668E-4)*tln*LOG(satratln)**3 + 8.0213604053330654E-7*tln**2*LOG(satratln)**3 + &563 &(-8.5797885383051337E-10)*tln**3*LOG(satratln)**3 + (1.0243693899717402E-3*LOG(satratln)**3)/x_n + &564 &(-1.7248695296299649E-2)*LOG(rhoaln) + 1.1294004162437157E-2*tln*LOG(rhoaln) + &565 &(-1.2283640163189278E-4)*tln**2*LOG(rhoaln) + 2.7391732258259009E-7*tln**3*LOG(rhoaln) + &566 &(6.8505583974029602E-2*LOG(rhoaln))/x_n +2.9750968179523635E-1*LOG(satratln)*LOG(rhoaln) + &567 &(-3.6681154503992296E-3)*tln*LOG(satratln)*LOG(rhoaln) + 1.0636473034653114E-5*tln**2*LOG(satratln)*LOG(rhoaln)+ &568 &5.8687098466515866E-9*tln**3*LOG(satratln)*LOG(rhoaln) + (-5.2028866094191509E-3*LOG(satratln)*LOG(rhoaln))/x_n+ &569 &7.6971988880587231E-4*LOG(satratln)**2*LOG(rhoaln) - 2.4605575820433763E-5*tln*LOG(satratln)**2*LOG(rhoaln) + &570 &2.3818484400893008E-7*tln**2*LOG(satratln)**2*LOG(rhoaln) + &571 &(-8.8474102392445200E-10)*tln**3*LOG(satratln)**2*LOG(rhoaln) + &572 &(-1.6640566678168968E-4*LOG(satratln)**2*LOG(rhoaln))/x_n - 7.7390093776705471E-2*LOG(rhoaln)**2 + &573 &5.8220163188828482E-4*tln*LOG(rhoaln)**2 + 1.2291679321523287E-6*tln**2*LOG(rhoaln)**2 + &574 &(-7.4690997508075749E-9)*tln**3*LOG(rhoaln)**2 + (-5.6357941220497648E-3*LOG(rhoaln)**2)/x_n + &575 &(-4.7170109625089768E-3)*LOG(satratln)*LOG(rhoaln)**2 + 6.9828868534370193E-5*tln*LOG(satratln)*LOG(rhoaln)**2 + &576 &(-3.1738912157036403E-7)*tln**2*LOG(satratln)*LOG(rhoaln)**2 + &577 &2.3975538706787416E-10*tln**3*LOG(satratln)*LOG(rhoaln)**2 + &578 &(4.2304213386288567E-4*LOG(satratln)*LOG(rhoaln)**2)/x_n + 1.3696520973423231E-3*LOG(rhoaln)**3 + &579 &(-1.6863387574788199E-5)*tln*LOG(rhoaln)**3 + 2.7959499278844516E-8*tln**2*LOG(rhoaln)**3 + &580 &3.9423927013227455E-11*tln**3*LOG(rhoaln)**3 + (8.6136359966337272E-5*LOG(rhoaln)**3)/x_n556 - 1.2080330016937095E-1/x_n + 1.1902674923928015E-3*LOG(satratln) - 1.9211358507172177E-2*tln*LOG(satratln) + & 557 2.4648094311204255E-4*tln**2*LOG(satratln) - 7.5641448594711666E-7*tln**3*LOG(satratln) + & 558 (-2.0668639384228818E-02*LOG(satratln))/x_n - 3.7593072011595188E-2*LOG(satratln)**2 + & 559 9.0993182774415718E-4 *tln*LOG(satratln)**2 + & 560 (-9.5698412164297149E-6)*tln**2*LOG(satratln)**2 + 3.7163166416110421E-8*tln**3*LOG(satratln)**2 + & 561 (1.1026579525210847E-2*LOG(satratln)**2)/x_n + 1.1530844115561925E-2 *LOG(satratln)**3 + & 562 (-1.8083253906466668E-4)*tln*LOG(satratln)**3 + 8.0213604053330654E-7*tln**2*LOG(satratln)**3 + & 563 (-8.5797885383051337E-10)*tln**3*LOG(satratln)**3 + (1.0243693899717402E-3*LOG(satratln)**3)/x_n + & 564 (-1.7248695296299649E-2)*LOG(rhoaln) + 1.1294004162437157E-2*tln*LOG(rhoaln) + & 565 (-1.2283640163189278E-4)*tln**2*LOG(rhoaln) + 2.7391732258259009E-7*tln**3*LOG(rhoaln) + & 566 (6.8505583974029602E-2*LOG(rhoaln))/x_n +2.9750968179523635E-1*LOG(satratln)*LOG(rhoaln) + & 567 (-3.6681154503992296E-3)*tln*LOG(satratln)*LOG(rhoaln) + 1.0636473034653114E-5*tln**2*LOG(satratln)*LOG(rhoaln)+ & 568 5.8687098466515866E-9*tln**3*LOG(satratln)*LOG(rhoaln) + (-5.2028866094191509E-3*LOG(satratln)*LOG(rhoaln))/x_n+ & 569 7.6971988880587231E-4*LOG(satratln)**2*LOG(rhoaln) - 2.4605575820433763E-5*tln*LOG(satratln)**2*LOG(rhoaln) + & 570 2.3818484400893008E-7*tln**2*LOG(satratln)**2*LOG(rhoaln) + & 571 (-8.8474102392445200E-10)*tln**3*LOG(satratln)**2*LOG(rhoaln) + & 572 (-1.6640566678168968E-4*LOG(satratln)**2*LOG(rhoaln))/x_n - 7.7390093776705471E-2*LOG(rhoaln)**2 + & 573 5.8220163188828482E-4*tln*LOG(rhoaln)**2 + 1.2291679321523287E-6*tln**2*LOG(rhoaln)**2 + & 574 (-7.4690997508075749E-9)*tln**3*LOG(rhoaln)**2 + (-5.6357941220497648E-3*LOG(rhoaln)**2)/x_n + & 575 (-4.7170109625089768E-3)*LOG(satratln)*LOG(rhoaln)**2 + 6.9828868534370193E-5*tln*LOG(satratln)*LOG(rhoaln)**2 + & 576 (-3.1738912157036403E-7)*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & 577 2.3975538706787416E-10*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & 578 (4.2304213386288567E-4*LOG(satratln)*LOG(rhoaln)**2)/x_n + 1.3696520973423231E-3*LOG(rhoaln)**3 + & 579 (-1.6863387574788199E-5)*tln*LOG(rhoaln)**3 + 2.7959499278844516E-8*tln**2*LOG(rhoaln)**3 + & 580 3.9423927013227455E-11*tln**3*LOG(rhoaln)**3 + (8.6136359966337272E-5*LOG(rhoaln)**3)/x_n 581 581 ntot_n=EXP(ntot_n) 582 582 … … 603 603 604 604 kinrhotresi = 5.3742280876674478E1 - 6.6837931590012266E-3 *LOG(satratli)**(-2) & 605 &- 1.0142598385422842E-01 * LOG(satratli)**(-1) - 6.4170597272606873E+00 * LOG(satratli) &606 &- 6.4315798914824518E-01 * LOG(satratli)**2 - 2.4428391714772721E-02 * LOG(satratli)**3 &607 &- 3.5356658734539019E-04 * LOG(satratli)**4 + 2.5400015099140506E-05 * tli * LOG(satratli)**(-2) &608 &- 2.7928900816637790E-04 * tli * LOG(satratli)**(-1) + 4.4108573484923690E-02 * tli * LOG(satratli) &609 &+ 6.3943789012475532E-03 * tli * LOG(satratli)**(2) + 2.3164296174966580E-04 * tli * LOG(satratli)**(3) &610 &+ 3.0372070669934950E-06 * tli * LOG(satratli)**4 + 3.8255873977423475E-06 * tli**2 * LOG(satratli)**(-1) &611 &- 1.2344793083561629E-04 * tli**2 * LOG(satratli) - 1.7959048869810192E-05 * tli**2 * LOG(satratli)**(2) &612 &- 3.2165622558722767E-07 * tli**2 * LOG(satratli)**3 - 4.7136923780988659E-09 * tli**3 * LOG(satratli)**(-1) &613 &+ 1.1873317184482216E-07 * tli**3 * LOG(satratli) + 1.5685860354866621E-08 * tli**3 * LOG(satratli)**2 &614 &- 1.4329645891059557E+04 * tli**(-1) + 1.3842599842575321E-01 * tli &615 &- 4.1376265912842938E-04 * tli**(2) + 3.9147639775826004E-07 * tli**3605 - 1.0142598385422842E-01 * LOG(satratli)**(-1) - 6.4170597272606873E+00 * LOG(satratli) & 606 - 6.4315798914824518E-01 * LOG(satratli)**2 - 2.4428391714772721E-02 * LOG(satratli)**3 & 607 - 3.5356658734539019E-04 * LOG(satratli)**4 + 2.5400015099140506E-05 * tli * LOG(satratli)**(-2) & 608 - 2.7928900816637790E-04 * tli * LOG(satratli)**(-1) + 4.4108573484923690E-02 * tli * LOG(satratli) & 609 + 6.3943789012475532E-03 * tli * LOG(satratli)**(2) + 2.3164296174966580E-04 * tli * LOG(satratli)**(3) & 610 + 3.0372070669934950E-06 * tli * LOG(satratli)**4 + 3.8255873977423475E-06 * tli**2 * LOG(satratli)**(-1) & 611 - 1.2344793083561629E-04 * tli**2 * LOG(satratli) - 1.7959048869810192E-05 * tli**2 * LOG(satratli)**(2) & 612 - 3.2165622558722767E-07 * tli**2 * LOG(satratli)**3 - 4.7136923780988659E-09 * tli**3 * LOG(satratli)**(-1) & 613 + 1.1873317184482216E-07 * tli**3 * LOG(satratli) + 1.5685860354866621E-08 * tli**3 * LOG(satratli)**2 & 614 - 1.4329645891059557E+04 * tli**(-1) + 1.3842599842575321E-01 * tli & 615 - 4.1376265912842938E-04 * tli**(2) + 3.9147639775826004E-07 * tli**3 616 616 617 617 kinrhotresi=EXP(kinrhotresi) !1/cm3 … … 621 621 IF (kinetic_i) THEN 622 622 jnuc_i1=1.0E6*(0.3E-9 + 0.487E-9)**2.*SQRT(8.*RPI*RKBOL*(1./mH2SO4mol+1./mH2SO4mol))* & 623 & SQRT(tli)*rhoali !1/cm3s623 SQRT(tli)*rhoali !1/cm3s 624 624 ntot_i=1. !set to 1 625 625 na_i=1. … … 628 628 ELSE 629 629 jnuc_i1 = 3.0108954259038608E+01+tli*6.1176722090512577E+01+(tli**2)*8.7240333618891663E-01+(tli**3)* & 630 &(-4.6191788649375719E-03)+(tli**(-1))*8.3537059107024481E-01 + &631 &(1.5028549216690628E+01+tli*(-1.9310989753720623E-01)+(tli**2)*8.0155514634860480E-04+(tli**3)* &632 &(-1.0832730707799128E-06)+(tli**(-1))*1.7577660457989019)*(LOG(satratli)**(-2)) + &633 &(-2.0487870170216488E-01 + tli * 1.3263949252910405E-03 + (tli**2) *(-8.4195688402450274E-06) + &634 &(tli**3)*1.6154895940993287E-08 + (tli**(-1))*3.8734212545203874E+01) * (LOG(satratli)**(-2)*LOG(rhoali)) + &635 &(1.4955918863858371 + tli * 9.2290004245522454E+01 + (tli**2) *(-8.9006965195392618E-01) + &636 &(tli**3) * 2.2319123411013099E-03 + (tli**(-1)) * 4.0180079996840852E-03) * &637 &(LOG(satratli)**(-1) * LOG(rhoali)**(-1)) + &638 &(7.9018031228561085 + tli *(-1.1649433968658949E+01) + (tli**2) * 1.1400827854910951E-01 + &639 &(tli**3) *(-3.1941526492127755E-04) + (tli**(-1)) *(-3.7662115740271446E-01)) * (LOG(satratli)**(-1)) + &640 &(1.5725237111225979E+02 + tli *(-1.0051649979836277) + (tli**2) * 1.1866484014507624E-03 + &641 &(tli**3) * 7.3557614998540389E-06 + (tli**(-1)) * 2.6270197023115189) * (LOG(satratli)**(-1) * LOG(rhoali)) + &642 &(-1.6973840122470968E+01 + tli * 1.1258423691432135E-01 + (tli**2) *(-2.9850139351463793E-04) + (tli**3) * &643 &1.4301286324827064E-07 + (tli**(-1)) * 1.3163389235253725E+01) * (LOG(satratli)**(-1) * LOG(rhoali)**2) + &644 &(-1.0399591631839757 + tli * 2.7022055588257691E-03 + (tli**2) *(-2.1507467231330936E-06) + (tli**3) * &645 &3.8059489037584171E-10 + (tli**(-1)) * 1.5000492788553410E+02) * (LOG(satratli)**(-1) * LOG(rhoali)**3) + &646 &(1.2250990965305315 + tli * 3.0495946490079444E+01 + (tli**2) * 2.1051563135187106E+01 + (tli**3) * &647 &(-8.2200682916580878E-02) + (tli**(-1)) * 2.9965871386685029E-02) * (LOG(rhoali)**(-2)) + &648 &(4.8281605955680433 + tli * 1.7346551710836445E+02 + (tli**2) *(-1.0113602140796010E+01) + (tli**3) * &649 &3.7482518458685089E-02 + (tli**(-1)) *(-1.4449998158558205E-01)) * (LOG(rhoali)**(-1)) + &650 &(2.3399230964451237E+02 + tli *(-2.3099267235261948E+01) + (tli**2) * 8.0122962140916354E-02 + &651 &(tli**3) * 6.1542576994557088E-05 + (tli**(-1)) * 5.3718413254843007) * (LOG(rhoali)) + &652 &(1.0299715519499360E+02 + tli *(-6.4663357203364136E-02) + (tli**2) *(-2.0487150565050316E-03) + &653 &(tli**3) * 8.7935289055530897E-07 + (tli**(-1)) * 3.6013204601215229E+01) * (LOG(rhoali)**2) + &654 &(-3.5452115439584042 + tli * 1.7083445731159330E-02 + (tli**2) *(-1.2552625290862626E-05) + (tli**3) * &655 &1.2968447449182847E-09 + (tli**(-1)) * 1.5748687512056560E+02) * (LOG(rhoali)**3) + &656 &(2.2338490119517975 + tli * 1.0229410216045540E+02 + (tli**2) *(-3.2103611955174052) + (tli**3) * &657 &1.3397152304977591E-02 + (tli**(-1)) *(-2.4155187776460030E-02)) * (LOG(satratli)* LOG(rhoali)**(-2)) + &658 &(3.7592282990713963 + tli *(-1.5257988769009816E+02) + (tli**2) * 2.6113805420558802 + (tli**3) * &659 &(-9.0380721653694363E-03) + (tli**(-1)) *(-1.3974197138171082E-01)) * (LOG(satratli)* LOG(rhoali)**(-1)) + &660 &(1.8293600730573988E+01 + tli * 1.8344728606002992E+01 + (tli**2) *(-4.0063363221106751E-01) + (tli**3) &661 &* 1.4842749371258522E-03 + (tli**(-1)) * 1.1848846003282287) * (LOG(satratli)) + &662 &(-1.7634531623032314E+02 + tli * 4.9011762441271278 + (tli**2) *(-1.3195821562746339E-02) + (tli**3) * &663 &(-2.8668619526430859E-05) + (tli**(-1)) *(-2.9823396976393551E-01)) * (LOG(satratli)* LOG(rhoali)) + &664 &(-3.2944043694275727E+01 + tli * 1.2517571921051887E-01 + (tli**2) * 8.3239769771186714E-05 + (tli**3) * &665 &2.8191859341519507E-07 + (tli**(-1)) *(-2.7352880736682319E+01)) * (LOG(satratli)* LOG(rhoali)**2) + &666 &(-1.1451811137553243 + tli * 2.0625997485732494E-03 + (tli**2) *(-3.4225389469233624E-06) + (tli**3) * &667 &4.4437613496984567E-10 + (tli**(-1)) * 1.8666644332606754E+02) * (LOG(satratli)* LOG(rhoali)**3) + &668 &(3.2270897099493567E+01 + tli * 7.7898447327513687E-01 + (tli**2) *(-6.5662738484679626E-03) + (tli**3) * &669 &3.7899330796456790E-06 + (tli**(-1)) * 7.1106427501756542E-01) * (LOG(satratli)**2 * LOG(rhoali)**(-1)) + &670 &(-2.8901906781697811E+01 + tli *(-1.5356398793054860) + (tli**2) * 1.9267271774384788E-02 + (tli**3) * &671 &(-5.3886270475516162E-05) + (tli**(-1)) * 5.0490415975693426E-01) * (LOG(satratli)**2) + &672 &(3.3365683645733924E+01 + tli *(-3.6114561564894537E-01) + (tli**2) * 9.2977354471929262E-04 + (tli**3) * &673 &1.9549769069511355E-07 + (tli**(-1)) *(-8.8865930095112855)) * (LOG(satratli)**2 * LOG(rhoali)) + &674 &(2.4592563042806375 + tli *(-8.3227071743101084E-03) + (tli**2) * 8.2563338043447783E-06 + (tli**3) * &675 &(-8.4374976698593496E-09) + (tli**(-1)) *(-2.0938173949893473E+02)) * (LOG(satratli)**2 * LOG(rhoali)**2) + &676 &(4.4099823444352317E+01 + tli * 2.5915665826835252 + (tli**2) *(-1.6449091819482634E-02) + (tli**3) * &677 &2.6797249816144721E-05 + (tli**(-1)) * 5.5045672663909995E-01)* satratli630 (-4.6191788649375719E-03)+(tli**(-1))*8.3537059107024481E-01 + & 631 (1.5028549216690628E+01+tli*(-1.9310989753720623E-01)+(tli**2)*8.0155514634860480E-04+(tli**3)* & 632 (-1.0832730707799128E-06)+(tli**(-1))*1.7577660457989019)*(LOG(satratli)**(-2)) + & 633 (-2.0487870170216488E-01 + tli * 1.3263949252910405E-03 + (tli**2) *(-8.4195688402450274E-06) + & 634 (tli**3)*1.6154895940993287E-08 + (tli**(-1))*3.8734212545203874E+01) * (LOG(satratli)**(-2)*LOG(rhoali)) + & 635 (1.4955918863858371 + tli * 9.2290004245522454E+01 + (tli**2) *(-8.9006965195392618E-01) + & 636 (tli**3) * 2.2319123411013099E-03 + (tli**(-1)) * 4.0180079996840852E-03) * & 637 (LOG(satratli)**(-1) * LOG(rhoali)**(-1)) + & 638 (7.9018031228561085 + tli *(-1.1649433968658949E+01) + (tli**2) * 1.1400827854910951E-01 + & 639 (tli**3) *(-3.1941526492127755E-04) + (tli**(-1)) *(-3.7662115740271446E-01)) * (LOG(satratli)**(-1)) + & 640 (1.5725237111225979E+02 + tli *(-1.0051649979836277) + (tli**2) * 1.1866484014507624E-03 + & 641 (tli**3) * 7.3557614998540389E-06 + (tli**(-1)) * 2.6270197023115189) * (LOG(satratli)**(-1) * LOG(rhoali)) + & 642 (-1.6973840122470968E+01 + tli * 1.1258423691432135E-01 + (tli**2) *(-2.9850139351463793E-04) + (tli**3) * & 643 1.4301286324827064E-07 + (tli**(-1)) * 1.3163389235253725E+01) * (LOG(satratli)**(-1) * LOG(rhoali)**2) + & 644 (-1.0399591631839757 + tli * 2.7022055588257691E-03 + (tli**2) *(-2.1507467231330936E-06) + (tli**3) * & 645 3.8059489037584171E-10 + (tli**(-1)) * 1.5000492788553410E+02) * (LOG(satratli)**(-1) * LOG(rhoali)**3) + & 646 (1.2250990965305315 + tli * 3.0495946490079444E+01 + (tli**2) * 2.1051563135187106E+01 + (tli**3) * & 647 (-8.2200682916580878E-02) + (tli**(-1)) * 2.9965871386685029E-02) * (LOG(rhoali)**(-2)) + & 648 (4.8281605955680433 + tli * 1.7346551710836445E+02 + (tli**2) *(-1.0113602140796010E+01) + (tli**3) * & 649 3.7482518458685089E-02 + (tli**(-1)) *(-1.4449998158558205E-01)) * (LOG(rhoali)**(-1)) + & 650 (2.3399230964451237E+02 + tli *(-2.3099267235261948E+01) + (tli**2) * 8.0122962140916354E-02 + & 651 (tli**3) * 6.1542576994557088E-05 + (tli**(-1)) * 5.3718413254843007) * (LOG(rhoali)) + & 652 (1.0299715519499360E+02 + tli *(-6.4663357203364136E-02) + (tli**2) *(-2.0487150565050316E-03) + & 653 (tli**3) * 8.7935289055530897E-07 + (tli**(-1)) * 3.6013204601215229E+01) * (LOG(rhoali)**2) + & 654 (-3.5452115439584042 + tli * 1.7083445731159330E-02 + (tli**2) *(-1.2552625290862626E-05) + (tli**3) * & 655 1.2968447449182847E-09 + (tli**(-1)) * 1.5748687512056560E+02) * (LOG(rhoali)**3) + & 656 (2.2338490119517975 + tli * 1.0229410216045540E+02 + (tli**2) *(-3.2103611955174052) + (tli**3) * & 657 1.3397152304977591E-02 + (tli**(-1)) *(-2.4155187776460030E-02)) * (LOG(satratli)* LOG(rhoali)**(-2)) + & 658 (3.7592282990713963 + tli *(-1.5257988769009816E+02) + (tli**2) * 2.6113805420558802 + (tli**3) * & 659 (-9.0380721653694363E-03) + (tli**(-1)) *(-1.3974197138171082E-01)) * (LOG(satratli)* LOG(rhoali)**(-1)) + & 660 (1.8293600730573988E+01 + tli * 1.8344728606002992E+01 + (tli**2) *(-4.0063363221106751E-01) + (tli**3) & 661 * 1.4842749371258522E-03 + (tli**(-1)) * 1.1848846003282287) * (LOG(satratli)) + & 662 (-1.7634531623032314E+02 + tli * 4.9011762441271278 + (tli**2) *(-1.3195821562746339E-02) + (tli**3) * & 663 (-2.8668619526430859E-05) + (tli**(-1)) *(-2.9823396976393551E-01)) * (LOG(satratli)* LOG(rhoali)) + & 664 (-3.2944043694275727E+01 + tli * 1.2517571921051887E-01 + (tli**2) * 8.3239769771186714E-05 + (tli**3) * & 665 2.8191859341519507E-07 + (tli**(-1)) *(-2.7352880736682319E+01)) * (LOG(satratli)* LOG(rhoali)**2) + & 666 (-1.1451811137553243 + tli * 2.0625997485732494E-03 + (tli**2) *(-3.4225389469233624E-06) + (tli**3) * & 667 4.4437613496984567E-10 + (tli**(-1)) * 1.8666644332606754E+02) * (LOG(satratli)* LOG(rhoali)**3) + & 668 (3.2270897099493567E+01 + tli * 7.7898447327513687E-01 + (tli**2) *(-6.5662738484679626E-03) + (tli**3) * & 669 3.7899330796456790E-06 + (tli**(-1)) * 7.1106427501756542E-01) * (LOG(satratli)**2 * LOG(rhoali)**(-1)) + & 670 (-2.8901906781697811E+01 + tli *(-1.5356398793054860) + (tli**2) * 1.9267271774384788E-02 + (tli**3) * & 671 (-5.3886270475516162E-05) + (tli**(-1)) * 5.0490415975693426E-01) * (LOG(satratli)**2) + & 672 (3.3365683645733924E+01 + tli *(-3.6114561564894537E-01) + (tli**2) * 9.2977354471929262E-04 + (tli**3) * & 673 1.9549769069511355E-07 + (tli**(-1)) *(-8.8865930095112855)) * (LOG(satratli)**2 * LOG(rhoali)) + & 674 (2.4592563042806375 + tli *(-8.3227071743101084E-03) + (tli**2) * 8.2563338043447783E-06 + (tli**3) * & 675 (-8.4374976698593496E-09) + (tli**(-1)) *(-2.0938173949893473E+02)) * (LOG(satratli)**2 * LOG(rhoali)**2) + & 676 (4.4099823444352317E+01 + tli * 2.5915665826835252 + (tli**2) *(-1.6449091819482634E-02) + (tli**3) * & 677 2.6797249816144721E-05 + (tli**(-1)) * 5.5045672663909995E-01)* satratli 678 678 jnuc_i1=EXP(jnuc_i1) 679 679 680 680 ntot_i = ABS((-4.8324296064013375E+04 + tli * 5.0469120697428906E+02 + (tli**2) *(-1.1528940488496042E+00) + & 681 &(tli**(-1)) *(-8.6892744676239192E+02) + (tli**(3)) * 4.0030302028120469E-04) + &682 &(-6.7259105232039847E+03 + tli * 1.9197488157452008E+02 + (tli**2) *(-1.3602976930126354E+00) + &683 &(tli**(-1)) *(-1.1212637938360332E+02) + (tli**(3)) * 2.8515597265933207E-03) * &684 &LOG(satratli)**(-2) * LOG(rhoali)**(-2) + &685 &(2.6216455217763342E+02 + tli *(-2.3687553252750821E+00) + (tli**2) * 7.4074554767517521E-03 + &686 &(tli**(-1)) *(-1.9213956820114927E+03) + (tli**(3)) *(-9.3839114856129453E-06)) * LOG(satratli)**(-2) + &687 &(3.9652478944137344E+00 + tli * 1.2469375098256536E-02 + (tli**2) *(-9.9837754694045633E-05) + (tli**(-1)) * &688 &(-5.1919499210175138E+02) + (tli**(3)) * 1.6489001324583862E-07) * LOG(satratli)**(-2) * LOG(rhoali) + &689 &(2.4975714429096206E+02 + tli * 1.7107594562445172E+02 + (tli**2) *(-7.8988711365135289E-01) + (tli**(-1)) * &690 &(-2.2243599782483177E+01) + (tli**(3)) *(-1.6291523004095427E-04)) * LOG(satratli)**(-1) * LOG(rhoali)**(-2) +&691 &(-8.9270715592533611E+02 + tli * 1.2053538883338946E+02 + (tli**2) *(-1.5490408828541018E+00) + (tli**(-1)) * &692 &(-1.1243275579419826E+01) + (tli**(3)) * 4.8053105606904655E-03) * LOG(satratli)**(-1) * LOG(rhoali)**(-1) + &693 &(7.6426441642091631E+03 + tli *(-7.1785462414656578E+01) + (tli**2) * 2.3851864923199523E-01 + (tli**(-1)) * &694 &8.5591775688708395E+01 + (tli**(3)) *(-3.7000473243342858E-04)) * LOG(satratli)**(-1) + &695 &(-5.1516826398607911E+01 + tli * 9.1385720811460558E-01 + (tli**2) *(-3.5477100262158974E-03) + &696 &(tli**(-1)) * 2.7545544507625586E+03 + (tli**(3)) * 5.4708262093640928E-06) * LOG(satratli)**(-1) * LOG(rhoali) + &697 &(-3.0386767129196176E+02 + tli *(-1.1033438883583569E+04) + (tli**2) * 8.1296859732896067E+01 + (tli**(-1)) * &698 &1.2625883141097162E+01 + (tli**(3)) *(-1.2728497822219101E-01)) * LOG(rhoali)**(-2) + &699 &(-3.3763494256461472E+03 + tli * 3.1916579136391006E+03 + (tli**2) *(-2.7234339474441143E+01) + (tli**(-1)) * &700 &(-2.1897653262707397E+01) + (tli**(3)) * 5.1788505812259071E-02) * LOG(rhoali)**(-1) + &701 &(-1.8817843873687068E+03 + tli * 4.3038072285882070E+00 + (tli**2) * 6.6244087689671860E-03 + (tli**(-1)) * &702 &(-2.7133073605696295E+03) + (tli**(3)) *(-1.7951557394285043E-05)) * LOG(rhoali) + &703 &(-1.7668827539244447E+02 + tli * 4.8160932330629913E-01 + (tli**2) *(-6.3133007671100293E-04) + (tli**(-1)) * &704 &2.5631774669873157E+04 + (tli**(3)) * 4.1534484127873519E-07) * LOG(rhoali)**(2) + &705 &(-1.6661835889222382E+03 + tli * 1.3708900504682877E+03 + (tli**2) *(-1.7919060052198969E+01) + (tli**(-1)) * &706 &(-3.5145029804436405E+01) + (tli**(3)) * 5.1047240947371224E-02) * LOG(satratli)* LOG(rhoali)**(-2) + &707 &(1.0843549363030939E+04 + tli *(-7.3557073636139577E+01) + (tli**2) * 1.2054625131778862E+00 + (tli**(-1)) * &708 &1.9358737917864391E+02 + (tli**(3)) *(-4.2871620775911338E-03)) * LOG(satratli)* LOG(rhoali)**(-1) + &709 &(-2.4269802549752835E+03 + tli * 1.1348265061941714E+01 + (tli**2) *(-5.0430423939495157E-02) + (tli**(-1)) * &710 &2.3709874548950634E+03 + (tli**(3)) * 1.4091851828620244E-04) * LOG(satratli) + &711 &(5.2745372575251588E+02 + tli *(-2.6080675912627314E+00) + (tli**2) * 5.6902218056670145E-03 + (tli**(-1)) * &712 &(-3.2149319482897838E+04) + (tli**(3)) *(-5.4121996056745853E-06)) * LOG(satratli)* LOG(rhoali) + &713 &(-1.6401959518360403E+01 + tli * 2.4322962162439640E-01 + (tli**2) * 1.1744366627725344E-03 + (tli**(-1)) * &714 &(-8.2694427518413195E+03) + (tli**(3)) *(-5.0028379203873102E-06))* LOG(satratli)**(2) + &715 &(-2.7556572017167782E+03 + tli * 4.9293344495058264E+01 + (tli**2) *(-2.6503456520676050E-01) + (tli**(-1)) * &716 &1.2130698030982167E+03 + (tli**(3)) * 4.3530610668042957E-04)* LOG(satratli)**2 * LOG(rhoali)**(-1) + &717 &(-6.3419182228959192E+00 + tli * 4.0636212834605827E-02 + (tli**2) *(-1.0450112687842742E-04) + (tli**(-1)) * &718 &3.1035882189759656E+02 + (tli**(3)) * 9.4328418657873500E-08)* LOG(satratli)**(-3) + &719 &(3.0189213304689042E+03 + tli *(-2.3804654203861684E+01) + (tli**2) * 6.8113013411972942E-02 + (tli**(-1)) * &720 &6.3112071081188913E+02 + (tli**(3)) *(-9.4460854261685723E-05))* (satratli) * LOG(rhoali) + &721 &(1.1924791930673702E+04 + tli *(-1.1973824959206000E+02) + (tli**2) * 1.6888713097971020E-01 + (tli**(-1)) * &722 &1.8735938211539585E+02 + (tli**(3)) * 5.0974564680442852E-04)* (satratli) + &723 &(3.6409071302482083E+01 + tli * 1.7919859306449623E-01 + (tli**2) *(-1.0020116255895206E-03) + (tli**(-1)) * &724 &(-8.3521083354432303E+03) + (tli**(3)) * 1.5879900546795635E-06)* satratli * LOG(rhoali)**(2))681 (tli**(-1)) *(-8.6892744676239192E+02) + (tli**(3)) * 4.0030302028120469E-04) + & 682 (-6.7259105232039847E+03 + tli * 1.9197488157452008E+02 + (tli**2) *(-1.3602976930126354E+00) + & 683 (tli**(-1)) *(-1.1212637938360332E+02) + (tli**(3)) * 2.8515597265933207E-03) * & 684 LOG(satratli)**(-2) * LOG(rhoali)**(-2) + & 685 (2.6216455217763342E+02 + tli *(-2.3687553252750821E+00) + (tli**2) * 7.4074554767517521E-03 + & 686 (tli**(-1)) *(-1.9213956820114927E+03) + (tli**(3)) *(-9.3839114856129453E-06)) * LOG(satratli)**(-2) + & 687 (3.9652478944137344E+00 + tli * 1.2469375098256536E-02 + (tli**2) *(-9.9837754694045633E-05) + (tli**(-1)) * & 688 (-5.1919499210175138E+02) + (tli**(3)) * 1.6489001324583862E-07) * LOG(satratli)**(-2) * LOG(rhoali) + & 689 (2.4975714429096206E+02 + tli * 1.7107594562445172E+02 + (tli**2) *(-7.8988711365135289E-01) + (tli**(-1)) * & 690 (-2.2243599782483177E+01) + (tli**(3)) *(-1.6291523004095427E-04)) * LOG(satratli)**(-1) * LOG(rhoali)**(-2) +& 691 (-8.9270715592533611E+02 + tli * 1.2053538883338946E+02 + (tli**2) *(-1.5490408828541018E+00) + (tli**(-1)) * & 692 (-1.1243275579419826E+01) + (tli**(3)) * 4.8053105606904655E-03) * LOG(satratli)**(-1) * LOG(rhoali)**(-1) + & 693 (7.6426441642091631E+03 + tli *(-7.1785462414656578E+01) + (tli**2) * 2.3851864923199523E-01 + (tli**(-1)) * & 694 8.5591775688708395E+01 + (tli**(3)) *(-3.7000473243342858E-04)) * LOG(satratli)**(-1) + & 695 (-5.1516826398607911E+01 + tli * 9.1385720811460558E-01 + (tli**2) *(-3.5477100262158974E-03) + & 696 (tli**(-1)) * 2.7545544507625586E+03 + (tli**(3)) * 5.4708262093640928E-06) * LOG(satratli)**(-1) * LOG(rhoali) + & 697 (-3.0386767129196176E+02 + tli *(-1.1033438883583569E+04) + (tli**2) * 8.1296859732896067E+01 + (tli**(-1)) * & 698 1.2625883141097162E+01 + (tli**(3)) *(-1.2728497822219101E-01)) * LOG(rhoali)**(-2) + & 699 (-3.3763494256461472E+03 + tli * 3.1916579136391006E+03 + (tli**2) *(-2.7234339474441143E+01) + (tli**(-1)) * & 700 (-2.1897653262707397E+01) + (tli**(3)) * 5.1788505812259071E-02) * LOG(rhoali)**(-1) + & 701 (-1.8817843873687068E+03 + tli * 4.3038072285882070E+00 + (tli**2) * 6.6244087689671860E-03 + (tli**(-1)) * & 702 (-2.7133073605696295E+03) + (tli**(3)) *(-1.7951557394285043E-05)) * LOG(rhoali) + & 703 (-1.7668827539244447E+02 + tli * 4.8160932330629913E-01 + (tli**2) *(-6.3133007671100293E-04) + (tli**(-1)) * & 704 2.5631774669873157E+04 + (tli**(3)) * 4.1534484127873519E-07) * LOG(rhoali)**(2) + & 705 (-1.6661835889222382E+03 + tli * 1.3708900504682877E+03 + (tli**2) *(-1.7919060052198969E+01) + (tli**(-1)) * & 706 (-3.5145029804436405E+01) + (tli**(3)) * 5.1047240947371224E-02) * LOG(satratli)* LOG(rhoali)**(-2) + & 707 (1.0843549363030939E+04 + tli *(-7.3557073636139577E+01) + (tli**2) * 1.2054625131778862E+00 + (tli**(-1)) * & 708 1.9358737917864391E+02 + (tli**(3)) *(-4.2871620775911338E-03)) * LOG(satratli)* LOG(rhoali)**(-1) + & 709 (-2.4269802549752835E+03 + tli * 1.1348265061941714E+01 + (tli**2) *(-5.0430423939495157E-02) + (tli**(-1)) * & 710 2.3709874548950634E+03 + (tli**(3)) * 1.4091851828620244E-04) * LOG(satratli) + & 711 (5.2745372575251588E+02 + tli *(-2.6080675912627314E+00) + (tli**2) * 5.6902218056670145E-03 + (tli**(-1)) * & 712 (-3.2149319482897838E+04) + (tli**(3)) *(-5.4121996056745853E-06)) * LOG(satratli)* LOG(rhoali) + & 713 (-1.6401959518360403E+01 + tli * 2.4322962162439640E-01 + (tli**2) * 1.1744366627725344E-03 + (tli**(-1)) * & 714 (-8.2694427518413195E+03) + (tli**(3)) *(-5.0028379203873102E-06))* LOG(satratli)**(2) + & 715 (-2.7556572017167782E+03 + tli * 4.9293344495058264E+01 + (tli**2) *(-2.6503456520676050E-01) + (tli**(-1)) * & 716 1.2130698030982167E+03 + (tli**(3)) * 4.3530610668042957E-04)* LOG(satratli)**2 * LOG(rhoali)**(-1) + & 717 (-6.3419182228959192E+00 + tli * 4.0636212834605827E-02 + (tli**2) *(-1.0450112687842742E-04) + (tli**(-1)) * & 718 3.1035882189759656E+02 + (tli**(3)) * 9.4328418657873500E-08)* LOG(satratli)**(-3) + & 719 (3.0189213304689042E+03 + tli *(-2.3804654203861684E+01) + (tli**2) * 6.8113013411972942E-02 + (tli**(-1)) * & 720 6.3112071081188913E+02 + (tli**(3)) *(-9.4460854261685723E-05))* (satratli) * LOG(rhoali) + & 721 (1.1924791930673702E+04 + tli *(-1.1973824959206000E+02) + (tli**2) * 1.6888713097971020E-01 + (tli**(-1)) * & 722 1.8735938211539585E+02 + (tli**(3)) * 5.0974564680442852E-04)* (satratli) + & 723 (3.6409071302482083E+01 + tli * 1.7919859306449623E-01 + (tli**2) *(-1.0020116255895206E-03) + (tli**(-1)) * & 724 (-8.3521083354432303E+03) + (tli**(3)) * 1.5879900546795635E-06)* satratli * LOG(rhoali)**(2)) 725 725 726 726 rc_i = (-3.6318550637865524E-08 + tli * 2.1740704135789128E-09 + (tli**2) * & 727 &(-8.5521429066506161E-12) + (tli**3) *(-9.3538647454573390E-15)) + &728 &(2.1366936839394922E-08 + tli *(-2.4087168827395623E-10) + (tli**2) * 8.7969869277074319E-13 + &729 &(tli**3) *(-1.0294466881303291E-15))* LOG(satratli)**(-2) * LOG(rhoali)**(-1) + &730 &(-7.7804007761164303E-10 + tli * 1.0327058173517932E-11 + (tli**2) *(-4.2557697639692428E-14) + &731 &(tli**3) * 5.4082507061618662E-17)* LOG(satratli)**(-2) + &732 &(3.2628927397420860E-12 + tli *(-7.6475692919751066E-14) + (tli**2) * 4.1985816845259788E-16 + &733 &(tli**3) *(-6.2281395889592719E-19))* LOG(satratli)**(-2) * LOG(rhoali) + &734 &(2.0442205540818555E-09 + tli * 4.0441858911249830E-08 + (tli**2) *(-3.3423487629482825E-10) + &735 &(tli**3) * 6.8000404742985678E-13)* LOG(satratli)**(-1) * LOG(rhoali)**(-2) + &736 &(1.8381489183824627E-08 + tli *(-8.9853322951518919E-09) + (tli**2) * 7.5888799566036185E-11 + &737 &(tli**3) *(-1.5823457864755549E-13))* LOG(satratli)**(-1) * LOG(rhoali)**(-1) + &738 &(1.1795760639695057E-07 + tli *(-8.1046722896375875E-10) + (tli**2) * 9.1868604369041857E-14 + &739 &(tli**3) * 4.7882428237444610E-15)* LOG(satratli)**(-1) + &740 &(-4.4028846582545952E-09 + tli * 4.6541269232626618E-11 + (tli**2) *(-1.1939929984285194E-13) + &741 &(tli**3) * 2.3602037016614437E-17)* LOG(satratli)**(-1) * LOG(rhoali) + &742 &(2.7885056884209128E-11 + tli *(-4.5167129624119121E-13) + (tli**2) * 1.6558404997394422E-15 + &743 &(tli**3) *(-1.2037336621218054E-18))* LOG(satratli)**(-1) * LOG(rhoali)**2 + &744 &(-2.3719627171699983E-09 + tli *(-1.5260127909292053E-07) + (tli**2) * 1.7177017944754134E-09 + &745 &(tli**3) *(-4.7031737537526395E-12))* LOG(rhoali)**(-2) + &746 &(-5.6946433724699646E-09 + tli * 8.4629788237081735E-09 + (tli**2) *(-1.7674135187061521E-10) + &747 &(tli**3) * 6.6236547903091862E-13)* LOG(rhoali)**(-1) + &748 &(-2.2808617930606012E-08 + tli * 1.4773376696847775E-10 + (tli**2) *(-1.3076953119957355E-13) + &749 &(tli**3) * 2.3625301497914000E-16)* LOG(rhoali) + &750 &(1.4014269939947841E-10 + tli *(-2.3675117757377632E-12) + (tli**2) * 5.1514033966707879E-15 + &751 &(tli**3) *(-4.8864233454747856E-18))* LOG(rhoali)**2 + &752 &(6.5464943868885886E-11 + tli * 1.6494354816942769E-08 + (tli**2) *(-1.7480097393483653E-10) + &753 &(tli**3) * 4.7460075628523984E-13)* LOG(satratli)* LOG(rhoali)**(-2) + &754 &(8.4737893183927871E-09 + tli *(-6.0243327445597118E-09) + (tli**2) * 5.8766070529814883E-11 + &755 &(tli**3) *(-1.4926748560042018E-13))* LOG(satratli)* LOG(rhoali)**(-1) + &756 &(1.0761964135701397E-07 + tli *(-1.0142496009071148E-09) + (tli**2) * 2.1337312466519190E-12 + &757 &(tli**3) * 1.6376014957685404E-15)* LOG(satratli) + &758 &(-3.5621571395968670E-09 + tli * 4.1175339587760905E-11 + (tli**2) *(-1.3535372357998504E-13) + &759 &(tli**3) * 8.9334219536920720E-17)* LOG(satratli)* LOG(rhoali) + &760 &(2.0700482083136289E-11 + tli *(-3.9238944562717421E-13) + (tli**2) * 1.5850961422040196E-15 + &761 &(tli**3) *(-1.5336775610911665E-18))* LOG(satratli)* LOG(rhoali)**2 + &762 &(1.8524255464416206E-09 + tli *(-2.1959816152743264E-11) + (tli**2) *(-6.4478119501677012E-14) + &763 &(tli**3) * 5.5135243833766056E-16)* LOG(satratli)**2 * LOG(rhoali)**(-1) + &764 &(1.9349488650922679E-09 + tli *(-2.2647295919976428E-11) + (tli**2) * 9.2917479748268751E-14 + &765 &(tli**3) *(-1.2741959892173170E-16))* LOG(satratli)**2 + &766 &(2.1484978031650972E-11 + tli *(-9.3976642475838013E-14) + (tli**2) *(-4.8892738002751923E-16) + &767 &(tli**3) * 1.4676120441783832E-18)* LOG(satratli)**2 * LOG(rhoali) + &768 &(6.7565715216420310E-13 + tli *(-3.5421162549480807E-15) + (tli**2) *(-3.4201196868693569E-18) + &769 &(tli**3) * 2.2260187650412392E-20)* LOG(satratli)**3 * LOG(rhoali)727 (-8.5521429066506161E-12) + (tli**3) *(-9.3538647454573390E-15)) + & 728 (2.1366936839394922E-08 + tli *(-2.4087168827395623E-10) + (tli**2) * 8.7969869277074319E-13 + & 729 (tli**3) *(-1.0294466881303291E-15))* LOG(satratli)**(-2) * LOG(rhoali)**(-1) + & 730 (-7.7804007761164303E-10 + tli * 1.0327058173517932E-11 + (tli**2) *(-4.2557697639692428E-14) + & 731 (tli**3) * 5.4082507061618662E-17)* LOG(satratli)**(-2) + & 732 (3.2628927397420860E-12 + tli *(-7.6475692919751066E-14) + (tli**2) * 4.1985816845259788E-16 + & 733 (tli**3) *(-6.2281395889592719E-19))* LOG(satratli)**(-2) * LOG(rhoali) + & 734 (2.0442205540818555E-09 + tli * 4.0441858911249830E-08 + (tli**2) *(-3.3423487629482825E-10) + & 735 (tli**3) * 6.8000404742985678E-13)* LOG(satratli)**(-1) * LOG(rhoali)**(-2) + & 736 (1.8381489183824627E-08 + tli *(-8.9853322951518919E-09) + (tli**2) * 7.5888799566036185E-11 + & 737 (tli**3) *(-1.5823457864755549E-13))* LOG(satratli)**(-1) * LOG(rhoali)**(-1) + & 738 (1.1795760639695057E-07 + tli *(-8.1046722896375875E-10) + (tli**2) * 9.1868604369041857E-14 + & 739 (tli**3) * 4.7882428237444610E-15)* LOG(satratli)**(-1) + & 740 (-4.4028846582545952E-09 + tli * 4.6541269232626618E-11 + (tli**2) *(-1.1939929984285194E-13) + & 741 (tli**3) * 2.3602037016614437E-17)* LOG(satratli)**(-1) * LOG(rhoali) + & 742 (2.7885056884209128E-11 + tli *(-4.5167129624119121E-13) + (tli**2) * 1.6558404997394422E-15 + & 743 (tli**3) *(-1.2037336621218054E-18))* LOG(satratli)**(-1) * LOG(rhoali)**2 + & 744 (-2.3719627171699983E-09 + tli *(-1.5260127909292053E-07) + (tli**2) * 1.7177017944754134E-09 + & 745 (tli**3) *(-4.7031737537526395E-12))* LOG(rhoali)**(-2) + & 746 (-5.6946433724699646E-09 + tli * 8.4629788237081735E-09 + (tli**2) *(-1.7674135187061521E-10) + & 747 (tli**3) * 6.6236547903091862E-13)* LOG(rhoali)**(-1) + & 748 (-2.2808617930606012E-08 + tli * 1.4773376696847775E-10 + (tli**2) *(-1.3076953119957355E-13) + & 749 (tli**3) * 2.3625301497914000E-16)* LOG(rhoali) + & 750 (1.4014269939947841E-10 + tli *(-2.3675117757377632E-12) + (tli**2) * 5.1514033966707879E-15 + & 751 (tli**3) *(-4.8864233454747856E-18))* LOG(rhoali)**2 + & 752 (6.5464943868885886E-11 + tli * 1.6494354816942769E-08 + (tli**2) *(-1.7480097393483653E-10) + & 753 (tli**3) * 4.7460075628523984E-13)* LOG(satratli)* LOG(rhoali)**(-2) + & 754 (8.4737893183927871E-09 + tli *(-6.0243327445597118E-09) + (tli**2) * 5.8766070529814883E-11 + & 755 (tli**3) *(-1.4926748560042018E-13))* LOG(satratli)* LOG(rhoali)**(-1) + & 756 (1.0761964135701397E-07 + tli *(-1.0142496009071148E-09) + (tli**2) * 2.1337312466519190E-12 + & 757 (tli**3) * 1.6376014957685404E-15)* LOG(satratli) + & 758 (-3.5621571395968670E-09 + tli * 4.1175339587760905E-11 + (tli**2) *(-1.3535372357998504E-13) + & 759 (tli**3) * 8.9334219536920720E-17)* LOG(satratli)* LOG(rhoali) + & 760 (2.0700482083136289E-11 + tli *(-3.9238944562717421E-13) + (tli**2) * 1.5850961422040196E-15 + & 761 (tli**3) *(-1.5336775610911665E-18))* LOG(satratli)* LOG(rhoali)**2 + & 762 (1.8524255464416206E-09 + tli *(-2.1959816152743264E-11) + (tli**2) *(-6.4478119501677012E-14) + & 763 (tli**3) * 5.5135243833766056E-16)* LOG(satratli)**2 * LOG(rhoali)**(-1) + & 764 (1.9349488650922679E-09 + tli *(-2.2647295919976428E-11) + (tli**2) * 9.2917479748268751E-14 + & 765 (tli**3) *(-1.2741959892173170E-16))* LOG(satratli)**2 + & 766 (2.1484978031650972E-11 + tli *(-9.3976642475838013E-14) + (tli**2) *(-4.8892738002751923E-16) + & 767 (tli**3) * 1.4676120441783832E-18)* LOG(satratli)**2 * LOG(rhoali) + & 768 (6.7565715216420310E-13 + tli *(-3.5421162549480807E-15) + (tli**2) *(-3.4201196868693569E-18) + & 769 (tli**3) * 2.2260187650412392E-20)* LOG(satratli)**3 * LOG(rhoali) 770 770 771 771 na_i=x_i*ntot_i -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/so2_to_h2so4.F90
r5082 r5087 81 81 ! SO2 (molec/cm3): convert from kg/kgA 82 82 rrak1 = tr_seri(ilon,ilev,id_SO2_strat) & 83 &*pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mSO2mol83 *pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mSO2mol 84 84 85 85 IF (rrak1 >= 0.0) THEN -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratdistrib.F90
r4601 r5087 41 41 alt=altLMDz(k)+float(i_int)*(altLMDz(k+1)-altLMDz(k))/float(n_int_alt) 42 42 f_lay_emiss(k)=f_lay_emiss(k)+1./(sqrt(2.*RPI)*sigma_alt)* & 43 & exp(-0.5*((alt-altemiss)/sigma_alt)**2.)* &44 &(altLMDz(k+1)-altLMDz(k))/float(n_int_alt)43 exp(-0.5*((alt-altemiss)/sigma_alt)**2.)* & 44 (altLMDz(k+1)-altLMDz(k))/float(n_int_alt) 45 45 ENDDO 46 46 f_lay_sum=f_lay_sum+f_lay_emiss(k) … … 54 54 DO k=1, klev 55 55 f_lay_emiss(k)=max(min(altemiss+sigma_alt,altLMDz(k+1))- & 56 &max(altemiss-sigma_alt, &57 &altLMDz(k)),0.)/(2.*sigma_alt)56 max(altemiss-sigma_alt, & 57 altLMDz(k)),0.)/(2.*sigma_alt) 58 58 f_lay_sum=f_lay_sum+f_lay_emiss(k) 59 59 ENDDO -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratemit.F90
r5082 r5087 54 54 IF (is_mpi_root) THEN 55 55 WRITE(*,*) 'IN STRATEMIT: date from phys_cal_mod=',year_cur,'-',& 56 &mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o56 mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o 57 57 ENDIF 58 58 … … 75 75 76 76 IF ( xlat(i)>=latmin-dlat_loc .AND. & 77 &xlat(i)<latmax+dlat_loc .AND. &78 &xlon(i)>=lonmin-dlon .AND. &79 &xlon(i)<lonmax+dlon ) THEN77 xlat(i)<latmax+dlat_loc .AND. & 78 xlon(i)>=lonmin-dlon .AND. & 79 xlon(i)<lonmax+dlon ) THEN 80 80 ! 81 81 WRITE(*,*) 'coordinates of volcanic injection point=',& 82 &xlat(i),xlon(i),day_cur,mth_cur,year_cur82 xlat(i),xlon(i),day_cur,mth_cur,year_cur 83 83 WRITE(*,*) 'DD m_emiss_vol_daily=', & 84 &m_emiss_vol_daily84 m_emiss_vol_daily 85 85 86 86 !compute altLMDz … … 97 97 IF (flag_emit==3) then 98 98 theta=(sin(theta_max/180.*RPI)-sin(theta_min/180.*RPI))/ & 99 &(sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI))99 (sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI)) 100 100 ELSE 101 101 theta=1. … … 108 108 ! stretch emission over stretchlong period 109 109 emission=m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k)/stretchlong/ & 110 &(86400.-pdt)*theta110 (86400.-pdt)*theta 111 111 112 112 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: emission avant/apres', & 113 &'alt= ', altLMDz(k), &114 &'flh2o= ',flh2o, &115 &'id_speac= ',id_spec,id_species_total, &116 &'emission= ',emission, &117 &'pdtphys= ',pdtphys, &118 &'rapport m_emiss/m_air*f_lay= ', m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k), &119 &'stretchlong= ', stretchlong, &120 &'theta= ', theta113 'alt= ', altLMDz(k), & 114 'flh2o= ',flh2o, & 115 'id_speac= ',id_spec,id_species_total, & 116 'emission= ',emission, & 117 'pdtphys= ',pdtphys, & 118 'rapport m_emiss/m_air*f_lay= ', m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k), & 119 'stretchlong= ', stretchlong, & 120 'theta= ', theta 121 121 122 122 IF(emission < 1.E-34) emission = 0.0 … … 124 124 IF (flh2o==0) THEN 125 125 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: tr_ser avant/apres',& 126 &'i= ',i,'k= ',k, 'flh2o= ',flh2o, &127 &tr_seri(i,k,id_spec), &128 &tr_seri(i,k,id_spec)+emission*pdtphys126 'i= ',i,'k= ',k, 'flh2o= ',flh2o, & 127 tr_seri(i,k,id_spec), & 128 tr_seri(i,k,id_spec)+emission*pdtphys 129 129 130 130 tr_seri(i,k,id_spec)=tr_seri(i,k,id_spec)+emission*pdtphys … … 141 141 142 142 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: ',& 143 &'i= ',i,'k= ',k, 'flh2o= ',flh2o, &144 &'emission= ',emission, &145 &'d_q_emiss(i,k)= ',d_q_emiss(i,k)143 'i= ',i,'k= ',k, 'flh2o= ',flh2o, & 144 'emission= ',emission, & 145 'd_q_emiss(i,k)= ',d_q_emiss(i,k) 146 146 147 147 IF(d_q_emiss(i,k) > 1.E34) THEN -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90
r5086 r5087 84 84 ! factor for converting dry to wet radius 85 85 f_r_wet(ilon,ilev) = (dens_aer_dry/(DENSO4(ilon,ilev)*1.e3)/ & 86 &(R2SO4(ilon,ilev)*1.e-2))**third86 (R2SO4(ilon,ilev)*1.e-2))**third 87 87 ! *** End of H2SO4-H2O flat surface *** 88 88 … … 121 121 ! wet radius (m) 122 122 radwet = 1.e-2*RRSI(IK)*(dens_aer_dry/(denso4ik*1.e3)/ & 123 &(r2so4ik*1.e-2))**third123 (r2so4ik*1.e-2))**third 124 124 fkelvin=exp( 2.*1.e-3*surtens*1.e-6*mvh2o / (radwet*rgas*temp) ) 125 125 pph2okel=pph2ogas(ilon,ilev) / fkelvin … … 133 133 ! factor for converting dry to wet radius 134 134 f_r_wetB(ilon,ilev,IK) = (dens_aer_dry/(DENSO4B(ilon,ilev,IK)*1.e3)/ & 135 &(R2SO4B(ilon,ilev,IK)*1.e-2))**third135 (R2SO4B(ilon,ilev,IK)*1.e-2))**third 136 136 ! 137 137 ! print*,'R,Rwet(m),kelvin,h2so4(%),ro=',RRSI(ik),radwet,fkelvin, & … … 565 565 ! H2SO4 activity 566 566 DATA X/ & 567 & 0.0,0.25,0.78,1.437,2.19,3.07,4.03,5.04,6.08 &568 & ,7.13,8.18,14.33,18.59,28.59,39.17,49.49 &569 & ,102.4,157.8,215.7,276.9,341.6,409.8,481.5,556.6 &570 & ,635.5,719.,808.,902.,1000.,1103.,1211.,1322.,1437.,1555. &571 & ,1677.,1800.,1926.,2054.,2183.,2312.,2442.,2572.,2701.,2829. &572 & ,2955.,3080.,3203.,3325.,3446.,3564.,3681.,3796.,3910.,4022. &573 & ,4134.,4351.,4564.,4771.,4974.,5171.,5364.,5551.,5732.,5908. &574 & ,6079.,6244.,6404.,6559.,6709.,6854.,6994.,7131.,7264.,7393. &575 & ,7520.,7821.,8105.,8373.,8627.,8867.,9093.,9308.,9511.,9703. &576 & ,9885.,10060.,10225.,10535.,10819.,11079.,11318.,11537. &577 & ,11740.,12097.,12407.,12676.,12915.,13126.,13564.,13910. &578 & ,14191.,14423.,14617.,14786.,10568.,15299.,15491.,15654. &579 &,15811./567 0.0,0.25,0.78,1.437,2.19,3.07,4.03,5.04,6.08 & 568 ,7.13,8.18,14.33,18.59,28.59,39.17,49.49 & 569 ,102.4,157.8,215.7,276.9,341.6,409.8,481.5,556.6 & 570 ,635.5,719.,808.,902.,1000.,1103.,1211.,1322.,1437.,1555. & 571 ,1677.,1800.,1926.,2054.,2183.,2312.,2442.,2572.,2701.,2829. & 572 ,2955.,3080.,3203.,3325.,3446.,3564.,3681.,3796.,3910.,4022. & 573 ,4134.,4351.,4564.,4771.,4974.,5171.,5364.,5551.,5732.,5908. & 574 ,6079.,6244.,6404.,6559.,6709.,6854.,6994.,7131.,7264.,7393. & 575 ,7520.,7821.,8105.,8373.,8627.,8867.,9093.,9308.,9511.,9703. & 576 ,9885.,10060.,10225.,10535.,10819.,11079.,11318.,11537. & 577 ,11740.,12097.,12407.,12676.,12915.,13126.,13564.,13910. & 578 ,14191.,14423.,14617.,14786.,10568.,15299.,15491.,15654. & 579 ,15811./ 580 580 ! H2SO4 weight fraction (percent) 581 581 DATA XC/ & 582 & 100.0,99.982,99.963,99.945,99.927,99.908,99.890,99.872 &583 & ,99.853,99.835,99.817,99.725,99.634,99.452,99.270 &584 & ,99.090,98.196,97.319,96.457,95.610,94.777,93.959,93.156 &585 & ,92.365,91.588,90.824,90.073,89.334,88.607,87.892,87.188 &586 & ,86.495,85.814,85.143,84.482,83.832,83.191,82.560,81.939 &587 & ,81.327,80.724,80.130,79.545,78.968,78.399,77.839,77.286 &588 & ,76.741,76.204,75.675,75.152,74.637,74.129,73.628,73.133 &589 & ,72.164,71.220,70.300,69.404,68.530,67.678,66.847,66.037 &590 & ,65.245,64.472,63.718,62.981,62.261,61.557,60.868,60.195 &591 & ,59.537,58.893,58.263,57.646,56.159,54.747,53.405,52.126 &592 & ,50.908,49.745,48.634,47.572,46.555,45.580,44.646,43.749 &593 & ,42.059,40.495,39.043,37.691,36.430,35.251,33.107,31.209 &594 & ,29.517,27.999,26.629,23.728,21.397,19.482,17.882,16.525 &595 &,15.360,13.461,11.980,10.792,9.819,8.932/582 100.0,99.982,99.963,99.945,99.927,99.908,99.890,99.872 & 583 ,99.853,99.835,99.817,99.725,99.634,99.452,99.270 & 584 ,99.090,98.196,97.319,96.457,95.610,94.777,93.959,93.156 & 585 ,92.365,91.588,90.824,90.073,89.334,88.607,87.892,87.188 & 586 ,86.495,85.814,85.143,84.482,83.832,83.191,82.560,81.939 & 587 ,81.327,80.724,80.130,79.545,78.968,78.399,77.839,77.286 & 588 ,76.741,76.204,75.675,75.152,74.637,74.129,73.628,73.133 & 589 ,72.164,71.220,70.300,69.404,68.530,67.678,66.847,66.037 & 590 ,65.245,64.472,63.718,62.981,62.261,61.557,60.868,60.195 & 591 ,59.537,58.893,58.263,57.646,56.159,54.747,53.405,52.126 & 592 ,50.908,49.745,48.634,47.572,46.555,45.580,44.646,43.749 & 593 ,42.059,40.495,39.043,37.691,36.430,35.251,33.107,31.209 & 594 ,29.517,27.999,26.629,23.728,21.397,19.482,17.882,16.525 & 595 ,15.360,13.461,11.980,10.792,9.819,8.932/ 596 596 597 597 DO I=1,klon … … 776 776 real, intent(in) :: T 777 777 real, parameter :: & 778 &b1=1.01325e5, &779 &b2=11.5, &780 &b3=1.0156e4, &781 &b4=0.38/545., &782 &tref=360.15778 b1=1.01325e5, & 779 b2=11.5, & 780 b3=1.0156e4, & 781 b4=0.38/545., & 782 tref=360.15 783 783 784 784 ! saturation vapor pressure ( N/m2 = Pa = kg/(m.s2) ) 785 785 psh2so4_out=b1*exp( -b2 +b3*( 1./tref-1./T & 786 & +b4*(1.+log(tref/T)-tref/T) ) )786 +b4*(1.+log(tref/T)-tref/T) ) ) 787 787 788 788 return … … 818 818 ! saturation vapor pressure (N/m2 = 1 Pa = 1 kg/(m·s2)) 819 819 psh2o_out=exp( 77.34491296 -7235.424651/T & 820 &-8.2*log(T) + 5.7133e-3*T )820 -8.2*log(T) + 5.7133e-3*T ) 821 821 else 822 822 ! Tabazadeh et al., 1997, parameterization for 185<T<260 … … 825 825 ; 826 826 psh2o_out=18.452406985 -3505.1578807/T & 827 &-330918.55082/(T*T) &828 &+12725068.262/(T*T*T)827 -330918.55082/(T*T) & 828 +12725068.262/(T*T*T) 829 829 ! in Pa 830 830 psh2o_out=100.*exp(psh2o_out) … … 844 844 real, intent(in) :: T, so4mfrac 845 845 real, parameter :: & 846 &a1= 0.7681724,&847 &a2= 2.184714, &848 &a3= 7.163002, &849 &a4=-44.31447, &850 &a5= 88.74606, &851 &a6=-75.73729, &852 &a7= 23.43228846 a1= 0.7681724,& 847 a2= 2.184714, & 848 a3= 7.163002, & 849 a4=-44.31447, & 850 a5= 88.74606, & 851 a6=-75.73729, & 852 a7= 23.43228 853 853 real, parameter :: & 854 &b1= 1.808225e-3, &855 &b2=-9.294656e-3, &856 &b3=-3.742148e-2, &857 &b4= 2.565321e-1, &858 &b5=-5.362872e-1, &859 &b6= 4.857736e-1, &860 &b7=-1.629592e-1854 b1= 1.808225e-3, & 855 b2=-9.294656e-3, & 856 b3=-3.742148e-2, & 857 b4= 2.565321e-1, & 858 b5=-5.362872e-1, & 859 b6= 4.857736e-1, & 860 b7=-1.629592e-1 861 861 real, parameter :: & 862 &c1=-3.478524e-6, &863 &c2= 1.335867e-5, &864 &c3= 5.195706e-5, &865 &c4=-3.717636e-4, &866 &c5= 7.990811e-4, &867 &c6=-7.458060e-4, &868 &c7= 2.581390e-4862 c1=-3.478524e-6, & 863 c2= 1.335867e-5, & 864 c3= 5.195706e-5, & 865 c4=-3.717636e-4, & 866 c5= 7.990811e-4, & 867 c6=-7.458060e-4, & 868 c7= 2.581390e-4 869 869 real :: a,b,c,so4m2,so4m3,so4m4,so4m5,so4m6 870 870 … … 876 876 877 877 a=+a1+a2*so4mfrac+a3*so4m2+a4*so4m3 & 878 &+a5*so4m4+a6*so4m5+a7*so4m6878 +a5*so4m4+a6*so4m5+a7*so4m6 879 879 b=+b1+b2*so4mfrac+b3*so4m2+b4*so4m3 & 880 &+b5*so4m4+b6*so4m5+b7*so4m6880 +b5*so4m4+b6*so4m5+b7*so4m6 881 881 c=+c1+c2*so4mfrac+c3*so4m2+c4*so4m3 & 882 &+c5*so4m4+c6*so4m5+c7*so4m6882 +c5*so4m4+c6*so4m5+c7*so4m6 883 883 density_out=(a+b*T+c*T*T) ! units are gm/cm**3 884 884 … … 895 895 real :: a,b,so4mfrac,so4m2,so4m3,so4m4,so4m5,so4sig 896 896 real, parameter :: & 897 &a1= 0.11864, &898 &a2=-0.11651, &899 &a3= 0.76852, &900 &a4=-2.40909, &901 &a5= 2.95434, &902 &a6=-1.25852897 a1= 0.11864, & 898 a2=-0.11651, & 899 a3= 0.76852, & 900 a4=-2.40909, & 901 a5= 2.95434, & 902 a6=-1.25852 903 903 real, parameter :: & 904 &b1=-1.5709e-4, &905 &b2= 4.0102e-4, &906 &b3=-2.3995e-3, &907 &b4= 7.611235e-3, &908 &b5=-9.37386e-3, &909 &b6= 3.89722e-3904 b1=-1.5709e-4, & 905 b2= 4.0102e-4, & 906 b3=-2.3995e-3, & 907 b4= 7.611235e-3, & 908 b5=-9.37386e-3, & 909 b6= 3.89722e-3 910 910 real, parameter :: convfac=1.e3 ! convert from newton/m to dyne/cm 911 911 real, parameter :: Mw=18.01528, Ma=98.079 … … 957 957 if(aw <= 0.05 .and. aw > 0.) then 958 958 y1=12.372089320*aw**(-0.16125516114) & 959 &-30.490657554*aw -2.1133114241959 -30.490657554*aw -2.1133114241 960 960 y2=13.455394705*aw**(-0.19213122550) & 961 &-34.285174607*aw -1.7620073078961 -34.285174607*aw -1.7620073078 962 962 else if(aw <= 0.85 .and. aw > 0.05) then 963 963 y1=11.820654354*aw**(-0.20786404244) & 964 &-4.8073063730*aw -5.1727540348964 -4.8073063730*aw -5.1727540348 965 965 y2=12.891938068*aw**(-0.23233847708) & 966 &-6.4261237757*aw -4.9005471319966 -6.4261237757*aw -4.9005471319 967 967 else 968 968 y1=-180.06541028*aw**(-0.38601102592) & 969 &-93.317846778*aw +273.88132245969 -93.317846778*aw +273.88132245 970 970 y2=-176.95814097*aw**(-0.36257048154) & 971 &-90.469744201*aw +267.45509988971 -90.469744201*aw +267.45509988 972 972 end if 973 973 ! h2so4 molality (m=moles of h2so4 (solute)/ kg of h2o(solvent)) … … 1029 1029 real, intent(in) :: T, ws 1030 1030 real, dimension(22),parameter :: x=(/ & 1031 &2.393284E-02,-4.359335E-05,7.961181E-08,0.0,-0.198716351, &1032 &1.39564574E-03,-2.020633E-06,0.51684706,-3.0539E-03,4.505475E-06, &1033 &-0.30119511,1.840408E-03,-2.7221253742E-06,-0.11331674116, &1034 &8.47763E-04,-1.22336185E-06,0.3455282,-2.2111E-03,3.503768245E-06, &1035 &-0.2315332,1.60074E-03,-2.5827835E-06/)1031 2.393284E-02,-4.359335E-05,7.961181E-08,0.0,-0.198716351, & 1032 1.39564574E-03,-2.020633E-06,0.51684706,-3.0539E-03,4.505475E-06, & 1033 -0.30119511,1.840408E-03,-2.7221253742E-06,-0.11331674116, & 1034 8.47763E-04,-1.22336185E-06,0.3455282,-2.2111E-03,3.503768245E-06, & 1035 -0.2315332,1.60074E-03,-2.5827835E-06/) 1036 1036 1037 1037 real :: w -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/traccoag_mod.F90
r5082 r5087 11 11 12 12 USE phys_local_var_mod, ONLY: mdw, R2SO4, DENSO4, f_r_wet, surf_PM25_sulf, & 13 &budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, &14 &R2SO4B, DENSO4B, f_r_wetB, sulfmmr, SAD_sulfate, sulfmmr_mode, nd_mode13 budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, & 14 R2SO4B, DENSO4B, f_r_wetB, sulfmmr, SAD_sulfate, sulfmmr_mode, nd_mode 15 15 16 16 USE dimphy … … 340 340 !and are dry at T = 20 deg. C and 50 perc. humidity 341 341 surf_PM25_sulf(i)=surf_PM25_sulf(i)+tr_seri(i,1,it+nbtr_sulgas) & 342 &*132./98.*dens_aer_dry*4./3.*RPI*(mdw(it)/2.)**3 &343 &*pplay(i,1)/t_seri(i,1)/RD*1.e9342 *132./98.*dens_aer_dry*4./3.*RPI*(mdw(it)/2.)**3 & 343 *pplay(i,1)/t_seri(i,1)/RD*1.e9 344 344 ENDIF 345 345 ENDDO … … 366 366 ! equivalent to condensed H2SO4 mmr= H2SO4 kg / kgA in bin it 367 367 sulfmmr_mode(i,j,it) = tr_seri(i,j,it+nbtr_sulgas) & ! [DRY part/kgA in bin it] 368 & *(4./3.)*RPI*(mdw(it)/2.)**3. & ! [mdw: dry diameter in m]369 & *dens_aer_dry ! [dry aerosol mass density in kg/m3]368 *(4./3.)*RPI*(mdw(it)/2.)**3. & ! [mdw: dry diameter in m] 369 *dens_aer_dry ! [dry aerosol mass density in kg/m3] 370 370 371 371 ! sulfmmr: Sulfate aerosol concentration (dry mass mixing ratio) … … 375 375 ! nd_mode: particle concentration in different modes (DRY part/m3) 376 376 nd_mode(i,j,it) = tr_seri(i,j,it+nbtr_sulgas) & ! [DRY part/kgA in bin it] 377 &*pplay(i,j)/t_seri(i,j)/RD ! [air mass concentration in kg air /m3A]377 *pplay(i,j)/t_seri(i,j)/RD ! [air mass concentration in kg air /m3A] 378 378 379 379 IF(flag_new_strat_compo) THEN 380 380 ! SAD_sulfate: SAD WET sulfate aerosols (cm2/cm3) 381 381 SAD_sulfate(i,j) = SAD_sulfate(i,j) + nd_mode(i,j,it) & ! [DRY part/m3A (in bin it)] 382 & *4.*RPI*( mdw(it)*f_r_wetB(i,j,it)/2. )**2. & ! [WET SA of part it in m2]383 &*1.e-2 ! conversion from m2/m3 to cm2/cm3A382 *4.*RPI*( mdw(it)*f_r_wetB(i,j,it)/2. )**2. & ! [WET SA of part it in m2] 383 *1.e-2 ! conversion from m2/m3 to cm2/cm3A 384 384 ELSE 385 385 ! SAD_sulfate: SAD WET sulfate aerosols (cm2/cm3) 386 386 SAD_sulfate(i,j) = SAD_sulfate(i,j) + nd_mode(i,j,it) & ! [DRY part/m3A (in bin it)] 387 & *4.*RPI*( mdw(it)*f_r_wet(i,j)/2. )**2. & ! [WET SA of part it in m2]388 &*1.e-2 ! conversion from m2/m3 to cm2/cm3A387 *4.*RPI*( mdw(it)*f_r_wet(i,j)/2. )**2. & ! [WET SA of part it in m2] 388 *1.e-2 ! conversion from m2/m3 to cm2/cm3A 389 389 ENDIF 390 390 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.
