Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found
Select Git revision
  • 840-unit-test-testtimeline-fails
  • 875-wendland-c6-missing-neighbour-contributions
  • 887-code-does-not-compile-with-parmetis-installed-locally-but-without-metis
  • CubeTest
  • FS_Del
  • GEARRT_Iliev1
  • GEARRT_Iliev3
  • GEARRT_Iliev4
  • GEARRT_Iliev5
  • GEARRT_Iliev5-fixed-nr-subcycles
  • GEARRT_Iliev7
  • GEARRT_Iliev_static
  • GEARRT_Ivanova
  • GEARRT_fixed_nr_subcycles
  • GEARRT_injection_tests_Iliev0
  • GPU_swift
  • GrackleCoolingUpdates2
  • Lambda-T-table
  • MAGMA2
  • MAGMA2_matthieu
  • MHD_FS
  • MHD_FS_TESTs
  • MHD_FS_VP_AdvectGauge
  • MHD_Orestis
  • MHD_canvas
  • MHD_canvas_RF_128
  • MHD_canvas_RF_growth_rate
  • MHD_canvas_RobertsFlow
  • MHD_canvas_SPH_errors
  • MHD_canvas_matthieu
  • MHD_canvas_nickishch
  • MHD_canvas_nickishch_Lorentz_force_test
  • MHD_canvas_nickishch_track_everything
  • MHD_canvas_sid
  • OAK/CPAW_updates
  • OAK/LoopAdvectionTest
  • OAK/adaptive_divv
  • OAK/kinetic_dedner
  • REMIX_cosmo
  • RT_dualc
  • RT_recombination_radiation
  • RT_test_mladen
  • SIDM
  • SIDM_wKDSDK
  • SNdust
  • SPHM1RT_CosmologicalStromgrenSphere
  • SPHM1RT_bincheck
  • SPHM1RT_smoothedRT
  • TangoSIDM
  • TestPropagation3D
  • Test_fixedhProb
  • activate_fewer_comms
  • active_h_max_optimization
  • adaptive_softening_Lieuwe
  • add_2p5D
  • add_black_holes_checks
  • adding_sidm_to_master
  • agn_crksph
  • agn_crksph_subtask_speedup
  • amd-optimization
  • arm_vec
  • automatic_tasks
  • better_ray_RNG
  • black_holes_accreted_angular_momenta_from_gas
  • burkert-potential
  • c11
  • c11_atomics_copy
  • cancel_all_sorts
  • cell_exchange_improvements
  • cell_types
  • cherry-pick-cd1c39e0
  • comm_tasks_are_special
  • conduction_velocities
  • cpp-fixes
  • cuda_test
  • darwin/adaptive_softening
  • darwin/gear_chemistry_fluxes
  • darwin/gear_mechanical_feedback
  • darwin/gear_preSN_feedback
  • darwin/gear_radiation
  • darwin/simulations
  • darwin/sink_formation_proba
  • darwin/sink_mpi
  • darwin/sink_mpi_physics
  • dead-time-stats
  • derijcke_cooling
  • dev_cms
  • do-not-activate-empty-star-pairs
  • domain_zoom_nometis
  • drift_flag_debug_check
  • driven_turbulence
  • driven_turbulence_forcings
  • engineering
  • eos_updates
  • evrard_disc
  • expand_fof_2022
  • explict_bkg_cdim
  • fewer_gpart_comms
  • fewer_star_comms
  • fewer_timestep_comms_no_empty_pairs
  • v0.0
  • v0.1
  • v0.1.0-pre
  • v0.2.0
  • v0.3.0
  • v0.4.0
  • v0.5.0
  • v0.6.0
  • v0.7.0
  • v0.8.0
  • v0.8.1
  • v0.8.2
  • v0.8.3
  • v0.8.4
  • v0.8.5
  • v0.9.0
  • v1.0.0
  • v2025.01
  • v2025.04
119 results

Target

Select target project
  • dc-oman1/swiftsim
  • swift/swiftsim
  • pdraper/swiftsim
  • tkchan/swiftsim
  • dc-turn5/swiftsim
