Skip to content
Snippets Groups Projects
Commit 454c2954 authored by Yves Revaz's avatar Yves Revaz Committed by Matthieu Schaller
Browse files

Grackle cooling primordial abundances

parent 5314ccde
No related branches found
No related tags found
7 merge requests!1997Yet another master into zoom buffer branch merge (with formatting all done!),!1994Draft: Testing master -> zoom merge,!1987Update zoom_merge with master updates (after wrangling immense conflicts),!1982Update zoom merge branch with latest master developments,!1956Rename space_getsid to space_getsid_and_swap_cells() to try to prevent...,!1945Another master->zoom_merge update,!1926Grackle cooling primordial abundances
Hide whitespace changes
Inline Side-by-side
Showing
with 880 additions and 21 deletions
examples/Cooling/CoolingWithPrimordialElements/CoolingWithPrimordialElements.yml 0 → 100644
+ 108
− 0
View file @ 454c2954
# Define some meta-data about the simulation
MetaData:
run_name: CosmoBox00
# Define the system of units to use internally.
InternalUnitSystem:
UnitMass_in_cgs: 1.98841e43 # 10^10 M_sun in grams
UnitLength_in_cgs: 3.08567758e24 # Mpc in centimeters
UnitVelocity_in_cgs: 1e5 # km/s in centimeters per second
UnitCurrent_in_cgs: 1 # Amperes
UnitTemp_in_cgs: 1 # Kelvin
# Cosmological parameters
Cosmology:
h: 0.67321 # Reduced qHubble constant
a_begin: 0.01282051 # Initial scale-factor of the simulation (z=77)
a_end: 1.0 # Final scale factor of the simulation
Omega_cdm: 0.26641139 # Cold Dark Matter density parameter (value without neutrinos : 0.26499)
Omega_lambda: 0.6842 # Dark-energy density parameter
Omega_b: 0.049389 # Baryon density parameter
# Parameters governing the time integration
TimeIntegration:
dt_min: 1e-10 # The minimal time-step size of the simulation (in internal units).
dt_max: 1e-2 # The maximal time-step size of the simulation (in internal units).
# Parameters governing the snapshots
Snapshots:
basename: snap/snapshot # Common part of the name of output files
output_list_on: 1
output_list: ./output_list.txt
compression: 4
# Parameters governing the conserved quantities statistics
Statistics:
delta_time: 1.01
scale_factor_first: 0.02
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025 # Constant dimensionless multiplier for time integration.
MAC: adaptive # Use the geometric opening angle condition
theta_cr: 0.7 # Opening angle (Multipole acceptance criterion)
epsilon_fmm: 0.001 # Adaptive opening angle
mesh_side_length: 256
comoving_DM_softening: 0.0066 # Comoving softening for DM
max_physical_DM_softening: 0.00174 # Physical softening for DM
comoving_baryon_softening: 0.0066 # Comoving softening for DM
max_physical_baryon_softening: 0.00174 # Physical softening for DM
# Parameters for the hydrodynamics scheme
SPH:
resolution_eta: 1.2348 # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
h_min_ratio: 0.01 # Minimal smoothing length in units of softening.
h_max: 0.5 # Maximal smoothing length in co-moving internal units.
CFL_condition: 0.2 # Courant-Friedrich-Levy condition for time integration.
minimal_temperature: 10.0 # (internal units)
# Cooling with Grackle 3.0
GrackleCooling:
cloudy_table: CloudyData_UVB=HM2012.h5 # Name of the Cloudy Table (available on the grackle bitbucket repository)
with_UV_background: 0 # Enable or not the UV background
redshift: -1 # Redshift to use (-1 means time based redshift)
with_metal_cooling: 0 # Enable or not the metal cooling
max_steps: 10000 # (optional) Max number of step when computing the initial composition
convergence_limit: 1e-2 # (optional) Convergence threshold (relative) for initial composition
thermal_time_myr: 5 # (optional) Time (in Myr) for adiabatic cooling after a feedback event.
self_shielding_method: 0 # (optional) Grackle (1->3 for Grackle's ones, 0 for none and -1 for GEAR)
self_shielding_threshold_atom_per_cm3: 0.007 # Required only with GEAR's self shielding. Density threshold of the self shielding
HydrogenFractionByMass : 0.76 # Hydrogen fraction by mass (default is 0.76)
use_radiative_transfer : 0 # Arrays of ionization and heating rates are provided
H2_three_body_rate : 1 # (optional) H2 formation three body rate (see grackle documentation)
H2_cie_cooling : 0 # Enable/disable H2 collision-induced emission cooling from Ripamonti & Abel (2004)
cmb_temperature_floor : 1 # Enable/disable an effective CMB temperature floor
initial_nHII_to_nH_ratio: 2.393761E-004 # initial nHII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeI_to_nH_ratio: 7.900000E-002 # initial nHeI to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeII_to_nH_ratio: 5.692753E-020 # initial nHeII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeIII_to_nH_ratio: 0 # initial nHeIII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nDI_to_nH_ratio: 2.529538E-005 # initial nDI to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nDII_to_nH_ratio: 4.284496E-009 # initial nDII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHM_to_nH_ratio: 6.509693E-012 # initial nHM to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nH2I_to_nH_ratio: 1.891128E-006 # initial nH2I to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nH2II_to_nH_ratio: 4.153181E-014 # initial nH2II to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHDI_to_nH_ratio: 3.356904E-010 # initial nHDI to nH ratio (number density ratio). Value is ignored if set to -1.
Scheduler:
max_top_level_cells: 12
#cell_split_size: 200
Restarts:
onexit: 1
delta_hours: 6.0
max_run_time: 71.5 # Three days minus fergie time
#resubmit_on_exit: 1
#resubmit_command: ./resub.sh
# Parameters related to the initial conditions
InitialConditions:
file_name: homogeneousCube_32.hdf5
periodic: 1
cleanup_h_factors: 0 # Remove the h-factors inherited from Gadget
cleanup_velocity_factors: 0 # Remove the sqrt(a) factor in the velocities inherited from Gadget
examples/Cooling/CoolingWithPrimordialElements/README 0 → 100644
+ 11
− 0
View file @ 454c2954
Runs a uniform cosmological box with imposed abundances for a set of
elements entering the Grackle cooling.
The simulation is run with the Grackle thermochemistry, down to redhift 0
where the resulting abundances are compared with predictions.
To run, configure swift with:
./configure --with-cooling=grackle_3 --with-grackle=$GRACKLE_ROOT
Note : the script getPrimordialAbundances.py helps in the setting of the
initial abundances.
examples/Cooling/CoolingWithPrimordialElements/faure2024.dat 0 → 100644
+ 193
− 0
View file @ 454c2954
! Fractional density (n(X)/nH)
! 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Redshift ctime Trad Tn nb nH e- H He D H2 HD D2 H+ D+ He+ H2+ HD+ D2+ HeH+ H3+ H2D+ HD2+ D3+ H-
3.401000E+003 0.000000E+000 9.269426E+003 9.269426E+003 8.151816E+003 7.547801E+003 1.080025E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.000000E+000 2.530000E-005 8.000000E-002 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000
3.400876E+003 3.168879E+000 9.269089E+003 9.269088E+003 8.150927E+003 7.546978E+003 1.080024E+000 1.796857E-010 1.705857E-006 4.546047E-015 6.941758E-036 3.211260E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.554308E-028 0.000000E+000 0.000000E+000 1.268827E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.306460E-028
3.400867E+003 3.485767E+000 9.269064E+003 9.269063E+003 8.150861E+003 7.546917E+003 1.080024E+000 1.796933E-010 1.705992E-006 4.546240E-015 6.942382E-036 3.211549E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.554494E-028 0.000000E+000 0.000000E+000 1.268924E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.306636E-028
3.400858E+003 3.834344E+000 9.269039E+003 9.269038E+003 8.150795E+003 7.546856E+003 1.080024E+000 1.797009E-010 1.706126E-006 4.546433E-015 6.943006E-036 3.211838E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.554679E-028 0.000000E+000 0.000000E+000 1.269021E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.306812E-028
3.400842E+003 4.217778E+000 9.268995E+003 9.268994E+003 8.150680E+003 7.546749E+003 1.080024E+000 1.797142E-010 1.706362E-006 4.546770E-015 6.944097E-036 3.212344E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.555003E-028 0.000000E+000 0.000000E+000 1.269191E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.307120E-028
3.400825E+003 4.639556E+000 9.268948E+003 9.268947E+003 8.150556E+003 7.546635E+003 1.080024E+000 1.797285E-010 1.706613E-006 4.547131E-015 6.945265E-036 3.212885E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.555350E-028 0.000000E+000 0.000000E+000 1.269372E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.307450E-028
3.400808E+003 5.103511E+000 9.268902E+003 9.268901E+003 8.150433E+003 7.546521E+003 1.080024E+000 1.797427E-010 1.706865E-006 4.547491E-015 6.946434E-036 3.213426E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.555697E-028 0.000000E+000 0.000000E+000 1.269554E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.307780E-028
3.400790E+003 5.613863E+000 9.268855E+003 9.268854E+003 8.150310E+003 7.546406E+003 1.080024E+000 1.797570E-010 1.707117E-006 4.547852E-015 6.947602E-036 3.213968E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.556044E-028 0.000000E+000 0.000000E+000 1.269735E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.308109E-028
3.400771E+003 6.175249E+000 9.268802E+003 9.268801E+003 8.150170E+003 7.546277E+003 1.080024E+000 1.797732E-010 1.707403E-006 4.548261E-015 6.948928E-036 3.214582E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.556437E-028 0.000000E+000 0.000000E+000 1.269941E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.308483E-028
3.400745E+003 6.792774E+000 9.268730E+003 9.268729E+003 8.149981E+003 7.546102E+003 1.080024E+000 1.797950E-010 1.707789E-006 4.548813E-015 6.950717E-036 3.215410E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.556968E-028 0.000000E+000 0.000000E+000 1.270219E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.308988E-028
3.400723E+003 7.472051E+000 9.268670E+003 9.268669E+003 8.149822E+003 7.545955E+003 1.080024E+000 1.798133E-010 1.708113E-006 4.549277E-015 6.952221E-036 3.216107E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.557414E-028 0.000000E+000 0.000000E+000 1.270453E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.309412E-028
3.400696E+003 8.219256E+000 9.268598E+003 9.268597E+003 8.149633E+003 7.545779E+003 1.080024E+000 1.798353E-010 1.708500E-006 4.549832E-015 6.954020E-036 3.216940E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.557948E-028 0.000000E+000 0.000000E+000 1.270732E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.309919E-028
3.400664E+003 9.041182E+000 9.268511E+003 9.268510E+003 8.149402E+003 7.545566E+003 1.080024E+000 1.798620E-010 1.708972E-006 4.550508E-015 6.956208E-036 3.217954E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.558597E-028 0.000000E+000 0.000000E+000 1.271072E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.310537E-028
3.400633E+003 9.945300E+000 9.268424E+003 9.268423E+003 8.149174E+003 7.545355E+003 1.080024E+000 1.798883E-010 1.709437E-006 4.551174E-015 6.958368E-036 3.218954E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.559238E-028 0.000000E+000 0.000000E+000 1.271407E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.311145E-028
3.400592E+003 1.093983E+001 9.268314E+003 9.268313E+003 8.148882E+003 7.545084E+003 1.080024E+000 1.799221E-010 1.710035E-006 4.552029E-015 6.961143E-036 3.220240E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.560060E-028 0.000000E+000 0.000000E+000 1.271838E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.311927E-028
3.400549E+003 1.203381E+001 9.268196E+003 9.268195E+003 8.148573E+003 7.544798E+003 1.080024E+000 1.799579E-010 1.710668E-006 4.552934E-015 6.964078E-036 3.221599E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.560931E-028 0.000000E+000 0.000000E+000 1.272294E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.312755E-028
3.400507E+003 1.323719E+001 9.268081E+003 9.268080E+003 8.148268E+003 7.544516E+003 1.080024E+000 1.799931E-010 1.711291E-006 4.553826E-015 6.966971E-036 3.222939E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.561788E-028 0.000000E+000 0.000000E+000 1.272743E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.313570E-028
3.400462E+003 1.456091E+001 9.267959E+003 9.267958E+003 8.147946E+003 7.544218E+003 1.080024E+000 1.800305E-010 1.711952E-006 4.554771E-015 6.970036E-036 3.224359E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.562696E-028 0.000000E+000 0.000000E+000 1.273219E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.314433E-028
3.400409E+003 1.601701E+001 9.267814E+003 9.267813E+003 8.147565E+003 7.543865E+003 1.080024E+000 1.800746E-010 1.712733E-006 4.555888E-015 6.973662E-036 3.226039E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.563771E-028 0.000000E+000 0.000000E+000 1.273782E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.315454E-028
3.400349E+003 1.761871E+001 9.267652E+003 9.267651E+003 8.147137E+003 7.543469E+003 1.080024E+000 1.801242E-010 1.713610E-006 4.557142E-015 6.977733E-036 3.227925E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.564977E-028 0.000000E+000 0.000000E+000 1.274414E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.316600E-028
3.400284E+003 1.938058E+001 9.267475E+003 9.267474E+003 8.146670E+003 7.543036E+003 1.080024E+000 1.801783E-010 1.714568E-006 4.558511E-015 6.982181E-036 3.229985E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999829E-002 4.566294E-028 0.000000E+000 0.000000E+000 1.275105E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.317852E-028
3.400213E+003 2.131863E+001 9.267281E+003 9.267280E+003 8.146160E+003 7.542564E+003 1.080024E+000 1.802375E-010 1.715616E-006 4.560008E-015 6.987044E-036 3.232237E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.567733E-028 0.000000E+000 0.000000E+000 1.275859E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.319220E-028
3.400135E+003 2.345050E+001 9.267068E+003 9.267067E+003 8.145596E+003 7.542042E+003 1.080024E+000 1.803029E-010 1.716774E-006 4.561662E-015 6.992421E-036 3.234728E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.569323E-028 0.000000E+000 0.000000E+000 1.276694E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.320732E-028
3.400048E+003 2.579555E+001 9.266831E+003 9.266830E+003 8.144972E+003 7.541464E+003 1.080024E+000 1.803753E-010 1.718057E-006 4.563495E-015 6.998379E-036 3.237488E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.571085E-028 0.000000E+000 0.000000E+000 1.277618E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.322407E-028
3.399950E+003 2.837510E+001 9.266563E+003 9.266562E+003 8.144265E+003 7.540809E+003 1.080024E+000 1.804574E-010 1.719512E-006 4.565571E-015 7.005134E-036 3.240617E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.573082E-028 0.000000E+000 0.000000E+000 1.278666E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.324304E-028
3.399845E+003 3.121261E+001 9.266278E+003 9.266277E+003 8.143515E+003 7.540115E+003 1.080024E+000 1.805445E-010 1.721057E-006 4.567777E-015 7.012313E-036 3.243942E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.575203E-028 0.000000E+000 0.000000E+000 1.279780E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.326320E-028
3.399730E+003 3.433387E+001 9.265964E+003 9.265964E+003 8.142687E+003 7.539349E+003 1.080024E+000 1.806407E-010 1.722763E-006 4.570210E-015 7.020235E-036 3.247612E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.577542E-028 0.000000E+000 0.000000E+000 1.281009E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.328543E-028
3.399586E+003 3.776726E+001 9.265571E+003 9.265570E+003 8.141651E+003 7.538389E+003 1.080024E+000 1.807613E-010 1.724903E-006 4.573260E-015 7.030175E-036 3.252216E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999828E-002 4.580475E-028 0.000000E+000 0.000000E+000 1.282550E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.331331E-028
3.399460E+003 4.154399E+001 9.265230E+003 9.265229E+003 8.140750E+003 7.537555E+003 1.080024E+000 1.808661E-010 1.726765E-006 4.575913E-015 7.038829E-036 3.256225E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999827E-002 4.583026E-028 0.000000E+000 0.000000E+000 1.283892E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.333756E-028
3.399310E+003 4.569838E+001 9.264820E+003 9.264819E+003 8.139670E+003 7.536555E+003 1.080024E+000 1.809921E-010 1.729002E-006 4.579100E-015 7.049228E-036 3.261041E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999827E-002 4.586090E-028 0.000000E+000 0.000000E+000 1.285504E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.336668E-028
3.399134E+003 5.026822E+001 9.264341E+003 9.264340E+003 8.138409E+003 7.535387E+003 1.080024E+000 1.811391E-010 1.731616E-006 4.582820E-015 7.061378E-036 3.266670E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999827E-002 4.589667E-028 0.000000E+000 0.000000E+000 1.287387E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.340067E-028
3.398944E+003 5.529505E+001 9.263823E+003 9.263823E+003 8.137044E+003 7.534124E+003 1.080024E+000 1.812985E-010 1.734451E-006 4.586851E-015 7.074558E-036 3.272775E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999827E-002 4.593543E-028 0.000000E+000 0.000000E+000 1.289429E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.343751E-028
3.398754E+003 6.082455E+001 9.263304E+003 9.263303E+003 8.135675E+003 7.532856E+003 1.080024E+000 1.814586E-010 1.737302E-006 4.590902E-015 7.087816E-036 3.278916E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999826E-002 4.597437E-028 0.000000E+000 0.000000E+000 1.291482E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.347453E-028
3.398528E+003 6.690701E+001 9.262688E+003 9.262687E+003 8.134052E+003 7.531353E+003 1.080024E+000 1.816485E-010 1.740686E-006 4.595707E-015 7.103558E-036 3.286208E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999826E-002 4.602057E-028 0.000000E+000 0.000000E+000 1.293920E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.351844E-028
3.398284E+003 7.359771E+001 9.262023E+003 9.262022E+003 8.132301E+003 7.529732E+003 1.080024E+000 1.818536E-010 1.744345E-006 4.600897E-015 7.120582E-036 3.294094E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999826E-002 4.607047E-028 0.000000E+000 0.000000E+000 1.296555E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.356586E-028
3.398001E+003 8.095748E+001 9.261253E+003 9.261252E+003 8.130273E+003 7.527854E+003 1.080024E+000 1.820918E-010 1.748597E-006 4.606922E-015 7.140372E-036 3.303261E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999825E-002 4.612839E-028 0.000000E+000 0.000000E+000 1.299617E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.362091E-028
3.397713E+003 8.905322E+001 9.260468E+003 9.260467E+003 8.128205E+003 7.525940E+003 1.080024E+000 1.823348E-010 1.752942E-006 4.613072E-015 7.160601E-036 3.312632E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999825E-002 4.618750E-028 0.000000E+000 0.000000E+000 1.302746E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.367709E-028
3.397387E+003 9.795855E+001 9.259577E+003 9.259577E+003 8.125861E+003 7.523769E+003 1.080024E+000 1.826110E-010 1.757883E-006 4.620058E-015 7.183620E-036 3.323295E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999824E-002 4.625466E-028 0.000000E+000 0.000000E+000 1.306305E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.374092E-028
3.397027E+003 1.077544E+002 9.258597E+003 9.258596E+003 8.123279E+003 7.521379E+003 1.080024E+000 1.829156E-010 1.763343E-006 4.627766E-015 7.209061E-036 3.335080E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999824E-002 4.632876E-028 0.000000E+000 0.000000E+000 1.310236E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.381134E-028
3.396623E+003 1.185298E+002 9.257497E+003 9.257496E+003 8.120386E+003 7.518699E+003 1.080024E+000 1.832579E-010 1.769487E-006 4.636426E-015 7.237702E-036 3.348348E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999823E-002 4.641200E-028 0.000000E+000 0.000000E+000 1.314659E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.389045E-028
3.396173E+003 1.303828E+002 9.256269E+003 9.256268E+003 8.117154E+003 7.515708E+003 1.080024E+000 1.836411E-010 1.776375E-006 4.646119E-015 7.269834E-036 3.363232E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999822E-002 4.650517E-028 0.000000E+000 0.000000E+000 1.319618E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.397900E-028
3.395667E+003 1.434211E+002 9.254892E+003 9.254891E+003 8.113531E+003 7.512353E+003 1.080024E+000 1.840719E-010 1.784135E-006 4.657019E-015 7.306053E-036 3.380010E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999822E-002 4.660994E-028 0.000000E+000 0.000000E+000 1.325204E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.407856E-028
3.395176E+003 1.577632E+002 9.253553E+003 9.253552E+003 8.110012E+003 7.509094E+003 1.080024E+000 1.844916E-010 1.791711E-006 4.667638E-015 7.341432E-036 3.396399E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999821E-002 4.671200E-028 0.000000E+000 0.000000E+000 1.330656E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.417556E-028
3.394607E+003 1.735395E+002 9.252001E+003 9.252000E+003 8.105931E+003 7.505316E+003 1.080023E+000 1.849797E-010 1.800539E-006 4.679987E-015 7.382689E-036 3.415511E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999820E-002 4.683068E-028 0.000000E+000 0.000000E+000 1.337010E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.428834E-028
3.393944E+003 1.908935E+002 9.250193E+003 9.250193E+003 8.101181E+003 7.500918E+003 1.080023E+000 1.855501E-010 1.810880E-006 4.694416E-015 7.431051E-036 3.437914E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999819E-002 4.696934E-028 0.000000E+000 0.000000E+000 1.344450E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.442011E-028
3.393270E+003 2.099828E+002 9.248357E+003 9.248356E+003 8.096357E+003 7.496451E+003 1.080023E+000 1.861315E-010 1.821450E-006 4.709126E-015 7.480525E-036 3.460833E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999818E-002 4.711069E-028 0.000000E+000 0.000000E+000 1.352055E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.455444E-028
3.392492E+003 2.309811E+002 9.246238E+003 9.246237E+003 8.090794E+003 7.491300E+003 1.080023E+000 1.868049E-010 1.833727E-006 4.726163E-015 7.538043E-036 3.487479E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999817E-002 4.727439E-028 0.000000E+000 0.000000E+000 1.360886E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.470999E-028
3.391652E+003 2.540792E+002 9.243949E+003 9.243948E+003 8.084785E+003 7.485737E+003 1.080023E+000 1.875355E-010 1.847092E-006 4.744649E-015 7.600719E-036 3.516514E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999815E-002 4.745199E-028 0.000000E+000 0.000000E+000 1.370497E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.487876E-028
3.390717E+003 2.794872E+002 9.241399E+003 9.241398E+003 8.078097E+003 7.479545E+003 1.080023E+000 1.883531E-010 1.862099E-006 4.765334E-015 7.671178E-036 3.549155E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999814E-002 4.765071E-028 0.000000E+000 0.000000E+000 1.381288E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.506759E-028
3.389698E+003 3.074359E+002 9.238623E+003 9.238622E+003 8.070818E+003 7.472805E+003 1.080023E+000 1.892481E-010 1.878592E-006 4.787977E-015 7.748699E-036 3.585068E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999812E-002 4.786822E-028 0.000000E+000 0.000000E+000 1.393144E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.527426E-028
3.388580E+003 3.381795E+002 9.235576E+003 9.235575E+003 8.062837E+003 7.465415E+003 1.080023E+000 1.902355E-010 1.896865E-006 4.812958E-015 7.834710E-036 3.624915E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999810E-002 4.810816E-028 0.000000E+000 0.000000E+000 1.406277E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.550225E-028
3.387340E+003 3.719974E+002 9.232196E+003 9.232195E+003 8.053987E+003 7.457221E+003 1.080023E+000 1.913380E-010 1.917366E-006 4.840852E-015 7.931343E-036 3.669684E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999808E-002 4.837605E-028 0.000000E+000 0.000000E+000 1.421006E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.575677E-028
3.385987E+003 4.091972E+002 9.228508E+003 9.228507E+003 8.044338E+003 7.448286E+003 1.080023E+000 1.925493E-010 1.940005E-006 4.871498E-015 8.038232E-036 3.719204E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999806E-002 4.867032E-028 0.000000E+000 0.000000E+000 1.437268E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.603635E-028
3.384476E+003 4.501169E+002 9.224391E+003 9.224390E+003 8.033576E+003 7.438322E+003 1.080023E+000 1.939116E-010 1.965613E-006 4.905963E-015 8.159354E-036 3.775320E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999803E-002 4.900123E-028 0.000000E+000 0.000000E+000 1.455655E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.635073E-028
3.382854E+003 4.951286E+002 9.219968E+003 9.219967E+003 8.022028E+003 7.427629E+003 1.080023E+000 1.953871E-010 1.993524E-006 4.943293E-015 8.291628E-036 3.836604E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999801E-002 4.935959E-028 0.000000E+000 0.000000E+000 1.475689E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.669116E-028
3.381056E+003 5.446414E+002 9.215068E+003 9.215067E+003 8.009243E+003 7.415792E+003 1.080023E+000 1.970369E-010 2.024948E-006 4.985034E-015 8.440868E-036 3.905750E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999798E-002 4.976022E-028 0.000000E+000 0.000000E+000 1.498236E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.707174E-028
3.379078E+003 5.991056E+002 9.209678E+003 9.209677E+003 7.995198E+003 7.402788E+003 1.080023E+000 1.988697E-010 2.060123E-006 5.031404E-015 8.608317E-036 3.983334E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999794E-002 5.020520E-028 0.000000E+000 0.000000E+000 1.523463E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.749441E-028
3.376897E+003 6.590161E+002 9.203734E+003 9.203733E+003 7.979728E+003 7.388464E+003 1.080023E+000 2.009134E-010 2.099674E-006 5.083110E-015 8.797093E-036 4.070802E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999790E-002 5.070128E-028 0.000000E+000 0.000000E+000 1.551814E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.796559E-028
3.374528E+003 7.249177E+002 9.197276E+003 9.197275E+003 7.962943E+003 7.372922E+003 1.080023E+000 2.031608E-010 2.143566E-006 5.139968E-015 9.007193E-036 4.168153E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999786E-002 5.124668E-028 0.000000E+000 0.000000E+000 1.583262E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.848357E-028
3.371891E+003 7.974095E+002 9.190088E+003 9.190087E+003 7.944286E+003 7.355648E+003 1.080023E+000 2.056959E-010 2.193581E-006 5.204107E-015 9.247356E-036 4.279438E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999781E-002 5.186176E-028 0.000000E+000 0.000000E+000 1.619075E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.906768E-028
3.369041E+003 8.771505E+002 9.182322E+003 9.182321E+003 7.924164E+003 7.337017E+003 1.080023E+000 2.084750E-010 2.249020E-006 5.274418E-015 9.514490E-036 4.403225E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999775E-002 5.253586E-028 0.000000E+000 0.000000E+000 1.658748E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.970777E-028
3.365875E+003 9.648655E+002 9.173692E+003 9.173692E+003 7.901843E+003 7.316350E+003 1.080023E+000 2.116131E-010 2.312387E-006 5.353812E-015 9.820992E-036 4.545260E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999769E-002 5.329682E-028 0.000000E+000 0.000000E+000 1.704062E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.043027E-028
3.362424E+003 1.061352E+003 9.164287E+003 9.164286E+003 7.877564E+003 7.293869E+003 1.080023E+000 2.150944E-010 2.383633E-006 5.441889E-015 1.016705E-035 4.705634E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999762E-002 5.414074E-028 0.000000E+000 0.000000E+000 1.754972E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.123143E-028
3.358634E+003 1.167487E+003 9.153957E+003 9.153956E+003 7.850956E+003 7.269233E+003 1.080023E+000 2.189926E-010 2.464592E-006 5.540512E-015 1.056212E-035 4.888726E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999754E-002 5.508539E-028 0.000000E+000 0.000000E+000 1.812776E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.212810E-028
3.354466E+003 1.284236E+003 9.142597E+003 9.142596E+003 7.821763E+003 7.242203E+003 1.080023E+000 2.233721E-010 2.557032E-006 5.651314E-015 1.101551E-035 5.098861E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999744E-002 5.614630E-028 0.000000E+000 0.000000E+000 1.878717E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.313496E-028
3.349924E+003 1.412660E+003 9.130217E+003 9.130217E+003 7.790032E+003 7.212824E+003 1.080023E+000 2.282579E-010 2.662010E-006 5.774925E-015 1.153328E-035 5.338847E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999734E-002 5.732940E-028 0.000000E+000 0.000000E+000 1.953527E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.425759E-028
3.344917E+003 1.553926E+003 9.116573E+003 9.116572E+003 7.755159E+003 7.180534E+003 1.080023E+000 2.337833E-010 2.783072E-006 5.914719E-015 1.213407E-035 5.617330E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999722E-002 5.866680E-028 0.000000E+000 0.000000E+000 2.039708E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.552639E-028
3.339472E+003 1.709318E+003 9.101731E+003 9.101730E+003 7.717344E+003 7.145522E+003 1.080022E+000 2.399658E-010 2.921468E-006 6.071136E-015 1.282554E-035 5.937868E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999708E-002 6.016254E-028 0.000000E+000 0.000000E+000 2.138110E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.694510E-028
3.333475E+003 1.880250E+003 9.085387E+003 9.085386E+003 7.675845E+003 7.107098E+003 1.080022E+000 2.469888E-010 3.082426E-006 6.248817E-015 1.363574E-035 6.313475E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999692E-002 6.186076E-028 0.000000E+000 0.000000E+000 2.252409E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.855547E-028
3.326913E+003 2.068275E+003 9.067501E+003 9.067500E+003 7.630600E+003 7.065205E+003 1.080022E+000 2.549421E-010 3.269501E-006 6.450035E-015 1.458517E-035 6.753661E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999673E-002 6.378291E-028 0.000000E+000 0.000000E+000 2.385065E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.037768E-028
3.319747E+003 2.275102E+003 9.047971E+003 9.047970E+003 7.581401E+003 7.019652E+003 1.080022E+000 2.639580E-010 3.487699E-006 6.678138E-015 1.570258E-035 7.271781E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999651E-002 6.596062E-028 0.000000E+000 0.000000E+000 2.539556E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.244151E-028
3.311914E+003 2.502613E+003 9.026622E+003 9.026621E+003 7.527862E+003 6.970079E+003 1.080022E+000 2.742283E-010 3.744148E-006 6.937975E-015 1.702897E-035 7.886861E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999626E-002 6.843977E-028 0.000000E+000 0.000000E+000 2.720826E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.479017E-028
3.303385E+003 2.752874E+003 9.003377E+003 9.003376E+003 7.469855E+003 6.916370E+003 1.080021E+000 2.859268E-010 4.046446E-006 7.233948E-015 1.860963E-035 8.619929E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999595E-002 7.126185E-028 0.000000E+000 0.000000E+000 2.934120E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.746263E-028
3.294034E+003 3.028161E+003 8.977890E+003 8.977889E+003 7.406596E+003 6.857799E+003 1.080021E+000 2.994054E-010 4.408109E-006 7.574957E-015 2.052350E-035 9.507644E-040 0.000000E+000 1.000000E+000 2.530000E-005 7.999559E-002 7.451109E-028 0.000000E+000 0.000000E+000 3.188798E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.053820E-028
3.283807E+003 3.330977E+003 8.950017E+003 8.950016E+003 7.337827E+003 6.794126E+003 1.080020E+000 3.149724E-010 4.843492E-006 7.968801E-015 2.285819E-035 1.059069E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999516E-002 7.826100E-028 0.000000E+000 0.000000E+000 3.494735E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.408579E-028
3.272705E+003 3.664075E+003 8.919758E+003 8.919757E+003 7.263652E+003 6.725446E+003 1.080020E+000 3.329145E-010 5.368646E-006 8.422737E-015 2.571554E-035 1.191639E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999463E-002 8.257973E-028 0.000000E+000 0.000000E+000 3.862900E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.816908E-028
3.260624E+003 4.030483E+003 8.886831E+003 8.886830E+003 7.183508E+003 6.651241E+003 1.080019E+000 3.537590E-010 6.009879E-006 8.950102E-015 2.926061E-035 1.356143E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999399E-002 8.759305E-028 0.000000E+000 0.000000E+000 4.311323E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.290586E-028
3.247426E+003 4.433531E+003 8.850859E+003 8.850858E+003 7.096629E+003 6.570799E+003 1.080018E+000 3.782288E-010 6.804944E-006 9.569189E-015 3.373417E-035 1.563769E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999320E-002 9.347344E-028 0.000000E+000 0.000000E+000 4.865826E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.845752E-028
3.233108E+003 4.876884E+003 8.811836E+003 8.811835E+003 7.003178E+003 6.484272E+003 1.080018E+000 4.069526E-010 7.795943E-006 1.029590E-014 3.941925E-035 1.827675E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999220E-002 1.003703E-027 0.000000E+000 0.000000E+000 5.554975E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 9.496292E-028
3.217536E+003 5.364572E+003 8.769394E+003 8.769393E+003 6.902471E+003 6.391027E+003 1.080016E+000 4.410151E-010 9.051038E-006 1.115768E-014 4.677445E-035 2.169181E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.999095E-002 1.085419E-027 0.000000E+000 0.000000E+000 6.425057E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.026629E-027
3.200676E+003 5.901030E+003 8.723442E+003 8.723441E+003 6.794530E+003 6.291084E+003 1.080015E+000 4.815571E-010 1.065653E-005 1.218339E-014 5.640617E-035 2.616493E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.998934E-002 1.182598E-027 0.000000E+000 0.000000E+000 7.534329E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.118089E-027
3.182373E+003 6.491133E+003 8.673558E+003 8.673557E+003 6.678635E+003 6.183776E+003 1.080013E+000 5.303788E-010 1.274882E-005 1.341859E-014 6.928636E-035 3.214819E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.998725E-002 1.299526E-027 0.000000E+000 0.000000E+000 8.974783E-024 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.227986E-027
3.162576E+003 7.140246E+003 8.619600E+003 8.619599E+003 6.554767E+003 6.069087E+003 1.080010E+000 5.895455E-010 1.551381E-005 1.491550E-014 8.679808E-035 4.028523E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.998449E-002 1.441119E-027 0.000000E+000 0.000000E+000 1.087111E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.360854E-027
3.141198E+003 7.854270E+003 8.561335E+003 8.561334E+003 6.422741E+003 5.946843E+003 1.080006E+000 6.618974E-010 1.923109E-005 1.674600E-014 1.110883E-034 5.157545E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.998077E-002 1.614142E-027 0.000000E+000 0.000000E+000 1.341024E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.522920E-027
3.118088E+003 8.639697E+003 8.498350E+003 8.498349E+003 6.282027E+003 5.816555E+003 1.080001E+000 7.515217E-010 2.433956E-005 1.901350E-014 1.456383E-034 6.763981E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.997566E-002 1.828341E-027 0.000000E+000 0.000000E+000 1.688471E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.723130E-027
3.093267E+003 9.503667E+003 8.430700E+003 8.430699E+003 6.133196E+003 5.678752E+003 1.079994E+000 8.632373E-010 3.147266E-005 2.183991E-014 1.957434E-034 9.094524E-039 0.000000E+000 1.000000E+000 2.530000E-005 7.996853E-002 2.095214E-027 0.000000E+000 0.000000E+000 2.171437E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.971966E-027
3.066538E+003 1.045403E+004 8.357849E+003 8.357848E+003 5.975571E+003 5.532806E+003 1.079984E+000 1.004794E-009 4.170518E-005 2.542129E-014 2.706796E-034 1.258144E-038 0.000000E+000 1.000000E+000 2.530000E-005 7.995829E-002 2.433284E-027 0.000000E+000 0.000000E+000 2.860991E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.286296E-027
3.037873E+003 1.149944E+004 8.279723E+003 8.279722E+003 5.809561E+003 5.379097E+003 1.079969E+000 1.186138E-009 5.671800E-005 3.000929E-014 3.857919E-034 1.794016E-038 0.000000E+000 1.000000E+000 2.530000E-005 7.994328E-002 2.866375E-027 0.000000E+000 0.000000E+000 3.867710E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.687652E-027
3.006971E+003 1.264938E+004 8.195500E+003 8.195499E+003 5.634070E+003 5.216609E+003 1.079946E+000 1.423677E-009 7.953418E-005 3.601904E-014 5.698576E-034 2.651294E-038 0.000000E+000 1.000000E+000 2.530000E-005 7.992047E-002 3.433866E-027 0.000000E+000 0.000000E+000 5.389900E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.211546E-027
2.974293E+003 1.391432E+004 8.106437E+003 8.106435E+003 5.452377E+003 5.048379E+003 1.079911E+000 1.734197E-009 1.145930E-004 4.387519E-014 8.689447E-034 4.045009E-038 0.000000E+000 1.000000E+000 2.530000E-005 7.988541E-002 4.176272E-027 0.000000E+000 0.000000E+000 7.716660E-023 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.893875E-027
2.939323E+003 1.530575E+004 8.011125E+003 8.011124E+003 5.262309E+003 4.872394E+003 1.079854E+000 2.152673E-009 1.709143E-004 5.446265E-014 1.379965E-033 6.427667E-038 0.000000E+000 1.000000E+000 2.530000E-005 7.982909E-002 5.178073E-027 0.000000E+000 0.000000E+000 1.143478E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.809855E-027
2.902049E+003 1.683633E+004 7.909535E+003 7.909534E+003 5.064642E+003 4.689374E+003 1.079761E+000 2.726529E-009 2.644566E-004 6.898120E-014 2.288874E-033 1.066816E-037 0.000000E+000 1.000000E+000 2.530000E-005 7.973554E-002 6.554510E-027 0.000000E+000 0.000000E+000 1.757715E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.060762E-027
2.862514E+003 1.851996E+004 7.801783E+003 7.801782E+003 4.860461E+003 4.500322E+003 1.079600E+000 3.527527E-009 4.252025E-004 8.924646E-014 3.974719E-033 1.853877E-037 0.000000E+000 1.000000E+000 2.530000E-005 7.957480E-002 8.481044E-027 0.000000E+000 0.000000E+000 2.807653E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.799136E-027
2.820812E+003 2.037195E+004 7.688123E+003 7.688122E+003 4.651112E+003 4.306484E+003 1.079314E+000 4.665929E-009 7.110243E-004 1.180481E-013 7.240458E-033 3.379666E-037 0.000000E+000 1.000000E+000 2.530000E-005 7.928898E-002 1.122926E-026 0.000000E+000 0.000000E+000 4.664929E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.025806E-026
2.776511E+003 2.240915E+004 7.567380E+003 7.567379E+003 4.435397E+003 4.106753E+003 1.078780E+000 6.339030E-009 1.245274E-003 1.603776E-013 1.397458E-032 6.528472E-037 0.000000E+000 1.000000E+000 2.530000E-005 7.875473E-002 1.528807E-026 0.000000E+000 0.000000E+000 8.119388E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.385278E-026
2.730003E+003 2.465007E+004 7.440623E+003 7.440622E+003 4.216225E+003 3.903821E+003 1.077752E+000 8.838186E-009 2.273604E-003 2.236063E-013 2.852758E-032 1.333932E-036 0.000000E+000 1.000000E+000 2.530000E-005 7.772640E-002 2.138915E-026 0.000000E+000 0.000000E+000 1.473816E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.918858E-026
2.681693E+003 2.711507E+004 7.308954E+003 7.308953E+003 3.996332E+003 3.700221E+003 1.075736E+000 1.263021E-008 4.289434E-003 3.195446E-013 6.142425E-032 2.874977E-036 0.000000E+000 1.000000E+000 2.530000E-005 7.571057E-002 3.072010E-026 0.000000E+000 0.000000E+000 2.766039E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.721933E-026
2.630516E+003 2.982658E+004 7.169472E+003 7.169471E+003 3.771875E+003 3.492395E+003 1.071645E+000 1.867746E-008 8.380230E-003 4.725407E-013 1.423594E-031 6.670353E-036 0.000000E+000 1.000000E+000 2.530000E-005 7.161977E-002 4.574927E-026 0.000000E+000 0.000000E+000 5.379452E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.987738E-026
2.577675E+003 3.280924E+004 7.025454E+003 7.025453E+003 3.549107E+003 3.286133E+003 1.063811E+000 2.835652E-008 1.621403E-002 7.174219E-013 3.489977E-031 1.637148E-035 0.000000E+000 1.000000E+000 2.530000E-005 6.378597E-002 7.010289E-026 0.000000E+000 0.000000E+000 1.037132E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.978013E-026
2.522514E+003 3.609016E+004 6.875111E+003 6.875110E+003 3.326098E+003 3.079648E+003 1.050384E+000 4.449115E-008 2.964092E-002 1.125631E-012 9.174778E-031 4.309316E-035 0.000000E+000 1.000000E+000 2.530000E-005 5.035908E-002 1.113169E-025 0.000000E+000 0.000000E+000 1.891555E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 9.213721E-026
2.465474E+003 3.969918E+004 6.719649E+003 6.719649E+003 3.105530E+003 2.875423E+003 1.032626E+000 7.221282E-008 4.739956E-002 1.826998E-012 2.593958E-030 1.220022E-034 0.000000E+000 9.999999E-001 2.530000E-005 3.260044E-002 1.834268E-025 0.000000E+000 0.000000E+000 3.022347E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.463181E-025
2.406345E+003 4.366909E+004 6.558493E+003 6.558492E+003 2.887407E+003 2.673462E+003 1.016357E+000 1.226510E-007 6.366810E-002 3.103111E-012 8.105726E-030 3.818032E-034 0.000000E+000 9.999999E-001 2.530000E-005 1.633190E-002 3.174686E-025 0.000000E+000 0.000000E+000 4.063849E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.435470E-025
2.345928E+003 4.803600E+004 6.393828E+003 6.393827E+003 2.675338E+003 2.477107E+003 1.006504E+000 2.184024E-007 7.352125E-002 5.525706E-012 2.816851E-029 1.328930E-033 0.000000E+000 9.999998E-001 2.529999E-005 6.478749E-003 5.784070E-025 0.000000E+000 0.000000E+000 4.707913E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.278973E-025
2.283783E+003 5.283960E+004 6.224451E+003 6.224449E+003 2.468306E+003 2.285415E+003 1.002162E+000 4.111166E-007 7.786308E-002 1.040170E-011 1.108248E-028 5.237495E-033 0.000000E+000 9.999996E-001 2.529999E-005 2.136923E-003 1.119470E-024 0.000000E+000 0.000000E+000 5.015197E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.996486E-025
2.220005E+003 5.812356E+004 6.050625E+003 6.050623E+003 2.267235E+003 2.099243E+003 1.000636E+000 8.198028E-007 7.938820E-002 2.074279E-011 4.955900E-028 2.346501E-032 0.000000E+000 9.999992E-001 2.529998E-005 6.118041E-004 2.308267E-024 0.000000E+000 0.000000E+000 5.159509E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.588982E-024
2.155714E+003 6.393592E+004 5.875399E+003 5.875397E+003 2.075907E+003 1.922091E+003 1.000182E+000 1.717212E-006 7.984194E-002 4.345315E-011 2.472590E-027 1.173052E-031 0.000000E+000 9.999983E-001 2.529996E-005 1.580624E-004 5.028965E-024 0.000000E+000 0.000000E+000 5.253913E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.324095E-024
2.090388E+003 7.032951E+004 5.697353E+003 5.697352E+003 1.892846E+003 1.752594E+003 1.000058E+000 3.819751E-006 7.996355E-002 9.667670E-011 1.409523E-026 6.701507E-031 0.000000E+000 9.999962E-001 2.529990E-005 3.645421E-005 1.171721E-023 0.000000E+000 0.000000E+000 5.349618E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.396908E-024
2.024039E+003 7.736246E+004 5.516517E+003 5.516516E+003 1.718268E+003 1.590951E+003 1.000024E+000 9.083787E-006 7.999256E-002 2.300160E-010 9.320662E-026 4.441800E-030 0.000000E+000 9.999909E-001 2.529977E-005 7.435585E-006 2.942456E-023 0.000000E+000 0.000000E+000 5.466781E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.762688E-023
1.957371E+003 8.509871E+004 5.334815E+003 5.334813E+003 1.554010E+003 1.438865E+003 1.000004E+000 2.304835E-005 7.999865E-002 5.842124E-010 7.122032E-025 3.402529E-029 0.000000E+000 9.999770E-001 2.529942E-005 1.347633E-006 7.952322E-023 0.000000E+000 0.000000E+000 5.614358E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.489950E-023
1.890718E+003 9.360858E+004 5.153152E+003 5.153150E+003 1.400602E+003 1.296824E+003 9.999630E-001 6.253200E-005 7.999978E-002 1.587891E-009 6.322436E-024 3.028616E-028 0.000000E+000 9.999375E-001 2.529841E-005 2.162309E-007 2.319867E-022 0.000000E+000 0.000000E+000 5.799232E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.225429E-022
1.823730E+003 1.029694E+005 4.970576E+003 4.970574E+003 1.256945E+003 1.163810E+003 9.998415E-001 1.838010E-004 7.999997E-002 4.678031E-009 6.715386E-023 3.226145E-027 0.000000E+000 9.998162E-001 2.529532E-005 2.998809E-008 7.413397E-022 0.000000E+000 0.000000E+000 6.030835E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.631876E-022
1.757110E+003 1.132664E+005 4.789003E+003 4.789001E+003 1.124168E+003 1.040872E+003 9.994419E-001 5.834377E-004 8.000000E-002 1.487696E-008 8.482586E-022 4.087718E-026 0.000000E+000 9.994166E-001 2.528512E-005 3.614538E-009 2.587754E-021 0.000000E+000 0.000000E+000 6.316241E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.165232E-021
1.690769E+003 1.245930E+005 4.608192E+003 4.608190E+003 1.001585E+003 9.273713E+002 9.980078E-001 2.017416E-003 8.000000E-002 5.149792E-008 1.300238E-020 6.286518E-025 0.000000E+000 9.979826E-001 2.524850E-005 3.720185E-010 9.960912E-021 0.000000E+000 0.000000E+000 6.663902E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.080555E-021
1.624984E+003 1.370523E+005 4.428893E+003 4.428891E+003 8.891632E+002 8.232800E+002 9.924409E-001 7.584178E-003 8.000000E-002 1.937150E-007 2.413797E-019 1.171141E-023 0.000000E+000 9.924158E-001 2.510628E-005 3.254678E-011 4.211770E-020 0.000000E+000 0.000000E+000 7.066382E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.552399E-020
1.560044E+003 1.507576E+005 4.251900E+003 4.251900E+003 7.867648E+002 7.284689E+002 9.697415E-001 3.028306E-002 8.000000E-002 7.736873E-007 5.189459E-018 2.526862E-022 0.000000E+000 9.697169E-001 2.452631E-005 2.438196E-012 1.887841E-019 0.000000E+000 0.000000E+000 7.437139E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.186790E-020
1.495959E+003 1.658333E+005 4.077237E+003 4.077246E+003 6.937353E+002 6.423324E+002 8.833419E-001 1.166804E-001 8.000000E-002 2.979604E-006 1.070255E-016 5.227758E-021 0.000000E+000 8.833196E-001 2.232040E-005 1.621813E-013 7.745461E-019 0.000000E+000 0.000000E+000 7.379167E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.227260E-019
1.432889E+003 1.824166E+005 3.905338E+003 3.905377E+003 6.096377E+002 5.644661E+002 6.610883E-001 3.389284E-001 8.000000E-002 8.637399E-006 1.296519E-015 6.344764E-020 0.000000E+000 6.610716E-001 1.666260E-005 1.079821E-014 2.006082E-018 0.000000E+000 0.000000E+000 6.089954E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.989448E-019
1.371360E+003 2.006583E+005 3.737641E+003 3.737762E+003 5.344273E+002 4.948285E+002 3.744022E-001 6.256072E-001 8.000000E-002 1.589593E-005 6.546740E-015 3.207561E-019 0.000000E+000 3.743928E-001 9.404068E-006 8.024331E-016 2.547514E-018 0.000000E+000 0.000000E+000 3.852133E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.401022E-019
1.311036E+003 2.207241E+005 3.573229E+003 3.573582E+003 4.669588E+002 4.323592E+002 1.696684E-001 8.303359E-001 8.000000E-002 2.105021E-005 1.764202E-014 8.662553E-019 0.000000E+000 1.696641E-001 4.249786E-006 6.410424E-017 1.903765E-018 0.000000E+000 0.000000E+000 1.978523E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.382740E-019
1.252051E+003 2.427966E+005 3.412465E+003 3.413494E+003 4.067246E+002 3.765880E+002 6.731842E-002 9.326833E-001 8.000000E-002 2.361679E-005 3.462626E-014 1.706246E-018 0.000000E+000 6.731674E-002 1.683212E-006 5.880040E-018 1.080424E-018 0.000000E+000 0.000000E+000 9.040086E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.579608E-019
1.194477E+003 2.670762E+005 3.255548E+003 3.258428E+003 3.531574E+002 3.269899E+002 2.531832E-002 9.746823E-001 8.000000E-002 2.466765E-005 5.742875E-014 2.842619E-018 0.000000E+000 2.531768E-002 6.323499E-007 8.766468E-019 5.556832E-019 0.000000E+000 0.000000E+000 3.983841E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.550635E-020
1.138933E+003 2.937838E+005 3.104163E+003 3.110772E+003 3.061466E+002 2.834625E+002 1.010285E-002 9.898974E-001 8.000000E-002 2.504788E-005 8.332902E-014 4.145769E-018 0.000000E+000 1.010260E-002 2.521219E-007 3.036711E-019 3.027660E-019 0.000000E+000 0.000000E+000 1.894570E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.821402E-020
1.085222E+003 3.231622E+005 2.957773E+003 2.968475E+003 2.648440E+002 2.452202E+002 4.799868E-003 9.952003E-001 8.000000E-002 2.518030E-005 1.092894E-013 5.466184E-018 0.000000E+000 4.799748E-003 1.196955E-007 1.873528E-019 2.001991E-019 0.000000E+000 0.000000E+000 1.093417E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.443178E-020
1.032744E+003 3.554784E+005 2.814743E+003 2.827338E+003 2.282506E+002 2.113382E+002 2.819167E-003 9.971809E-001 8.000000E-002 2.522975E-005 1.355824E-013 6.821360E-018 0.000000E+000 2.819096E-003 7.024961E-008 1.459420E-019 1.687948E-019 0.000000E+000 0.000000E+000 7.988044E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 9.172137E-021
9.821738E+002 3.910263E+005 2.676915E+003 2.689644E+003 1.963358E+002 1.817882E+002 1.957161E-003 9.980429E-001 8.000000E-002 2.525127E-005 1.638260E-013 8.294937E-018 0.000000E+000 1.957112E-003 4.873113E-008 1.250446E-019 1.736497E-019 0.000000E+000 0.000000E+000 7.057700E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.939794E-021
9.339302E+002 4.301289E+005 2.545427E+003 2.557763E+003 1.688021E+002 1.562946E+002 1.504991E-003 9.984950E-001 8.000000E-002 2.526256E-005 1.960920E-013 9.995393E-018 0.000000E+000 1.504954E-003 3.744152E-008 1.121874E-019 2.042375E-019 0.000000E+000 0.000000E+000 7.060247E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.862608E-021
8.870018E+002 4.731418E+005 2.417523E+003 2.429306E+003 1.446132E+002 1.338980E+002 1.225188E-003 9.987748E-001 8.000000E-002 2.526955E-005 2.361168E-013 1.212204E-017 0.000000E+000 1.225157E-003 3.045331E-008 1.031439E-019 2.656682E-019 0.000000E+000 0.000000E+000 7.699139E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.301280E-021
8.416509E+002 5.204560E+005 2.293920E+003 2.305038E+003 1.235465E+002 1.143922E+002 1.036264E-003 9.989638E-001 8.000000E-002 2.527427E-005 2.886628E-013 1.493235E-017 0.000000E+000 1.036238E-003 2.573277E-008 9.638274E-020 3.757184E-019 0.000000E+000 0.000000E+000 8.992751E-022 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.038734E-021
7.986964E+002 5.725016E+005 2.176847E+003 2.187198E+003 1.055795E+002 9.775650E+001 9.025311E-004 9.990975E-001 8.000000E-002 2.527761E-005 3.605374E-013 1.879652E-017 0.000000E+000 9.025087E-004 2.238935E-008 9.118621E-020 5.685059E-019 0.000000E+000 0.000000E+000 1.110228E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.981913E-021
7.572431E+002 6.297517E+005 2.063866E+003 2.073321E+003 8.997886E+001 8.331181E+001 8.012929E-004 9.991987E-001 8.000000E-002 2.528014E-005 4.672246E-013 2.455759E-017 0.000000E+000 8.012730E-004 1.985642E-008 8.697504E-020 9.237185E-019 0.000000E+000 0.000000E+000 1.446941E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.089102E-021
7.176872E+002 6.927269E+005 1.956057E+003 1.964482E+003 7.660200E+001 7.092612E+001 7.228051E-004 9.992772E-001 8.000000E-002 2.528211E-005 6.345095E-013 3.362845E-017 0.000000E+000 7.227872E-004 1.789088E-008 8.351591E-020 1.602406E-018 0.000000E+000 0.000000E+000 1.981419E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.350934E-021
6.795804E+002 7.619996E+005 1.852196E+003 1.859434E+003 6.503648E+001 6.021756E+001 6.596547E-004 9.993404E-001 8.000000E-002 2.528369E-005 9.181211E-013 4.907423E-017 0.000000E+000 6.596384E-004 1.630759E-008 8.059096E-020 2.987794E-018 0.000000E+000 0.000000E+000 2.859385E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.782234E-021
6.431924E+002 8.381995E+005 1.753021E+003 1.758905E+003 5.513879E+001 5.105324E+001 6.081762E-004 9.993918E-001 8.000000E-002 2.528498E-005 1.429905E-012 7.708265E-017 0.000000E+000 6.081612E-004 1.501511E-008 7.810134E-020 5.973003E-018 0.000000E+000 0.000000E+000 4.340157E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.411292E-021
6.089863E+002 9.220195E+005 1.659792E+003 1.664169E+003 4.680122E+001 4.333344E+001 5.660495E-004 9.994340E-001 8.000000E-002 2.528604E-005 2.400340E-012 1.304801E-016 0.000000E+000 5.660355E-004 1.395561E-008 7.598628E-020 1.267937E-017 0.000000E+000 0.000000E+000 6.885307E-021 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.274095E-021
5.757396E+002 1.014221E+006 1.569178E+003 1.571825E+003 3.954693E+001 3.661667E+001 5.298575E-004 9.994702E-001 8.000000E-002 2.528696E-005 4.453827E-012 2.441502E-016 0.000000E+000 5.298444E-004 1.304350E-008 7.410845E-020 2.943669E-017 0.000000E+000 0.000000E+000 1.162341E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.481751E-021
5.445589E+002 1.115644E+006 1.484195E+003 1.484935E+003 3.346332E+001 3.098383E+001 4.994679E-004 9.995005E-001 8.000000E-002 2.528772E-005 8.979111E-012 4.962371E-016 0.000000E+000 4.994556E-004 1.227574E-008 7.248523E-020 7.288551E-017 0.000000E+000 0.000000E+000 2.055583E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.011345E-020
5.145064E+002 1.227208E+006 1.402287E+003 1.400879E+003 2.822324E+001 2.613202E+001 4.729392E-004 9.995271E-001 8.000000E-002 2.528840E-005 2.009058E-011 1.119241E-015 0.000000E+000 4.729276E-004 1.160355E-008 7.103154E-020 1.981782E-016 0.000000E+000 0.000000E+000 3.878911E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.239724E-020
4.860155E+002 1.349929E+006 1.324635E+003 1.320855E+003 2.378947E+001 2.202678E+001 4.499443E-004 9.995501E-001 8.000000E-002 2.528898E-005 4.941042E-011 2.773665E-015 0.000000E+000 4.499333E-004 1.101886E-008 6.974258E-020 5.859542E-016 0.000000E+000 0.000000E+000 7.762859E-020 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.559031E-020
4.591644E+002 1.484922E+006 1.251453E+003 1.245085E+003 2.006038E+001 1.857399E+001 4.299473E-004 9.995701E-001 8.000000E-002 2.528949E-005 1.329528E-010 7.515672E-015 0.000000E+000 4.299368E-004 1.050831E-008 6.859909E-020 1.878781E-015 0.000000E+000 0.000000E+000 1.645004E-019 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.009704E-020
4.334508E+002 1.633414E+006 1.181370E+003 1.172150E+003 1.687539E+001 1.562499E+001 4.121395E-004 9.995879E-001 8.000000E-002 2.528995E-005 3.989955E-010 2.269635E-014 0.000000E+000 4.121295E-004 1.005145E-008 6.756275E-020 6.690440E-015 0.000000E+000 0.000000E+000 3.750779E-019 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.669675E-020
4.090984E+002 1.796755E+006 1.114998E+003 1.102685E+003 1.418789E+001 1.313662E+001 3.963512E-004 9.996037E-001 8.000000E-002 2.529036E-005 1.324580E-009 7.574860E-014 0.000000E+000 3.963416E-004 9.644087E-009 6.662959E-020 2.617282E-014 0.000000E+000 0.000000E+000 9.164692E-019 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.649314E-020
3.857217E+002 1.976431E+006 1.051285E+003 1.035592E+003 1.189205E+001 1.101090E+001 3.820777E-004 9.996179E-001 8.000000E-002 2.529073E-005 4.927859E-009 2.832806E-013 0.000000E+000 3.820685E-004 9.273328E-009 6.577438E-020 1.127801E-013 0.000000E+000 0.000000E+000 2.443255E-018 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.166126E-020
3.639938E+002 2.174074E+006 9.920651E+002 9.728189E+002 9.993472E+000 9.252998E+000 3.695091E-004 9.996305E-001 8.000000E-002 2.529105E-005 1.866458E-008 1.084869E-012 0.000000E+000 3.695002E-004 8.944286E-009 6.501238E-020 4.530293E-013 0.000000E+000 0.000000E+000 6.908488E-018 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.500241E-020
3.432870E+002 2.391481E+006 9.356287E+002 9.125825E+002 8.383134E+000 7.761980E+000 3.580978E-004 9.996418E-001 8.000000E-002 2.529135E-005 5.989336E-008 3.608318E-012 0.000000E+000 3.580892E-004 8.642818E-009 6.431350E-020 1.165447E-012 0.000000E+000 0.000000E+000 2.126364E-017 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.128468E-019
3.236451E+002 2.630629E+006 8.820948E+002 8.550316E+002 7.024925E+000 6.504408E+000 3.477379E-004 9.996520E-001 8.000000E-002 2.529163E-005 1.255368E-007 7.993816E-012 0.000000E+000 3.477295E-004 8.366207E-009 6.367348E-020 1.472356E-012 0.000000E+000 0.000000E+000 7.042222E-017 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.760637E-019
3.050449E+002 2.893692E+006 8.313998E+002 8.001292E+002 5.882007E+000 5.446176E+000 3.383067E-004 9.996613E-001 8.000000E-002 2.529188E-005 1.886700E-007 1.260394E-011 0.000000E+000 3.382986E-004 8.111268E-009 6.308658E-020 1.362583E-012 0.000000E+000 0.000000E+000 2.399975E-016 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.853146E-019
2.874388E+002 3.183061E+006 7.834146E+002 7.477748E+002 4.921196E+000 4.556558E+000 3.296891E-004 9.996698E-001 8.000000E-002 2.529211E-005 2.405386E-007 1.674265E-011 0.000000E+000 3.296812E-004 7.874972E-009 6.254708E-020 1.190105E-012 0.000000E+000 0.000000E+000 7.231293E-016 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.813245E-019
2.707080E+002 3.501368E+006 7.378146E+002 6.976579E+002 4.110906E+000 3.806306E+000 3.217534E-004 9.996777E-001 8.000000E-002 2.529233E-005 2.822013E-007 2.043435E-011 0.000000E+000 3.217458E-004 7.653745E-009 6.204788E-020 1.030719E-012 0.000000E+000 0.000000E+000 1.472924E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.496359E-019
2.549640E+002 3.851504E+006 6.949044E+002 6.501654E+002 3.434558E+000 3.180073E+000 3.144902E-004 9.996849E-001 8.000000E-002 2.529253E-005 3.152023E-007 2.375973E-011 0.000000E+000 3.144828E-004 7.447373E-009 6.158931E-020 8.914602E-013 0.000000E+000 0.000000E+000 1.948183E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.565098E-018
2.400542E+002 4.236655E+006 6.542678E+002 6.048959E+002 2.866569E+000 2.654169E+000 3.077765E-004 9.996915E-001 8.000000E-002 2.529272E-005 3.414619E-007 2.685360E-011 0.000000E+000 3.077692E-004 7.252407E-009 6.116435E-020 7.694722E-013 0.000000E+000 0.000000E+000 2.130842E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.028962E-018
2.260149E+002 4.660320E+006 6.160037E+002 5.620211E+002 2.392465E+000 2.215194E+000 3.015853E-004 9.996977E-001 8.000000E-002 2.529290E-005 3.622012E-007 2.979956E-011 0.000000E+000 3.015782E-004 7.068093E-009 6.077187E-020 6.634237E-013 0.000000E+000 0.000000E+000 2.224080E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.155043E-018
2.128577E+002 5.126352E+006 5.801436E+002 5.216440E+002 1.998491E+000 1.850412E+000 2.958844E-004 9.997034E-001 8.000000E-002 2.529307E-005 3.784893E-007 3.267646E-011 0.000000E+000 2.958775E-004 6.893548E-009 6.041027E-020 5.717113E-013 0.000000E+000 0.000000E+000 2.304485E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.312013E-017
2.003365E+002 5.638987E+006 5.460171E+002 4.830739E+002 1.666152E+000 1.542697E+000 2.905378E-004 9.997087E-001 8.000000E-002 2.529324E-005 3.914559E-007 3.561552E-011 0.000000E+000 2.905311E-004 6.724588E-009 6.007128E-020 4.912252E-013 0.000000E+000 0.000000E+000 2.390321E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.979033E-017
1.884700E+002 6.202886E+006 5.136751E+002 4.464301E+002 1.387271E+000 1.284481E+000 2.855303E-004 9.997137E-001 8.000000E-002 2.529340E-005 4.017047E-007 3.869226E-011 0.000000E+000 2.855237E-004 6.560588E-009 5.975419E-020 4.209836E-013 0.000000E+000 0.000000E+000 2.484928E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.219384E-017
1.773194E+002 6.823175E+006 4.832842E+002 4.119633E+002 1.155323E+000 1.069719E+000 2.808671E-004 9.997183E-001 8.000000E-002 2.529356E-005 4.097224E-007 4.196705E-011 0.000000E+000 2.808607E-004 6.401657E-009 5.945953E-020 3.602798E-013 0.000000E+000 0.000000E+000 2.588695E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.861310E-016
1.668922E+002 7.505492E+006 4.548648E+002 3.797551E+002 9.632579E-001 8.918849E-001 2.765339E-004 9.997226E-001 8.000000E-002 2.529371E-005 4.159676E-007 4.550673E-011 0.000000E+000 2.765276E-004 6.247353E-009 5.918651E-020 3.081267E-013 0.000000E+000 0.000000E+000 2.701793E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.097017E-016
1.569445E+002 8.256041E+006 4.277524E+002 3.491054E+002 8.010747E-001 7.417188E-001 2.724155E-004 9.997267E-001 8.000000E-002 2.529386E-005 4.209474E-007 4.947789E-011 0.000000E+000 2.724094E-004 6.093447E-009 5.892799E-020 2.625034E-013 0.000000E+000 0.000000E+000 2.827500E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.520678E-015
1.475456E+002 9.081646E+006 4.021356E+002 3.202765E+002 6.656001E-001 6.162822E-001 2.685295E-004 9.997306E-001 8.000000E-002 2.529401E-005 4.249571E-007 5.395526E-011 0.000000E+000 2.685235E-004 5.940310E-009 5.868515E-020 2.230842E-013 0.000000E+000 0.000000E+000 2.965979E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.926504E-015
1.388552E+002 9.989810E+006 3.784499E+002 2.937939E+002 5.547801E-001 5.136734E-001 2.649327E-004 9.997342E-001 8.000000E-002 2.529415E-005 4.284069E-007 5.894256E-011 0.000000E+000 2.649269E-004 5.790258E-009 5.846153E-020 1.898079E-013 0.000000E+000 0.000000E+000 3.115063E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.694196E-014
1.305486E+002 1.098879E+007 3.558103E+002 2.686983E+002 4.610535E-001 4.268915E-001 2.614835E-004 9.997377E-001 8.000000E-002 2.529430E-005 4.325177E-007 6.475503E-011 0.000000E+000 2.614778E-004 5.637300E-009 5.824833E-020 1.608325E-013 0.000000E+000 0.000000E+000 3.280948E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.466064E-014
1.227552E+002 1.208767E+007 3.345694E+002 2.454073E+002 3.833136E-001 3.549118E-001 2.582296E-004 9.997409E-001 8.000000E-002 2.529445E-005 4.410285E-007 7.159025E-011 0.000000E+000 2.582241E-004 5.483251E-009 5.804850E-020 1.361446E-013 0.000000E+000 0.000000E+000 3.462336E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.668644E-013
1.153897E+002 1.329644E+007 3.144947E+002 2.236816E+002 3.183725E-001 2.947826E-001 2.551309E-004 9.997439E-001 8.000000E-002 2.529459E-005 4.694552E-007 8.034811E-011 0.000000E+000 2.551256E-004 5.326054E-009 5.785958E-020 1.150133E-013 0.000000E+000 0.000000E+000 3.662241E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.154340E-012
1.084853E+002 1.462608E+007 2.956768E+002 2.036265E+002 2.645741E-001 2.449704E-001 2.521974E-004 9.997467E-001 8.000000E-002 2.529474E-005 5.714643E-007 9.380365E-011 0.000000E+000 2.521923E-004 5.166091E-009 5.768216E-020 9.712431E-014 0.000000E+000 0.000000E+000 3.880899E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.182390E-012
1.019356E+002 1.608869E+007 2.778255E+002 1.849347E+002 2.194886E-001 2.032256E-001 2.493805E-004 9.997490E-001 8.000000E-002 2.529488E-005 8.216663E-007 1.187613E-010 0.000000E+000 2.493755E-004 5.000532E-009 5.751333E-020 8.184546E-014 0.000000E+000 0.000000E+000 4.123221E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.160158E-012
9.581458E+001 1.769756E+007 2.611427E+002 1.678096E+002 1.822757E-001 1.687701E-001 2.467118E-004 9.997510E-001 8.000000E-002 2.529501E-005 1.142950E-006 1.590901E-010 0.000000E+000 2.467070E-004 4.831034E-009 5.735486E-020 6.902839E-014 0.000000E+000 0.000000E+000 4.388152E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.910460E-012
9.006415E+001 1.946731E+007 2.454699E+002 1.520695E+002 1.513874E-001 1.401704E-001 2.441674E-004 9.997530E-001 8.000000E-002 2.529513E-005 1.436475E-006 2.115081E-010 0.000000E+000 2.441628E-004 4.656084E-009 5.720515E-020 5.824676E-014 0.000000E+000 0.000000E+000 4.679584E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 8.190582E-012
8.463566E+001 2.141405E+007 2.306745E+002 1.375633E+002 1.256302E-001 1.163217E-001 2.417261E-004 9.997549E-001 8.000000E-002 2.529525E-005 1.684471E-006 2.715331E-010 0.000000E+000 2.417217E-004 4.474129E-009 5.706287E-020 4.916109E-014 0.000000E+000 0.000000E+000 5.002302E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 7.316823E-012
7.951202E+001 2.355545E+007 2.167100E+002 1.242255E+002 1.041675E-001 9.644927E-002 2.393804E-004 9.997568E-001 8.000000E-002 2.529538E-005 1.891128E-006 3.356904E-010 0.000000E+000 2.393761E-004 4.284496E-009 5.692753E-020 4.153181E-014 0.000000E+000 0.000000E+000 5.360465E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 6.509693E-012
7.467543E+001 2.591099E+007 2.035279E+002 1.119863E+002 8.629126E-002 7.989760E-002 2.371226E-004 9.997588E-001 8.000000E-002 2.529551E-005 2.062866E-006 4.013810E-010 0.000000E+000 2.371185E-004 4.086547E-009 5.679862E-020 3.514648E-014 0.000000E+000 0.000000E+000 5.759118E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.787095E-012
7.017007E+001 2.850209E+007 1.912485E+002 1.009243E+002 7.159611E-002 6.629127E-002 2.349752E-004 9.997606E-001 8.000000E-002 2.529565E-005 2.203600E-006 4.655855E-010 0.000000E+000 2.349713E-004 3.882664E-009 5.667735E-020 2.988679E-014 0.000000E+000 0.000000E+000 6.197841E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 5.151647E-012
6.592292E+001 3.135230E+007 1.796729E+002 9.082379E+001 5.936669E-002 5.496798E-002 2.329057E-004 9.997625E-001 8.000000E-002 2.529580E-005 2.320266E-006 5.273060E-010 0.000000E+000 2.329020E-004 3.670394E-009 5.656180E-020 2.551285E-014 0.000000E+000 0.000000E+000 6.687499E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.587058E-012
6.194158E+001 3.448753E+007 1.688218E+002 8.166895E+001 4.924706E-002 4.559816E-002 2.309199E-004 9.997643E-001 8.000000E-002 2.529596E-005 2.416364E-006 5.848059E-010 0.000000E+000 2.309164E-004 3.451026E-009 5.645221E-020 2.190560E-014 0.000000E+000 0.000000E+000 7.232558E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 4.089009E-012
5.810253E+001 3.793629E+007 1.583585E+002 7.315137E+001 4.064608E-002 3.763446E-002 2.289569E-004 9.997661E-001 8.000000E-002 2.529614E-005 2.497607E-006 6.385807E-010 0.000000E+000 2.289537E-004 3.218247E-009 5.634519E-020 1.885758E-014 0.000000E+000 0.000000E+000 7.860068E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.638098E-012
5.463936E+001 4.172992E+007 1.489196E+002 6.574776E+001 3.380262E-002 3.129806E-002 2.271399E-004 9.997677E-001 8.000000E-002 2.529633E-005 2.561867E-006 6.846069E-010 0.000000E+000 2.271369E-004 2.988244E-009 5.624736E-020 1.645156E-014 0.000000E+000 0.000000E+000 8.537545E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 3.256642E-012
5.130393E+001 4.590291E+007 1.398289E+002 5.888496E+001 2.798242E-002 2.590910E-002 2.253428E-004 9.997694E-001 8.000000E-002 2.529653E-005 2.616166E-006 7.258462E-010 0.000000E+000 2.253400E-004 2.747058E-009 5.615181E-020 1.442451E-014 0.000000E+000 0.000000E+000 9.320518E-015 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.912471E-012
4.821007E+001 5.049320E+007 1.313965E+002 5.276826E+001 2.321917E-002 2.149878E-002 2.236289E-004 9.997711E-001 8.000000E-002 2.529673E-005 2.660290E-006 7.607670E-010 0.000000E+000 2.236263E-004 2.504933E-009 5.606187E-020 1.278221E-014 0.000000E+000 0.000000E+000 1.019681E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.614047E-012
4.522805E+001 5.554252E+007 1.232691E+002 4.711226E+001 1.917154E-002 1.775105E-002 2.219285E-004 9.997727E-001 8.000000E-002 2.529695E-005 2.697497E-006 7.909465E-010 0.000000E+000 2.219263E-004 2.254244E-009 5.597383E-020 1.140009E-014 0.000000E+000 0.000000E+000 1.122239E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.345813E-012
4.251272E+001 6.109677E+007 1.158684E+002 4.217719E+001 1.592172E-002 1.474202E-002 2.203335E-004 9.997742E-001 8.000000E-002 2.529717E-005 2.727125E-006 8.151859E-010 0.000000E+000 2.203315E-004 2.011385E-009 5.589237E-020 1.029823E-014 0.000000E+000 0.000000E+000 1.236195E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 2.118536E-012
3.986658E+001 6.720645E+007 1.086564E+002 3.757461E+001 1.312986E-002 1.215703E-002 2.187309E-004 9.997758E-001 8.000000E-002 2.529740E-005 2.752354E-006 8.356743E-010 0.000000E+000 2.187291E-004 1.762374E-009 5.581165E-020 9.353267E-015 0.000000E+000 0.000000E+000 1.372642E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.913042E-012
3.747910E+001 7.392709E+007 1.021493E+002 3.360472E+001 1.090940E-002 1.010109E-002 2.172388E-004 9.997772E-001 8.000000E-002 2.529762E-005 2.772234E-006 8.514375E-010 0.000000E+000 2.172373E-004 1.529050E-009 5.573755E-020 8.596120E-015 0.000000E+000 0.000000E+000 1.524847E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.741611E-012
3.517334E+001 8.131980E+007 9.586494E+001 2.994199E+001 9.017251E-003 8.349132E-003 2.157511E-004 9.997787E-001 8.000000E-002 2.529784E-005 2.789008E-006 8.642026E-010 0.000000E+000 2.157498E-004 1.298583E-009 5.566472E-020 7.937007E-015 0.000000E+000 0.000000E+000 1.707692E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.589122E-012
3.302470E+001 8.945178E+007 9.000884E+001 2.668604E+001 7.463631E-003 6.910626E-003 2.143184E-004 9.997801E-001 8.000000E-002 2.529804E-005 2.802639E-006 8.739627E-010 0.000000E+000 2.143173E-004 1.082663E-009 5.559560E-020 7.373232E-015 0.000000E+000 0.000000E+000 1.921788E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.459148E-012
3.100436E+001 9.839696E+007 8.450238E+001 2.376766E+001 6.175918E-003 5.718324E-003 2.129253E-004 9.997814E-001 8.000000E-002 2.529824E-005 2.813813E-006 8.813342E-010 0.000000E+000 2.129244E-004 8.825753E-010 5.552938E-020 6.873859E-015 0.000000E+000 0.000000E+000 2.176804E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.348265E-012
2.913074E+001 1.082367E+008 7.939585E+001 2.118924E+001 5.122572E-003 4.743024E-003 2.115888E-004 9.997828E-001 8.000000E-002 2.529841E-005 2.822848E-006 8.867003E-010 0.000000E+000 2.115881E-004 7.037638E-010 5.546679E-020 6.424560E-015 0.000000E+000 0.000000E+000 2.478812E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.256027E-012
2.734339E+001 1.190603E+008 7.452442E+001 1.884777E+001 4.236338E-003 3.922453E-003 2.102690E-004 9.997841E-001 8.000000E-002 2.529857E-005 2.830358E-006 8.906108E-010 0.000000E+000 2.102685E-004 5.435595E-010 5.540591E-020 5.995520E-015 0.000000E+000 0.000000E+000 2.850692E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.178497E-012
2.568486E+001 1.309664E+008 7.000409E+001 1.678121E+001 3.511277E-003 3.251114E-003 2.090005E-004 9.997853E-001 8.000000E-002 2.529870E-005 2.836422E-006 8.932981E-010 0.000000E+000 2.090001E-004 4.079754E-010 5.534829E-020 5.584660E-015 0.000000E+000 0.000000E+000 3.301291E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.116866E-012
2.411109E+001 1.440630E+008 6.571478E+001 1.491720E+001 2.904584E-003 2.689374E-003 2.077532E-004 9.997866E-001 8.000000E-002 2.529881E-005 2.841425E-006 8.951240E-010 0.000000E+000 2.077529E-004 2.944540E-010 5.529250E-020 5.173227E-015 0.000000E+000 0.000000E+000 3.866340E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.069046E-012
2.262709E+001 1.584693E+008 6.167012E+001 1.324807E+001 2.400597E-003 2.222729E-003 2.065333E-004 9.997878E-001 8.000000E-002 2.529890E-005 2.845519E-006 8.963068E-010 0.000000E+000 2.065331E-004 2.036036E-010 5.523878E-020 4.754201E-015 0.000000E+000 0.000000E+000 4.581304E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.035351E-012
2.122592E+001 1.743162E+008 5.785125E+001 1.175282E+001 1.981678E-003 1.834848E-003 2.053381E-004 9.997890E-001 8.000000E-002 2.529897E-005 2.848865E-006 8.970395E-010 0.000000E+000 2.053379E-004 1.339654E-010 5.518699E-020 4.322880E-015 0.000000E+000 0.000000E+000 5.501588E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.016136E-012
1.995168E+001 1.917478E+008 5.437831E+001 1.046252E+001 1.645780E-003 1.523839E-003 2.042100E-004 9.997901E-001 8.000000E-002 2.529902E-005 2.851499E-006 8.974576E-010 0.000000E+000 2.042099E-004 8.486382E-011 5.513888E-020 3.894012E-015 0.000000E+000 0.000000E+000 6.654599E-014 0.000000E+000 0.000000E+000 0.000000E+000 0.000000E+000 1.012244E-012
examples/Cooling/CoolingWithPrimordialElements/getICs.sh 0 → 100755
+ 2
− 0
View file @ 454c2954
#!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/homogeneousCube_32.hdf5
examples/Cooling/CoolingWithPrimordialElements/getPrimordialAbundances.py 0 → 100755
+ 62
− 0
View file @ 454c2954
#!/usr/bin/env python3
"""
This script is an example of how to compute the
fractional abundances of different elements with
respect to Hydrogen.
For comparison, see
for example Dalgarno 2005
"""
# some parameters (Planck 2018)
Omega0 = 0.315
Omega_b = 0.049
h0 = 0.674
T0 = 2.73
fH = 0.76
z0 = 300
z = 77
# Temperature
T = T0 / (1 + 200) * (1 + z) ** 2 # Anninos & Norman 1996 eq.11
# all abundances per unit of nH
nHII = 1.2e-5 * Omega0 ** 0.5 / (h0 * Omega_b) # Anninos & Norman 1996 eq. 9
nHeI = (1 - fH) / fH / 4
nHeII = 0 # assume neutral He
nHeIII = 0 # assume neutral He
nHM = 2e-11 * (1 + z) ** (1.76) * nHII # Anninos & Norman 1996 eq. 12
nH2I = (
2e-20 * fH * Omega0 ** (3 / 2.0) / (h0 * Omega_b) * (1 + z0) ** (5.1)
) # Anninos & Norman 1996 eq. 15
nH2II = 1e-17 * (1 + z) ** (3.6) * nHII # Anninos & Norman 1996 eq. 13
nDI = 4e-5 # Bromm 20002 from Galli & Palla 1998
nDII = 1e-8 # Bromm 20002
nHDI = 1e-3 * nH2I # Bromm 20002
print("T : %g K" % T)
print("initial_nHII_to_nH_ratio : %g" % nHII)
print("initial_nHeI_to_nH_ratio : %g" % nHeI)
print("initial_nHeII_to_nH_ratio : %g" % nHeII)
print("initial_nHeIII_to_nH_ratio : %g" % nHeIII)
print("initial_nDI_to_nH_ratio : %g" % nDI)
print("initial_nDII_to_nH_ratio : %g" % nDII)
print("initial_nHM_to_nH_ratio : %g" % nHM)
print("initial_nH2I_to_nH_ratio : %g" % nH2I)
print("initial_nH2II_to_nH_ratio : %g" % nH2II)
print("initial_nHDI_to_nH_ratio : %g" % nHDI)
examples/Cooling/CoolingWithPrimordialElements/local.mplstyle 0 → 100644
+ 28
− 0
View file @ 454c2954
# Line Properties
lines.linewidth: 2
# Font options
font.size: 8
font.family: STIXGeneral
mathtext.fontset: stix
# LaTeX options
text.usetex: False
# Legend options
#legend.frameon: False
#legend.fontsize: 8
# Figure options for publications
figure.dpi: 300
figure.figsize: 3.321, 3.0 # Correct width for MNRAS
# Histogram options
hist.bins: auto
# Ticks inside plots; more space devoted to plot.
xtick.direction: in
ytick.direction: in
# Always plot ticks on top of data
axes.axisbelow: False
examples/Cooling/CoolingWithPrimordialElements/output_list.txt 0 → 100644
+ 28
− 0
View file @ 454c2954
# Redshift
77
70
60
50
40
30
20
19
18
17
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
examples/Cooling/CoolingWithPrimordialElements/plotAbundances.py 0 → 100755
+ 185
− 0
View file @ 454c2954
#!/usr/bin/env python3
###############################################################################
# This file is part of SWIFT.
# Copyright (c) 2023 Yves Revaz (yves.revaz@epfl.ch)
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
##############################################################################
import matplotlib.pyplot as plt
plt.style.use("local.mplstyle")
import numpy as np
from glob import glob
import h5py
import sys
filenames = sys.argv[1:]
n = len(filenames)
rs = np.zeros(n)
nHIs = np.zeros(n)
nHIIs = np.zeros(n)
nHeIs = np.zeros(n)
nHeIIs = np.zeros(n)
nHeIIIs = np.zeros(n)
nHMs = np.zeros(n)
nH2Is = np.zeros(n)
nH2IIs = np.zeros(n)
nDIs = np.zeros(n)
nDIIs = np.zeros(n)
nHDIs = np.zeros(n)
for i, filename in enumerate(filenames):
with h5py.File(filename, "r") as f:
XHI = f["PartType0/HI"][:]
XHII = f["PartType0/HII"][:]
XHeI = f["PartType0/HeI"][:]
XHeII = f["PartType0/HeII"][:]
XHeIII = f["PartType0/HeIII"][:]
XHM = f["PartType0/HM"][:]
XH2I = f["PartType0/H2I"][:]
XH2II = f["PartType0/H2II"][:]
XDI = f["PartType0/DI"][:]
XDII = f["PartType0/DII"][:]
XHDI = f["PartType0/HDI"][:]
redshift = f["Header"].attrs["Redshift"][0]
rs[i] = redshift
# compute hydrogen mass fraction
XH = (
XHI.mean()
+ XHII.mean()
+ XHM.mean()
+ XH2I.mean()
+ XH2II.mean()
+ XDI.mean()
+ XDII.mean()
+ XHDI.mean()
)
# make it equal to the hydrogen number density (do not consider density on purpose)
nH = XH
# all abundances per unit of nH
nHI = XHI / nH
nHII = XHII / nH
nHeI = XHeI / nH / 4
nHeII = XHeII / nH / 4
nHeIII = XHeIII / nH / 4
nHM = XHM / nH
nH2I = XH2I / nH / 2
nH2II = XH2II / nH / 2
nDI = XDI / nH / 2
nDII = XDII / nH / 2
nHDI = XHDI / nH / 3
nHIs[i] = nHI.mean()
nHIIs[i] = nHII.mean()
nHeIs[i] = nHeI.mean()
nHeIIs[i] = nHeII.mean()
nHeIIIs[i] = nHeIII.mean()
nHMs[i] = nHM.mean()
nH2Is[i] = nH2I.mean()
nH2IIs[i] = nH2II.mean()
nDIs[i] = nDI.mean()
nDIIs[i] = nDII.mean()
nHDIs[i] = nHDI.mean()
# data from Faure, A., Hily-Blant, Pineau des Forêts, Flower 2024
# The last lines have been removed on purpose. Their predictions
# are not reliable below z=20.
# Note: the "Redshift" field = z+1
with open("faure2024.dat", "r") as f:
f.readline()
f.readline()
header = str.split(f.readline())
lines = f.readlines()
data = np.array(list(map(str.split, lines)), float)
dic = {x: i for i, x in enumerate(header)}
##########################################
# plot
##########################################
plt.figure()
plt.subplots_adjust(left=0.17, bottom=0.12, right=0.97, top=0.97)
ax = plt.gca()
rs = rs + 1
pHI = ax.plot(rs, nHIs, ms=2, lw=1, label="HI")
pHII = ax.plot(rs, nHIIs, ms=2, lw=1, label="HII")
pHeI = ax.plot(rs, nHeIs, ms=2, lw=1, label="HeI")
pHeII = ax.plot(rs, nHeIIs, ms=2, lw=1, label="HeII")
pHeIII = ax.plot(rs, nHeIIIs, ms=2, lw=1, label="HeIII")
pHM = ax.plot(rs, nHMs, ms=2, lw=1, label="HM")
pH2I = ax.plot(rs, nH2Is, ms=2, lw=1, label="H2I")
pH2II = ax.plot(rs, nH2IIs, ms=2, lw=1, label="H2II")
pDI = ax.plot(rs, nDIs, ms=2, lw=1, label="DI")
pDII = ax.plot(rs, nDIIs, ms=2, lw=1, label="DII")
pHDI = ax.plot(rs, nHDIs, ms=2, lw=1, label="HDI")
# add data
ax.plot(data[:, dic["Redshift"]], data[:, dic["H"]], "--", color=pHI[0].get_color())
ax.plot(
data[:, dic["Redshift"]], data[:, dic["H+"]], "--", color=pHII[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["He"]], "--", color=pHeI[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["He+"]], "--", color=pHeII[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["H-"]], "--", color=pHM[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["H2"]], "--", color=pH2I[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["H2+"]], "--", color=pH2II[0].get_color()
)
ax.plot(data[:, dic["Redshift"]], data[:, dic["D"]], "--", color=pDI[0].get_color())
ax.plot(
data[:, dic["Redshift"]], data[:, dic["D+"]], "--", color=pDII[0].get_color()
)
ax.plot(
data[:, dic["Redshift"]], data[:, dic["HD"]], "--", color=pHDI[0].get_color()
)
ax.set_xlabel(r"$\rm{1+z}$")
ax.set_xlim(100, 0.9)
ax.set_ylabel(r"$\rm{Fractional abundances}$")
ax.set_ylim(1e-21, 2)
ax.loglog()
ax.legend(loc="lower right")
plt.show()
examples/Cooling/CoolingWithPrimordialElements/run.sh 0 → 100755
+ 28
− 0
View file @ 454c2954
#!/bin/bash
# Generate the initial conditions if they are not present.
if [ ! -e homogeneousCube_32.hdf5 ]
then
echo "Fetching initial conditions..."
./getICS.sh
fi
# Get the Grackle cooling table
if [ ! -e CloudyData_UVB=HM2012.h5 ]
then
echo "Fetching the Cloudy tables required by Grackle..."
../getGrackleCoolingTable.sh
fi
# Create output directory
rm -rf snap
mkdir snap
# Run SWIFT
../../../swift --hydro --cooling --cosmology --threads=14 CoolingWithPrimordialElements.yml
# Check energy conservation and cooling rate
python3 plotAbundances.py snap/snapshot*.hdf5
examples/Cooling/getGrackleCoolingTable.sh
+ 2
− 1
View file @ 454c2954
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/CloudyData_UVB=HM2012.h5 wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/CloudyData_UVB=HM2012.h5
wget https://obswww.unige.ch/~lhausamm/links/CloudyData_UVB=HM2012_shielded.h5 wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/CloudyData_UVB=HM2012_shielded.h5
examples/parameter_example.yml
+ 12
− 0
View file @ 454c2954
...@@ -555,6 +555,18 @@ GrackleCooling: ...@@ -555,6 +555,18 @@ GrackleCooling:
H2_on_dust: 0 # Flag to enable H2 formation on dust grains H2_on_dust: 0 # Flag to enable H2 formation on dust grains
local_dust_to_gas_ratio : -1 # The ratio of total dust mass to gas mass in the local Universe (-1 to use the Grackle default value). local_dust_to_gas_ratio : -1 # The ratio of total dust mass to gas mass in the local Universe (-1 to use the Grackle default value).
cmb_temperature_floor : 1 # Enable/disable an effective CMB temperature floor cmb_temperature_floor : 1 # Enable/disable an effective CMB temperature floor
initial_nHII_to_nH_ratio: -1 # initial nHII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeI_to_nH_ratio: -1 # initial nHeI to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeII_to_nH_ratio: -1 # initial nHeII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHeIII_to_nH_ratio: -1 # initial nHeIII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nDI_to_nH_ratio: -1 # initial nDI to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nDII_to_nH_ratio: -1 # initial nDII to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHM_to_nH_ratio: -1 # initial nHM to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nH2I_to_nH_ratio: -1 # initial nH2I to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nH2II_to_nH_ratio: -1 # initial nH2II to nH ratio (number density ratio). Value is ignored if set to -1.
initial_nHDI_to_nH_ratio: -1 # initial nHDI to nH ratio (number density ratio). Value is ignored if set to -1.
......
src/cooling/grackle/cooling.c
+ 161
− 17
View file @ 454c2954
...@@ -207,32 +207,157 @@ void cooling_first_init_part(const struct phys_const* phys_const, ...@@ -207,32 +207,157 @@ void cooling_first_init_part(const struct phys_const* phys_const,
#if COOLING_GRACKLE_MODE >= 1 #if COOLING_GRACKLE_MODE >= 1
gr_float zero = 1.e-20; gr_float zero = 1.e-20;
/* NOTE: at this stage, we assume neutral gas /* NOTE: if the ratio with respect to hydrogen is given, use it. Instead,
* a better determination will be done in cooling_post_init_part */ * we assume neutral gas.
* A better determination of the abundances can be done in
* cooling_post_init_part */
/* Compute nH (formally divided by the gas density and assuming the proton
* mass to be one) */
double nH = grackle_data->HydrogenFractionByMass;
/* Compute nHe (formally divided by the gas density and assuming the proton
* mass to be one) */
double nHe = (1.f - grackle_data->HydrogenFractionByMass) / 4.f;
/* Electron density */
double ne = zero;
/* Sum of mass fraction (to test its consistency) */
double Xtot = 0.f;
/* primordial chemistry >= 1 */ /* primordial chemistry >= 1 */
/* assume neutral gas */
xp->cooling_data.HI_frac = grackle_data->HydrogenFractionByMass; /* Hydrogen I */
xp->cooling_data.HII_frac = zero; if (cooling->initial_nHII_to_nH_ratio >= 0.f)
xp->cooling_data.HeI_frac = 1. - grackle_data->HydrogenFractionByMass; xp->cooling_data.HI_frac = nH * (1.f - cooling->initial_nHII_to_nH_ratio);
xp->cooling_data.HeII_frac = zero; else
xp->cooling_data.HeIII_frac = zero; xp->cooling_data.HI_frac = nH;
xp->cooling_data.e_frac = zero;
Xtot += xp->cooling_data.HI_frac;
/* Hydrogen II */
if (cooling->initial_nHII_to_nH_ratio >= 0.f) {
double nHII = nH * cooling->initial_nHII_to_nH_ratio;
xp->cooling_data.HII_frac = nHII;
ne += nHII;
} else
xp->cooling_data.HII_frac = zero;
Xtot += xp->cooling_data.HII_frac;
/* Helium I */
if (cooling->initial_nHeI_to_nH_ratio >= 0.f) {
double nHeI = nH * cooling->initial_nHeI_to_nH_ratio;
xp->cooling_data.HeI_frac = nHeI * 4.f;
} else
xp->cooling_data.HeI_frac = nHe * 4.f;
Xtot += xp->cooling_data.HeI_frac;
/* Helium II */
if (cooling->initial_nHeII_to_nH_ratio >= 0.f) {
double nHeII = nH * cooling->initial_nHeII_to_nH_ratio;
xp->cooling_data.HeII_frac = nHeII * 4.f;
ne += nHeII;
} else
xp->cooling_data.HeII_frac = zero;
Xtot += xp->cooling_data.HeII_frac;
/* Helium III */
if (cooling->initial_nHeIII_to_nH_ratio >= 0.f) {
double nHeIII = nH * cooling->initial_nHeIII_to_nH_ratio;
xp->cooling_data.HeIII_frac = nHeIII * 4.f;
ne += 2.f * nHeIII;
} else
xp->cooling_data.HeIII_frac = zero;
Xtot += xp->cooling_data.HeIII_frac;
/* electron mass fraction (multiplied by the proton mass (Grackle convention)
*/
xp->cooling_data.e_frac = ne;
#endif // MODE >= 1 #endif // MODE >= 1
#if COOLING_GRACKLE_MODE >= 2 #if COOLING_GRACKLE_MODE >= 2
/* primordial chemistry >= 2 */ /* primordial chemistry >= 2 */
xp->cooling_data.HM_frac = zero;
xp->cooling_data.H2I_frac = zero; /* Hydrogen- */
xp->cooling_data.H2II_frac = zero; if (cooling->initial_nHM_to_nH_ratio >= 0.f) {
double nHM = nH * cooling->initial_nHM_to_nH_ratio;
xp->cooling_data.HM_frac = nHM;
ne -= nHM;
} else
xp->cooling_data.HM_frac = zero;
Xtot += xp->cooling_data.HM_frac;
/* H2I */
if (cooling->initial_nH2I_to_nH_ratio >= 0.f) {
double nH2I = nH * cooling->initial_nH2I_to_nH_ratio;
xp->cooling_data.H2I_frac = nH2I * 2.f;
} else
xp->cooling_data.H2I_frac = zero;
Xtot += xp->cooling_data.H2I_frac;
/* H2II */
if (cooling->initial_nH2II_to_nH_ratio >= 0.f) {
double nH2II = nH * cooling->initial_nH2II_to_nH_ratio;
xp->cooling_data.H2II_frac = nH2II * 2.f;
ne += nH2II;
} else
xp->cooling_data.H2II_frac = zero;
Xtot += xp->cooling_data.H2II_frac;
/* electron mass fraction (multiplied by the proton mass (Grackle convention)
*/
xp->cooling_data.e_frac = ne;
#endif // MODE >= 2 #endif // MODE >= 2
#if COOLING_GRACKLE_MODE >= 3 #if COOLING_GRACKLE_MODE >= 3
/* primordial chemistry >= 3 */ /* primordial chemistry >= 3 */
xp->cooling_data.DI_frac = grackle_data->DeuteriumToHydrogenRatio *
grackle_data->HydrogenFractionByMass; /* Deuterium I */
xp->cooling_data.DII_frac = zero; if (cooling->initial_nDI_to_nH_ratio >= 0.f) {
xp->cooling_data.HDI_frac = zero; double nDI = nH * cooling->initial_nDI_to_nH_ratio;
xp->cooling_data.DI_frac = nDI * 2.f;
} else
xp->cooling_data.DI_frac = grackle_data->DeuteriumToHydrogenRatio *
grackle_data->HydrogenFractionByMass;
Xtot += xp->cooling_data.DI_frac;
/* Deuterium II */
if (cooling->initial_nDII_to_nH_ratio >= 0.f) {
double nDII = nH * cooling->initial_nDII_to_nH_ratio;
xp->cooling_data.DII_frac = nDII * 2.f;
ne += nDII;
} else
xp->cooling_data.DII_frac = zero;
Xtot += xp->cooling_data.DII_frac;
/* HD I */
if (cooling->initial_nHDI_to_nH_ratio >= 0.f) {
double nHDI = nH * cooling->initial_nHDI_to_nH_ratio;
xp->cooling_data.HDI_frac = nHDI * 3.f;
} else
xp->cooling_data.HDI_frac = zero;
Xtot += xp->cooling_data.HDI_frac;
/* electron mass fraction (multiplied by the proton mass (Grackle convention)
*/
xp->cooling_data.e_frac = ne;
if (fabs(Xtot - 1.0) > 1e-3)
error("Got total mass fraction of gas = %.6g", Xtot);
#endif // MODE >= 3 #endif // MODE >= 3
} }
...@@ -311,7 +436,26 @@ void cooling_print_backend(const struct cooling_function_data* cooling) { ...@@ -311,7 +436,26 @@ void cooling_print_backend(const struct cooling_function_data* cooling) {
cooling->RT_HeII_ionization_rate); cooling->RT_HeII_ionization_rate);
message("grackle_chemistry_data.RT_H2_dissociation_rate = %g", message("grackle_chemistry_data.RT_H2_dissociation_rate = %g",
cooling->RT_H2_dissociation_rate); cooling->RT_H2_dissociation_rate);
message("cooling.initial_nHII_to_nH_ratio= %g",
cooling->initial_nHII_to_nH_ratio);
message("cooling.initial_nHeI_to_nH_ratio= %g",
cooling->initial_nHeI_to_nH_ratio);
message("cooling.initial_nHeII_to_nH_ratio= %g",
cooling->initial_nHeII_to_nH_ratio);
message("cooling.initial_nHeIII_to_nH_ratio= %g",
cooling->initial_nHeIII_to_nH_ratio);
message("cooling.initial_nDI_to_nH_ratio= %g",
cooling->initial_nDI_to_nH_ratio);
message("cooling.initial_nDII_to_nH_ratio= %g",
cooling->initial_nDII_to_nH_ratio);
message("cooling.initial_nHM_to_nH_ratio= %g",
cooling->initial_nHM_to_nH_ratio);
message("cooling.initial_nH2I_to_nH_ratio= %g",
cooling->initial_nH2I_to_nH_ratio);
message("cooling.initial_nH2II_to_nH_ratio= %g",
cooling->initial_nH2II_to_nH_ratio);
message("cooling.initial_nHDI_to_nH_ratio= %g",
cooling->initial_nHDI_to_nH_ratio);
message("Units:"); message("Units:");
message("\tComoving = %i", cooling->units.comoving_coordinates); message("\tComoving = %i", cooling->units.comoving_coordinates);
message("\tLength = %g", cooling->units.length_units); message("\tLength = %g", cooling->units.length_units);
......
src/cooling/grackle/cooling_io.h
+ 30
− 3
View file @ 454c2954
...@@ -207,6 +207,36 @@ __attribute__((always_inline)) INLINE static void cooling_read_parameters( ...@@ -207,6 +207,36 @@ __attribute__((always_inline)) INLINE static void cooling_read_parameters(
cooling->HydrogenFractionByMass = parser_get_opt_param_double( cooling->HydrogenFractionByMass = parser_get_opt_param_double(
parameter_file, "GrackleCooling:HydrogenFractionByMass", 0.76); parameter_file, "GrackleCooling:HydrogenFractionByMass", 0.76);
cooling->initial_nHII_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHII_to_nH_ratio", -1);
cooling->initial_nHeI_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHeI_to_nH_ratio", -1);
cooling->initial_nHeII_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHeII_to_nH_ratio", -1);
cooling->initial_nHeIII_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHeIII_to_nH_ratio", -1);
cooling->initial_nDI_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nDI_to_nH_ratio", -1);
cooling->initial_nDII_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nDII_to_nH_ratio", -1);
cooling->initial_nHM_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHM_to_nH_ratio", -1);
cooling->initial_nH2I_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nH2I_to_nH_ratio", -1);
cooling->initial_nH2II_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nH2II_to_nH_ratio", -1);
cooling->initial_nHDI_to_nH_ratio = parser_get_opt_param_double(
parameter_file, "GrackleCooling:initial_nHDI_to_nH_ratio", -1);
/* Self shielding */ /* Self shielding */
cooling->self_shielding_method = parser_get_opt_param_int( cooling->self_shielding_method = parser_get_opt_param_int(
parameter_file, "GrackleCooling:self_shielding_method", 0); parameter_file, "GrackleCooling:self_shielding_method", 0);
...@@ -216,9 +246,6 @@ __attribute__((always_inline)) INLINE static void cooling_read_parameters( ...@@ -216,9 +246,6 @@ __attribute__((always_inline)) INLINE static void cooling_read_parameters(
parameter_file, "GrackleCooling:self_shielding_threshold_atom_per_cm3"); parameter_file, "GrackleCooling:self_shielding_threshold_atom_per_cm3");
} }
cooling->HydrogenFractionByMass = parser_get_opt_param_double(
parameter_file, "GrackleCooling:HydrogenFractionByMass", 0.76);
/* Initial step convergence */ /* Initial step convergence */
cooling->max_step = parser_get_opt_param_int( cooling->max_step = parser_get_opt_param_int(
parameter_file, "GrackleCooling:max_steps", 10000); parameter_file, "GrackleCooling:max_steps", 10000);
......
src/cooling/grackle/cooling_properties.h
+ 30
− 0
View file @ 454c2954
...@@ -103,6 +103,36 @@ struct cooling_function_data { ...@@ -103,6 +103,36 @@ struct cooling_function_data {
/*! Hydrogen fraction by mass */ /*! Hydrogen fraction by mass */
float HydrogenFractionByMass; float HydrogenFractionByMass;
/*! initial nHII to nH ratio (number density ratio) */
float initial_nHII_to_nH_ratio;
/*! initial nHeI to nH ratio (number density ratio) */
float initial_nHeI_to_nH_ratio;
/*! initial nHeII to nH ratio (number density ratio) */
float initial_nHeII_to_nH_ratio;
/*! initial nHeIII to nH ratio (number density ratio) */
float initial_nHeIII_to_nH_ratio;
/*! initial nDI to nH ratio (number density ratio) */
float initial_nDI_to_nH_ratio;
/*! initial nDII to nH ratio (number density ratio) */
float initial_nDII_to_nH_ratio;
/*! initial nHM to nH ratio (number density ratio) */
float initial_nHM_to_nH_ratio;
/*! initial nH2I to nH ratio (number density ratio) */
float initial_nH2I_to_nH_ratio;
/*! initial nH2II to nH ratio (number density ratio) */
float initial_nH2II_to_nH_ratio;
/*! initial nHDI to nH ratio (number density ratio) */
float initial_nHDI_to_nH_ratio;
/*! Self shielding method (1 -> 3 for grackle's ones, 0 for none and -1 for /*! Self shielding method (1 -> 3 for grackle's ones, 0 for none and -1 for
* GEAR) */ * GEAR) */
int self_shielding_method; int self_shielding_method;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment