Changeset 3798 for LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer
- Timestamp:
- Jan 11, 2021, 11:24:08 PM (4 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 13 edited
-
. (modified) (1 prop)
-
libf/phylmd/StratAer/aer_sedimnt.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/calcaerosolstrato_rrtm.F90 (modified) (2 diffs)
-
libf/phylmd/StratAer/coagulate.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/cond_evap_tstep_mod.F90 (modified) (2 diffs)
-
libf/phylmd/StratAer/interp_sulf_input.F90 (modified) (2 diffs)
-
libf/phylmd/StratAer/micphy_tstep.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/miecalc_aer.F90 (modified) (6 diffs)
-
libf/phylmd/StratAer/nucleation_tstep_mod.F90 (modified) (7 diffs)
-
libf/phylmd/StratAer/ocs_to_so2.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/so2_to_h2so4.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/sulfate_aer_mod.F90 (modified) (1 diff)
-
libf/phylmd/StratAer/traccoag_mod.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
-
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/aer_sedimnt.F90
r3605 r3798 19 19 USE phys_local_var_mod, ONLY: mdw, budg_sed_part, DENSO4, f_r_wet, vsed_aer 20 20 USE dimphy, ONLY : klon,klev 21 USE infotrac 21 USE infotrac_phy 22 22 USE aerophys 23 23 USE YOMCST -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/calcaerosolstrato_rrtm.F90
r3605 r3798 4 4 SUBROUTINE calcaerosolstrato_rrtm(pplay,t_seri,paprs,debut) 5 5 6 USE infotrac, ONLY : nbtr7 6 USE phys_state_var_mod, ONLY: tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, tau_aero_lw_rrtm 8 7 USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_550, tau_strat_1020, stratomask … … 118 117 zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG !thickness of layer in m 119 118 tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz 120 tau_strat_1020(i,k)=tau_strat_wave(i,k, 6)/zdz119 tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz 121 120 ENDDO 122 121 ENDDO -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/coagulate.F90
r3605 r3798 25 25 USE dimphy, ONLY : klon,klev 26 26 USE aerophys 27 USE infotrac 27 USE infotrac_phy 28 28 USE phys_local_var_mod, ONLY: DENSO4, f_r_wet 29 29 USE YOMCST -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
r3605 r3798 20 20 21 21 USE aerophys 22 USE infotrac 22 USE infotrac_phy 23 23 USE YOMCST, ONLY : RPI 24 24 … … 141 141 142 142 USE aerophys 143 USE infotrac 143 USE infotrac_phy 144 144 USE YOMCST, ONLY : RPI 145 145 -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/interp_sulf_input.F90
r3605 r3798 16 16 USE dimphy 17 17 USE phys_cal_mod 18 USE infotrac 18 USE infotrac_phy 19 19 USE aerophys 20 20 USE YOMCST … … 34 34 INTEGER n_lon ! number of longitudes in the input data 35 35 INTEGER, SAVE :: n_lev ! number of levels in the input data 36 !$OMP THREADPRIVATE(n_lev) 36 37 INTEGER n_mth ! number of months in the input data 37 38 INTEGER, SAVE :: mth_pre -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/micphy_tstep.F90
r3605 r3798 7 7 USE dimphy, ONLY : klon,klev 8 8 USE aerophys 9 USE infotrac 9 USE infotrac_phy 10 10 USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, budg_h2so4_to_part, R2SO4, DENSO4, f_r_wet 11 11 USE nucleation_tstep_mod -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/miecalc_aer.F90
r3605 r3798 1 !2 ! $Id$3 !4 1 SUBROUTINE MIECALC_AER(tau_strat, piz_strat, cg_strat, tau_strat_wave, tau_lw_abs_rrtm, paprs, debut) 5 2 … … 21 18 USE aerophys 22 19 USE aero_mod 23 USE infotrac , ONLY : nbtr, nbtr_bin, nbtr_sulgas, id_SO2_strat20 USE infotrac_phy, ONLY : nbtr, nbtr_bin, nbtr_sulgas, id_SO2_strat 24 21 USE dimphy 25 22 USE YOMCST , ONLY : RG, RPI … … 157 154 REAL ss_g(nbands_sw_rrtm+nbands_lw_rrtm+nwave_sw+nwave_lw) 158 155 156 !-- fichier h2so4_0.75_300.00_hummel_1988_p_q.dat 157 ! -- wavenumber (cm-1), wavelength (um), n_r, n_i 159 158 INTEGER, PARAMETER :: nb_lambda_h2so4=62 160 159 REAL, DIMENSION (nb_lambda_h2so4,4) :: ref_ind 161 !-- fichier h2so4_0.75_300.00_hummel_1988_p_q.dat 162 ! -- wavenumber (cm-1), wavelength (um), n_r, n_i 163 DATA ref_ind / & 164 200.000, 50.0000, 2.01000, 6.5000E-01, & 165 250.000, 40.0000, 1.94000, 6.3000E-01, & 166 285.714, 35.0000, 1.72000, 5.2000E-01, & 167 333.333, 30.0000, 1.73000, 2.9000E-01, & 168 358.423, 27.9000, 1.78000, 2.5000E-01, & 169 400.000, 25.0000, 1.84000, 2.4000E-01, & 170 444.444, 22.5000, 1.82000, 2.9000E-01, & 171 469.484, 21.3000, 1.79000, 2.5000E-01, & 172 500.000, 20.0000, 1.81000, 2.3000E-01, & 173 540.541, 18.5000, 1.92700, 3.0200E-01, & 174 555.556, 18.0000, 1.95000, 4.1000E-01, & 175 581.395, 17.2000, 1.72400, 5.9000E-01, & 176 609.756, 16.4000, 1.52000, 4.1400E-01, & 177 666.667, 15.0000, 1.59000, 2.1100E-01, & 178 675.676, 14.8000, 1.61000, 2.0500E-01, & 179 714.286, 14.0000, 1.64000, 1.9500E-01, & 180 769.231, 13.0000, 1.69000, 1.9500E-01, & 181 800.000, 12.5000, 1.74000, 1.9800E-01, & 182 869.565, 11.5000, 1.89000, 3.7400E-01, & 183 909.091, 11.0000, 1.67000, 4.8500E-01, & 184 944.198, 10.5910, 1.72000, 3.4000E-01, & 185 1000.000, 10.0000, 1.89000, 4.5500E-01, & 186 1020.408, 9.8000, 1.91000, 6.8000E-01, & 187 1052.632, 9.5000, 1.67000, 7.5000E-01, & 188 1086.957, 9.2000, 1.60000, 5.8600E-01, & 189 1111.111, 9.0000, 1.65000, 6.3300E-01, & 190 1149.425, 8.7000, 1.53000, 7.7200E-01, & 191 1176.471, 8.5000, 1.37000, 7.5500E-01, & 192 1219.512, 8.2000, 1.20000, 6.4500E-01, & 193 1265.823, 7.9000, 1.14000, 4.8800E-01, & 194 1388.889, 7.2000, 1.21000, 1.7600E-01, & 195 1538.462, 6.5000, 1.37000, 1.2800E-01, & 196 1612.903, 6.2000, 1.42400, 1.6500E-01, & 197 1666.667, 6.0000, 1.42500, 1.9500E-01, & 198 1818.182, 5.5000, 1.33700, 1.8300E-01, & 199 2000.000, 5.0000, 1.36000, 1.2100E-01, & 200 2222.222, 4.5000, 1.38500, 1.2000E-01, & 201 2500.000, 4.0000, 1.39800, 1.2600E-01, & 202 2666.667, 3.7500, 1.39600, 1.3100E-01, & 203 2857.143, 3.5000, 1.37600, 1.5800E-01, & 204 2948.113, 3.3920, 1.35200, 1.5900E-01, & 205 3125.000, 3.2000, 1.31100, 1.3500E-01, & 206 3333.333, 3.0000, 1.29300, 9.5500E-02, & 207 3703.704, 2.7000, 1.30300, 5.7000E-03, & 208 4000.000, 2.5000, 1.34400, 3.7600E-03, & 209 4444.444, 2.2500, 1.37000, 1.8000E-03, & 210 5000.000, 2.0000, 1.38400, 1.2600E-03, & 211 5555.556, 1.8000, 1.39000, 5.5000E-04, & 212 6510.417, 1.5360, 1.40300, 1.3700E-04, & 213 7692.308, 1.3000, 1.41000, 1.0000E-05, & 214 9433.962, 1.0600, 1.42000, 1.5000E-06, & 215 11627.907, 0.8600, 1.42500, 1.7900E-07, & 216 14409.222, 0.6940, 1.42800, 1.9900E-08, & 217 15797.788, 0.6330, 1.42900, 1.4700E-08, & 218 18181.818, 0.5500, 1.43000, 1.0000E-08, & 219 19417.476, 0.5150, 1.43100, 1.0000E-08, & 220 20491.803, 0.4880, 1.43200, 1.0000E-08, & 221 25000.000, 0.4000, 1.44000, 1.0000E-08, & 222 29673.591, 0.3370, 1.45900, 1.0000E-08, & 223 33333.333, 0.3000, 1.46900, 1.0000E-08, & 224 40000.000, 0.2500, 1.48400, 1.0000E-08, & 225 50000.000, 0.2000, 1.49800, 1.0000E-08 / 160 226 161 !--------------------------------------------------------- 227 162 228 163 IF (debut) THEN 229 164 165 ref_ind = RESHAPE( (/ & 166 200.000, 50.0000, 2.01000, 6.5000E-01, & 167 250.000, 40.0000, 1.94000, 6.3000E-01, & 168 285.714, 35.0000, 1.72000, 5.2000E-01, & 169 333.333, 30.0000, 1.73000, 2.9000E-01, & 170 358.423, 27.9000, 1.78000, 2.5000E-01, & 171 400.000, 25.0000, 1.84000, 2.4000E-01, & 172 444.444, 22.5000, 1.82000, 2.9000E-01, & 173 469.484, 21.3000, 1.79000, 2.5000E-01, & 174 500.000, 20.0000, 1.81000, 2.3000E-01, & 175 540.541, 18.5000, 1.92700, 3.0200E-01, & 176 555.556, 18.0000, 1.95000, 4.1000E-01, & 177 581.395, 17.2000, 1.72400, 5.9000E-01, & 178 609.756, 16.4000, 1.52000, 4.1400E-01, & 179 666.667, 15.0000, 1.59000, 2.1100E-01, & 180 675.676, 14.8000, 1.61000, 2.0500E-01, & 181 714.286, 14.0000, 1.64000, 1.9500E-01, & 182 769.231, 13.0000, 1.69000, 1.9500E-01, & 183 800.000, 12.5000, 1.74000, 1.9800E-01, & 184 869.565, 11.5000, 1.89000, 3.7400E-01, & 185 909.091, 11.0000, 1.67000, 4.8500E-01, & 186 944.198, 10.5910, 1.72000, 3.4000E-01, & 187 1000.000, 10.0000, 1.89000, 4.5500E-01, & 188 1020.408, 9.8000, 1.91000, 6.8000E-01, & 189 1052.632, 9.5000, 1.67000, 7.5000E-01, & 190 1086.957, 9.2000, 1.60000, 5.8600E-01, & 191 1111.111, 9.0000, 1.65000, 6.3300E-01, & 192 1149.425, 8.7000, 1.53000, 7.7200E-01, & 193 1176.471, 8.5000, 1.37000, 7.5500E-01, & 194 1219.512, 8.2000, 1.20000, 6.4500E-01, & 195 1265.823, 7.9000, 1.14000, 4.8800E-01, & 196 1388.889, 7.2000, 1.21000, 1.7600E-01, & 197 1538.462, 6.5000, 1.37000, 1.2800E-01, & 198 1612.903, 6.2000, 1.42400, 1.6500E-01, & 199 1666.667, 6.0000, 1.42500, 1.9500E-01, & 200 1818.182, 5.5000, 1.33700, 1.8300E-01, & 201 2000.000, 5.0000, 1.36000, 1.2100E-01, & 202 2222.222, 4.5000, 1.38500, 1.2000E-01, & 203 2500.000, 4.0000, 1.39800, 1.2600E-01, & 204 2666.667, 3.7500, 1.39600, 1.3100E-01, & 205 2857.143, 3.5000, 1.37600, 1.5800E-01, & 206 2948.113, 3.3920, 1.35200, 1.5900E-01, & 207 3125.000, 3.2000, 1.31100, 1.3500E-01, & 208 3333.333, 3.0000, 1.29300, 9.5500E-02, & 209 3703.704, 2.7000, 1.30300, 5.7000E-03, & 210 4000.000, 2.5000, 1.34400, 3.7600E-03, & 211 4444.444, 2.2500, 1.37000, 1.8000E-03, & 212 5000.000, 2.0000, 1.38400, 1.2600E-03, & 213 5555.556, 1.8000, 1.39000, 5.5000E-04, & 214 6510.417, 1.5360, 1.40300, 1.3700E-04, & 215 7692.308, 1.3000, 1.41000, 1.0000E-05, & 216 9433.962, 1.0600, 1.42000, 1.5000E-06, & 217 11627.907, 0.8600, 1.42500, 1.7900E-07, & 218 14409.222, 0.6940, 1.42800, 1.9900E-08, & 219 15797.788, 0.6330, 1.42900, 1.4700E-08, & 220 18181.818, 0.5500, 1.43000, 1.0000E-08, & 221 19417.476, 0.5150, 1.43100, 1.0000E-08, & 222 20491.803, 0.4880, 1.43200, 1.0000E-08, & 223 25000.000, 0.4000, 1.44000, 1.0000E-08, & 224 29673.591, 0.3370, 1.45900, 1.0000E-08, & 225 33333.333, 0.3000, 1.46900, 1.0000E-08, & 226 40000.000, 0.2500, 1.48400, 1.0000E-08, & 227 50000.000, 0.2000, 1.49800, 1.0000E-08 /), (/nb_lambda_h2so4,4/), order=(/2,1/) ) 228 230 229 !--initialising dry diameters to geometrically spaced mass/volume (see Jacobson 1994) 231 232 233 234 235 236 237 238 239 240 241 242 PRINT *,'init mdw=', mdw230 mdw(1)=mdwmin 231 IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio 232 mdw(2)=mdw(1)*2.**(1./3.) 233 DO it=3, nbtr_bin 234 mdw(it)=mdw(it-1)*V_rat**(1./3.) 235 ENDDO 236 ELSE 237 DO it=2, nbtr_bin 238 mdw(it)=mdw(it-1)*V_rat**(1./3.) 239 ENDDO 240 ENDIF 241 WRITE(lunout,*) 'init mdw=', mdw 243 242 244 243 !--compute particle radius for a composition of 75% H2SO4 / 25% H2O at T=293K … … 286 285 287 286 IF (refr_ind_interpol) THEN 288 287 289 288 ilambda_max=ref_ind(1,2)/1.e6 !--in m 290 289 ilambda_min=ref_ind(nb_lambda_h2so4,2)/1.e6 !--in m … … 381 380 Nmax=INT(x+4*x**(1./3.)+2.)+1 382 381 ELSE 383 PRINT *,'x out of bound, x=', x382 WRITE(lunout,*) 'x out of bound, x=', x 384 383 STOP 385 384 ENDIF … … 461 460 omegatot=omegatot+r**2*Q_ext*omega*number 462 461 gtot =gtot+r**2*Q_sca*g*number 463 462 464 463 ENDDO !---bin 464 465 465 !------------------------------------------------------------------ 466 466 -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/nucleation_tstep_mod.F90
r3605 r3798 9 9 10 10 USE aerophys 11 USE infotrac 11 USE infotrac_phy 12 12 USE YOMCST, ONLY : RPI, RD, RMD, RKBOL, RNAVO 13 13 … … 73 73 74 74 USE aerophys 75 USE infotrac 75 USE infotrac_phy 76 76 77 77 IMPLICIT NONE … … 206 206 & + 0.0000504021689382576*t*LOG(rhoa) 207 207 208 jnuc= 0.1430901615568665 + 2.219563673425199*t - 0.02739106114964264*t**2 + 209 & 0.00007228107239317088*t**3 + 5.91822263375044/x + 208 jnuc= 0.1430901615568665 + 2.219563673425199*t - 0.02739106114964264*t**2 + & 209 & 0.00007228107239317088*t**3 + 5.91822263375044/x + & 210 210 & 0.1174886643003278*LOG(rh) + 0.4625315047693772*t*LOG(rh) - & 211 & 0.01180591129059253*t**2*LOG(rh) + 211 & 0.01180591129059253*t**2*LOG(rh) + & 212 212 & 0.0000404196487152575*t**3*LOG(rh) + (15.79628615047088*LOG(rh))/x - & 213 & 0.215553951893509*LOG(rh)**2 - 0.0810269192332194*t*LOG(rh)**2 + 213 & 0.215553951893509*LOG(rh)**2 - 0.0810269192332194*t*LOG(rh)**2 + & 214 214 & 0.001435808434184642*t**2*LOG(rh)**2 - & 215 215 & 4.775796947178588E-6*t**3*LOG(rh)**2 - & 216 & (2.912974063702185*LOG(rh)**2)/x - 3.588557942822751*LOG(rh)**3 + 216 & (2.912974063702185*LOG(rh)**2)/x - 3.588557942822751*LOG(rh)**3 + & 217 217 & 0.04950795302831703*t*LOG(rh)**3 - & 218 & 0.0002138195118737068*t**2*LOG(rh)**3 + 218 & 0.0002138195118737068*t**2*LOG(rh)**3 + & 219 219 & 3.108005107949533E-7*t**3*LOG(rh)**3 - & 220 & (0.02933332747098296*LOG(rh)**3)/x + 220 & (0.02933332747098296*LOG(rh)**3)/x + & 221 221 & 1.145983818561277*LOG(rhoa) - & 222 & 0.6007956227856778*t*LOG(rhoa) + 222 & 0.6007956227856778*t*LOG(rhoa) + & 223 223 & 0.00864244733283759*t**2*LOG(rhoa) - & 224 224 & 0.00002289467254710888*t**3*LOG(rhoa) - & 225 & (8.44984513869014*LOG(rhoa))/x + 226 & 2.158548369286559*LOG(rh)*LOG(rhoa) + 225 & (8.44984513869014*LOG(rhoa))/x + & 226 & 2.158548369286559*LOG(rh)*LOG(rhoa) + & 227 227 & 0.0808121412840917*t*LOG(rh)*LOG(rhoa) - & 228 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 + 229 & 4.019572560156515E-7*t**3*LOG(rh)*LOG(rhoa) + & 230 & (0.7213255852557236*LOG(rh)*LOG(rhoa))/x + & 231 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) + 232 & 0.01601062035325362*t*LOG(rh)**2*LOG(rhoa) + & 233 & 0.00003771238979714162*t**2*LOG(rh)**2*LOG(rhoa) + & 234 234 & 3.217942606371182E-8*t**3*LOG(rh)**2*LOG(rhoa) - & 235 & (0.01132550810022116*LOG(rh)**2*LOG(rhoa))/x + 235 & (0.01132550810022116*LOG(rh)**2*LOG(rhoa))/x + & 236 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 + 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 240 & (0.7118597859976135*LOG(rhoa)**2)/x - & 241 & 1.056105824379897*LOG(rh)*LOG(rhoa)**2 + 241 & 1.056105824379897*LOG(rh)*LOG(rhoa)**2 + & 242 242 & 0.00903377584628419*t*LOG(rh)*LOG(rhoa)**2 - & 243 & 0.00001984167387090606*t**2*LOG(rh)*LOG(rhoa)**2 + 243 & 0.00001984167387090606*t**2*LOG(rh)*LOG(rhoa)**2 + & 244 244 & 2.460478196482179E-8*t**3*LOG(rh)*LOG(rhoa)**2 - & 245 245 & (0.05790872906645181*LOG(rh)*LOG(rhoa)**2)/x - & 246 & 0.1487119673397459*LOG(rhoa)**3 + 246 & 0.1487119673397459*LOG(rhoa)**3 + & 247 247 & 0.002835082097822667*t*LOG(rhoa)**3 - & 248 & 9.24618825471694E-6*t**2*LOG(rhoa)**3 + 248 & 9.24618825471694E-6*t**2*LOG(rhoa)**3 + & 249 249 & 5.004267665960894E-9*t**3*LOG(rhoa)**3 - & 250 250 & (0.01270805101481648*LOG(rhoa)**3)/x … … 252 252 253 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) + 254 & 0.1017165718716887/x - 0.002050640345231486*LOG(rh) - 0.007585041382707174*t*LOG(rh) + & 255 255 & 0.0001926539658089536*t**2*LOG(rh) - 6.70429719683894E-7*t**3*LOG(rh) - & 256 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 + 257 & 0.00001547571354871789*t**2*LOG(rh)**2 + 5.666608424980593E-8*t**3*LOG(rh)**2 + & 258 258 & (0.03384437400744206*LOG(rh)**2)/x + 0.04743226764572505*LOG(rh)**3 - & 259 259 & 0.0006251042204583412*t*LOG(rh)**3 + 2.650663328519478E-6*t**2*LOG(rh)**3 - & 260 260 & 3.674710848763778E-9*t**3*LOG(rh)**3 - (0.0002672510825259393*LOG(rh)**3)/x - & 261 261 & 0.01252108546759328*LOG(rhoa) + 0.005806550506277202*t*LOG(rhoa) - & 262 & 0.0001016735312443444*t**2*LOG(rhoa) + 2.881946187214505E-7*t**3*LOG(rhoa) + 262 & 0.0001016735312443444*t**2*LOG(rhoa) + 2.881946187214505E-7*t**3*LOG(rhoa) + & 263 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) + 264 & 0.0006723156277391984*t*LOG(rh)*LOG(rhoa) + 2.602884877659698E-6*t**2*LOG(rh)*LOG(rhoa) + & 265 265 & 1.194163699688297E-8*t**3*LOG(rh)*LOG(rhoa) - (0.00851515345806281*LOG(rh)*LOG(rhoa))/x - & 266 266 & 0.01837488495738111*LOG(rh)**2*LOG(rhoa) + 0.0001720723574407498*t*LOG(rh)**2*LOG(rhoa) - & 267 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 + 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 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 + 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 273 & 2.534598655067518E-7*t**2*LOG(rh)*LOG(rhoa)**2 - & 274 & 3.635186504599571E-10*t**3*LOG(rh)*LOG(rhoa)**2 + 274 & 3.635186504599571E-10*t**3*LOG(rh)*LOG(rhoa)**2 + & 275 275 & (0.0006100650851863252*LOG(rh)*LOG(rhoa)**2)/x + 0.0003201836700403512*LOG(rhoa)**3 - & 276 276 & 0.0000174761713262546*t*LOG(rhoa)**3 + 6.065037668052182E-8*t**2*LOG(rhoa)**3 - & … … 527 527 rc_n=0.3E-9 528 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 +&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 571 & (-6.9425267104126316E-3*LOG(rhoaln)**3)/x_n 572 572 jnuc_n=EXP(jnuc_n) 573 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 +&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 599 & 3.9423927013227455E-11*tln**3*LOG(rhoaln)**3 + (8.6136359966337272E-5*LOG(rhoaln)**3)/x_n 600 600 ntot_n=EXP(ntot_n) … … 634 634 635 635 IF (tln .LT. 185. .AND. tln .GT. 155.) THEN 636 rhoatres=1.1788859232398459E+5 - 1.0244255702550814E+4*satratln + 636 rhoatres=1.1788859232398459E+5 - 1.0244255702550814E+4*satratln + & 637 637 & 4.6815029684321962E+3*satratln**2 -1.6755952338499657E+2*tln 638 638 ENDIF … … 672 672 rc_i=0.487E-9 673 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) * &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 722 & 2.6797249816144721E-05 + (tli**(-1)) * 5.5045672663909995E-01)* satratli 723 723 jnuc_i1=EXP(jnuc_i1) 724 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+ &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 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))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 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 +&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 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 +&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 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 +&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 784 & (tli**3) * 4.7882428237444610E-15)* LOG(satratli)**(-1) + & 785 & (-4.4028846582545952E-09 + tli * 4.6541269232626618E-11 + (tli**2) * -1.1939929984285194E-13 +&785 & (-4.4028846582545952E-09 + tli * 4.6541269232626618E-11 + (tli**2) *(-1.1939929984285194E-13) + & 786 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 +&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 792 & (tli**3) * 6.6236547903091862E-13)* LOG(rhoali)**(-1) + & 793 & (-2.2808617930606012E-08 + tli * 1.4773376696847775E-10 + (tli**2) * -1.3076953119957355E-13 +&793 & (-2.2808617930606012E-08 + tli * 1.4773376696847775E-10 + (tli**2) *(-1.3076953119957355E-13) + & 794 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 +&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 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 +&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 802 & (tli**3) * 1.6376014957685404E-15)* LOG(satratli) + & 803 & (-3.5621571395968670E-09 + tli * 4.1175339587760905E-11 + (tli**2) * -1.3535372357998504E-13 +&803 & (-3.5621571395968670E-09 + tli * 4.1175339587760905E-11 + (tli**2) *(-1.3535372357998504E-13) + & 804 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+ &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 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+ &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 812 & (tli**3) * 1.4676120441783832E-18)* LOG(satratli)**2 * LOG(rhoali) + & 813 & (6.7565715216420310E-13 + tli * -3.5421162549480807E-15 + (tli**2) * -3.4201196868693569E-18+ &813 & (6.7565715216420310E-13 + tli *(-3.5421162549480807E-15) + (tli**2) *(-3.4201196868693569E-18) + & 814 814 & (tli**3) * 2.2260187650412392E-20)* LOG(satratli)**3 * LOG(rhoali) 815 815 -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/ocs_to_so2.F90
r3605 r3798 6 6 USE dimphy, ONLY : klon,klev 7 7 USE aerophys 8 USE infotrac 8 USE infotrac_phy 9 9 USE YOMCST, ONLY : RG 10 10 USE phys_local_var_mod, ONLY : OCS_lifetime, budg_3D_ocs_to_so2, budg_ocs_to_so2 -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/so2_to_h2so4.F90
r3605 r3798 6 6 USE dimphy, ONLY : klon,klev 7 7 USE aerophys 8 USE infotrac 8 USE infotrac_phy 9 9 USE YOMCST, ONLY : RG 10 10 USE phys_local_var_mod, ONLY : SO2_lifetime, budg_3D_so2_to_h2so4, budg_so2_to_h2so4 -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/sulfate_aer_mod.F90
r3605 r3798 50 50 ! 51 51 SAVE INSTEP,F,XC,YC,XC1,XC16,YC1,YC28 52 !$OMP THREADPRIVATE(INSTEP,F,XC,YC,XC1,XC16,YC1,YC28) 52 53 53 54 ! convert pplay (in Pa) to PMB (in mb) -
LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/traccoag_mod.F90
r3605 r3798 17 17 18 18 USE dimphy 19 USE infotrac 19 USE infotrac_phy 20 20 USE aerophys 21 21 USE geometry_mod, ONLY : cell_area, boundslat
Note: See TracChangeset
for help on using the changeset viewer.