5 results
Select Git revision
  • CubeTest
  • GPU_swift
  • TangoSIDM
  • active_h_max_optimization
  • arm_vec
  • c11
  • c11_atomics_copy
  • comm_tasks_are_special
  • cuda_test
  • domain_zoom_nometis
  • drift_flag_debug_check
  • driven_turbulence
  • engineering
  • evrard_disc
  • expand_fof
  • fix_sink_timestep
  • fixed_hSIDM
  • fof_snapshots
  • gear_metal_diffusion
  • generic_cache
  • genetic_partitioning2
  • gizmo
  • gizmo_entropy_switch
  • gizmo_mfv_entropy
  • hashmap_mesh
  • isotropic_feedback
  • ivanova-testing
  • jsw/6dfof
  • kahip
  • lean_gparts
  • load-balance-testing
  • locked_hydro
  • logger_read_history
  • logger_read_history2
  • logger_write_hdf5
  • mass_dependent_h_max
  • master
  • mpi-one-thread
  • mpi-packed-parts
  • mpi-send-subparts
  • mpi-send-subparts-vector
  • mpi-subparts-vector-grav
  • mpi-testsome
  • mpi-threads
  • mpi_force_checks
  • numa_awareness
  • onesided-mpi-rdma
  • onesided-mpi-recv-cache
  • onesided-mpi-recv-window
  • onesided-mpi-single-recv-window
  • origin-master
  • parallel_exchange_cells
  • paranoid
  • phantom
  • planetary
  • planetary_boundary
  • queue-timers
  • queue-timers-clean
  • rdma-only
  • rdma-only-multiple-sends
  • rdma-only-subcopies
  • rdma-only-subparts
  • rdma-only-subparts-update
  • rdma-only-subparts-update-flamingo
  • rdma-only-subparts-update-flamingo-cellids
  • rdma-only-subparts-update-keep
  • rdma-only-subparts-update-keep-update
  • rdma-only-subsends
  • reweight-fitted-costs
  • reweight-scaled-costs
  • rgb-engineering
  • rt-gas-interactions
  • rt-ghost2-and-thermochemistry
  • scheduler_determinism
  • search-window-tests
  • signal-handler-dump
  • simba-stellar-feedback
  • sink_formation2
  • sink_merger
  • sink_merger2
  • skeleton
  • smarter_sends
  • snipes_data
  • spiral_potential
  • subgrid_SF_threshold
  • subsends
  • swift-rdma
  • swift_zoom_support
  • sync-send
  • thread-dump-extract-waiters
  • threadpool_rmapper
  • traphic
  • variable_hSIDM
  • whe-nu-bg-cosmo
  • when_to_proxy
  • yb-bhdev
  • yb-sndev
  • yb-sndev-dev
  • yb-varsndt-isotropic
  • yb-vi-gastrack
  • v0.0
  • v0.1
  • v0.1.0-pre
  • v0.2.0
  • v0.3.0
  • v0.4.0
  • v0.5.0
  • v0.6.0
  • v0.7.0
  • v0.8.0
  • v0.8.1
  • v0.8.2
  • v0.8.3
  • v0.8.4
  • v0.8.5
  • v0.9.0
116 results
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 645 additions and 15 deletions
examples/Cooling/CoolingRedshiftDependence/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/gravity_glassCube_32.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/gravity_glassCube_32.hdf5
examples/Cooling/CoolingSedovBlast_3D/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_64.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_64.hdf5
examples/Cooling/CoolingWithPrimordialElements/CoolingWithPrimordialElements.yml 0 → 100644
View file @ ed042bc2
# 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.
maximal_density_Hpcm3: -1 # Maximal density (in hydrogen atoms/cm^3) for cooling. Higher densities are floored to this value to ensure grackle works properly when interpolating beyond the cloudy_table maximal density. A value < 0 deactivates this parameter.
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
View file @ ed042bc2
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
View file @ ed042bc2
! 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
View file @ ed042bc2
#!/bin/bash
wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/homogeneousCube_32.hdf5
examples/Cooling/CoolingWithPrimordialElements/getPrimordialAbundances.py 0 → 100755
View file @ ed042bc2
#!/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
View file @ ed042bc2
# 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
View file @ ed042bc2
# 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
View file @ ed042bc2
#!/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
View file @ ed042bc2
#!/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/FeedbackEvent_3D/feedback.yml
View file @ ed042bc2
...@@ -54,6 +54,7 @@ GrackleCooling: ...@@ -54,6 +54,7 @@ GrackleCooling:
max_steps: 1000 max_steps: 1000
convergence_limit: 1e-2 convergence_limit: 1e-2
thermal_time_myr: 5 # (optional) Time (in Myr) for adiabatic cooling after a feedback event. thermal_time_myr: 5 # (optional) Time (in Myr) for adiabatic cooling after a feedback event.
maximal_density_Hpcm3: -1 # Maximal density (in hydrogen atoms/cm^3) for cooling. Higher densities are floored to this value to ensure grackle works properly when interpolating beyond the cloudy_table maximal density. A value < 0 deactivates this parameter.
GearChemistry: GearChemistry:
initial_metallicity: 0.01295 initial_metallicity: 0.01295
......
examples/Cooling/FeedbackEvent_3D/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_32.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_32.hdf5
examples/Cooling/getEagleCoolingTable.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/EAGLE/coolingtables.tar.gz wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/EAGLE/coolingtables.tar.gz
tar -xf coolingtables.tar.gz tar -xf coolingtables.tar.gz
examples/Cooling/getGrackleCoolingTable.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/CloudyData_UVB=HM2012.h5 wget https://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 https://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/CloudyData_UVB=HM2012_shielded.h5
examples/Cooling/getPS2020CoolingTables.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/COLIBRE/UV_dust1_CR1_G1_shield1.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/COLIBRE/UV_dust1_CR1_G1_shield1.hdf5
examples/Cosmology/ComovingSodShock_2D/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassPlane_128.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassPlane_128.hdf5
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassPlane_48.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassPlane_48.hdf5
examples/Cosmology/ComovingSodShock_3D/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_64.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_64.hdf5
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_32.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/glassCube_32.hdf5
examples/Cosmology/ConstantCosmoVolume/getGlass.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/gravity_glassCube_32.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/gravity_glassCube_32.hdf5
examples/EAGLE_DMO_low_z/EAGLE_DMO_100/getIC.sh
View file @ ed042bc2
#!/bin/bash #!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/EAGLE_DMO_low_z/EAGLE_DMO_ICs_100.hdf5 wget https://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/EAGLE_DMO_low_z/EAGLE_DMO_ICs_100.hdf5