{
"cells": [
{
"cell_type": "markdown",
"id": "7968866b-24b8-4ed6-bc13-6035310cb4f5",
"metadata": {},
"source": [
"[](https://github.com/open-atmos/PyPartMC/blob/main/examples/chamber.ipynb) \n",
"[](https://colab.research.google.com/github/open-atmos/PyPartMC/blob/main/examples/chamber.ipynb) \n",
"[](https://mybinder.org/v2/gh/open-atmos/PyPartMC.git/main?urlpath=lab/tree/examples/chamber.ipynb) \n",
"[](https://jupyterhub.arm.gov/hub/user-redirect/git-pull?repo=https%3A//github.com/open-atmos/PyPartMC&branch=main&urlPath=) (requires [logging in with ARM account](https://www.arm.gov/capabilities/computing-resources) and directing Jupyter to a notebook within the cloned repo)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "159edeb4",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:47.303697Z",
"start_time": "2024-11-26T04:57:47.300966Z"
}
},
"outputs": [],
"source": [
"# This file is a part of PyPartMC licensed under the GNU General Public License v3\n",
"# Copyright (C) 2023 University of Illinois Urbana-Champaign\n",
"# Authors:\n",
"# - https://github.com/compdyn/partmc/graphs/contributors\n",
"# - https://github.com/open-atmos/PyPartMC/graphs/contributors"
]
},
{
"cell_type": "markdown",
"id": "491fc802",
"metadata": {},
"source": [
"#### Chamber Simulation Example from Barrel Study, refer to [Tian el al., 2017 (Aerosol Science and Technology)](https://doi.org/10.1080/02786826.2017.1311988)"
]
},
{
"cell_type": "markdown",
"id": "f4df6740",
"metadata": {},
"source": [
"#### Equations and parameters in this notebook\n",
"1. **Wall-loss treatment**\n",
"* Diffusional wall-loss\n",
" ##### $\\alpha^{D}_{\\mu}$ = $\\frac{D(R_{me,\\mu})A_D}{\\delta_DV}$\n",
"\n",
"* Sedimental wall-loss\n",
" ##### $\\alpha^{S}_{\\mu}$ = $\\frac{4\\pi \\rho R^{3}_{m, \\mu}D(R_{me,\\mu})A_S}{3kTV}$\n",
"\n",
"* Parameters:\n",
" - $D(R_{me,\\mu})$ $\\rm (m^2s^{-1})$: Diffusion coefficient for particle $\\mu$\n",
" - $A_D$ $\\rm (m^2)$: Diffusional depostion area\n",
" - Variable name: \"area_diffuse\"\n",
" - $A_S$ $\\rm (m^2)$: Sedimentation area\n",
" - Variable name: \"area_sedi\"\n",
" - V $\\rm (m^3)$: Volume of the chamber\n",
" - Variable name: \"chamber_vol\"\n",
" - $\\delta_D$ $(\\rm m)$: Diffusional boundary layer thickness\n",
" \n",
" - $\\delta_D$ = $k_D$ $\\left(\\frac{D}{D_0}\\right)^a$\n",
" - $k_D$: Chamber-specific parameter.\n",
" - Variable name: \"prefactor_BL\"\n",
" - Can be determined for each experimental setup\n",
" - $a$: 0.25, Theoretically determined coefficient, [Fuchs, N. A. (1964). Mechanics of Aerosols.Pergamon, New York].\n",
" - Variable name: \"exponent_BL\"\n",
" - It should be in the range of (1, 3).\n",
"\n",
"2. **Fractal Treatment**\n",
"\n",
" ##### $ N$ = $\\frac{1}{f}$ $\\left(\\frac{R_{\\rm geo}}{R_0}\\right)^{d_{\\rm f}}$\n",
"\n",
"* Parameters:\n",
" - $N$: Number of monomers in a fractal-like agglomerate\n",
" - $R_{\\rm geo}$ $(\\rm m)$: Particle geometric radius\n",
" - $f$: Volume filling factor\n",
" - Variable name: \"aero_data.vol_fill_factor\"\n",
" - For spherical primary particles, this can occupy as much as 74% of the available volume\n",
" - $R_0$ $(\\rm m)$: Radius of primary particles\n",
" - Variable name: \"aero_data.prime_radius\"\n",
" - $d_{\\rm f}$: Fractal dimension\n",
" - Variable name: \"aero_data.frac_dim\"\n",
"\n",
"* Note on Fractal Parameters:\n",
" - $f$, $R_0$ $(\\rm m)$: Can be estimated by SEM images\n",
" - $d_{\\rm f}$: Needs to be specified\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "4f8359c2",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:47.331171Z",
"start_time": "2024-11-26T04:57:47.325296Z"
}
},
"outputs": [],
"source": [
"import sys\n",
"import os\n",
"if 'google.colab' in sys.modules:\n",
" !pip --quiet install open-atmos-jupyter-utils\n",
" from open_atmos_jupyter_utils import pip_install_on_colab\n",
" pip_install_on_colab('PyPartMC')\n",
"elif 'JUPYTER_IMAGE' in os.environ and '.arm.gov' in os.environ['JUPYTER_IMAGE']:\n",
" !pip --quiet install PyPartMC open_atmos_jupyter_utils\n",
" _pypartmc_path = !pip show PyPartMC | fgrep Location | cut -f2 -d' '\n",
" sys.path.extend(_pypartmc_path if _pypartmc_path[0] not in sys.path else [])"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "b494ea6e",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.195912Z",
"start_time": "2024-11-26T04:57:47.390641Z"
}
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from open_atmos_jupyter_utils import show_plot\n",
"import PyPartMC as ppmc\n",
"from PyPartMC import si"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "aa2c0567-2b75-47ba-bf6d-d5a0c22b5b83",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.214894Z",
"start_time": "2024-11-26T04:57:48.213477Z"
}
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 4,
"id": "df8d1a84",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.224403Z",
"start_time": "2024-11-26T04:57:48.221378Z"
}
},
"outputs": [],
"source": [
"gas_data = ppmc.GasData((\"H2SO4\",\"HNO3\",\"HCl\",\"NH3\",\"NO\",\"NO2\", \"NO3\",\n",
" \"N2O5\", \"HONO\", \"HNO4\", \"O3\", \"O1D\", \"O3P\", \"OH\",\n",
" \"HO2\", \"H2O2\", \"CO\", \"SO2\", \"CH4\", \"C2H6\", \"CH3O2\", \n",
" \"ETHP\", \"HCHO\", \"CH3OH\", \"ANOL\", \"CH3OOH\", \"ETHOOH\",\n",
" \"ALD2\", \"HCOOH\", \"RCOOH\", \"C2O3\", \"PAN\", \"ARO1\", \"ARO2\",\n",
" \"ALK1\", \"OLE1\", \"API1\", \"API2\", \"LIM1\", \"LIM2\", \"PAR\", \"AONE\",\n",
" \"MGLY\", \"ETH\", \"OLET\", \"OLEI\", \"TOL\", \"XYL\", \"CRES\", \"TO2\",\n",
" \"CRO\", \"OPEN\", \"ONIT\", \"ROOH\", \"RO2\", \"ANO2\", \"NAP\", \"XO2\",\n",
" \"XPAR\", \"ISOP\", \"ISOPRD\", \"ISOPP\", \"ISOPN\", \"ISOPO2\", \"API\",\n",
" \"LIM\", \"DMS\", \"MSA\", \"DMSO\", \"DMSO2\", \"CH3SO2H\", \"CH3SCH2OO\", \n",
" \"CH3SO2\", \"CH3SO3\", \"CH3SO2OO\", \"CH3SO2CH2OO\", \"SULFHOX\"\n",
" ))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "a2d7bad8",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.232295Z",
"start_time": "2024-11-26T04:57:48.230550Z"
}
},
"outputs": [],
"source": [
"env_state = ppmc.EnvState(\n",
" {\n",
" \"rel_humidity\": 0.00,\n",
" \"latitude\": 0,\n",
" \"longitude\": 0,\n",
" \"altitude\": 0 * si.m,\n",
" \"start_time\": 21600 * si.s,\n",
" \"start_day\": 200,\n",
" }\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "815283a0",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.245408Z",
"start_time": "2024-11-26T04:57:48.237835Z"
}
},
"outputs": [],
"source": [
"aero_data = ppmc.AeroData(\n",
" (\n",
" # density ions in soln (1) molecular weight kappa (1)\n",
" # | | | |\n",
" {\"SO4\": [1760 * si.kg / si.m**3, 1, 96.0 * si.g / si.mol, 0.00]},\n",
" {\"NO3\": [1800 * si.kg / si.m**3, 1, 62.0 * si.g / si.mol, 0.00]},\n",
" {\"Cl\": [2200 * si.kg / si.m**3, 1, 35.5 * si.g / si.mol, 0.00]},\n",
" {\"NH4\": [1760 * si.kg / si.m**3, 1, 18.0 * si.g / si.mol, 0.00]},\n",
" {\"MSA\": [1800 * si.kg / si.m**3, 0, 95.0 * si.g / si.mol, 0.53]},\n",
" {\"ARO1\": [1400 * si.kg / si.m**3, 0, 150.0 * si.g / si.mol, 0.10]},\n",
" {\"ARO2\": [1400 * si.kg / si.m**3, 0, 150.0 * si.g / si.mol, 0.10]},\n",
" {\"ALK1\": [1400 * si.kg / si.m**3, 0, 140.0 * si.g / si.mol, 0.10]},\n",
" {\"OLE1\": [1400 * si.kg / si.m**3, 0, 140.0 * si.g / si.mol, 0.10]},\n",
" {\"API1\": [1400 * si.kg / si.m**3, 0, 184.0 * si.g / si.mol, 0.10]},\n",
" {\"API2\": [1400 * si.kg / si.m**3, 0, 184.0 * si.g / si.mol, 0.10]},\n",
" {\"LIM1\": [1400 * si.kg / si.m**3, 0, 200.0 * si.g / si.mol, 0.10]},\n",
" {\"LIM2\": [1400 * si.kg / si.m**3, 0, 200.0 * si.g / si.mol, 0.10]},\n",
" {\"CO3\": [2600 * si.kg / si.m**3, 1, 60.0 * si.g / si.mol, 0.00]},\n",
" {\"Na\": [2200 * si.kg / si.m**3, 1, 23.0 * si.g / si.mol, 0.00]},\n",
" {\"Ca\": [2600 * si.kg / si.m**3, 1, 40.0 * si.g / si.mol, 0.00]},\n",
" {\"OIN\": [2600 * si.kg / si.m**3, 0, 1.0 * si.g / si.mol, 0.10]},\n",
" {\"OC\": [1566 * si.kg / si.m**3, 0, 1.0 * si.g / si.mol, 0.001]},\n",
" {\"BC\": [1700 * si.kg / si.m**3, 0, 1.0 * si.g / si.mol, 0.00]},\n",
" {\"H2O\": [1000 * si.kg / si.m**3, 0, 18.0 * si.g / si.mol, 0.00]},\n",
" )\n",
")\n",
"\n",
"#Adjust the following variables for fractal treatment: \n",
"#For spherical particles, comment out the variables.\n",
"#aero_data.frac_dim: Should be in the range of (1,3). \n",
"#aero_data.vol_fill_factor: Set this as 1/(percent of the available volume).\n",
"aero_data.frac_dim = 2.3\n",
"aero_data.vol_fill_factor = 1.43\n",
"aero_data.prime_radius = 4.5e-8"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "c822823d",
"metadata": {
"ExecuteTime": {
"end_time": "2024-11-26T04:57:48.257813Z",
"start_time": "2024-11-26T04:57:48.254697Z"
}
},
"outputs": [],
"source": [
"gas_state = ppmc.GasState(gas_data)\n",
"\n",
"input_gas_state = (\n",
" {\"NO\": [0.1E+00]},\n",
")\n",
"\n",
"gas_state.mix_rats = input_gas_state"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "58706963",
"metadata": {
"jupyter": {
"is_executing": true
}
},
"outputs": [],
"source": [
"times = [0 * si.s]\n",
"back_gas = [{\"time\": times},\n",
" {\"rate\": [0 / si.s]},\n",
" {\"NO\": [0.1E+00]},\n",
" ]\n",
"\n",
"gas_emit_times = [0]\n",
"gas_emit_rates = np.zeros(len(gas_emit_times))\n",
"SO2 = [4.234E-09]\n",
"\n",
"emit_gas = [\n",
" {\"time\": gas_emit_times},\n",
" {\"rate\": list(gas_emit_rates)},\n",
" {\"SO2\": SO2},\n",
"]\n",
"\n",
"AERO_DIST_BACKGROUND = {\n",
" \"back_small\": {\n",
" \"mass_frac\": [{\"SO4\": [1]}, {\"OC\": [1.375]}, {\"NH4\": [0.375]}],\n",
" \"diam_type\": \"geometric\",\n",
" \"mode_type\": \"log_normal\",\n",
" \"num_conc\": 0 / si.m**3,\n",
" \"geom_mean_diam\": 0.02 * si.um,\n",
" \"log10_geom_std_dev\": 0.161,\n",
" },\n",
"}\n",
"\n",
"AERO_DIST_EMIT = {\n",
" \"gasoline\": {\n",
" \"mass_frac\": [{\"OC\": [0.8]}, {\"BC\": [0.2]}],\n",
" \"diam_type\": \"geometric\",\n",
" \"mode_type\": \"log_normal\",\n",
" \"num_conc\": 0 / si.m**3,\n",
" \"geom_mean_diam\": 5e-8 * si.m,\n",
" \"log10_geom_std_dev\": 0.24,\n",
" },\n",
"}"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "c6a96b7d",
"metadata": {},
"outputs": [],
"source": [
"time_timeseries = [0]\n",
"pressure_timeseries = list(np.ones(1) * 1e5)\n",
"temp_timeseries = [293.4]\n",
"height_timeseries = [1000]"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "920e41e3",
"metadata": {},
"outputs": [],
"source": [
"scenario = ppmc.Scenario(\n",
" gas_data,\n",
" aero_data,\n",
" {\n",
" \"temp_profile\": [{\"time\": time_timeseries}, {\"temp\": temp_timeseries}],\n",
" \"pressure_profile\": [\n",
" {\"time\": time_timeseries},\n",
" {\"pressure\": pressure_timeseries},\n",
" ],\n",
" \"height_profile\": [{\"time\": time_timeseries}, {\"height\": height_timeseries}],\n",
" \"gas_emissions\": emit_gas,\n",
" \"gas_background\": back_gas,\n",
" \"aero_emissions\": [\n",
" {\"time\": [0 * si.s]},\n",
" {\"rate\": [0 / si.s]},\n",
" {\"dist\": [[AERO_DIST_EMIT]]},\n",
" ],\n",
" \"aero_background\": [\n",
" {\"time\": [0 * si.s]},\n",
" {\"rate\": [2.293e-5 / si.s]},\n",
" {\"dist\": [[AERO_DIST_BACKGROUND]]},\n",
" ],\n",
" \n",
" #Adjust the following variables for wall-loss treatment.\n",
" \n",
" \"loss_function\": \"chamber\",\n",
" \"chamber_vol\": 0.2093,\n",
" \"area_diffuse\": 1.988,\n",
" \"area_sedi\": 0.2463,\n",
" \"prefactor_BL\": 0.06, \n",
" \"exponent_BL\": 0.26, \n",
" },\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "6722ba83",
"metadata": {},
"outputs": [],
"source": [
"T_INITIAL = 0.0\n",
"scenario.init_env_state(env_state, T_INITIAL)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "9781ca2f",
"metadata": {},
"outputs": [],
"source": [
"diams = [1.33567e-08, 1.38477e-08, 1.43478e-08, 1.48479e-08, 1.53971e-08,\n",
" 1.59972e-08, 1.65481e-08, 1.71464e-08, 1.77975e-08, 1.84467e-08,\n",
" 1.91468e-08, 1.98469e-08, 2.05470e-08, 2.12962e-08, 2.20964e-08,\n",
" 2.28965e-08, 2.36966e-08, 2.45459e-08, 2.54460e-08, 2.63953e-08,\n",
" 2.73954e-08, 2.83956e-08, 2.94449e-08, 3.05450e-08, 3.16452e-08,\n",
" 3.27945e-08, 3.39947e-08, 3.52440e-08, 3.65442e-08, 3.78444e-08,\n",
" 3.92428e-08, 4.06940e-08, 4.21433e-08, 4.36927e-08, 4.52929e-08,\n",
" 4.69423e-08, 4.86917e-08, 5.04920e-08, 5.23414e-08, 5.42417e-08,\n",
" 5.62402e-08, 5.83406e-08, 6.04409e-08, 6.26394e-08, 6.49398e-08,\n",
" 6.72893e-08, 6.97388e-08, 7.23374e-08, 7.50379e-08, 7.77383e-08,\n",
" 8.05369e-08, 8.35356e-08, 8.66361e-08, 8.97857e-08, 9.30354e-08,\n",
" 9.64341e-08, 9.99838e-08, 1.03633e-07, 1.07432e-07, 1.11382e-07,\n",
" 1.15481e-07, 1.19731e-07, 1.24080e-07, 1.28629e-07, 1.33378e-07,\n",
" 1.38228e-07, 1.43276e-07, 1.48526e-07, 1.53975e-07, 1.59625e-07,\n",
" 1.65473e-07, 1.71572e-07, 1.77871e-07, 1.84370e-07, 1.91119e-07,\n",
" 1.98118e-07, 2.05367e-07, 2.12866e-07, 2.20664e-07, 2.28763e-07,\n",
" 2.37112e-07, 2.45810e-07, 2.54858e-07, 2.64207e-07, 2.73856e-07,\n",
" 2.83854e-07, 2.94252e-07, 3.05050e-07, 3.16249e-07, 3.27847e-07,\n",
" 3.39845e-07, 3.52293e-07, 3.65191e-07, 3.78539e-07, 3.92387e-07,\n",
" 4.06784e-07, 4.21732e-07, 4.37180e-07, 4.53177e-07, 4.69774e-07,\n",
" 4.86971e-07, 5.04818e-07, 5.23316e-07, 5.42462e-07, 5.62309e-07,\n",
" 5.82905e-07, 6.04302e-07, 6.26449e-07, 6.49395e-07, 6.73191e-07,\n",
" 6.97837e-07, 7.23383e-07, 7.49878e-07, 7.77374e-07, 8.05870e-07,\n",
" 8.35396e-07]\n",
"\n",
"concs = [1.03e+07, 1.04e+07, 4.39e+06, 4.44e+06, 8.03e+06, 7.53e+06,\n",
" 3.77e+06, 9.10e+06, 9.77e+06, 1.32e+07, 1.54e+07, 2.26e+07,\n",
" 2.14e+07, 3.20e+07, 2.68e+07, 3.36e+07, 4.62e+07, 6.90e+07,\n",
" 8.63e+07, 9.12e+07, 1.11e+08, 1.50e+08, 1.57e+08, 1.96e+08,\n",
" 2.44e+08, 2.83e+08, 3.48e+08, 3.99e+08, 4.47e+08, 5.42e+08,\n",
" 5.91e+08, 6.37e+08, 7.77e+08, 8.68e+08, 9.50e+08, 1.13e+09,\n",
" 1.25e+09, 1.36e+09, 1.47e+09, 1.66e+09, 1.90e+09, 1.94e+09,\n",
" 2.14e+09, 2.38e+09, 2.40e+09, 2.58e+09, 2.81e+09, 2.87e+09,\n",
" 2.96e+09, 3.13e+09, 3.27e+09, 3.53e+09, 3.51e+09, 3.60e+09,\n",
" 3.81e+09, 3.90e+09, 3.87e+09, 3.89e+09, 3.91e+09, 3.94e+09,\n",
" 3.94e+09, 3.76e+09, 3.76e+09, 3.74e+09, 3.53e+09, 3.38e+09,\n",
" 3.34e+09, 3.19e+09, 3.00e+09, 2.81e+09, 2.56e+09, 2.42e+09,\n",
" 2.18e+09, 1.96e+09, 1.74e+09, 1.52e+09, 1.36e+09, 1.15e+09,\n",
" 9.86e+08, 8.00e+08, 6.91e+08, 5.82e+08, 4.58e+08, 3.65e+08,\n",
" 2.80e+08, 2.08e+08, 1.64e+08, 1.25e+08, 8.43e+07, 6.77e+07,\n",
" 6.23e+07, 3.54e+07, 3.59e+07, 3.04e+07, 1.90e+07, 1.31e+07,\n",
" 1.23e+07, 1.26e+07, 9.70e+06, 1.12e+07, 1.25e+07, 6.67e+06,\n",
" 8.67e+06, 8.51e+06, 6.23e+06, 5.12e+06, 5.63e+06, 5.36e+06,\n",
" 4.29e+06, 5.55e+06, 4.83e+06, 4.47e+06, 2.90e+06, 2.11e+06,\n",
" 3.86e+06]\n",
"\n",
"# Update \"mass_frac:\" with the desired aerosol type\n",
"# If \"mode_type\" is \"sampled\", use the data input in \"size_dist\"\n",
"# If \"mode_type\" is \"log_normal\", replace \"size_dist\" with the following:\n",
"# \"num_conc\": your_number concentration_data / ppmc.si.m**3,\n",
"# \"geom_mean_diam\": your_geomean_diam_data* ppmc.si.m,\n",
"# \"log10_geom_std_dev\": your_log10_geometric_std_dev_data,\n",
"\n",
"AERO_DIST_INIT = [\n",
" {\n",
" \"test_mode\": {\n",
" \"mass_frac\": [{\"SO4\": [96]}, {\"NH4\": [36]}],\n",
" \"diam_type\": \"geometric\",\n",
" \"mode_type\": \"sampled\",\n",
" \"size_dist\": [\n",
" {\"diam\": diams},\n",
" {\"num_conc\": concs},\n",
" ],\n",
" }\n",
" }\n",
"]\n",
"\n",
"aero_dist_init = ppmc.AeroDist(aero_data, AERO_DIST_INIT)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "d8d9c8fd",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"33273"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"run_part_opt = ppmc.RunPartOpt(\n",
" {\n",
" \"output_prefix\": \"chamber\",\n",
" \"do_coagulation\": True,\n",
" \"coag_kernel\": \"brown\",\n",
" \"t_max\": 20160 * si.s,\n",
" \"del_t\": 60 * si.s,\n",
" }\n",
")\n",
"\n",
"N_PART = 100000\n",
"aero_state = ppmc.AeroState(aero_data, N_PART, 'nummass_source')\n",
"aero_state.dist_sample(\n",
" aero_dist_init,\n",
" sample_prop=1.0,\n",
" create_time=0.0,\n",
" allow_doubling=True,\n",
" allow_halving=True,\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "5c8aef2b",
"metadata": {},
"outputs": [],
"source": [
"camp_core = ppmc.CampCore()\n",
"photolysis = ppmc.Photolysis()"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "0f7c29fe",
"metadata": {},
"outputs": [],
"source": [
"N_STEPS = int(run_part_opt.t_max / run_part_opt.del_t)\n",
"num_conc = np.zeros(N_STEPS + 1)\n",
"num_conc[0] = aero_state.total_num_conc\n",
"mass_conc = np.zeros(N_STEPS + 1)\n",
"mass_conc[0] = aero_state.total_mass_conc\n",
"time = np.zeros(N_STEPS + 1)\n",
"gas_mix_rat = np.zeros((N_STEPS + 1, gas_state.n_spec))\n",
"gas_mix_rat[0, :] = gas_state.mix_rats\n",
"\n",
"height = np.zeros((N_STEPS + 1))\n",
"temperature = np.zeros((N_STEPS + 1))\n",
"rh = np.zeros((N_STEPS + 1))\n",
"\n",
"height[0] = env_state.height\n",
"temperature[0] = env_state.temp\n",
"rh[0] = env_state.rh\n",
"\n",
"diam_grid = ppmc.BinGrid(30, \"log\", 1e-9, 1e-6)\n",
"dists = []\n",
"dry_diameters = aero_state.dry_diameters\n",
"num_concs = aero_state.num_concs\n",
"dists.append(ppmc.histogram_1d(diam_grid, dry_diameters, num_concs))\n",
"dist_times = [0]\n",
"\n",
"last_output_time = 0.\n",
"last_progress_time = 0.\n",
"i_output = 1\n",
"\n",
"for i_time in range(1,N_STEPS + 1):\n",
" (last_output_time, last_progress_time, i_output) = ppmc.run_part_timestep(\n",
" scenario,\n",
" env_state,\n",
" aero_data,\n",
" aero_state,\n",
" gas_data,\n",
" gas_state,\n",
" run_part_opt,\n",
" camp_core,\n",
" photolysis,\n",
" i_time,\n",
" T_INITIAL,\n",
" last_output_time,\n",
" last_progress_time,\n",
" i_output\n",
" )\n",
"\n",
" num_conc[i_time] = aero_state.total_num_conc\n",
" mass_conc[i_time] = aero_state.total_mass_conc\n",
" time[i_time] = env_state.elapsed_time\n",
" gas_mix_rat[i_time, :] = gas_state.mix_rats\n",
" height[i_time] = env_state.height\n",
" temperature[i_time] = env_state.temp\n",
" rh[i_time] = env_state.rh\n",
" if np.mod(i_time * run_part_opt.del_t, 420.0) == 0:\n",
" dry_diameters = aero_state.dry_diameters\n",
" mobility_diameters = aero_state.mobility_diameters(env_state)\n",
" num_concs = aero_state.num_concs\n",
" dists.append(ppmc.histogram_1d(diam_grid, dry_diameters, num_concs))\n",
" dist_times.append(env_state.elapsed_time)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "82c0cebb",
"metadata": {},
"outputs": [],
"source": [
"plt.rcParams.update({'font.size': 9})\n",
"plt.rcParams.update({'figure.figsize': (3.08,2.5)})\n",
"plt.rcParams.update({\"axes.grid\" : True})"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "9ab2aeed-c8f3-4683-82f0-2837d7191372",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Writing ref_aero_time.txt\n"
]
}
],
"source": [
"%%file ref_aero_time.txt\n",
"0\t1.30898E+11\t3.25738E-07\n",
"420\t1.19515E+11\t3.21076E-07\n",
"840\t1.09495E+11\t3.18544E-07\n",
"1260\t1.01068E+11\t3.1369E-07\n",
"1680\t94134900000\t3.1467E-07\n",
"2100\t87897200000\t3.10829E-07\n",
"2520\t82551400000\t3.08678E-07\n",
"2940\t77647400000\t3.04676E-07\n",
"3360\t73623800000\t3.01067E-07\n",
"3780\t67930500000\t2.95461E-07\n",
"4200\t65888100000\t2.90071E-07\n",
"4620\t62582400000\t2.84073E-07\n",
"5040\t59765600000\t2.80337E-07\n",
"5460\t57719900000\t2.8295E-07\n",
"5880\t54980200000\t2.83708E-07\n",
"6300\t52876900000\t2.78871E-07\n",
"6720\t51163400000\t2.76975E-07\n",
"7140\t49071700000\t2.75866E-07\n",
"7560\t46814100000\t2.71413E-07\n",
"7980\t45264200000\t2.67989E-07\n",
"8400\t43484200000\t2.6377E-07\n",
"8820\t42400400000\t2.64128E-07\n",
"9240\t40871900000\t2.58219E-07\n",
"9660\t39407200000\t2.53514E-07\n",
"10080\t38229100000\t2.52278E-07\n",
"10500\t37014700000\t2.4769E-07\n",
"10920\t35644000000\t2.42154E-07\n",
"11340\t34410900000\t2.40518E-07\n",
"11760\t33681200000\t2.43037E-07\n",
"12180\t32030300000\t2.40009E-07\n",
"12600\t31212500000\t2.38165E-07\n",
"13020\t30291900000\t2.34896E-07\n",
"13440\t29387600000\t2.30355E-07\n",
"13860\t28567600000\t2.30407E-07\n",
"14280\t27855500000\t2.28884E-07\n",
"14700\t26997800000\t2.23884E-07\n",
"15120\t26348700000\t2.21466E-07\n",
"15540\t25458700000\t2.15623E-07\n",
"15960\t24840600000\t2.14156E-07\n",
"16380\t24219400000\t2.09616E-07\n",
"16800\t23710900000\t2.05986E-07\n",
"17220\t23038300000\t2.05518E-07\n",
"17640\t22625600000\t2.02444E-07\n",
"18060\t22069200000\t2.03548E-07\n",
"18480\t21447500000\t1.99564E-07\n",
"18900\t20930500000\t1.98804E-07\n",
"19320\t20511200000\t1.98153E-07\n",
"19740\t20020700000\t1.95232E-07\n",
"20160\t19622600000\t1.90199E-07"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "72077383-7d47-4265-ae8c-a1dcd10db293",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Writing ref_aero_size_num.txt\n"
]
}
],
"source": [
"%%file ref_aero_size_num.txt\n",
"1.36E-08\t6.60E+08\t6.53E+07\t1.50E+08\t2.61E+08\t6.53E+07\t6.53E+07\t1.96E+08\t9.09E+07\t1.31E+08\t1.31E+08\t6.53E+07\t0\t6.53E+07\t6.53E+07\t0\t1.96E+08\t6.53E+07\t1.31E+08\t6.53E+07\t0\t0\t1.96E+08\t0\t6.58E+06\t0\t0\t0\t1.96E+08\t1.31E+08\t0\t1.31E+08\t0\t6.53E+07\t6.53E+07\t6.53E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"1.41E-08\t6.76E+08\t4.34E+08\t1.24E+08\t6.21E+07\t0\t6.21E+07\t1.24E+08\t5.58E+07\t0\t0\t0\t0\t0\t1.24E+08\t6.21E+07\t1.24E+08\t6.21E+07\t0\t0\t0\t0\t0\t6.21E+07\t5.58E+07\t0\t0\t0\t6.21E+07\t6.21E+07\t0\t0\t0\t0\t6.21E+07\t0\t0\t0\t0\t6.21E+07\t0\t2.94E+07\t6.21E+07\t6.21E+07\t6.21E+07\t0\t6.21E+07\t0\t0\t6.21E+07\n",
"1.46E-08\t2.95E+08\t1.18E+08\t1.77E+08\t4.87E+07\t0\t5.90E+07\t5.90E+07\t1.77E+08\t0\t5.90E+07\t0\t0\t0\t5.90E+07\t5.90E+07\t1.18E+08\t0\t5.90E+07\t0\t0\t0\t0\t0\t0\t4.87E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t5.90E+07\t0\t0\t5.90E+07\t0\t1.49E+08\t5.90E+07\t0\t5.90E+07\t0\t0\t5.90E+07\t0\t0\n",
"1.51E-08\t2.82E+08\t3.38E+08\t2.82E+08\t9.82E+06\t5.63E+07\t0\t8.81E+06\t0\t1.13E+08\t0\t0\t0\t0\t0\t0\t0\t1.13E+08\t0\t0\t0\t0\t1.13E+08\t5.63E+07\t0\t9.82E+06\t0\t0\t5.63E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t5.63E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"1.57E-08\t4.83E+08\t3.76E+08\t1.32E+08\t5.37E+07\t5.37E+07\t1.07E+08\t9.90E+07\t0\t1.07E+08\t1.07E+08\t0\t0\t0\t5.37E+07\t0\t0\t0\t0\t0\t0\t5.37E+07\t0\t0\t0\t0\t0\t0\t5.37E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t2.50E+07\t0\t0\t5.37E+07\t0\t0\t0\t0\n",
"1.63E-08\t5.12E+08\t3.59E+08\t1.30E+08\t5.12E+07\t5.12E+07\t0\t1.54E+08\t5.12E+07\t0\t0\t0\t0\t0\t5.12E+07\t0\t0\t0\t0\t0\t0\t0\t0\t5.12E+07\t0\t5.12E+07\t0\t0\t5.12E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t5.12E+07\t7.87E+07\t5.12E+07\t0\t0\t5.12E+07\t0\t0\t0\n",
"1.68E-08\t2.45E+08\t2.94E+08\t1.47E+08\t0\t0\t1.47E+08\t1.47E+08\t0\t4.89E+07\t0\t0\t0\t0\t0\t9.79E+07\t4.89E+07\t4.89E+07\t4.89E+07\t0\t4.89E+07\t0\t0\t0\t0\t0\t0\t4.89E+07\t4.89E+07\t4.89E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t4.89E+07\t0\t0\t0\t0\t0\t0\t0\n",
"1.75E-08\t5.62E+08\t4.33E+08\t3.40E+08\t4.68E+07\t9.37E+07\t4.68E+07\t4.68E+07\t0\t0\t4.68E+07\t9.37E+07\t0\t0\t4.68E+07\t0\t4.68E+07\t0\t0\t0\t0\t0\t1.16E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"1.81E-08\t6.28E+08\t6.62E+08\t3.48E+08\t1.35E+08\t8.97E+07\t1.35E+08\t0\t8.97E+07\t4.48E+07\t8.97E+07\t8.97E+07\t0\t0\t0\t0\t0\t0\t4.48E+07\t4.48E+07\t4.48E+07\t0\t1.23E+08\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t4.48E+07\t0\t0\t0\t4.48E+07\n",
"1.88E-08\t8.17E+08\t6.02E+08\t2.15E+08\t1.29E+08\t8.60E+07\t0\t8.60E+07\t4.30E+07\t0\t0\t0\t4.30E+07\t0\t4.30E+07\t4.30E+07\t4.30E+07\t0\t0\t0\t0\t8.60E+07\t4.30E+07\t0\t0\t0\t0\t0\t0\t4.30E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t4.30E+07\t0\t0\t0\t0\t0\t0\t0\t0\n",
"1.95E-08\t9.90E+08\t7.42E+08\t3.30E+08\t4.12E+07\t0\t8.25E+07\t1.24E+08\t8.25E+07\t0\t0\t4.12E+07\t0\t0\t0\t4.12E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t4.12E+07\t4.12E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.02E-08\t1.50E+09\t8.31E+08\t3.17E+08\t1.58E+08\t1.98E+08\t7.92E+07\t7.92E+07\t7.92E+07\t3.96E+07\t0\t0\t3.96E+07\t3.96E+07\t0\t0\t0\t0\t0\t0\t0\t0\t3.96E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.09E-08\t1.38E+09\t8.92E+08\t7.89E+08\t2.26E+08\t1.25E+08\t1.36E+08\t1.72E+08\t1.01E+08\t1.03E+08\t0\t3.22E+07\t0\t1.11E+08\t7.33E+07\t3.52E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t3.77E+07\t3.79E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t3.80E+07\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.17E-08\t2.00E+09\t1.74E+09\t8.05E+08\t3.11E+08\t3.68E+08\t4.79E+07\t1.19E+08\t9.36E+07\t0\t3.46E+07\t0\t6.49E+07\t0\t0\t3.59E+07\t0\t3.33E+07\t0\t3.41E+07\t3.47E+07\t0\t0\t0\t0\t0\t7.21E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.25E-08\t1.74E+09\t1.43E+09\t7.44E+08\t6.24E+08\t2.68E+08\t1.31E+08\t1.79E+08\t1.48E+08\t4.94E+07\t6.06E+07\t5.76E+07\t0\t2.61E+07\t1.33E+08\t9.22E+07\t3.34E+07\t0\t2.85E+07\t3.03E+07\t6.78E+07\t0\t3.48E+07\t3.18E+07\t0\t0\t0\t0\t0\t3.50E+07\t3.52E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.33E-08\t2.25E+09\t1.53E+09\t7.73E+08\t2.85E+08\t2.63E+08\t2.47E+08\t1.21E+08\t6.02E+07\t1.40E+08\t0\t4.33E+07\t0\t2.41E+07\t5.73E+07\t5.72E+07\t2.69E+07\t0\t0\t2.95E+07\t0\t0\t2.80E+07\t0\t0\t2.17E+06\t0\t3.18E+07\t0\t3.38E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.41E-08\t3.02E+09\t2.50E+09\t1.73E+09\t1.13E+09\t6.91E+08\t4.27E+08\t4.20E+08\t1.84E+08\t2.53E+08\t4.79E+07\t7.33E+07\t0\t5.26E+07\t1.78E+08\t0\t1.86E+07\t5.82E+07\t0\t0\t0\t0\t3.07E+07\t0\t2.44E+07\t2.77E+07\t6.32E+07\t3.15E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t3.25E+07\n",
"2.50E-08\t4.41E+09\t2.75E+09\t2.08E+09\t9.45E+08\t5.30E+08\t3.82E+08\t2.80E+08\t1.46E+08\t1.63E+08\t0\t2.12E+07\t1.18E+08\t3.96E+07\t7.59E+07\t7.44E+07\t1.12E+08\t8.10E+07\t7.80E+07\t2.26E+07\t0\t2.09E+07\t0\t0\t0\t0\t5.79E+07\t0\t0\t2.77E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.59E-08\t5.43E+09\t3.04E+09\t2.39E+09\t1.44E+09\t6.82E+08\t6.77E+08\t5.30E+08\t3.84E+08\t2.02E+08\t2.10E+07\t8.61E+07\t1.85E+08\t9.72E+06\t2.09E+08\t1.43E+08\t603192\t3.29E+07\t3.89E+07\t4.13E+07\t1.17E+07\t7.47E+07\t7.45E+07\t2.15E+07\t2.41E+07\t0\t2.37E+07\t0\t2.26E+07\t2.75E+07\t0\t0\t2.89E+07\t0\t0\t0\t0\t0\t0\t0\t0\t2.98E+07\t0\t0\t0\t0\t0\t0\t0\t0\n",
"2.69E-08\t5.65E+09\t4.55E+09\t2.25E+09\t1.67E+09\t7.85E+08\t9.82E+08\t7.37E+08\t4.46E+08\t2.11E+08\t2.32E+08\t2.84E+08\t6.96E+07\t1.55E+08\t6.18E+07\t5.11E+07\t7.31E+07\t5.01E+07\t1.76E+08\t2.97E+07\t0\t7.49E+06\t6.60E+06\t9.73E+07\t1.44E+07\t0\t2.31E+07\t7.65E+07\t2.00E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t2.87E+07\t0\t0\t0\t2.88E+07\t0\t0\t0\t0\n",
"2.79E-08\t7.11E+09\t5.12E+09\t2.72E+09\t1.87E+09\t1.27E+09\t1.23E+09\t7.09E+08\t5.77E+08\t3.66E+08\t1.40E+08\t4.46E+08\t2.85E+08\t1.07E+08\t1.06E+08\t1.27E+08\t1.76E+08\t1.48E+08\t9.51E+07\t8.91E+07\t5.07E+07\t8.31E+07\t0\t0\t2.39E+07\t3.60E+07\t0\t6.69E+07\t1.18E+07\t0\t7.85E+07\t0\t0\t0\t0\t0\t2.78E+07\t0\t0\t0\t0\t0\t2.80E+07\t0\t0\t0\t2.63E+07\t0\t0\t0\n",
"2.89E-08\t9.50E+09\t6.29E+09\t4.07E+09\t2.62E+09\t1.68E+09\t1.31E+09\t1.32E+09\t9.67E+08\t6.81E+08\t1.57E+08\t4.06E+08\t2.93E+08\t3.71E+08\t3.05E+08\t0\t1.72E+08\t5.23E+07\t3.38E+07\t2.20E+08\t9.46E+06\t0\t1.69E+07\t4.98E+07\t1.07E+08\t4.33E+07\t7.88E+07\t8.85E+06\t0\t6.89E+06\t0\t1.86E+07\t1.87E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"3.00E-08\t9.83E+09\t8.83E+09\t4.24E+09\t3.69E+09\t2.48E+09\t1.49E+09\t1.63E+09\t1.17E+09\t9.70E+08\t1.02E+08\t5.24E+08\t2.63E+08\t2.97E+08\t2.48E+08\t2.54E+08\t2.31E+08\t6.91E+07\t1.33E+08\t8.09E+07\t1.82E+08\t1.55E+08\t1.25E+08\t5.57E+07\t9.31E+07\t1.18E+08\t2.03E+07\t0\t2.22E+07\t3.01E+07\t1.24E+07\t0\t0\t1.28E+07\t0\t1.44E+07\t2.42E+07\t2.39E+07\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n",
"3.11E-08\t1.28E+10\t9.66E+09\t5.80E+09\t4.39E+09\t2.63E+09\t1.95E+09\t2.07E+09\t1.24E+09\t1.25E+09\t1.48E+08\t4.21E+08\t6.44E+08\t4.63E+08\t2.50E+08\t4.54E+07\t2.23E+08\t2.56E+08\t2.31E+08\t1.94E+08\t1.38E+08\t6.06E+07\t1.75E+08\t2.01E+08\t1.65E+08\t1.21E+07\t7.43E+07\t1.21E+08\t1.51E+08\t1.91E+07\t0\t8.22E+07\t3.22E+06\t1.07E+07\t3.89E+07\t0\t0\t0\t0\t4.77E+07\t0\t0\t0\t0\t2.09E+07\t0\t0\t4.82E+07\t1.81E+07\t0\n",
"3.22E-08\t1.57E+10\t1.11E+10\t6.88E+09\t5.18E+09\t3.33E+09\t3.08E+09\t2.17E+09\t1.89E+09\t1.42E+09\t6.70E+08\t8.75E+08\t6.33E+08\t6.19E+08\t5.31E+08\t9.20E+08\t1.17E+08\t3.93E+08\t4.58E+08\t3.31E+08\t1.73E+08\t1.64E+08\t2.54E+07\t2.31E+08\t6.30E+07\t3.78E+06\t1.01E+08\t2.01E+08\t0\t1.30E+07\t0\t0\t0\t0\t0\t2.23E+07\t0\t1.15E+06\t3.13E+07\t4.24E+07\t1.66E+07\t0\t0\t2.18E+07\t0\t0\t0\t0\t0\t0\n",
"3.34E-08\t1.81E+10\t1.37E+10\t8.71E+09\t5.70E+09\t4.79E+09\t3.67E+09\t2.91E+09\t2.48E+09\t1.76E+09\t7.01E+08\t1.47E+09\t8.27E+08\t5.86E+08\t6.25E+08\t6.32E+08\t4.21E+08\t5.68E+08\t3.83E+08\t2.63E+08\t1.75E+08\t1.27E+08\t3.52E+08\t1.77E+08\t1.31E+08\t1.01E+08\t1.49E+08\t1.24E+08\t3.87E+07\t5.74E+06\t4.03E+07\t5.85E+07\t2.84E+07\t0\t2.30E+06\t4.28E+07\t0\t0\t4.63E+07\t0\t0\t2.71E+07\t0\t0\t9.04E+06\t0\t1.30E+07\t0\t0\t1.54E+07\n",
"3.46E-08\t2.22E+10\t1.54E+10\t1.03E+10\t8.71E+09\t6.25E+09\t4.72E+09\t3.54E+09\t2.29E+09\t2.10E+09\t8.96E+08\t1.41E+09\t1.30E+09\t6.49E+08\t8.84E+08\t8.07E+08\t7.42E+08\t3.54E+08\t7.43E+08\t5.07E+08\t4.71E+08\t2.63E+08\t9.41E+07\t2.79E+08\t4.34E+08\t5.40E+07\t9.26E+07\t5.30E+07\t2.42E+07\t0\t2.95E+07\t0\t7.42E+07\t0\t1.85E+07\t6.07E+06\t0\t0\t0\t2.25E+07\t0\t5.19E+07\t4.18E+07\t0\t2.35E+07\t4.02E+06\t0\t0\t2.82E+07\t0\n",
"3.59E-08\t2.53E+10\t1.84E+10\t1.24E+10\t8.72E+09\t6.01E+09\t5.66E+09\t4.55E+09\t3.45E+09\t2.61E+09\t1.57E+09\t1.82E+09\t1.84E+09\t1.06E+09\t8.47E+08\t9.17E+08\t6.10E+08\t6.04E+08\t7.61E+08\t4.04E+08\t5.01E+08\t4.88E+08\t1.94E+08\t3.48E+08\t2.72E+08\t3.49E+08\t2.32E+08\t2.29E+07\t2.33E+08\t1.79E+08\t1.63E+07\t0\t8.28E+06\t5.94E+06\t0\t0\t1.08E+08\t1.83E+07\t7.18E+07\t2.58E+07\t0\t105038\t3.43E+07\t0\t1.81E+07\t1.70E+07\t6.72E+06\t7.56E+07\t6.23E+07\t0\n",
"3.72E-08\t2.95E+10\t2.16E+10\t1.40E+10\t1.08E+10\t9.02E+09\t6.98E+09\t5.68E+09\t4.27E+09\t4.17E+09\t1.64E+09\t2.74E+09\t2.42E+09\t2.18E+09\t1.39E+09\t1.50E+09\t1.23E+09\t1.16E+09\t9.14E+08\t6.67E+08\t7.69E+08\t6.91E+08\t6.98E+08\t3.53E+08\t2.53E+08\t3.22E+08\t2.88E+08\t2.39E+08\t3.31E+08\t1.11E+08\t0\t4.21E+07\t5.33E+07\t1.32E+08\t2.95E+07\t3.55E+07\t1.30E+08\t1.43E+08\t4.85E+07\t6.72E+07\t8.07E+07\t2.79E+07\t3.89E+07\t0\t4.49E+06\t0\t6.48E+07\t0\t1.22E+07\t1.44E+07\n",
"3.85E-08\t3.44E+10\t2.58E+10\t1.85E+10\t1.37E+10\t1.13E+10\t8.69E+09\t6.28E+09\t5.24E+09\t4.26E+09\t2.27E+09\t3.37E+09\t2.87E+09\t1.83E+09\t2.12E+09\t1.46E+09\t1.64E+09\t1.43E+09\t1.12E+09\t9.26E+08\t1.15E+09\t8.15E+08\t8.54E+08\t7.59E+08\t5.20E+08\t3.63E+08\t2.05E+08\t4.35E+08\t3.43E+08\t3.10E+08\t2.48E+08\t2.97E+08\t2.11E+08\t1.23E+07\t1.20E+08\t7.61E+07\t0\t2.44E+07\t1.40E+08\t0\t0\t3.01E+07\t2.98E+07\t0\t93744.2\t2.87E+07\t1.67E+07\t0\t2.20E+07\t0\n",
"4.00E-08\t3.75E+10\t2.97E+10\t2.17E+10\t1.67E+10\t1.26E+10\t9.92E+09\t8.88E+09\t7.27E+09\t5.98E+09\t3.45E+09\t4.27E+09\t3.75E+09\t3.08E+09\t2.05E+09\t2.38E+09\t2.23E+09\t1.61E+09\t1.55E+09\t1.02E+09\t1.13E+09\t1.01E+09\t9.76E+08\t9.79E+08\t1.14E+09\t6.96E+08\t6.54E+08\t6.26E+08\t5.34E+08\t2.80E+08\t2.22E+08\t7.33E+07\t1.05E+08\t9.90E+07\t7.71E+06\t9.01E+07\t0\t1.90E+08\t1.77E+08\t7.13E+07\t1.68E+08\t1.71E+07\t0\t0\t0\t0\t7.79E+07\t0\t0\t2.20E+07\n",
"4.14E-08\t4.19E+10\t3.44E+10\t2.45E+10\t1.93E+10\t1.57E+10\t1.26E+10\t1.04E+10\t8.75E+09\t6.00E+09\t3.84E+09\t5.40E+09\t4.45E+09\t3.67E+09\t3.51E+09\t3.06E+09\t2.66E+09\t2.65E+09\t1.96E+09\t1.63E+09\t1.61E+09\t1.28E+09\t1.11E+09\t9.58E+08\t9.42E+08\t8.03E+08\t9.72E+08\t5.57E+08\t8.85E+08\t6.38E+08\t1.89E+08\t2.83E+08\t2.78E+08\t3.89E+07\t1.38E+08\t1.93E+08\t1.66E+08\t1.34E+08\t8.61E+07\t1.73E+08\t5.51E+07\t0\t1.22E+08\t0\t8.46E+07\t0\t0\t0\t1.25E+08\t6.91E+07\n",
"4.29E-08\t4.95E+10\t3.81E+10\t2.90E+10\t2.46E+10\t1.79E+10\t1.42E+10\t1.24E+10\t9.46E+09\t8.70E+09\t5.25E+09\t6.54E+09\t5.16E+09\t4.97E+09\t3.86E+09\t3.95E+09\t3.14E+09\t2.45E+09\t2.45E+09\t2.18E+09\t2.09E+09\t1.51E+09\t1.68E+09\t1.37E+09\t1.02E+09\t1.43E+09\t8.54E+08\t7.12E+08\t7.83E+08\t6.17E+08\t3.72E+08\t1.96E+08\t2.19E+08\t3.91E+08\t1.94E+08\t3.52E+08\t3.34E+07\t5.89E+07\t2.28E+08\t3.06E+08\t2.88E+08\t0\t1.48E+08\t1.65E+08\t1.05E+08\t0\t6.89E+07\t9.83E+07\t603336\t1.16E+08\n",
"4.45E-08\t5.56E+10\t4.24E+10\t3.39E+10\t2.72E+10\t2.25E+10\t1.82E+10\t1.42E+10\t1.20E+10\t1.13E+10\t5.89E+09\t7.24E+09\t6.13E+09\t5.54E+09\t4.81E+09\t4.51E+09\t3.73E+09\t3.01E+09\t3.02E+09\t2.33E+09\t2.32E+09\t2.34E+09\t2.00E+09\t1.83E+09\t1.62E+09\t1.56E+09\t1.36E+09\t1.31E+09\t1.40E+09\t8.08E+08\t5.63E+08\t5.19E+08\t5.82E+08\t3.41E+08\t4.73E+07\t2.62E+08\t2.49E+08\t1.65E+08\t2.15E+08\t9.28E+07\t1.88E+08\t0\t1.35E+08\t3.06E+07\t1.16E+08\t1.96E+08\t0\t6.73E+07\t1.77E+08\t1.29E+08\n",
"4.61E-08\t6.11E+10\t4.99E+10\t3.81E+10\t3.06E+10\t2.55E+10\t2.02E+10\t1.73E+10\t1.35E+10\t1.26E+10\t8.69E+09\t8.46E+09\t7.46E+09\t6.64E+09\t6.19E+09\t6.08E+09\t4.93E+09\t4.93E+09\t4.81E+09\t3.54E+09\t2.77E+09\t2.62E+09\t2.45E+09\t2.16E+09\t2.09E+09\t2.24E+09\t1.57E+09\t1.23E+09\t1.40E+09\t9.22E+08\t7.65E+08\t5.82E+08\t5.38E+08\t6.14E+08\t5.12E+08\t6.52E+08\t2.67E+08\t5.53E+08\t4.19E+08\t1.07E+08\t1.33E+08\t3.72E+08\t1.55E+08\t6.99E+07\t1.63E+08\t8.81E+07\t2.13E+08\t1.71E+08\t2.95E+08\t5.37E+06\n",
"4.78E-08\t7.09E+10\t5.67E+10\t4.39E+10\t3.78E+10\t2.94E+10\t2.28E+10\t2.18E+10\t1.58E+10\t1.35E+10\t1.04E+10\t1.10E+10\t9.13E+09\t8.20E+09\t7.48E+09\t6.45E+09\t6.44E+09\t5.36E+09\t5.07E+09\t4.50E+09\t3.83E+09\t3.47E+09\t2.75E+09\t3.41E+09\t2.75E+09\t2.25E+09\t2.46E+09\t1.74E+09\t1.66E+09\t1.61E+09\t1.10E+09\t7.97E+08\t8.82E+08\t9.13E+08\t3.78E+08\t6.08E+08\t7.46E+08\t5.62E+08\t5.14E+08\t5.02E+08\t4.88E+08\t3.59E+08\t3.98E+08\t9.88E+07\t3.20E+08\t3.04E+08\t1.70E+08\t4.41E+08\t2.60E+08\t1.14E+08\n",
"4.96E-08\t7.93E+10\t6.52E+10\t5.23E+10\t3.97E+10\t3.50E+10\t2.88E+10\t2.35E+10\t2.07E+10\t1.77E+10\t1.35E+10\t1.51E+10\t1.08E+10\t9.70E+09\t1.00E+10\t8.72E+09\t8.34E+09\t6.90E+09\t6.42E+09\t5.32E+09\t4.43E+09\t4.99E+09\t3.69E+09\t3.32E+09\t2.60E+09\t3.48E+09\t2.70E+09\t3.02E+09\t2.14E+09\t2.45E+09\t1.62E+09\t9.00E+08\t1.19E+09\t1.21E+09\t8.25E+08\t7.52E+08\t7.02E+08\t6.39E+08\t7.48E+08\t8.53E+08\t7.35E+08\t8.36E+08\t7.59E+08\t3.90E+08\t3.40E+08\t2.22E+08\t5.19E+08\t4.78E+08\t3.56E+08\t4.77E+08\n",
"5.14E-08\t8.69E+10\t7.17E+10\t5.89E+10\t4.52E+10\t4.03E+10\t3.26E+10\t2.86E+10\t2.35E+10\t2.12E+10\t1.64E+10\t1.67E+10\t1.54E+10\t1.23E+10\t1.24E+10\t1.23E+10\t9.53E+09\t9.11E+09\t8.05E+09\t7.20E+09\t5.48E+09\t5.28E+09\t4.61E+09\t4.30E+09\t4.18E+09\t3.83E+09\t3.14E+09\t3.45E+09\t2.82E+09\t2.98E+09\t2.10E+09\t1.87E+09\t1.75E+09\t1.29E+09\t1.70E+09\t9.82E+08\t1.36E+09\t1.32E+09\t1.14E+09\t1.26E+09\t1.03E+09\t8.08E+08\t7.55E+08\t9.22E+08\t8.92E+08\t9.37E+08\t7.31E+08\t4.17E+08\t4.18E+08\t4.70E+08\n",
"5.33E-08\t9.48E+10\t7.67E+10\t6.55E+10\t5.64E+10\t4.54E+10\t3.80E+10\t3.40E+10\t2.86E+10\t2.47E+10\t1.74E+10\t1.96E+10\t1.73E+10\t1.54E+10\t1.51E+10\t1.23E+10\t1.10E+10\t9.85E+09\t9.65E+09\t8.33E+09\t8.07E+09\t7.56E+09\t6.84E+09\t5.85E+09\t5.65E+09\t4.61E+09\t4.76E+09\t4.33E+09\t3.64E+09\t3.52E+09\t2.43E+09\t2.70E+09\t2.42E+09\t2.31E+09\t1.58E+09\t1.94E+09\t1.35E+09\t1.33E+09\t1.32E+09\t9.77E+08\t7.50E+08\t8.73E+08\t1.18E+09\t1.19E+09\t9.75E+08\t9.17E+08\t6.26E+08\t8.59E+08\t6.57E+08\t7.76E+08\n",
"5.52E-08\t1.05E+11\t8.91E+10\t7.02E+10\t6.23E+10\t5.09E+10\t4.29E+10\t3.75E+10\t3.32E+10\t2.90E+10\t2.35E+10\t2.21E+10\t1.94E+10\t1.78E+10\t1.72E+10\t1.51E+10\t1.40E+10\t1.34E+10\t1.11E+10\t9.86E+09\t9.37E+09\t8.46E+09\t8.33E+09\t7.18E+09\t6.75E+09\t5.16E+09\t5.37E+09\t5.30E+09\t4.75E+09\t4.00E+09\t3.76E+09\t2.62E+09\t2.93E+09\t2.47E+09\t2.53E+09\t2.58E+09\t2.06E+09\t1.75E+09\t1.68E+09\t1.45E+09\t1.24E+09\t1.41E+09\t1.28E+09\t1.23E+09\t1.20E+09\t1.42E+09\t1.17E+09\t9.55E+08\t7.83E+08\t1.13E+09\n",
"5.73E-08\t1.19E+11\t9.70E+10\t8.09E+10\t6.85E+10\t6.08E+10\t5.10E+10\t4.22E+10\t3.80E+10\t3.55E+10\t2.64E+10\t2.70E+10\t2.33E+10\t2.10E+10\t1.92E+10\t1.67E+10\t1.57E+10\t1.64E+10\t1.35E+10\t1.30E+10\t1.28E+10\t1.07E+10\t9.02E+09\t8.10E+09\t8.44E+09\t7.96E+09\t6.62E+09\t6.37E+09\t5.78E+09\t5.58E+09\t4.07E+09\t3.41E+09\t3.13E+09\t3.32E+09\t3.19E+09\t2.71E+09\t2.99E+09\t2.14E+09\t2.55E+09\t1.97E+09\t2.33E+09\t1.86E+09\t2.17E+09\t2.08E+09\t1.65E+09\t1.90E+09\t1.58E+09\t1.75E+09\t1.17E+09\t1.11E+09\n",
"5.94E-08\t1.26E+11\t1.07E+11\t8.95E+10\t7.92E+10\t6.51E+10\t5.74E+10\t4.89E+10\t4.27E+10\t3.97E+10\t3.06E+10\t3.10E+10\t2.71E+10\t2.28E+10\t2.20E+10\t2.14E+10\t1.91E+10\t1.63E+10\t1.47E+10\t1.40E+10\t1.35E+10\t1.25E+10\t1.26E+10\t1.05E+10\t1.04E+10\t9.26E+09\t7.79E+09\t7.46E+09\t7.26E+09\t6.43E+09\t5.90E+09\t4.25E+09\t4.71E+09\t4.50E+09\t3.26E+09\t3.76E+09\t3.48E+09\t2.90E+09\t3.35E+09\t2.39E+09\t2.64E+09\t2.83E+09\t2.34E+09\t2.16E+09\t2.20E+09\t1.73E+09\t1.92E+09\t2.25E+09\t1.92E+09\t1.37E+09\n",
"6.15E-08\t1.38E+11\t1.17E+11\t9.79E+10\t8.16E+10\t7.31E+10\t6.15E+10\t5.32E+10\t5.02E+10\t4.25E+10\t3.61E+10\t3.54E+10\t3.16E+10\t2.68E+10\t2.61E+10\t2.48E+10\t2.21E+10\t1.97E+10\t1.98E+10\t1.70E+10\t1.63E+10\t1.45E+10\t1.37E+10\t1.23E+10\t1.10E+10\t1.06E+10\t9.31E+09\t9.00E+09\t9.62E+09\t9.40E+09\t6.29E+09\t6.52E+09\t5.83E+09\t5.22E+09\t4.48E+09\t4.71E+09\t3.95E+09\t3.81E+09\t3.81E+09\t3.16E+09\t2.39E+09\t3.20E+09\t3.66E+09\t3.16E+09\t2.86E+09\t3.04E+09\t2.17E+09\t2.16E+09\t1.84E+09\t1.99E+09\n",
"6.38E-08\t1.52E+11\t1.28E+11\t1.04E+11\t9.06E+10\t7.71E+10\t7.10E+10\t6.18E+10\t5.33E+10\t5.01E+10\t4.15E+10\t3.95E+10\t3.70E+10\t3.34E+10\t2.99E+10\t2.69E+10\t2.41E+10\t2.27E+10\t2.12E+10\t1.81E+10\t1.84E+10\t1.54E+10\t1.59E+10\t1.65E+10\t1.40E+10\t1.28E+10\t1.19E+10\t1.11E+10\t1.03E+10\t9.97E+09\t6.76E+09\t7.54E+09\t6.83E+09\t6.35E+09\t6.14E+09\t5.73E+09\t5.14E+09\t5.08E+09\t4.90E+09\t3.75E+09\t4.39E+09\t4.19E+09\t3.50E+09\t3.87E+09\t3.06E+09\t3.51E+09\t3.60E+09\t3.10E+09\t2.45E+09\t2.48E+09\n",
"6.61E-08\t1.55E+11\t1.35E+11\t1.11E+11\t9.92E+10\t8.84E+10\t7.85E+10\t6.97E+10\t5.86E+10\t5.09E+10\t4.59E+10\t4.41E+10\t3.80E+10\t3.78E+10\t3.51E+10\t3.01E+10\t2.90E+10\t2.57E+10\t2.57E+10\t2.40E+10\t1.98E+10\t1.97E+10\t1.78E+10\t1.68E+10\t1.61E+10\t1.42E+10\t1.44E+10\t1.25E+10\t1.25E+10\t1.09E+10\t8.22E+09\t8.58E+09\t8.81E+09\t7.22E+09\t7.62E+09\t6.15E+09\t7.26E+09\t6.66E+09\t6.83E+09\t6.30E+09\t5.23E+09\t4.83E+09\t4.44E+09\t4.10E+09\t4.20E+09\t4.64E+09\t4.13E+09\t3.61E+09\t3.15E+09\t3.50E+09\n",
"6.85E-08\t1.66E+11\t1.46E+11\t1.28E+11\t1.09E+11\t9.68E+10\t8.52E+10\t7.52E+10\t6.83E+10\t6.04E+10\t5.31E+10\t5.11E+10\t4.48E+10\t4.25E+10\t4.06E+10\t3.59E+10\t3.37E+10\t3.08E+10\t2.95E+10\t2.63E+10\t2.35E+10\t2.24E+10\t2.18E+10\t1.96E+10\t1.72E+10\t1.54E+10\t1.56E+10\t1.51E+10\t1.49E+10\t1.45E+10\t1.15E+10\t1.03E+10\t9.68E+09\t8.72E+09\t9.50E+09\t7.98E+09\t8.06E+09\t7.31E+09\t6.68E+09\t7.05E+09\t6.33E+09\t6.63E+09\t5.60E+09\t5.33E+09\t5.29E+09\t4.95E+09\t4.60E+09\t5.08E+09\t4.55E+09\t4.85E+09\n",
"7.10E-08\t1.77E+11\t1.49E+11\t1.32E+11\t1.11E+11\t9.98E+10\t8.87E+10\t7.73E+10\t7.27E+10\t6.41E+10\t5.48E+10\t5.37E+10\t4.91E+10\t4.35E+10\t4.18E+10\t3.62E+10\t3.54E+10\t3.43E+10\t2.94E+10\t2.85E+10\t2.67E+10\t2.31E+10\t2.40E+10\t2.14E+10\t1.98E+10\t1.80E+10\t1.63E+10\t1.62E+10\t1.45E+10\t1.43E+10\t1.12E+10\t1.22E+10\t9.67E+09\t1.11E+10\t9.21E+09\t9.76E+09\t8.20E+09\t8.04E+09\t6.74E+09\t6.59E+09\t6.98E+09\t6.70E+09\t6.30E+09\t6.17E+09\t6.08E+09\t5.42E+09\t5.06E+09\t5.29E+09\t4.13E+09\t5.17E+09\n",
"7.37E-08\t1.80E+11\t1.55E+11\t1.38E+11\t1.22E+11\t1.07E+11\t9.28E+10\t8.03E+10\t7.54E+10\t6.55E+10\t5.85E+10\t5.67E+10\t5.14E+10\t4.69E+10\t4.48E+10\t3.88E+10\t3.73E+10\t3.48E+10\t3.06E+10\t2.96E+10\t2.74E+10\t2.50E+10\t2.43E+10\t2.42E+10\t2.19E+10\t2.12E+10\t1.97E+10\t1.86E+10\t1.57E+10\t1.52E+10\t1.16E+10\t1.31E+10\t1.20E+10\t1.04E+10\t9.84E+09\t8.52E+09\t9.27E+09\t9.42E+09\t8.63E+09\t7.10E+09\t6.96E+09\t7.07E+09\t6.66E+09\t6.98E+09\t6.45E+09\t6.59E+09\t4.60E+09\t5.43E+09\t5.13E+09\t5.07E+09\n",
"7.64E-08\t1.93E+11\t1.70E+11\t1.43E+11\t1.30E+11\t1.09E+11\t1.03E+11\t8.94E+10\t8.41E+10\t7.30E+10\t6.76E+10\t6.32E+10\t5.82E+10\t5.31E+10\t4.91E+10\t4.39E+10\t4.30E+10\t3.92E+10\t3.64E+10\t3.42E+10\t3.05E+10\t3.07E+10\t2.72E+10\t2.69E+10\t2.39E+10\t2.31E+10\t2.17E+10\t1.94E+10\t2.03E+10\t1.79E+10\t1.45E+10\t1.59E+10\t1.31E+10\t1.29E+10\t1.24E+10\t1.11E+10\t1.01E+10\t9.88E+09\t9.61E+09\t8.53E+09\t9.21E+09\t8.04E+09\t7.45E+09\t7.83E+09\t7.82E+09\t5.81E+09\t7.10E+09\t6.62E+09\t5.45E+09\t5.35E+09\n",
"7.91E-08\t2.04E+11\t1.81E+11\t1.54E+11\t1.34E+11\t1.23E+11\t1.12E+11\t9.97E+10\t8.81E+10\t8.14E+10\t7.19E+10\t6.71E+10\t6.27E+10\t5.92E+10\t5.46E+10\t4.92E+10\t4.65E+10\t4.40E+10\t3.82E+10\t3.90E+10\t3.55E+10\t3.23E+10\t3.04E+10\t2.85E+10\t2.76E+10\t2.60E+10\t2.36E+10\t2.28E+10\t1.95E+10\t2.04E+10\t1.83E+10\t1.65E+10\t1.51E+10\t1.39E+10\t1.42E+10\t1.21E+10\t1.28E+10\t1.19E+10\t1.19E+10\t1.00E+10\t1.11E+10\t9.97E+09\t8.76E+09\t8.23E+09\t8.29E+09\t8.38E+09\t8.83E+09\t7.51E+09\t7.10E+09\t7.80E+09\n",
"8.20E-08\t2.06E+11\t1.85E+11\t1.65E+11\t1.43E+11\t1.27E+11\t1.17E+11\t1.08E+11\t9.64E+10\t8.54E+10\t7.51E+10\t7.56E+10\t6.95E+10\t6.46E+10\t5.96E+10\t5.35E+10\t5.06E+10\t5.01E+10\t4.55E+10\t4.18E+10\t3.98E+10\t3.45E+10\t3.32E+10\t3.22E+10\t3.06E+10\t2.98E+10\t2.71E+10\t2.71E+10\t2.46E+10\t2.34E+10\t1.95E+10\t1.93E+10\t1.83E+10\t1.64E+10\t1.53E+10\t1.53E+10\t1.26E+10\t1.40E+10\t1.29E+10\t1.17E+10\t1.26E+10\t1.23E+10\t1.04E+10\t9.59E+09\t1.04E+10\t9.91E+09\t1.01E+10\t8.35E+09\t9.06E+09\t8.22E+09\n",
"8.51E-08\t2.23E+11\t2.00E+11\t1.75E+11\t1.55E+11\t1.34E+11\t1.22E+11\t1.14E+11\t1.02E+11\t9.23E+10\t8.31E+10\t8.15E+10\t7.18E+10\t6.75E+10\t6.70E+10\t6.07E+10\t5.68E+10\t5.30E+10\t5.06E+10\t4.67E+10\t4.24E+10\t3.98E+10\t3.84E+10\t3.33E+10\t3.44E+10\t3.41E+10\t3.08E+10\t2.87E+10\t2.73E+10\t2.46E+10\t2.22E+10\t2.17E+10\t2.00E+10\t2.05E+10\t1.89E+10\t1.81E+10\t1.50E+10\t1.56E+10\t1.53E+10\t1.53E+10\t1.36E+10\t1.31E+10\t1.27E+10\t1.14E+10\t1.18E+10\t1.14E+10\t1.13E+10\t9.17E+09\t9.69E+09\t9.47E+09\n",
"8.82E-08\t2.26E+11\t2.04E+11\t1.83E+11\t1.63E+11\t1.41E+11\t1.32E+11\t1.20E+11\t1.11E+11\t1.02E+11\t8.85E+10\t8.61E+10\t7.94E+10\t7.65E+10\t6.89E+10\t6.47E+10\t5.97E+10\t5.84E+10\t5.45E+10\t4.86E+10\t4.74E+10\t4.34E+10\t4.16E+10\t3.80E+10\t3.71E+10\t3.43E+10\t3.45E+10\t3.27E+10\t2.94E+10\t2.82E+10\t2.45E+10\t2.30E+10\t2.42E+10\t2.15E+10\t2.27E+10\t1.84E+10\t1.84E+10\t1.74E+10\t1.54E+10\t1.69E+10\t1.61E+10\t1.43E+10\t1.49E+10\t1.45E+10\t1.16E+10\t1.23E+10\t1.17E+10\t1.20E+10\t1.06E+10\t1.14E+10\n",
"9.14E-08\t2.33E+11\t2.05E+11\t1.88E+11\t1.71E+11\t1.52E+11\t1.40E+11\t1.27E+11\t1.10E+11\t1.08E+11\t9.41E+10\t9.10E+10\t8.01E+10\t7.99E+10\t7.62E+10\t7.14E+10\t6.40E+10\t6.29E+10\t6.09E+10\t5.31E+10\t5.23E+10\t4.56E+10\t4.57E+10\t4.26E+10\t4.18E+10\t3.92E+10\t3.69E+10\t3.40E+10\t3.28E+10\t3.12E+10\t2.82E+10\t2.77E+10\t2.58E+10\t2.45E+10\t2.37E+10\t2.16E+10\t2.05E+10\t1.94E+10\t1.94E+10\t1.93E+10\t1.82E+10\t1.68E+10\t1.60E+10\t1.48E+10\t1.53E+10\t1.43E+10\t1.45E+10\t1.28E+10\t1.24E+10\t1.17E+10\n",
"9.47E-08\t2.45E+11\t2.11E+11\t1.96E+11\t1.73E+11\t1.59E+11\t1.47E+11\t1.29E+11\t1.21E+11\t1.08E+11\t9.71E+10\t9.87E+10\t9.01E+10\t8.11E+10\t7.93E+10\t7.21E+10\t7.09E+10\t6.68E+10\t6.18E+10\t5.94E+10\t5.47E+10\t5.23E+10\t4.89E+10\t4.53E+10\t4.49E+10\t4.32E+10\t3.95E+10\t3.68E+10\t3.76E+10\t3.41E+10\t3.17E+10\t2.88E+10\t2.95E+10\t2.76E+10\t2.69E+10\t2.44E+10\t2.34E+10\t2.14E+10\t1.96E+10\t2.12E+10\t1.92E+10\t1.90E+10\t1.71E+10\t1.84E+10\t1.56E+10\t1.60E+10\t1.60E+10\t1.33E+10\t1.44E+10\t1.44E+10\n",
"9.82E-08\t2.48E+11\t2.28E+11\t2.04E+11\t1.85E+11\t1.68E+11\t1.56E+11\t1.43E+11\t1.29E+11\t1.20E+11\t1.10E+11\t1.04E+11\t9.70E+10\t8.96E+10\t8.27E+10\t8.05E+10\t7.53E+10\t7.39E+10\t6.76E+10\t6.13E+10\t5.93E+10\t5.71E+10\t5.47E+10\t5.34E+10\t4.87E+10\t4.53E+10\t4.42E+10\t3.89E+10\t4.02E+10\t3.79E+10\t3.48E+10\t3.35E+10\t3.00E+10\t3.11E+10\t2.87E+10\t2.65E+10\t2.55E+10\t2.30E+10\t2.27E+10\t2.06E+10\t2.20E+10\t2.24E+10\t1.95E+10\t1.95E+10\t1.93E+10\t1.73E+10\t1.78E+10\t1.74E+10\t1.64E+10\t1.44E+10\n",
"1.02E-07\t2.48E+11\t2.30E+11\t2.13E+11\t1.93E+11\t1.73E+11\t1.59E+11\t1.47E+11\t1.36E+11\t1.23E+11\t1.15E+11\t1.08E+11\t1.01E+11\t9.44E+10\t9.01E+10\t8.77E+10\t7.87E+10\t7.45E+10\t7.10E+10\t6.44E+10\t6.38E+10\t6.06E+10\t5.62E+10\t5.32E+10\t5.21E+10\t4.98E+10\t4.82E+10\t4.57E+10\t4.31E+10\t3.77E+10\t3.87E+10\t3.53E+10\t3.26E+10\t3.41E+10\t3.03E+10\t3.11E+10\t3.01E+10\t2.68E+10\t2.64E+10\t2.56E+10\t2.44E+10\t2.27E+10\t2.25E+10\t2.20E+10\t2.03E+10\t2.06E+10\t1.98E+10\t1.82E+10\t1.91E+10\t1.69E+10\n",
"1.06E-07\t2.49E+11\t2.32E+11\t2.03E+11\t1.92E+11\t1.78E+11\t1.66E+11\t1.57E+11\t1.38E+11\t1.35E+11\t1.16E+11\t1.10E+11\t1.06E+11\t9.88E+10\t9.52E+10\t9.51E+10\t8.29E+10\t7.99E+10\t7.70E+10\t6.97E+10\t6.75E+10\t6.24E+10\t5.78E+10\t5.83E+10\t5.53E+10\t5.17E+10\t4.94E+10\t4.76E+10\t4.80E+10\t4.38E+10\t4.01E+10\t3.64E+10\t3.77E+10\t3.51E+10\t3.38E+10\t3.04E+10\t3.16E+10\t3.00E+10\t2.75E+10\t2.76E+10\t2.48E+10\t2.46E+10\t2.58E+10\t2.42E+10\t2.18E+10\t2.34E+10\t2.27E+10\t2.00E+10\t1.91E+10\t1.76E+10\n",
"1.09E-07\t2.49E+11\t2.32E+11\t2.10E+11\t1.97E+11\t1.82E+11\t1.68E+11\t1.58E+11\t1.47E+11\t1.38E+11\t1.24E+11\t1.15E+11\t1.07E+11\t1.05E+11\t9.87E+10\t9.47E+10\t9.03E+10\t8.46E+10\t7.98E+10\t7.51E+10\t7.13E+10\t6.74E+10\t6.64E+10\t6.11E+10\t5.84E+10\t5.48E+10\t5.12E+10\t4.72E+10\t5.00E+10\t4.57E+10\t4.37E+10\t4.05E+10\t3.67E+10\t3.76E+10\t3.68E+10\t3.60E+10\t3.34E+10\t3.17E+10\t2.99E+10\t2.92E+10\t2.82E+10\t2.69E+10\t2.72E+10\t2.64E+10\t2.42E+10\t2.31E+10\t2.35E+10\t2.18E+10\t2.22E+10\t2.14E+10\n",
"1.13E-07\t2.51E+11\t2.32E+11\t2.18E+11\t1.99E+11\t1.89E+11\t1.68E+11\t1.67E+11\t1.49E+11\t1.43E+11\t1.27E+11\t1.20E+11\t1.19E+11\t1.10E+11\t1.03E+11\t9.75E+10\t9.13E+10\t8.86E+10\t8.11E+10\t7.90E+10\t7.63E+10\t7.22E+10\t7.16E+10\t6.37E+10\t6.36E+10\t5.95E+10\t5.77E+10\t5.35E+10\t5.17E+10\t4.97E+10\t4.69E+10\t4.14E+10\t4.60E+10\t4.03E+10\t3.94E+10\t3.76E+10\t3.75E+10\t3.60E+10\t3.38E+10\t3.19E+10\t3.03E+10\t3.03E+10\t3.09E+10\t2.86E+10\t2.81E+10\t2.76E+10\t2.47E+10\t2.47E+10\t2.63E+10\t2.38E+10\n",
"1.18E-07\t2.51E+11\t2.30E+11\t2.12E+11\t1.99E+11\t1.89E+11\t1.78E+11\t1.59E+11\t1.51E+11\t1.46E+11\t1.33E+11\t1.28E+11\t1.17E+11\t1.11E+11\t1.08E+11\t9.84E+10\t9.34E+10\t8.71E+10\t8.83E+10\t8.40E+10\t7.61E+10\t7.41E+10\t7.27E+10\t7.07E+10\t6.42E+10\t6.11E+10\t5.81E+10\t5.87E+10\t5.50E+10\t5.33E+10\t4.73E+10\t4.65E+10\t4.67E+10\t4.34E+10\t4.45E+10\t4.36E+10\t3.84E+10\t3.71E+10\t3.71E+10\t3.54E+10\t3.23E+10\t3.21E+10\t3.19E+10\t3.10E+10\t3.02E+10\t2.82E+10\t2.59E+10\t2.77E+10\t2.58E+10\t2.45E+10\n",
"1.22E-07\t2.43E+11\t2.28E+11\t2.18E+11\t2.07E+11\t1.90E+11\t1.81E+11\t1.69E+11\t1.58E+11\t1.51E+11\t1.35E+11\t1.28E+11\t1.20E+11\t1.13E+11\t1.08E+11\t1.04E+11\t9.66E+10\t9.72E+10\t9.01E+10\t7.99E+10\t7.94E+10\t8.12E+10\t7.62E+10\t7.19E+10\t6.97E+10\t6.58E+10\t6.29E+10\t5.86E+10\t5.50E+10\t5.66E+10\t5.39E+10\t5.03E+10\t5.00E+10\t4.48E+10\t4.57E+10\t4.37E+10\t4.26E+10\t4.00E+10\t3.70E+10\t3.57E+10\t3.47E+10\t3.65E+10\t3.40E+10\t3.49E+10\t3.14E+10\t3.19E+10\t2.93E+10\t3.05E+10\t2.78E+10\t2.85E+10\n",
"1.26E-07\t2.40E+11\t2.30E+11\t2.13E+11\t2.00E+11\t1.97E+11\t1.78E+11\t1.64E+11\t1.59E+11\t1.52E+11\t1.40E+11\t1.26E+11\t1.25E+11\t1.19E+11\t1.10E+11\t1.07E+11\t1.05E+11\t9.49E+10\t9.31E+10\t8.93E+10\t8.60E+10\t8.16E+10\t8.09E+10\t7.71E+10\t7.15E+10\t6.51E+10\t6.53E+10\t6.41E+10\t6.28E+10\t5.74E+10\t5.67E+10\t5.28E+10\t4.81E+10\t4.87E+10\t4.68E+10\t4.69E+10\t4.49E+10\t4.26E+10\t3.90E+10\t3.96E+10\t3.71E+10\t3.78E+10\t3.72E+10\t3.53E+10\t3.20E+10\t3.22E+10\t3.07E+10\t3.07E+10\t3.05E+10\t2.95E+10\n",
"1.31E-07\t2.38E+11\t2.24E+11\t2.09E+11\t1.99E+11\t1.91E+11\t1.77E+11\t1.67E+11\t1.55E+11\t1.48E+11\t1.41E+11\t1.34E+11\t1.27E+11\t1.18E+11\t1.16E+11\t1.11E+11\t1.06E+11\t9.71E+10\t9.57E+10\t9.06E+10\t8.83E+10\t8.60E+10\t8.08E+10\t7.65E+10\t7.38E+10\t7.21E+10\t6.92E+10\t6.50E+10\t6.25E+10\t6.07E+10\t5.66E+10\t5.95E+10\t5.43E+10\t5.25E+10\t5.01E+10\t4.83E+10\t4.57E+10\t4.45E+10\t4.46E+10\t4.31E+10\t3.99E+10\t4.12E+10\t3.68E+10\t3.85E+10\t3.70E+10\t3.55E+10\t3.51E+10\t3.08E+10\t3.20E+10\t3.06E+10\n",
"1.36E-07\t2.28E+11\t2.13E+11\t2.09E+11\t1.94E+11\t1.87E+11\t1.73E+11\t1.68E+11\t1.60E+11\t1.54E+11\t1.44E+11\t1.38E+11\t1.23E+11\t1.23E+11\t1.19E+11\t1.10E+11\t1.06E+11\t1.04E+11\t1.00E+11\t9.72E+10\t8.91E+10\t8.84E+10\t8.24E+10\t8.08E+10\t7.48E+10\t7.55E+10\t7.32E+10\t6.87E+10\t6.66E+10\t6.36E+10\t6.10E+10\t5.86E+10\t5.80E+10\t5.54E+10\t5.15E+10\t5.00E+10\t4.67E+10\t4.60E+10\t4.67E+10\t4.24E+10\t4.19E+10\t4.26E+10\t4.08E+10\t3.93E+10\t3.99E+10\t3.65E+10\t3.61E+10\t3.45E+10\t3.54E+10\t3.21E+10\n",
"1.41E-07\t2.17E+11\t2.07E+11\t2.02E+11\t1.96E+11\t1.84E+11\t1.75E+11\t1.68E+11\t1.62E+11\t1.54E+11\t1.44E+11\t1.37E+11\t1.31E+11\t1.22E+11\t1.20E+11\t1.15E+11\t1.07E+11\t1.04E+11\t1.00E+11\t9.42E+10\t9.00E+10\t9.04E+10\t8.51E+10\t8.06E+10\t8.15E+10\t7.47E+10\t7.24E+10\t6.66E+10\t6.88E+10\t6.67E+10\t6.52E+10\t6.02E+10\t5.98E+10\t5.58E+10\t5.37E+10\t5.66E+10\t5.17E+10\t5.04E+10\t4.95E+10\t4.70E+10\t4.43E+10\t4.45E+10\t4.35E+10\t4.31E+10\t4.06E+10\t4.04E+10\t3.66E+10\t3.92E+10\t3.51E+10\t3.51E+10\n",
"1.46E-07\t2.13E+11\t2.01E+11\t1.95E+11\t1.86E+11\t1.77E+11\t1.69E+11\t1.61E+11\t1.59E+11\t1.53E+11\t1.45E+11\t1.37E+11\t1.30E+11\t1.26E+11\t1.22E+11\t1.13E+11\t1.12E+11\t1.05E+11\t1.02E+11\t9.63E+10\t9.78E+10\t8.99E+10\t8.70E+10\t8.40E+10\t8.08E+10\t7.81E+10\t7.83E+10\t7.54E+10\t7.08E+10\t6.76E+10\t6.66E+10\t6.36E+10\t6.11E+10\t5.87E+10\t5.61E+10\t5.56E+10\t5.26E+10\t5.28E+10\t5.03E+10\t4.78E+10\t4.75E+10\t4.65E+10\t4.47E+10\t4.40E+10\t4.50E+10\t4.25E+10\t3.97E+10\t3.84E+10\t3.58E+10\t3.70E+10\n",
"1.51E-07\t2.04E+11\t1.92E+11\t1.90E+11\t1.84E+11\t1.73E+11\t1.61E+11\t1.60E+11\t1.56E+11\t1.44E+11\t1.41E+11\t1.33E+11\t1.29E+11\t1.26E+11\t1.21E+11\t1.13E+11\t1.11E+11\t1.06E+11\t9.93E+10\t9.81E+10\t9.36E+10\t9.15E+10\t8.96E+10\t8.74E+10\t8.28E+10\t8.14E+10\t7.72E+10\t7.33E+10\t6.82E+10\t6.90E+10\t6.84E+10\t6.64E+10\t6.59E+10\t6.25E+10\t5.98E+10\t5.68E+10\t5.40E+10\t5.45E+10\t5.10E+10\t5.05E+10\t4.72E+10\t4.87E+10\t4.59E+10\t4.56E+10\t4.44E+10\t4.20E+10\t4.27E+10\t4.21E+10\t3.90E+10\t3.78E+10\n",
"1.57E-07\t1.92E+11\t1.87E+11\t1.80E+11\t1.73E+11\t1.69E+11\t1.65E+11\t1.54E+11\t1.52E+11\t1.44E+11\t1.35E+11\t1.32E+11\t1.28E+11\t1.25E+11\t1.16E+11\t1.15E+11\t1.10E+11\t1.11E+11\t1.05E+11\t9.91E+10\t9.45E+10\t9.16E+10\t8.97E+10\t8.60E+10\t8.36E+10\t8.25E+10\t7.83E+10\t7.46E+10\t6.98E+10\t7.00E+10\t6.71E+10\t6.88E+10\t6.35E+10\t6.16E+10\t6.07E+10\t6.12E+10\t5.75E+10\t5.57E+10\t5.48E+10\t5.27E+10\t5.07E+10\t4.99E+10\t4.67E+10\t4.71E+10\t4.49E+10\t4.32E+10\t4.35E+10\t4.13E+10\t4.10E+10\t4.11E+10\n",
"1.63E-07\t1.80E+11\t1.74E+11\t1.77E+11\t1.68E+11\t1.64E+11\t1.58E+11\t1.49E+11\t1.48E+11\t1.45E+11\t1.32E+11\t1.32E+11\t1.24E+11\t1.19E+11\t1.19E+11\t1.12E+11\t1.13E+11\t1.10E+11\t1.01E+11\t9.85E+10\t9.80E+10\t9.10E+10\t9.05E+10\t8.69E+10\t8.49E+10\t8.05E+10\t7.88E+10\t7.56E+10\t7.30E+10\t7.25E+10\t6.96E+10\t6.93E+10\t6.69E+10\t6.57E+10\t6.16E+10\t5.92E+10\t5.85E+10\t5.77E+10\t5.49E+10\t5.33E+10\t5.12E+10\t5.18E+10\t4.97E+10\t4.93E+10\t4.78E+10\t4.73E+10\t4.28E+10\t4.25E+10\t4.13E+10\t4.25E+10\n",
"1.69E-07\t1.63E+11\t1.63E+11\t1.59E+11\t1.58E+11\t1.51E+11\t1.49E+11\t1.46E+11\t1.42E+11\t1.37E+11\t1.29E+11\t1.25E+11\t1.23E+11\t1.19E+11\t1.15E+11\t1.09E+11\t1.08E+11\t1.04E+11\t1.01E+11\t9.79E+10\t9.54E+10\t9.13E+10\t9.06E+10\t8.56E+10\t8.68E+10\t8.17E+10\t8.20E+10\t7.88E+10\t7.43E+10\t7.04E+10\t7.09E+10\t6.91E+10\t6.70E+10\t6.66E+10\t6.26E+10\t6.10E+10\t5.83E+10\t5.90E+10\t5.63E+10\t5.56E+10\t5.32E+10\t5.12E+10\t5.17E+10\t4.99E+10\t4.96E+10\t4.70E+10\t4.50E+10\t4.47E+10\t4.08E+10\t4.03E+10\n",
"1.75E-07\t1.55E+11\t1.50E+11\t1.50E+11\t1.48E+11\t1.46E+11\t1.40E+11\t1.37E+11\t1.36E+11\t1.26E+11\t1.28E+11\t1.25E+11\t1.18E+11\t1.17E+11\t1.16E+11\t1.13E+11\t1.07E+11\t1.03E+11\t1.02E+11\t9.63E+10\t9.38E+10\t9.07E+10\t8.85E+10\t8.62E+10\t8.12E+10\t8.06E+10\t7.86E+10\t8.12E+10\t7.39E+10\t7.25E+10\t7.11E+10\t7.18E+10\t6.64E+10\t6.58E+10\t6.30E+10\t6.25E+10\t6.17E+10\t5.72E+10\t5.55E+10\t5.78E+10\t5.62E+10\t5.23E+10\t5.06E+10\t5.08E+10\t4.98E+10\t4.61E+10\t4.75E+10\t4.43E+10\t4.20E+10\t4.28E+10\n",
"1.81E-07\t1.40E+11\t1.37E+11\t1.38E+11\t1.35E+11\t1.37E+11\t1.35E+11\t1.30E+11\t1.23E+11\t1.25E+11\t1.21E+11\t1.18E+11\t1.16E+11\t1.11E+11\t1.08E+11\t1.05E+11\t1.03E+11\t1.02E+11\t9.85E+10\t9.29E+10\t9.32E+10\t8.92E+10\t8.54E+10\t8.58E+10\t8.26E+10\t8.09E+10\t7.89E+10\t7.99E+10\t7.31E+10\t7.40E+10\t7.09E+10\t7.09E+10\t6.86E+10\t6.45E+10\t6.26E+10\t6.07E+10\t6.13E+10\t5.95E+10\t5.71E+10\t5.50E+10\t5.53E+10\t5.22E+10\t5.31E+10\t4.95E+10\t4.84E+10\t4.83E+10\t4.81E+10\t4.43E+10\t4.69E+10\t4.46E+10\n",
"1.88E-07\t1.25E+11\t1.26E+11\t1.26E+11\t1.24E+11\t1.26E+11\t1.22E+11\t1.22E+11\t1.19E+11\t1.21E+11\t1.15E+11\t1.14E+11\t1.12E+11\t1.05E+11\t1.07E+11\t1.04E+11\t1.02E+11\t9.85E+10\t9.40E+10\t9.22E+10\t9.13E+10\t8.45E+10\t8.44E+10\t8.32E+10\t8.09E+10\t7.97E+10\t7.55E+10\t7.59E+10\t7.56E+10\t7.58E+10\t6.80E+10\t6.93E+10\t6.71E+10\t6.76E+10\t6.23E+10\t6.25E+10\t6.17E+10\t6.10E+10\t5.93E+10\t5.62E+10\t5.60E+10\t5.42E+10\t5.36E+10\t5.11E+10\t5.13E+10\t4.76E+10\t4.71E+10\t4.48E+10\t4.50E+10\t4.74E+10\n",
"1.95E-07\t1.12E+11\t1.11E+11\t1.15E+11\t1.16E+11\t1.12E+11\t1.14E+11\t1.16E+11\t1.10E+11\t1.09E+11\t1.07E+11\t1.04E+11\t1.04E+11\t9.79E+10\t9.90E+10\t9.34E+10\t9.33E+10\t9.41E+10\t9.12E+10\t8.75E+10\t8.27E+10\t8.35E+10\t8.05E+10\t8.04E+10\t7.68E+10\t7.98E+10\t7.65E+10\t7.32E+10\t7.10E+10\t7.11E+10\t6.66E+10\t6.68E+10\t6.51E+10\t6.37E+10\t6.44E+10\t6.30E+10\t6.17E+10\t5.77E+10\t5.76E+10\t5.58E+10\t5.62E+10\t5.39E+10\t5.18E+10\t5.25E+10\t4.95E+10\t4.97E+10\t4.86E+10\t4.70E+10\t4.56E+10\t4.53E+10\n",
"2.02E-07\t9.77E+10\t1.01E+11\t1.02E+11\t1.02E+11\t1.03E+11\t1.04E+11\t1.00E+11\t1.04E+11\t1.00E+11\t9.95E+10\t9.61E+10\t9.73E+10\t9.40E+10\t9.20E+10\t8.84E+10\t9.08E+10\t8.94E+10\t9.00E+10\t8.29E+10\t8.26E+10\t7.90E+10\t8.06E+10\t7.88E+10\t7.54E+10\t7.29E+10\t7.30E+10\t7.10E+10\t6.89E+10\t6.81E+10\t6.80E+10\t6.54E+10\t6.45E+10\t6.10E+10\t6.09E+10\t6.07E+10\t5.70E+10\t5.80E+10\t5.64E+10\t5.26E+10\t5.59E+10\t5.27E+10\t5.25E+10\t4.93E+10\t4.94E+10\t4.84E+10\t4.52E+10\t4.67E+10\t4.59E+10\t4.32E+10\n",
"2.09E-07\t8.74E+10\t8.83E+10\t9.14E+10\t9.19E+10\t9.26E+10\t9.21E+10\t9.38E+10\t9.27E+10\t9.26E+10\t9.19E+10\t9.08E+10\t9.10E+10\t8.81E+10\t8.61E+10\t8.46E+10\t8.35E+10\t8.23E+10\t8.00E+10\t8.24E+10\t7.69E+10\t7.49E+10\t7.67E+10\t7.51E+10\t7.16E+10\t7.03E+10\t7.06E+10\t6.94E+10\t6.38E+10\t6.57E+10\t6.35E+10\t6.37E+10\t6.15E+10\t6.03E+10\t6.07E+10\t5.85E+10\t5.71E+10\t5.72E+10\t5.44E+10\t5.34E+10\t5.31E+10\t5.17E+10\t4.97E+10\t4.95E+10\t4.88E+10\t4.81E+10\t4.45E+10\t4.67E+10\t4.53E+10\t4.46E+10\n",
"2.17E-07\t7.35E+10\t7.61E+10\t8.09E+10\t7.94E+10\t8.34E+10\t8.41E+10\t8.45E+10\t8.35E+10\t8.05E+10\t8.05E+10\t8.09E+10\t8.16E+10\t8.12E+10\t8.06E+10\t7.94E+10\t7.77E+10\t7.56E+10\t7.36E+10\t7.26E+10\t7.42E+10\t7.22E+10\t7.03E+10\t7.01E+10\t6.76E+10\t6.73E+10\t6.61E+10\t6.60E+10\t6.34E+10\t6.44E+10\t6.40E+10\t6.02E+10\t5.97E+10\t5.97E+10\t5.70E+10\t5.36E+10\t5.67E+10\t5.19E+10\t5.20E+10\t5.11E+10\t5.11E+10\t4.98E+10\t4.84E+10\t4.90E+10\t4.65E+10\t4.52E+10\t4.48E+10\t4.57E+10\t4.27E+10\t4.40E+10\n",
"2.25E-07\t6.30E+10\t6.52E+10\t6.90E+10\t7.03E+10\t7.22E+10\t7.19E+10\t7.47E+10\t7.40E+10\t7.52E+10\t7.38E+10\t7.46E+10\t7.31E+10\t7.14E+10\t7.43E+10\t7.08E+10\t6.95E+10\t7.04E+10\t6.88E+10\t6.76E+10\t6.52E+10\t6.58E+10\t6.39E+10\t6.58E+10\t6.46E+10\t6.46E+10\t6.16E+10\t6.35E+10\t6.02E+10\t5.89E+10\t5.78E+10\t5.66E+10\t5.62E+10\t5.72E+10\t5.28E+10\t5.13E+10\t5.22E+10\t5.14E+10\t5.12E+10\t4.87E+10\t4.69E+10\t4.85E+10\t4.58E+10\t4.53E+10\t4.61E+10\t4.30E+10\t4.40E+10\t4.36E+10\t4.36E+10\t4.16E+10\n",
"2.33E-07\t5.14E+10\t5.74E+10\t5.74E+10\t5.96E+10\t6.14E+10\t6.49E+10\t6.36E+10\t6.33E+10\t6.54E+10\t6.63E+10\t6.37E+10\t6.58E+10\t6.52E+10\t6.35E+10\t6.47E+10\t6.53E+10\t6.58E+10\t6.31E+10\t6.18E+10\t6.36E+10\t6.21E+10\t6.04E+10\t6.09E+10\t5.85E+10\t5.72E+10\t5.77E+10\t5.60E+10\t5.56E+10\t5.53E+10\t5.45E+10\t5.33E+10\t5.28E+10\t5.08E+10\t5.04E+10\t5.10E+10\t5.08E+10\t4.95E+10\t4.91E+10\t4.75E+10\t4.72E+10\t4.73E+10\t4.49E+10\t4.41E+10\t4.29E+10\t4.18E+10\t4.11E+10\t3.93E+10\t3.93E+10\t3.95E+10\n",
"2.41E-07\t4.41E+10\t4.78E+10\t4.87E+10\t4.93E+10\t5.31E+10\t5.44E+10\t5.43E+10\t5.55E+10\t5.69E+10\t5.95E+10\t5.58E+10\t5.86E+10\t5.74E+10\t5.69E+10\t5.64E+10\t5.76E+10\t5.71E+10\t5.67E+10\t5.52E+10\t5.48E+10\t5.63E+10\t5.67E+10\t5.29E+10\t5.43E+10\t5.27E+10\t5.34E+10\t5.21E+10\t5.16E+10\t5.23E+10\t5.09E+10\t5.01E+10\t4.76E+10\t4.97E+10\t4.81E+10\t4.55E+10\t4.64E+10\t4.76E+10\t4.53E+10\t4.46E+10\t4.39E+10\t4.32E+10\t4.30E+10\t4.10E+10\t4.14E+10\t4.15E+10\t3.98E+10\t3.96E+10\t3.77E+10\t3.74E+10\n",
"2.50E-07\t3.71E+10\t3.87E+10\t4.12E+10\t4.11E+10\t4.41E+10\t4.65E+10\t4.83E+10\t4.66E+10\t4.96E+10\t5.11E+10\t4.94E+10\t4.92E+10\t4.94E+10\t4.96E+10\t5.11E+10\t5.01E+10\t4.78E+10\t5.08E+10\t5.14E+10\t4.96E+10\t4.98E+10\t4.99E+10\t4.95E+10\t4.84E+10\t4.90E+10\t4.61E+10\t4.58E+10\t4.61E+10\t4.62E+10\t4.79E+10\t4.62E+10\t4.58E+10\t4.33E+10\t4.46E+10\t4.22E+10\t4.21E+10\t4.11E+10\t3.98E+10\t4.10E+10\t4.07E+10\t3.85E+10\t3.81E+10\t3.79E+10\t3.72E+10\t3.69E+10\t3.69E+10\t3.81E+10\t3.65E+10\t3.60E+10\n",
"2.60E-07\t2.93E+10\t3.05E+10\t3.16E+10\t3.43E+10\t3.44E+10\t3.73E+10\t3.89E+10\t3.96E+10\t3.89E+10\t4.16E+10\t4.19E+10\t4.02E+10\t4.38E+10\t4.22E+10\t4.27E+10\t4.18E+10\t4.35E+10\t4.20E+10\t4.41E+10\t4.22E+10\t4.27E+10\t4.25E+10\t4.29E+10\t3.99E+10\t4.39E+10\t4.22E+10\t4.13E+10\t4.13E+10\t4.20E+10\t4.00E+10\t3.90E+10\t3.99E+10\t3.94E+10\t3.82E+10\t3.96E+10\t3.86E+10\t3.78E+10\t3.70E+10\t3.73E+10\t3.63E+10\t3.51E+10\t3.50E+10\t3.49E+10\t3.37E+10\t3.35E+10\t3.22E+10\t3.37E+10\t3.37E+10\t3.33E+10\n",
"2.69E-07\t2.34E+10\t2.43E+10\t2.68E+10\t2.60E+10\t3.06E+10\t2.99E+10\t3.16E+10\t3.31E+10\t3.45E+10\t3.45E+10\t3.49E+10\t3.45E+10\t3.51E+10\t3.65E+10\t3.54E+10\t3.66E+10\t3.75E+10\t3.70E+10\t3.78E+10\t3.66E+10\t3.70E+10\t3.79E+10\t3.72E+10\t3.82E+10\t3.78E+10\t3.72E+10\t3.67E+10\t3.68E+10\t3.69E+10\t3.48E+10\t3.51E+10\t3.44E+10\t3.43E+10\t3.67E+10\t3.58E+10\t3.47E+10\t3.51E+10\t3.30E+10\t3.40E+10\t3.38E+10\t3.18E+10\t3.15E+10\t3.11E+10\t3.08E+10\t3.04E+10\t3.04E+10\t3.07E+10\t3.13E+10\t2.84E+10\n",
"2.79E-07\t1.80E+10\t1.83E+10\t1.98E+10\t2.12E+10\t2.26E+10\t2.30E+10\t2.56E+10\t2.77E+10\t2.73E+10\t2.67E+10\t2.90E+10\t3.00E+10\t3.01E+10\t2.93E+10\t2.92E+10\t2.98E+10\t3.04E+10\t2.99E+10\t3.08E+10\t3.11E+10\t3.17E+10\t3.19E+10\t3.16E+10\t2.97E+10\t3.24E+10\t3.24E+10\t3.10E+10\t3.04E+10\t3.30E+10\t3.15E+10\t3.07E+10\t3.00E+10\t2.89E+10\t3.03E+10\t3.11E+10\t3.00E+10\t3.11E+10\t2.87E+10\t2.93E+10\t2.82E+10\t2.74E+10\t2.81E+10\t2.74E+10\t2.77E+10\t2.75E+10\t2.71E+10\t2.60E+10\t2.77E+10\t2.61E+10\n",
"2.89E-07\t1.33E+10\t1.48E+10\t1.71E+10\t1.63E+10\t1.79E+10\t1.94E+10\t1.97E+10\t2.09E+10\t2.23E+10\t2.24E+10\t2.21E+10\t2.36E+10\t2.50E+10\t2.34E+10\t2.46E+10\t2.39E+10\t2.48E+10\t2.47E+10\t2.62E+10\t2.68E+10\t2.53E+10\t2.70E+10\t2.68E+10\t2.70E+10\t2.65E+10\t2.68E+10\t2.53E+10\t2.53E+10\t2.68E+10\t2.58E+10\t2.66E+10\t2.72E+10\t2.62E+10\t2.70E+10\t2.56E+10\t2.60E+10\t2.77E+10\t2.59E+10\t2.57E+10\t2.50E+10\t2.39E+10\t2.39E+10\t2.52E+10\t2.46E+10\t2.30E+10\t2.47E+10\t2.46E+10\t2.31E+10\t2.33E+10\n",
"3.00E-07\t1.05E+10\t1.15E+10\t1.22E+10\t1.34E+10\t1.36E+10\t1.39E+10\t1.63E+10\t1.50E+10\t1.65E+10\t1.70E+10\t1.76E+10\t1.80E+10\t1.89E+10\t1.76E+10\t1.89E+10\t1.99E+10\t1.93E+10\t2.01E+10\t2.13E+10\t2.17E+10\t1.98E+10\t2.09E+10\t2.12E+10\t2.17E+10\t2.16E+10\t2.16E+10\t2.10E+10\t2.10E+10\t2.26E+10\t2.23E+10\t2.24E+10\t2.17E+10\t2.09E+10\t2.19E+10\t2.37E+10\t2.22E+10\t2.26E+10\t2.01E+10\t2.25E+10\t2.10E+10\t2.08E+10\t2.11E+10\t2.11E+10\t2.10E+10\t2.00E+10\t2.01E+10\t2.03E+10\t2.06E+10\t1.92E+10\n",
"3.11E-07\t7.95E+09\t9.03E+09\t9.09E+09\t9.62E+09\t1.04E+10\t1.07E+10\t1.13E+10\t1.21E+10\t1.28E+10\t1.33E+10\t1.38E+10\t1.42E+10\t1.57E+10\t1.46E+10\t1.51E+10\t1.49E+10\t1.52E+10\t1.59E+10\t1.64E+10\t1.70E+10\t1.63E+10\t1.74E+10\t1.81E+10\t1.66E+10\t1.80E+10\t1.82E+10\t1.83E+10\t1.81E+10\t1.70E+10\t1.86E+10\t1.82E+10\t1.88E+10\t1.91E+10\t1.81E+10\t1.81E+10\t1.84E+10\t1.86E+10\t1.81E+10\t1.77E+10\t1.85E+10\t1.71E+10\t1.85E+10\t1.74E+10\t1.74E+10\t1.69E+10\t1.60E+10\t1.73E+10\t1.78E+10\t1.68E+10\n",
"3.22E-07\t5.39E+09\t6.75E+09\t7.41E+09\t7.49E+09\t8.71E+09\t8.81E+09\t9.03E+09\t8.83E+09\t9.89E+09\t1.04E+10\t1.03E+10\t1.16E+10\t1.07E+10\t1.16E+10\t1.21E+10\t1.17E+10\t1.28E+10\t1.36E+10\t1.28E+10\t1.37E+10\t1.40E+10\t1.42E+10\t1.33E+10\t1.41E+10\t1.50E+10\t1.49E+10\t1.36E+10\t1.38E+10\t1.41E+10\t1.48E+10\t1.48E+10\t1.49E+10\t1.53E+10\t1.46E+10\t1.57E+10\t1.45E+10\t1.59E+10\t1.43E+10\t1.56E+10\t1.47E+10\t1.53E+10\t1.39E+10\t1.39E+10\t1.42E+10\t1.35E+10\t1.47E+10\t1.42E+10\t1.44E+10\t1.48E+10\n",
"3.34E-07\t4.33E+09\t5.00E+09\t5.53E+09\t5.65E+09\t6.14E+09\t6.81E+09\t7.10E+09\t6.65E+09\t7.80E+09\t7.96E+09\t8.33E+09\t8.64E+09\t8.36E+09\t7.87E+09\t9.38E+09\t8.90E+09\t9.33E+09\t9.55E+09\t1.07E+10\t1.11E+10\t1.18E+10\t1.11E+10\t1.06E+10\t1.13E+10\t1.15E+10\t1.16E+10\t1.08E+10\t1.10E+10\t1.16E+10\t1.20E+10\t1.15E+10\t1.20E+10\t1.16E+10\t1.24E+10\t1.24E+10\t1.24E+10\t1.22E+10\t1.27E+10\t1.22E+10\t1.22E+10\t1.24E+10\t1.17E+10\t1.16E+10\t1.24E+10\t1.10E+10\t1.22E+10\t1.22E+10\t1.16E+10\t1.12E+10\n",
"3.46E-07\t3.99E+09\t3.58E+09\t3.88E+09\t4.54E+09\t4.51E+09\t4.71E+09\t5.30E+09\t5.70E+09\t5.42E+09\t6.26E+09\t5.85E+09\t5.86E+09\t6.72E+09\t6.34E+09\t6.63E+09\t6.74E+09\t7.27E+09\t7.29E+09\t7.74E+09\t7.92E+09\t8.54E+09\t8.42E+09\t8.75E+09\t8.43E+09\t9.65E+09\t8.57E+09\t8.20E+09\t8.83E+09\t9.52E+09\t9.11E+09\t9.37E+09\t9.26E+09\t9.39E+09\t9.02E+09\t1.04E+10\t9.51E+09\t9.72E+09\t1.04E+10\t9.79E+09\t1.02E+10\t1.03E+10\t9.10E+09\t9.64E+09\t9.82E+09\t9.86E+09\t1.06E+10\t9.66E+09\t9.75E+09\t9.56E+09\n",
"3.59E-07\t2.26E+09\t2.90E+09\t2.99E+09\t3.82E+09\t3.96E+09\t3.75E+09\t3.73E+09\t3.75E+09\t4.10E+09\t4.25E+09\t4.35E+09\t4.03E+09\t5.05E+09\t4.79E+09\t4.71E+09\t5.09E+09\t5.51E+09\t6.23E+09\t6.37E+09\t6.79E+09\t6.26E+09\t6.62E+09\t6.78E+09\t7.55E+09\t6.51E+09\t7.00E+09\t6.13E+09\t6.29E+09\t6.26E+09\t6.99E+09\t7.49E+09\t7.54E+09\t7.46E+09\t7.99E+09\t7.31E+09\t7.61E+09\t7.26E+09\t8.19E+09\t8.28E+09\t7.97E+09\t8.30E+09\t7.24E+09\t7.60E+09\t7.71E+09\t8.17E+09\t8.06E+09\t8.20E+09\t7.75E+09\t7.67E+09\n",
"3.72E-07\t2.30E+09\t1.97E+09\t2.34E+09\t2.40E+09\t2.67E+09\t2.60E+09\t3.07E+09\t3.14E+09\t3.16E+09\t2.89E+09\t3.64E+09\t3.71E+09\t3.85E+09\t2.86E+09\t3.64E+09\t4.01E+09\t3.82E+09\t4.69E+09\t4.44E+09\t4.70E+09\t4.74E+09\t4.91E+09\t4.61E+09\t5.36E+09\t5.30E+09\t5.35E+09\t4.49E+09\t5.27E+09\t5.25E+09\t5.56E+09\t6.17E+09\t5.49E+09\t5.66E+09\t6.12E+09\t6.26E+09\t5.81E+09\t5.93E+09\t6.75E+09\t6.33E+09\t6.46E+09\t6.68E+09\t5.38E+09\t6.40E+09\t5.82E+09\t7.08E+09\t5.77E+09\t6.12E+09\t6.58E+09\t6.18E+09\n",
"3.85E-07\t1.95E+09\t1.44E+09\t1.96E+09\t2.07E+09\t2.23E+09\t2.40E+09\t2.03E+09\t2.68E+09\t2.73E+09\t2.76E+09\t2.95E+09\t2.53E+09\t2.64E+09\t2.77E+09\t2.35E+09\t2.87E+09\t3.25E+09\t3.17E+09\t3.39E+09\t3.56E+09\t3.41E+09\t3.85E+09\t3.71E+09\t3.95E+09\t3.97E+09\t3.97E+09\t3.67E+09\t3.32E+09\t3.74E+09\t4.76E+09\t4.43E+09\t4.12E+09\t4.39E+09\t4.98E+09\t4.67E+09\t5.10E+09\t4.48E+09\t4.08E+09\t4.70E+09\t4.59E+09\t4.35E+09\t4.49E+09\t4.71E+09\t5.20E+09\t5.01E+09\t4.62E+09\t5.27E+09\t5.21E+09\t4.89E+09\n",
"4.00E-07\t1.21E+09\t1.30E+09\t1.41E+09\t1.69E+09\t1.78E+09\t1.93E+09\t1.90E+09\t1.56E+09\t1.52E+09\t2.00E+09\t2.20E+09\t1.85E+09\t2.15E+09\t2.31E+09\t2.12E+09\t2.25E+09\t2.53E+09\t2.77E+09\t2.75E+09\t3.09E+09\t2.81E+09\t3.22E+09\t3.21E+09\t3.05E+09\t2.94E+09\t3.08E+09\t2.49E+09\t2.64E+09\t3.57E+09\t3.17E+09\t3.54E+09\t2.82E+09\t3.46E+09\t3.31E+09\t3.83E+09\t3.74E+09\t3.66E+09\t3.54E+09\t4.27E+09\t3.53E+09\t3.24E+09\t3.67E+09\t3.64E+09\t3.72E+09\t3.57E+09\t4.19E+09\t3.93E+09\t3.93E+09\t4.13E+09\n",
"4.14E-07\t8.33E+08\t1.12E+09\t1.16E+09\t1.53E+09\t1.14E+09\t1.52E+09\t1.06E+09\t1.45E+09\t1.80E+09\t1.47E+09\t1.73E+09\t1.40E+09\t2.05E+09\t2.08E+09\t1.34E+09\t1.82E+09\t1.71E+09\t1.94E+09\t1.72E+09\t1.95E+09\t1.64E+09\t2.32E+09\t2.30E+09\t2.18E+09\t2.46E+09\t2.40E+09\t2.36E+09\t2.08E+09\t2.17E+09\t2.48E+09\t2.47E+09\t2.52E+09\t2.96E+09\t2.81E+09\t2.85E+09\t3.20E+09\t2.73E+09\t3.32E+09\t3.03E+09\t2.86E+09\t2.52E+09\t3.05E+09\t2.78E+09\t3.26E+09\t2.85E+09\t2.59E+09\t2.75E+09\t2.71E+09\t2.49E+09\n",
"4.29E-07\t7.84E+08\t9.46E+08\t1.08E+09\t9.76E+08\t1.22E+09\t1.26E+09\t1.14E+09\t1.28E+09\t1.09E+09\t9.57E+08\t1.02E+09\t1.32E+09\t1.48E+09\t1.36E+09\t1.30E+09\t1.25E+09\t1.42E+09\t1.43E+09\t1.91E+09\t1.42E+09\t1.65E+09\t1.60E+09\t1.75E+09\t1.68E+09\t1.63E+09\t1.92E+09\t1.45E+09\t1.65E+09\t2.00E+09\t1.60E+09\t1.50E+09\t1.89E+09\t2.26E+09\t1.81E+09\t1.93E+09\t2.00E+09\t2.30E+09\t2.37E+09\t2.08E+09\t2.56E+09\t1.82E+09\t2.12E+09\t1.83E+09\t2.34E+09\t2.32E+09\t2.40E+09\t2.04E+09\t2.16E+09\t2.17E+09\n",
"4.45E-07\t8.06E+08\t6.92E+08\t6.18E+08\t9.80E+08\t8.33E+08\t9.62E+08\t8.34E+08\t1.19E+09\t9.21E+08\t1.18E+09\t9.16E+08\t7.98E+08\t7.63E+08\t7.01E+08\t1.34E+09\t8.86E+08\t1.14E+09\t1.43E+09\t1.05E+09\t1.08E+09\t1.23E+09\t1.23E+09\t1.11E+09\t1.40E+09\t1.46E+09\t1.22E+09\t1.34E+09\t1.28E+09\t1.22E+09\t1.40E+09\t1.34E+09\t1.61E+09\t1.49E+09\t1.78E+09\t1.61E+09\t1.51E+09\t1.58E+09\t1.30E+09\t1.92E+09\t1.67E+09\t1.42E+09\t1.60E+09\t1.44E+09\t1.49E+09\t1.55E+09\t1.97E+09\t1.90E+09\t1.61E+09\t1.90E+09\n",
"4.61E-07\t6.21E+08\t6.92E+08\t6.96E+08\t6.50E+08\t1.03E+09\t6.26E+08\t6.22E+08\t5.76E+08\t9.25E+08\t6.72E+08\t8.40E+08\t5.74E+08\t7.83E+08\t4.92E+08\t9.04E+08\t7.01E+08\t8.97E+08\t1.15E+09\t9.94E+08\t8.90E+08\t8.02E+08\t1.16E+09\t1.23E+09\t8.92E+08\t8.54E+08\t9.08E+08\t7.41E+08\t8.77E+08\t6.98E+08\t9.54E+08\t1.04E+09\t1.10E+09\t1.39E+09\t1.05E+09\t1.20E+09\t1.11E+09\t1.09E+09\t1.07E+09\t1.32E+09\t1.41E+09\t8.83E+08\t1.17E+09\t1.02E+09\t1.51E+09\t1.21E+09\t1.43E+09\t1.36E+09\t1.33E+09\t1.65E+09\n",
"4.78E-07\t7.15E+08\t6.17E+08\t6.18E+08\t5.33E+08\t6.61E+08\t6.49E+08\t8.38E+08\t6.37E+08\t8.92E+08\t5.45E+08\t6.24E+08\t4.05E+08\t5.44E+08\t6.98E+08\t7.40E+08\t9.07E+08\t6.47E+08\t8.06E+08\t7.65E+08\t6.61E+08\t7.77E+08\t9.33E+08\t9.23E+08\t7.33E+08\t7.90E+08\t9.01E+08\t5.27E+08\t7.31E+08\t8.62E+08\t7.11E+08\t1.06E+09\t1.03E+09\t7.43E+08\t1.12E+09\t8.95E+08\t1.09E+09\t9.51E+08\t9.67E+08\t1.24E+09\t9.54E+08\t8.06E+08\t9.55E+08\t9.87E+08\t9.27E+08\t8.56E+08\t1.11E+09\t1.03E+09\t1.16E+09\t9.16E+08\n",
"4.96E-07\t7.97E+08\t4.93E+08\t6.66E+08\t7.02E+08\t6.46E+08\t7.61E+08\t6.27E+08\t4.28E+08\t5.14E+08\t6.30E+08\t5.84E+08\t3.41E+08\t4.37E+08\t5.60E+08\t6.44E+08\t4.49E+08\t7.40E+08\t6.58E+08\t6.92E+08\t7.28E+08\t6.71E+08\t7.32E+08\t6.11E+08\t6.75E+08\t3.87E+08\t5.10E+08\t4.27E+08\t5.24E+08\t5.67E+08\t4.37E+08\t8.00E+08\t9.94E+08\t8.17E+08\t8.12E+08\t9.85E+08\t7.90E+08\t8.36E+08\t5.48E+08\t7.43E+08\t1.10E+09\t5.81E+08\t6.46E+08\t7.62E+08\t6.03E+08\t8.58E+08\t8.74E+08\t9.54E+08\t7.57E+08\t7.33E+08\n",
"5.14E-07\t4.27E+08\t6.43E+08\t5.35E+08\t5.64E+08\t6.65E+08\t4.71E+08\t5.15E+08\t6.67E+08\t4.24E+08\t6.01E+08\t3.46E+08\t3.34E+08\t2.06E+08\t6.62E+08\t4.51E+08\t5.37E+08\t5.35E+08\t6.77E+08\t5.54E+08\t4.69E+08\t5.25E+08\t5.59E+08\t4.24E+08\t6.43E+08\t4.88E+08\t4.02E+08\t6.23E+08\t5.79E+08\t7.53E+08\t6.94E+08\t7.53E+08\t6.53E+08\t5.79E+08\t6.99E+08\t6.11E+08\t7.63E+08\t4.09E+08\t5.83E+08\t4.64E+08\t4.86E+08\t5.99E+08\t6.40E+08\t4.51E+08\t5.79E+08\t7.09E+08\t6.87E+08\t5.62E+08\t3.83E+08\t6.65E+08\n",
"5.33E-07\t5.55E+08\t4.03E+08\t4.73E+08\t4.55E+08\t4.90E+08\t3.78E+08\t4.13E+08\t2.98E+08\t5.13E+08\t2.49E+08\t3.18E+08\t2.24E+08\t2.71E+08\t3.41E+08\t4.26E+08\t4.88E+08\t5.20E+08\t4.80E+08\t4.23E+08\t3.01E+08\t5.78E+08\t4.93E+08\t3.98E+08\t3.36E+08\t4.23E+08\t4.68E+08\t4.30E+08\t5.65E+08\t2.91E+08\t5.60E+08\t4.06E+08\t4.55E+08\t5.35E+08\t5.00E+08\t4.90E+08\t4.21E+08\t5.13E+08\t4.43E+08\t3.63E+08\t4.98E+08\t3.14E+08\t5.07E+08\t6.65E+08\t3.38E+08\t4.48E+08\t4.90E+08\t7.37E+08\t5.58E+08\t3.56E+08\n",
"5.52E-07\t5.46E+08\t3.00E+08\t3.14E+08\t2.91E+08\t3.00E+08\t3.45E+08\t3.18E+08\t4.50E+08\t2.32E+08\t3.59E+08\t3.36E+08\t2.13E+08\t2.45E+08\t3.41E+08\t4.82E+08\t2.82E+08\t4.32E+08\t4.68E+08\t4.14E+08\t4.77E+08\t3.36E+08\t3.91E+08\t4.59E+08\t5.68E+08\t2.77E+08\t4.36E+08\t3.22E+08\t4.77E+08\t3.41E+08\t5.00E+08\t3.82E+08\t6.50E+08\t3.18E+08\t3.04E+08\t4.09E+08\t3.91E+08\t3.63E+08\t3.68E+08\t2.91E+08\t2.50E+08\t3.41E+08\t4.90E+08\t3.86E+08\t2.23E+08\t5.00E+08\t4.32E+08\t3.59E+08\t4.68E+08\t6.00E+08\n",
"5.73E-07\t3.99E+08\t4.60E+08\t1.86E+08\t2.99E+08\t3.91E+08\t4.37E+08\t3.41E+08\t2.99E+08\t2.28E+08\t3.02E+08\t2.51E+08\t2.70E+08\t1.63E+08\t2.55E+08\t3.40E+08\t3.63E+08\t3.44E+08\t2.78E+08\t3.48E+08\t4.09E+08\t2.09E+08\t4.09E+08\t2.53E+08\t3.72E+08\t4.14E+08\t2.53E+08\t3.86E+08\t2.74E+08\t2.33E+08\t2.74E+08\t3.75E+08\t3.04E+08\t2.33E+08\t3.49E+08\t4.14E+08\t3.21E+08\t3.67E+08\t3.40E+08\t3.74E+08\t2.51E+08\t2.34E+08\t3.03E+08\t3.29E+08\t2.33E+08\t4.43E+08\t3.95E+08\t3.04E+08\t5.69E+08\t4.33E+08\n",
"5.94E-07\t3.27E+08\t3.90E+08\t2.32E+08\t2.78E+08\t2.36E+08\t2.83E+08\t3.92E+08\t3.96E+08\t3.07E+08\t3.03E+08\t2.32E+08\t1.38E+08\t1.89E+08\t3.00E+08\t3.70E+08\t3.54E+08\t2.32E+08\t4.69E+08\t3.47E+08\t3.03E+08\t3.07E+08\t3.98E+08\t1.13E+08\t2.12E+08\t2.32E+08\t1.64E+08\t3.31E+08\t2.29E+08\t3.54E+08\t3.07E+08\t3.27E+08\t2.11E+08\t2.12E+08\t3.07E+08\t1.42E+08\t3.07E+08\t2.36E+08\t2.60E+08\t1.60E+08\t9.44E+07\t1.87E+08\t2.98E+08\t1.38E+08\t3.54E+08\t3.25E+08\t2.56E+08\t4.00E+08\t2.51E+08\t1.89E+08\n",
"6.15E-07\t3.60E+08\t3.35E+08\t2.68E+08\t2.92E+08\t1.44E+08\t9.64E+07\t2.71E+08\t2.32E+08\t2.16E+08\t1.96E+08\t5.17E+07\t1.72E+08\t4.80E+07\t3.39E+08\t2.23E+08\t2.40E+08\t3.16E+08\t2.64E+08\t1.48E+08\t2.56E+08\t2.88E+08\t2.44E+08\t2.80E+08\t4.44E+08\t3.04E+08\t1.44E+08\t1.80E+08\t5.35E+08\t2.76E+08\t4.20E+08\t1.68E+08\t1.68E+08\t1.68E+08\t4.20E+08\t3.36E+08\t2.40E+08\t2.64E+08\t1.80E+08\t1.96E+08\t1.44E+08\t2.64E+08\t2.56E+08\t1.68E+08\t3.12E+08\t2.44E+08\t3.76E+08\t2.64E+08\t3.92E+08\t2.88E+08\n",
"6.38E-07\t3.43E+08\t2.82E+08\t3.19E+08\t2.95E+08\t1.72E+08\t2.45E+08\t1.47E+08\t1.84E+08\t1.47E+08\t2.94E+08\t2.46E+08\t1.96E+08\t7.37E+07\t1.96E+08\t1.96E+08\t2.95E+08\t2.21E+08\t2.95E+08\t3.21E+08\t2.33E+08\t2.95E+08\t7.74E+07\t3.33E+08\t1.12E+08\t1.12E+08\t1.72E+08\t2.82E+08\t2.95E+08\t2.12E+08\t2.82E+08\t1.96E+08\t2.46E+08\t2.74E+08\t1.59E+08\t3.68E+08\t2.47E+08\t1.72E+08\t2.08E+08\t2.21E+08\t7.37E+07\t2.46E+08\t1.84E+08\t2.72E+08\t2.21E+08\t2.21E+08\t2.33E+08\t1.76E+08\t2.70E+08\t7.37E+07\n",
"6.61E-07\t2.75E+08\t3.02E+08\t2.25E+08\t2.25E+08\t2.50E+08\t1.75E+08\t1.50E+08\t1.02E+08\t1.50E+08\t1.25E+08\t1.25E+08\t1.50E+08\t1.36E+07\t2.00E+08\t3.27E+08\t2.38E+08\t2.00E+08\t1.75E+08\t2.36E+08\t1.25E+08\t2.00E+08\t1.71E+08\t1.61E+08\t1.73E+08\t1.23E+08\t1.50E+08\t2.38E+08\t1.88E+08\t2.98E+08\t1.88E+08\t1.75E+08\t2.38E+08\t2.96E+08\t2.99E+08\t9.98E+07\t2.48E+08\t1.38E+08\t1.50E+08\t1.50E+08\t1.13E+08\t1.63E+08\t2.75E+08\t9.79E+07\t1.75E+08\t2.13E+08\t4.99E+07\t2.21E+08\t7.49E+07\t2.13E+08\n",
"6.85E-07\t3.55E+08\t3.76E+08\t2.26E+08\t1.78E+08\t1.01E+08\t1.01E+08\t2.52E+08\t2.24E+08\t7.61E+07\t1.27E+08\t2.03E+08\t2.54E+07\t1.13E+08\t2.79E+08\t3.51E+08\t1.89E+08\t1.27E+08\t1.78E+08\t6.22E+07\t2.26E+08\t1.76E+08\t5.07E+07\t1.89E+08\t1.27E+08\t1.76E+08\t7.61E+07\t2.14E+08\t2.63E+08\t2.77E+08\t1.13E+08\t2.77E+08\t1.38E+08\t2.26E+08\t7.61E+07\t1.78E+08\t2.01E+08\t1.38E+08\t1.01E+08\t5.07E+07\t6.22E+07\t1.36E+08\t2.54E+07\t1.27E+08\t2.03E+08\t2.12E+08\t1.01E+08\t1.50E+08\t1.52E+08\t8.76E+07\n",
"7.11E-07\t3.10E+08\t1.57E+08\t2.60E+08\t2.58E+07\t2.58E+08\t2.12E+08\t2.34E+08\t1.83E+08\t2.63E+08\t1.55E+08\t2.38E+08\t0\t8.27E+07\t1.09E+08\t3.17E+08\t2.06E+08\t2.17E+08\t2.58E+08\t7.74E+07\t5.90E+07\t2.08E+08\t2.32E+08\t2.06E+08\t1.55E+08\t1.31E+08\t5.16E+07\t7.74E+07\t7.94E+07\t5.16E+07\t1.55E+08\t1.05E+08\t2.58E+08\t1.57E+08\t1.29E+08\t1.81E+08\t7.94E+07\t1.09E+08\t1.81E+08\t5.16E+07\t1.55E+08\t1.57E+08\t2.06E+08\t3.10E+08\t1.55E+08\t1.83E+08\t1.81E+08\t1.11E+08\t1.03E+08\t1.60E+08\n",
"7.37E-07\t2.86E+08\t1.04E+08\t1.41E+08\t5.20E+07\t5.20E+07\t1.51E+08\t7.80E+07\t1.56E+08\t9.86E+07\t2.45E+08\t4.66E+07\t1.56E+08\t7.26E+07\t1.77E+08\t2.13E+08\t2.45E+08\t2.86E+08\t1.56E+08\t1.41E+08\t8.32E+07\t3.54E+08\t2.08E+08\t1.30E+08\t5.20E+07\t5.20E+07\t1.67E+08\t1.56E+08\t1.30E+08\t1.04E+08\t1.82E+08\t1.41E+08\t7.80E+07\t0\t1.56E+08\t1.04E+08\t7.80E+07\t7.26E+07\t5.20E+07\t1.30E+08\t5.20E+07\t8.86E+07\t1.82E+08\t1.04E+08\t2.86E+08\t1.30E+08\t2.34E+08\t2.39E+08\t8.86E+07\t1.35E+08\n",
"7.64E-07\t1.85E+08\t1.85E+08\t9.51E+07\t1.06E+08\t2.52E+08\t1.59E+08\t7.94E+07\t1.72E+08\t1.06E+08\t6.87E+07\t7.94E+07\t7.94E+07\t5.29E+07\t1.32E+08\t1.62E+08\t2.01E+08\t1.75E+08\t1.72E+08\t2.27E+08\t3.21E+08\t8.94E+07\t1.32E+08\t5.29E+07\t1.06E+08\t5.29E+07\t9.51E+07\t1.85E+08\t1.32E+08\t1.20E+08\t2.38E+08\t1.62E+08\t1.06E+08\t1.46E+08\t1.99E+08\t1.59E+08\t1.06E+08\t1.06E+08\t5.29E+07\t1.06E+08\t0\t1.22E+08\t7.94E+07\t1.59E+08\t1.32E+08\t5.29E+07\t1.59E+08\t1.22E+08\t1.75E+08\t5.59E+07\n",
"7.92E-07\t1.35E+08\t1.88E+08\t1.62E+08\t2.03E+08\t2.01E+08\t8.08E+07\t2.15E+08\t6.68E+07\t6.86E+07\t2.69E+07\t5.39E+07\t1.08E+08\t2.69E+07\t1.62E+08\t3.24E+08\t8.08E+07\t8.08E+07\t1.89E+08\t1.62E+08\t9.38E+07\t8.08E+07\t1.88E+08\t4.17E+07\t0\t9.56E+07\t5.39E+07\t8.08E+07\t1.76E+08\t1.75E+08\t1.22E+08\t2.01E+08\t1.76E+08\t6.68E+07\t6.68E+07\t5.39E+07\t2.69E+07\t1.62E+08\t8.08E+07\t0\t2.69E+07\t1.08E+08\t1.35E+08\t8.08E+07\t9.56E+07\t1.88E+08\t1.35E+08\t8.08E+07\t5.39E+07\t1.30E+07\n",
"8.21E-07\t2.47E+08\t2.74E+08\t1.51E+08\t1.36E+08\t2.74E+07\t1.92E+08\t1.37E+08\t1.92E+08\t1.49E+08\t5.48E+07\t5.48E+07\t5.48E+07\t0\t1.92E+08\t2.04E+08\t3.02E+08\t1.37E+08\t9.46E+07\t1.78E+08\t1.10E+08\t8.22E+07\t2.74E+07\t6.72E+07\t8.22E+07\t9.46E+07\t5.48E+07\t8.22E+07\t1.49E+08\t2.47E+08\t6.72E+07\t1.10E+08\t1.49E+08\t8.22E+07\t9.59E+07\t1.10E+08\t5.48E+07\t2.74E+07\t1.10E+08\t0\t2.74E+07\t8.22E+07\t1.64E+08\t5.48E+07\t1.22E+08\t5.48E+07\t5.48E+07\t1.23E+08\t1.64E+08\t1.64E+08"
]
},
{
"cell_type": "code",
"execution_count": 19,
"id": "8e1b89e0",
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "9c26425d47b04adcad673cf0ae24cc09",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(HTML(value=\"./tmptxihlvmw.pdf
\"), HTML(value…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.xlabel(\"Time (s)\")\n",
"plt.plot(time, num_conc, \"g\", label=\"num conc\")\n",
"plt.xticks(np.linspace(0, time[-1], 5))\n",
"plt.xlim([time[0],time[-1]])\n",
"plt.ylabel(r\"Number concentration ($\\#$ m$^{-3}$)\", color='g')\n",
"barrel = np.loadtxt('ref_aero_time.txt')\n",
"plt.plot(barrel[:,0],barrel[:,1],'k')\n",
"show_plot()"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "faa1de28",
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "77a425a42eb147d09bf9f2ce1e10495a",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(HTML(value=\"./tmpa4ve8hoy.pdf
\"), HTML(value…"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plot_times = [0,1,10,30]\n",
"barrel = np.loadtxt('ref_aero_size_num.txt')\n",
"diams_barrel = barrel[:,0]\n",
"\n",
"for plot_time in plot_times:\n",
" p = plt.plot(diam_grid.centers, np.array(dists[plot_time]) * np.log(10),\n",
" label=f'$t = {dist_times[plot_time] / 60:0.0f}$ m', lw=1)\n",
" plt.plot(diams_barrel, barrel[:,plot_time+1],'--',c=p[0].get_color(),lw=1)\n",
"plt.xscale(\"log\")\n",
"plt.xlabel(\"Dry diameter (m)\")\n",
"plt.ylabel(r\"Number concentration $n(D)$ ($\\#$ m$^{-3}$)\")\n",
"plt.ylim(bottom=0)\n",
"plt.legend()\n",
"plt.xlim([diam_grid.edges[0],diam_grid.edges[-1]])\n",
"show_plot()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "44475ce8",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.13"
}
},
"nbformat": 4,
"nbformat_minor": 5
}