diff --git a/.github/workflows/Downgrade.yml b/.github/workflows/Downgrade.yml index 0d6c0b132..a3c64a7a0 100644 --- a/.github/workflows/Downgrade.yml +++ b/.github/workflows/Downgrade.yml @@ -49,7 +49,7 @@ jobs: - uses: julia-actions/cache@v2 - uses: julia-actions/julia-downgrade-compat@v1 with: - skip: LinearAlgebra,Printf,SparseArrays,Test + skip: LinearAlgebra,Printf,SparseArrays,DelimitedFiles,Test projects: ., test - uses: julia-actions/julia-buildpkg@v1 env: diff --git a/NEWS.md b/NEWS.md index d1368fd0a..d75b3fd1d 100644 --- a/NEWS.md +++ b/NEWS.md @@ -5,6 +5,13 @@ DispersiveShallowWater.jl follows the interpretation of used in the Julia ecosystem. Notable changes will be documented in this file for human readability. +## Changes in the v0.8 lifecycle + +#### Added + +- Add function `data_dingemans` to load experimental data for the Dingemans experiment and expand the documentation + about the experiment ([#238]). + ## Changes when updating to v0.8 from v0.7.x #### Changed diff --git a/Project.toml b/Project.toml index 5cb6a0758..5805a4225 100644 --- a/Project.toml +++ b/Project.toml @@ -5,6 +5,7 @@ version = "0.8.1-DEV" [deps] BandedMatrices = "aae01518-5342-5314-be14-df237901396f" +DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab" DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e" FastBroadcast = "7034ab61-46d4-4ed7-9d0f-46aef9175898" ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210" @@ -27,6 +28,7 @@ TrixiBase = "9a0f1c46-06d5-4909-a5a3-ce25d3fa3284" [compat] BandedMatrices = "1.5" +DelimitedFiles = "1" DiffEqBase = "6.160" FastBroadcast = "0.3.5" ForwardDiff = "0.10.36, 1" diff --git a/data/Dingemans.csv b/data/Dingemans.csv new file mode 100644 index 000000000..ff6ca956d --- /dev/null +++ b/data/Dingemans.csv @@ -0,0 +1,1203 @@ +time,x1,x2,x3,x4,x5,x6 +10.000,0.79640557,0.7999320514,0.800267507,0.7999659904,0.800113309,0.8003631130000001 +10.050,0.79532754,0.8000288719,0.80031518,0.8000138195000001,0.800113309,0.8003631130000001 +10.100,0.7941515100000001,0.8000288719,0.800267507,0.7999659904,0.8000208649,0.8002645540000001 +10.150,0.7928284800000001,0.8000288719,0.80031518,0.8000138195000001,0.8000670871000001,0.800313833 +10.200,0.79130944,0.8000288719,0.8003628530000001,0.8000616485,0.8000670871000001,0.8002645540000001 +10.250,0.7899864000000001,0.8000288719,0.8004105260000001,0.7999659904,0.8000670871000001,0.800313833 +10.300,0.7889574,0.7999804616,0.8004581980000001,0.7999659904,0.800159532,0.8002645540000001 +10.350,0.7878794,0.8000772821000001,0.8004105260000001,0.8000616485,0.800113309,0.800313833 +10.400,0.7868013,0.800174103,0.8004105260000001,0.8000138195000001,0.800113309,0.8003631130000001 +10.450,0.7858703,0.800270923,0.8004105260000001,0.8000138195000001,0.800205754,0.800313833 +10.500,0.7852823,0.800222513,0.8003628530000001,0.8000138195000001,0.800113309,0.8002645540000001 +10.550,0.7851353000000001,0.800270923,0.8004105260000001,0.7999659904,0.800113309,0.8002645540000001 +10.600,0.7851353000000001,0.800270923,0.8004105260000001,0.7999659904,0.800159532,0.8002645540000001 +10.650,0.7853803,0.800270923,0.8004581980000001,0.8000616485,0.800159532,0.800313833 +10.700,0.7858213000000001,0.800319333,0.80031518,0.8000138195000001,0.800205754,0.800215274 +10.750,0.7870463000000001,0.800464564,0.8004105260000001,0.8000138195000001,0.800159532,0.8002645540000001 +10.800,0.7884184000000001,0.800464564,0.8003628530000001,0.8000138195000001,0.800159532,0.800313833 +10.850,0.79018241,0.8005129740000001,0.8003628530000001,0.8000616485,0.800159532,0.8002645540000001 +10.900,0.7920444600000001,0.8005129740000001,0.80031518,0.8000616485,0.800113309,0.8002645540000001 +10.950,0.79444552,0.8006582050000001,0.8003628530000001,0.8000616485,0.800159532,0.800313833 +11.000,0.7972385900000001,0.8006097950000001,0.80031518,0.8000616485,0.800113309,0.8002645540000001 +11.050,0.800129661,0.800755026,0.8003628530000001,0.8001094780000001,0.800159532,0.800313833 +11.100,0.8030207300000001,0.800706616,0.80031518,0.8001094780000001,0.800113309,0.8002645540000001 +11.150,0.8056178,0.800851846,0.80031518,0.8000616485,0.800159532,0.8002645540000001 +11.200,0.80875388,0.800851846,0.800267507,0.8001094780000001,0.800113309,0.800215274 +11.250,0.8116449,0.800948667,0.800267507,0.8000616485,0.800113309,0.8002645540000001 +11.300,0.8147810000000001,0.800948667,0.800267507,0.8000616485,0.800113309,0.8002645540000001 +11.350,0.8175251,0.800948667,0.80031518,0.8001094780000001,0.800159532,0.800313833 +11.400,0.8193871,0.800997077,0.80031518,0.8000616485,0.8000670871000001,0.800313833 +11.450,0.8212002,0.8010939,0.800267507,0.8000138195000001,0.800113309,0.8002645540000001 +11.500,0.8227682000000001,0.8010939,0.800267507,0.7999659904,0.8000670871000001,0.8002645540000001 +11.550,0.8235032,0.80104549,0.8003628530000001,0.7999659904,0.8000670871000001,0.8003631130000001 +11.600,0.8236993,0.800997077,0.80031518,0.8001094780000001,0.8000208649,0.8002645540000001 +11.650,0.8237483000000001,0.800948667,0.80031518,0.8000616485,0.8000670871000001,0.8002645540000001 +11.700,0.8233562000000001,0.800948667,0.80031518,0.8001094780000001,0.8000670871000001,0.800215274 +11.750,0.8221312000000001,0.800851846,0.8003628530000001,0.8000616485,0.7999746426000001,0.800313833 +11.800,0.8209552,0.800755026,0.8003628530000001,0.8001094780000001,0.8000670871000001,0.8002645540000001 +11.850,0.8193381000000001,0.8006097950000001,0.8003628530000001,0.8001573070000001,0.8000208649,0.800313833 +11.900,0.8173291,0.800367744,0.80031518,0.8001094780000001,0.7999746426000001,0.800313833 +11.950,0.815369,0.8000772821000001,0.80031518,0.8000616485,0.8000670871000001,0.8002645540000001 +12.000,0.8133600000000001,0.7998836410000001,0.800267507,0.8001094780000001,0.8000670871000001,0.800313833 +12.050,0.8111549,0.79959318,0.800267507,0.8001094780000001,0.8000670871000001,0.800313833 +12.100,0.80894988,0.7993027180000001,0.80031518,0.8002051360000001,0.800113309,0.800313833 +12.150,0.80694083,0.79891544,0.800219835,0.8001094780000001,0.8000208649,0.800313833 +12.200,0.80502978,0.79843133,0.800267507,0.8001094780000001,0.800113309,0.800412392 +12.250,0.80292273,0.79799564,0.800267507,0.8001573070000001,0.8000208649,0.8003631130000001 +12.300,0.8009626820000001,0.79751154,0.80031518,0.8001573070000001,0.800113309,0.800313833 +12.350,0.799100636,0.79697903,0.800267507,0.8001573070000001,0.800113309,0.800313833 +12.400,0.79728759,0.79644651,0.800267507,0.8001573070000001,0.800113309,0.800313833 +12.450,0.79537654,0.79586559,0.800267507,0.8000616485,0.800113309,0.8002645540000001 +12.500,0.7933184900000001,0.79533308,0.800267507,0.8001573070000001,0.800113309,0.8002645540000001 +12.550,0.79145644,0.79475215,0.80031518,0.8002529650000001,0.8000670871000001,0.8002645540000001 +12.600,0.79003541,0.79421964,0.8003628530000001,0.8001573070000001,0.800113309,0.8002645540000001 +12.650,0.7885654000000001,0.79363872,0.80031518,0.8001573070000001,0.8000670871000001,0.800313833 +12.700,0.7869973,0.79315462,0.8004105260000001,0.8001094780000001,0.8000670871000001,0.800215274 +12.750,0.7855763,0.79267051,0.80031518,0.8001094780000001,0.8000208649,0.8003631130000001 +12.800,0.7844003,0.79233164,0.8003628530000001,0.8001094780000001,0.8000670871000001,0.800313833 +12.850,0.7834692000000001,0.79194436,0.80031518,0.8001094780000001,0.7999746426000001,0.800313833 +12.900,0.7826362,0.79170231,0.8003628530000001,0.8000616485,0.7999746426000001,0.800313833 +12.950,0.7818032,0.79150867,0.800267507,0.8000138195000001,0.7999746426000001,0.800313833 +13.000,0.7813622,0.79146026,0.8003628530000001,0.8000616485,0.8000208649,0.800313833 +13.050,0.7811662,0.79146026,0.8003628530000001,0.8000138195000001,0.8000208649,0.800313833 +13.100,0.7812152,0.79175072,0.8004105260000001,0.8001094780000001,0.800113309,0.800313833 +13.150,0.7813622,0.792138,0.80031518,0.7999659904,0.8000208649,0.8002645540000001 +13.200,0.7816562,0.79267051,0.8003628530000001,0.8000616485,0.800113309,0.8003631130000001 +13.250,0.7818522,0.79334826,0.8003628530000001,0.8000616485,0.8000670871000001,0.8002645540000001 +13.300,0.7824892,0.79421964,0.8004105260000001,0.8001094780000001,0.800159532,0.8003631130000001 +13.350,0.7833222000000001,0.79523626,0.800267507,0.8001094780000001,0.8000670871000001,0.8002645540000001 +13.400,0.7846943000000001,0.79649492,0.8003628530000001,0.8000138195000001,0.8000670871000001,0.800313833 +13.450,0.7861153000000001,0.797802,0.800219835,0.8001094780000001,0.8000208649,0.8002645540000001 +13.500,0.7877814000000001,0.799447949,0.80031518,0.8000616485,0.8000670871000001,0.800313833 +13.550,0.7895944,0.80114231,0.800267507,0.8000138195000001,0.7999746426000001,0.8002645540000001 +13.600,0.79170145,0.8029819,0.8003628530000001,0.8001094780000001,0.8000670871000001,0.8002645540000001 +13.650,0.7937595000000001,0.80453103,0.80031518,0.8000138195000001,0.7999746426000001,0.800313833 +13.700,0.7958665500000001,0.80637062,0.8003628530000001,0.8001094780000001,0.8000670871000001,0.800313833 +13.750,0.79792461,0.80840385,0.80031518,0.8000616485,0.7999746426000001,0.800313833 +13.800,0.800374667,0.8106791000000001,0.8003628530000001,0.8001573070000001,0.800113309,0.8003631130000001 +13.850,0.80267772,0.8126639,0.800267507,0.8001094780000001,0.8000208649,0.8002645540000001 +13.900,0.80468678,0.8147456000000001,0.8004105260000001,0.8001094780000001,0.8000670871000001,0.800313833 +13.950,0.8064998200000001,0.8164884,0.80031518,0.8000616485,0.8000208649,0.800215274 +14.000,0.80850887,0.8182311,0.8004581980000001,0.8000616485,0.8000670871000001,0.8003631130000001 +14.050,0.8104199000000001,0.8194898,0.8003628530000001,0.8000138195000001,0.8000208649,0.800313833 +14.100,0.811988,0.8206032000000001,0.8003628530000001,0.8001094780000001,0.8000208649,0.800412392 +14.150,0.813605,0.8214262,0.800267507,0.8000138195000001,0.8000208649,0.800313833 +14.200,0.815516,0.8219103000000001,0.8003628530000001,0.8000616485,0.8000208649,0.800412392 +14.250,0.8168881,0.8217651,0.800267507,0.8000138195000001,0.7999284204,0.8003631130000001 +14.300,0.8182111000000001,0.8214746,0.8003628530000001,0.8001094780000001,0.8000208649,0.8003631130000001 +14.350,0.8192401,0.8207485,0.8003628530000001,0.8000138195000001,0.7999746426000001,0.800313833 +14.400,0.8203672000000001,0.8194898,0.8004105260000001,0.8000616485,0.8000208649,0.8003631130000001 +14.450,0.8210532,0.8178438,0.80031518,0.7999659904,0.7999746426000001,0.8002645540000001 +14.500,0.8212982000000001,0.815859,0.8004105260000001,0.8000138195000001,0.8000670871000001,0.800313833 +14.550,0.8212002,0.8134385000000001,0.80031518,0.7999659904,0.7999746426000001,0.800165994 +14.600,0.8209552,0.8110664000000001,0.8004581980000001,0.8001094780000001,0.8000670871000001,0.800313833 +14.650,0.8199262,0.80821021,0.8004105260000001,0.7999659904,0.7999746426000001,0.8002645540000001 +14.700,0.8189461,0.80545082,0.800553544,0.8000138195000001,0.800113309,0.8003631130000001 +14.750,0.8176721,0.80225574,0.8004581980000001,0.7999181613,0.8000208649,0.8002645540000001 +14.800,0.8161041,0.7992543080000001,0.800553544,0.8000616485,0.800113309,0.800313833 +14.850,0.8143400000000001,0.79605923,0.800505871,0.7999659904,0.8000670871000001,0.8002645540000001 +14.900,0.812674,0.79320303,0.8006488900000001,0.8000138195000001,0.800159532,0.800412392 +14.950,0.8108609,0.79034682,0.800553544,0.7999659904,0.8000208649,0.8002645540000001 +15.000,0.80870487,0.7878295000000001,0.8006965620000001,0.8000138195000001,0.800113309,0.800313833 +15.050,0.8064998200000001,0.7856026,0.8006488900000001,0.7999659904,0.8000208649,0.800215274 +15.100,0.8044417700000001,0.7836178,0.800744235,0.8000616485,0.800113309,0.800313833 +15.150,0.8018937100000001,0.7817782000000001,0.8006488900000001,0.8000138195000001,0.8000208649,0.800116715 +15.200,0.799835654,0.7803743000000001,0.800791908,0.8000616485,0.800113309,0.8002645540000001 +15.250,0.7976306000000001,0.7790672000000001,0.8006965620000001,0.7999181613,0.8000670871000001,0.800165994 +15.300,0.7954745400000001,0.7781474,0.8008395810000001,0.8000616485,0.800113309,0.8002645540000001 +15.350,0.7932204900000001,0.7774213000000001,0.8008395810000001,0.7999181613,0.800113309,0.800165994 +15.400,0.79170145,0.7773245000000001,0.800934926,0.8001094780000001,0.800113309,0.8002645540000001 +15.450,0.7898394000000001,0.7774213000000001,0.8008395810000001,0.7999659904,0.8000670871000001,0.800215274 +15.500,0.7881244000000001,0.7779538,0.800934926,0.8000616485,0.800113309,0.8002645540000001 +15.550,0.7863603,0.7785831000000001,0.8008872530000001,0.7999181613,0.8000670871000001,0.800215274 +15.600,0.7849393,0.7794545,0.800982599,0.8000138195000001,0.800113309,0.800313833 +15.650,0.7834202,0.7806164,0.800982599,0.7999659904,0.8000670871000001,0.8002645540000001 +15.700,0.7824892,0.7822139,0.8011256200000001,0.8000616485,0.8000670871000001,0.8002645540000001 +15.750,0.7817052,0.783763,0.80103027,0.7999181613,0.8000208649,0.800215274 +15.800,0.7813132,0.7856510000000001,0.8010779400000001,0.8000616485,0.800113309,0.8002645540000001 +15.850,0.7808722,0.7874422000000001,0.800982599,0.7999181613,0.7999746426000001,0.800165994 +15.900,0.7808232,0.789427,0.8010779400000001,0.8000616485,0.800113309,0.8003631130000001 +15.950,0.7811172000000001,0.79126662,0.8008872530000001,0.7999659904,0.7999746426000001,0.800215274 +16.000,0.7818032,0.79339667,0.80103027,0.8000616485,0.8000670871000001,0.8003631130000001 +16.050,0.7821952000000001,0.79552672,0.8008872530000001,0.7999659904,0.7999746426000001,0.8002645540000001 +16.100,0.7826852000000001,0.797802,0.8008872530000001,0.8000138195000001,0.800159532,0.800313833 +16.150,0.7832732,0.7998352310000001,0.8006965620000001,0.7999659904,0.8000208649,0.8002645540000001 +16.200,0.7842533,0.80191687,0.800744235,0.8000616485,0.800113309,0.8002645540000001 +16.250,0.7852333,0.80380487,0.8004581980000001,0.7999181613,0.7999746426000001,0.800165994 +16.300,0.7865073,0.80564446,0.8004105260000001,0.8000138195000001,0.800113309,0.8003631130000001 +16.350,0.7880754000000001,0.80719359,0.8001721620000001,0.7999659904,0.7999746426000001,0.800215274 +16.400,0.7898394000000001,0.80893636,0.8001721620000001,0.8000616485,0.800159532,0.8003631130000001 +16.450,0.7916524500000001,0.8104371,0.799838452,0.7999659904,0.8000208649,0.800215274 +16.500,0.7933674900000001,0.812083,0.79979078,0.8000616485,0.8000670871000001,0.800313833 +16.550,0.7950335300000001,0.8132449,0.799409397,0.7999181613,0.8000208649,0.8002645540000001 +16.600,0.79728759,0.8145519,0.79945707,0.8000616485,0.800113309,0.8003631130000001 +16.650,0.7995416470000001,0.8154233000000001,0.7990756880000001,0.7999181613,0.7999746426000001,0.800215274 +16.700,0.8017467,0.8166336000000001,0.79893267,0.8000616485,0.800113309,0.8002645540000001 +16.750,0.8036087500000001,0.8175534000000001,0.79855129,0.7999181613,0.8000208649,0.800215274 +16.800,0.8058138,0.8183764,0.79840827,0.8000616485,0.800113309,0.800313833 +16.850,0.8079208600000001,0.8188605,0.79797921,0.7999181613,0.8000670871000001,0.800165994 +16.900,0.8101259000000001,0.8194414000000001,0.7977408500000001,0.8000616485,0.800113309,0.8002645540000001 +16.950,0.8119390000000001,0.8194898,0.7973118,0.7999659904,0.8000208649,0.8000674351 +17.000,0.8139970000000001,0.8195866,0.7970734300000001,0.8000138195000001,0.800113309,0.8002645540000001 +17.050,0.8158591000000001,0.8192477,0.7966920500000001,0.7999181613,0.8000208649,0.800165994 +17.100,0.8175251,0.8187152000000001,0.79650136,0.8000138195000001,0.800113309,0.8002645540000001 +17.150,0.8187011000000001,0.817747,0.7960723000000001,0.7999659904,0.7999746426000001,0.800165994 +17.200,0.8198772000000001,0.8169241,0.79588161,0.8000616485,0.800113309,0.8002645540000001 +17.250,0.8202692,0.815617,0.7955479000000001,0.7999181613,0.8000208649,0.800165994 +17.300,0.8206612000000001,0.8141163,0.7954048800000001,0.8001094780000001,0.800113309,0.800313833 +17.350,0.8207102000000001,0.8122283,0.79507117,0.7999659904,0.8000208649,0.800165994 +17.400,0.8206612000000001,0.8102918,0.79507117,0.8001094780000001,0.800113309,0.8002645540000001 +17.450,0.8203182,0.80787133,0.79478514,0.8000616485,0.8000208649,0.800165994 +17.500,0.8198772000000001,0.80559605,0.79488048,0.8001573070000001,0.800159532,0.8002645540000001 +17.550,0.8190441,0.80288508,0.7946897900000001,0.8000616485,0.8000208649,0.800165994 +17.600,0.8180151,0.800367744,0.7948328100000001,0.8001573070000001,0.800159532,0.8002645540000001 +17.650,0.8167411,0.79755995,0.7948328100000001,0.8000616485,0.8000208649,0.800116715 +17.700,0.8152710000000001,0.79518785,0.79521419,0.8001573070000001,0.800113309,0.800313833 +17.750,0.8129190000000001,0.79271892,0.7954048800000001,0.8000616485,0.8000670871000001,0.800116715 +17.800,0.8108609,0.79068569,0.7959292800000001,0.8001094780000001,0.800159532,0.8002645540000001 +17.850,0.8086068700000001,0.7886525000000001,0.79626299,0.8000138195000001,0.8000208649,0.800116715 +17.900,0.80654882,0.7870065,0.7969780900000001,0.8001573070000001,0.800113309,0.8002645540000001 +17.950,0.8044417700000001,0.7852153000000001,0.79755016,0.8000616485,0.7999746426000001,0.800116715 +18.000,0.80267772,0.7838598,0.79855129,0.8001094780000001,0.8000670871000001,0.8002645540000001 +18.050,0.80047267,0.7825528,0.799361725,0.8000616485,0.8000208649,0.800165994 +18.100,0.79831662,0.781633,0.800553544,0.8002051360000001,0.8000670871000001,0.8002645540000001 +18.150,0.7959155600000001,0.78081,0.80160235,0.8000616485,0.7999746426000001,0.800215274 +18.200,0.7938085,0.7805679000000001,0.80298486,0.8001573070000001,0.800113309,0.8002645540000001 +18.250,0.79170145,0.7802775000000001,0.8040336600000001,0.8000616485,0.7999746426000001,0.800165994 +18.300,0.7899864000000001,0.7805195,0.8054161700000001,0.8002051360000001,0.800113309,0.8002645540000001 +18.350,0.7880754000000001,0.78081,0.80679868,0.8000138195000001,0.7999284204,0.800116715 +18.400,0.7866543,0.7813909,0.8085625700000001,0.8001094780000001,0.8000208649,0.8002645540000001 +18.450,0.7852823,0.7818750000000001,0.8100881000000001,0.8000616485,0.7999284204,0.800215274 +18.500,0.7843023,0.7827464000000001,0.8118043,0.8002051360000001,0.800113309,0.8002645540000001 +18.550,0.7831262,0.7833757,0.8130915000000001,0.8001094780000001,0.7999746426000001,0.800165994 +18.600,0.7822442000000001,0.7843924,0.8145693,0.8002529650000001,0.800113309,0.800313833 +18.650,0.7813622,0.7853606000000001,0.8157135000000001,0.8001094780000001,0.7999746426000001,0.800215274 +18.700,0.7813132,0.7867161,0.8170007,0.8002529650000001,0.8000208649,0.8002645540000001 +18.750,0.7812642000000001,0.7879263000000001,0.8177157,0.8001094780000001,0.8000208649,0.800165994 +18.800,0.7815092,0.7894754,0.8184785,0.8003486230000001,0.800113309,0.8002645540000001 +18.850,0.7815092,0.79087933,0.8186692000000001,0.8003486230000001,0.7999746426000001,0.800165994 +18.900,0.7818032,0.7926221,0.8189076000000001,0.8003964520000001,0.8000670871000001,0.800313833 +18.950,0.7822442000000001,0.79436487,0.8185739000000001,0.8002529650000001,0.8000208649,0.800116715 +19.000,0.7831262,0.79644651,0.8182878,0.8003486230000001,0.800113309,0.8002645540000001 +19.050,0.7838123,0.7982861,0.8174774,0.8003007940000001,0.8000208649,0.800165994 +19.100,0.7849393,0.800464564,0.8166669000000001,0.8004921100000001,0.8000670871000001,0.800313833 +19.150,0.7862623000000001,0.8024978,0.8152844,0.8003486230000001,0.8000208649,0.800165994 +19.200,0.7882714000000001,0.80491831,0.8139019000000001,0.8005399390000001,0.800113309,0.8002645540000001 +19.250,0.7897904,0.80714518,0.8120904,0.8004442810000001,0.7999746426000001,0.800165994 +19.300,0.7916034500000001,0.80966251,0.8104695000000001,0.8005877680000001,0.800113309,0.8003631130000001 +19.350,0.7931714900000001,0.8117442,0.8084195500000001,0.8004921100000001,0.7999284204,0.800215274 +19.400,0.7952785400000001,0.8140678,0.8064649700000001,0.8005877680000001,0.800113309,0.800313833 +19.450,0.79714059,0.8159558,0.80417667,0.8004442810000001,0.7999746426000001,0.800165994 +19.500,0.799345642,0.8178923,0.80222209,0.8006355970000001,0.800113309,0.800313833 +19.550,0.80135469,0.8191993000000001,0.7999814707,0.8005399390000001,0.7999746426000001,0.800165994 +19.600,0.8036577500000001,0.8206032000000001,0.79802689,0.8006834260000001,0.800113309,0.8002645540000001 +19.650,0.8057648000000001,0.8212326000000001,0.79588161,0.8005399390000001,0.7999746426000001,0.800165994 +19.700,0.8078718500000001,0.8216683,0.7939747,0.8006834260000001,0.800113309,0.800313833 +19.750,0.80929289,0.8214746,0.7919247700000001,0.8005399390000001,0.7999284204,0.800165994 +19.800,0.8113509000000001,0.8213294,0.79011321,0.8006834260000001,0.8000670871000001,0.8002645540000001 +19.850,0.813752,0.8205548,0.788111,0.8005399390000001,0.7999284204,0.800215274 +19.900,0.8161531000000001,0.8197803,0.7864424000000001,0.8005877680000001,0.8000208649,0.800313833 +19.950,0.8176231,0.81857,0.7846785000000001,0.8004921100000001,0.7999746426000001,0.800165994 +20.000,0.8187011000000001,0.8172629,0.7832007000000001,0.8004921100000001,0.8000670871000001,0.800313833 +20.050,0.8194361,0.8155202,0.7815798,0.8003964520000001,0.7999284204,0.800165994 +20.100,0.8197792,0.8137774000000001,0.7804356,0.8004442810000001,0.8000208649,0.8003631130000001 +20.150,0.8199752,0.8116957,0.7791008,0.8002529650000001,0.7999284204,0.800165994 +20.200,0.8208082000000001,0.80990457,0.7781950000000001,0.8003007940000001,0.8000208649,0.800313833 +20.250,0.8210532,0.80787133,0.7772892,0.8000138195000001,0.7999284204,0.800215274 +20.300,0.8209552,0.80588651,0.7768125,0.8000616485,0.8000670871000001,0.800313833 +20.350,0.8197302000000001,0.80356282,0.7764311,0.799774674,0.7999284204,0.800215274 +20.400,0.8187011000000001,0.80162641,0.7764788,0.799822503,0.800113309,0.8003631130000001 +20.450,0.8178681000000001,0.799447949,0.7764311,0.7994877,0.7999284204,0.800215274 +20.500,0.8169371000000001,0.79770518,0.7768602,0.7994877,0.800113309,0.8003631130000001 +20.550,0.815222,0.79576877,0.7773846000000001,0.799248554,0.7999746426000001,0.800165994 +20.600,0.8137030000000001,0.79407441,0.7785764000000001,0.7991528960000001,0.800159532,0.800313833 +20.650,0.8110079,0.79223482,0.7799112,0.79881809,0.7999746426000001,0.800165994 +20.700,0.80880288,0.7907341,0.7818181000000001,0.7987224300000001,0.800159532,0.800313833 +20.750,0.8066468200000001,0.7888461,0.7839634000000001,0.7983398,0.8000208649,0.800116715 +20.800,0.80468678,0.7873454000000001,0.7867284,0.79824414,0.800205754,0.800313833 +20.850,0.80277573,0.7856026,0.7897795000000001,0.79786151,0.8000670871000001,0.800215274 +20.900,0.8010606800000001,0.7843924,0.79354565,0.7978136800000001,0.8002519760000001,0.8002645540000001 +20.950,0.79865962,0.7830853000000001,0.7972164500000001,0.79743105,0.8000670871000001,0.800116715 +21.000,0.79650357,0.7824075,0.80160235,0.79728756,0.800205754,0.800313833 +21.050,0.7941025100000001,0.7815846000000001,0.8059405700000001,0.7968571000000001,0.8000670871000001,0.800116715 +21.100,0.7920444600000001,0.7812941,0.8108509,0.79680927,0.800205754,0.8002645540000001 +21.150,0.7897414,0.78081,0.8152844,0.79633098,0.8000670871000001,0.8000674351 +21.200,0.7880264,0.7807616000000001,0.819432,0.79628315,0.800205754,0.8002645540000001 +21.250,0.7865563000000001,0.7806164,0.8226737000000001,0.79590052,0.800113309,0.800165994 +21.300,0.7855763,0.7809552000000001,0.8261538,0.79585269,0.800205754,0.8002645540000001 +21.350,0.7842533,0.7811489,0.8284421000000001,0.79542223,0.800113309,0.8000674351 +21.400,0.7832242,0.7819718,0.8303014000000001,0.7954700600000001,0.8002519760000001,0.800313833 +21.450,0.7821952000000001,0.7827948,0.8309211000000001,0.79513526,0.800113309,0.800116715 +21.500,0.7818522,0.7841019,0.8310164,0.79518308,0.8003444200000001,0.800313833 +21.550,0.7813132,0.7851669,0.8298723,0.79489611,0.800113309,0.800116715 +21.600,0.7811172000000001,0.7869097,0.8284421000000001,0.7950396000000001,0.8003444200000001,0.800313833 +21.650,0.7808722,0.7885072000000001,0.8261062,0.79489611,0.800159532,0.8000674351 +21.700,0.7813132,0.79068569,0.8236748,0.79513526,0.8003906430000001,0.8002645540000001 +21.750,0.7816072000000001,0.79242846,0.8205284,0.79513526,0.800205754,0.8000674351 +21.800,0.7825382000000001,0.79465533,0.8177634,0.79556572,0.800436865,0.800313833 +21.850,0.7829792000000001,0.79654333,0.8146647,0.7957092,0.8002519760000001,0.8000674351 +21.900,0.7842533,0.79891544,0.8122334,0.79637881,0.800436865,0.800313833 +21.950,0.7850863,0.80114231,0.80942068,0.79680927,0.800298198,0.8000674351 +22.000,0.7865073,0.80370805,0.80713239,0.79786151,0.800436865,0.800313833 +22.050,0.7876343,0.80569287,0.8046534000000001,0.79857895,0.8003444200000001,0.8000674351 +22.100,0.7893004,0.80811339,0.80279416,0.8000616485,0.8005293090000001,0.800313833 +22.150,0.79091743,0.8100982,0.8006012170000001,0.8013052,0.8003444200000001,0.8000674351 +22.200,0.79326949,0.8123735000000001,0.7989803400000001,0.80312271,0.8005755320000001,0.800313833 +22.250,0.7952785400000001,0.8142131,0.7971211,0.80455758,0.8003906430000001,0.800116715 +22.300,0.7975326,0.8162463000000001,0.79573859,0.8068533800000001,0.8005293090000001,0.800313833 +22.350,0.799198638,0.8175534000000001,0.7940700500000001,0.80910134,0.8003444200000001,0.800165994 +22.400,0.80120769,0.8189573,0.79283055,0.8118754,0.800483087,0.800412392 +22.450,0.80311874,0.8196834000000001,0.7913050300000001,0.8142669,0.8002519760000001,0.800116715 +22.500,0.80547079,0.8207001,0.7902562200000001,0.8168975,0.800483087,0.8003631130000001 +22.550,0.8075288500000001,0.8208937000000001,0.7889691000000001,0.8188106,0.800298198,0.800116715 +22.600,0.8098809100000001,0.8211357,0.7881586,0.8206760000000001,0.800483087,0.800313833 +22.650,0.8115469000000001,0.8207001,0.7870145000000001,0.8214891000000001,0.800113309,0.800116715 +22.700,0.8135560000000001,0.8205548,0.7864424000000001,0.8222065000000001,0.8003444200000001,0.800313833 +22.750,0.815369,0.8197318,0.7853459,0.8217282,0.8000670871000001,0.8000181555 +22.800,0.8172311000000001,0.8190541,0.7847739,0.8211064,0.8002519760000001,0.800313833 +22.850,0.8183091,0.8176018,0.7837251000000001,0.8195759,0.7999746426000001,0.8000674351 +22.900,0.8195831,0.8163915,0.7832483,0.818141,0.800113309,0.800313833 +22.950,0.8201222,0.8144551000000001,0.7823425,0.8160844,0.7997897540000001,0.8000674351 +23.000,0.8209062,0.8126639,0.7821042,0.8144104000000001,0.799882198,0.800313833 +23.050,0.8208082000000001,0.8103403,0.7814368,0.8123059,0.799604865,0.8000674351 +23.100,0.8209552,0.80850067,0.7812937000000001,0.8109667,0.7996510870000001,0.800313833 +23.150,0.8206122,0.80608015,0.7806740000000001,0.8093404900000001,0.799373754,0.800116715 +23.200,0.8205142000000001,0.80409533,0.7807217000000001,0.80819259,0.799419976,0.800313833 +23.250,0.8189951000000001,0.80186846,0.7804356,0.80675772,0.79909642,0.8000674351 +23.300,0.8179661,0.8000288719,0.7809124000000001,0.80594462,0.7991888650000001,0.8003631130000001 +23.350,0.8168881,0.79775359,0.7811984000000001,0.8047489000000001,0.7988190900000001,0.800116715 +23.400,0.8158591000000001,0.79596241,0.7821519,0.80403146,0.7988190900000001,0.8003631130000001 +23.450,0.813801,0.79383236,0.7826763,0.80293139,0.7984955300000001,0.8000674351 +23.500,0.812037,0.79223482,0.7838681000000001,0.80230961,0.7985879800000001,0.800313833 +23.550,0.80914589,0.79020159,0.7850122,0.80120955,0.79817198,0.800116715 +23.600,0.8070388300000001,0.7887009,0.7871575000000001,0.8004921100000001,0.7982644200000001,0.8003631130000001 +23.650,0.8049317800000001,0.7867645000000001,0.7892551000000001,0.799296384,0.7978022,0.800165994 +23.700,0.8032167400000001,0.7856026,0.7919247700000001,0.7984832900000001,0.7979408600000001,0.800412392 +23.750,0.800766677,0.7841503000000001,0.7940700500000001,0.7970962500000001,0.79752486,0.800116715 +23.800,0.7990026330000001,0.7834242,0.79726412,0.7962353200000001,0.79761731,0.800412392 +23.850,0.7969935800000001,0.7824075,0.8004105260000001,0.79489611,0.79729375,0.800215274 +23.900,0.7950335300000001,0.7820687000000001,0.80436737,0.79413084,0.79743242,0.8004616720000001 +23.950,0.79224046,0.7811973000000001,0.8078474800000001,0.79283946,0.79706264,0.800165994 +24.000,0.79013341,0.7809068,0.8118043,0.7920263700000001,0.79724753,0.800510952 +24.050,0.7880264,0.7802775000000001,0.8151414,0.7907349800000001,0.79687775,0.800215274 +24.100,0.7868013,0.7803743000000001,0.8189552,0.7899219000000001,0.79706264,0.800510952 +24.150,0.7851353000000001,0.7803743000000001,0.8220063000000001,0.7886305,0.79687775,0.8002645540000001 +24.200,0.7843513000000001,0.7812457,0.8250573,0.7879131,0.7971550900000001,0.800609511 +24.250,0.7833712,0.7817782000000001,0.8268689,0.7867173000000001,0.7970164200000001,0.8002645540000001 +24.300,0.7830772,0.7830369,0.8285851000000001,0.7860956,0.7973862,0.8005602310000001 +24.350,0.7822442000000001,0.7838114,0.8288712,0.7848998,0.7973862,0.8002645540000001 +24.400,0.7817542000000001,0.7853606000000001,0.8289665,0.7843259,0.7978946400000001,0.8005602310000001 +24.450,0.7808722,0.7865708,0.82787,0.7832258000000001,0.79803331,0.800313833 +24.500,0.7810192,0.7886041,0.8267259,0.7828432000000001,0.79872664,0.8005602310000001 +24.550,0.7810192,0.79025,0.8245806,0.7819344,0.799050198,0.800313833 +24.600,0.7818522,0.79252528,0.8226260000000001,0.7816475,0.7999746426000001,0.8005602310000001 +24.650,0.7822932,0.79417123,0.8197657,0.7808344,0.800483087,0.800313833 +24.700,0.7833712,0.79649492,0.817382,0.7806909000000001,0.8015924200000001,0.8005602310000001 +24.750,0.7840573000000001,0.79838292,0.814474,0.7802126,0.8023319800000001,0.800313833 +24.800,0.7853803,0.800803436,0.8122334,0.7805952,0.80371864,0.800609511 +24.850,0.7862133,0.80273985,0.80942068,0.7807387,0.80464309,0.8002645540000001 +24.900,0.7876834,0.80520877,0.8071800600000001,0.7819344,0.80607598,0.8005602310000001 +24.950,0.7889574,0.807242,0.80460573,0.783178,0.80727775,0.8002645540000001 +25.000,0.79106443,0.80985616,0.80269882,0.7857608,0.8090342,0.8004616720000001 +25.050,0.7928284800000001,0.811841,0.8004105260000001,0.7891088,0.8103746000000001,0.800215274 +25.100,0.79498453,0.8142131,0.79864663,0.7945613100000001,0.8121773,0.8004616720000001 +25.150,0.79679758,0.8160043,0.79650136,0.80149652,0.8133791,0.800165994 +25.200,0.799198638,0.8178923,0.7950235,0.8113015,0.8149506000000001,0.800313833 +25.250,0.80101168,0.8190057000000001,0.7931642600000001,0.8229239,0.8159213000000001,0.8000674351 +25.300,0.8032167400000001,0.8203128000000001,0.7918771,0.8367465000000001,0.8170306,0.8002645540000001 +25.350,0.8049807800000001,0.8207001,0.7901608800000001,0.8504257000000001,0.8174466,0.7999688759000001 +25.400,0.80733284,0.8212326000000001,0.7891121000000001,0.8587001000000001,0.8179089,0.800165994 +25.450,0.80934189,0.8210389,0.7875865000000001,0.8573609000000001,0.8175853000000001,0.7999195963000001 +25.500,0.8117920000000001,0.8211842,0.7866331000000001,0.8479864,0.817308,0.8000674351 +25.550,0.8135070000000001,0.8205548,0.7852029,0.8339725,0.8161524,0.799771757 +25.600,0.815222,0.8201191000000001,0.7844401000000001,0.8202455000000001,0.8150893,0.799870317 +25.650,0.8162021,0.8188121,0.783296,0.80766647,0.8131018000000001,0.799574639 +25.700,0.8180641000000001,0.817747,0.7827239,0.799344213,0.8112066,0.7997224780000001 +25.750,0.8192891,0.8160043,0.7816751000000001,0.79561355,0.80871064,0.7994268 +25.800,0.8201712000000001,0.8144551000000001,0.7813891000000001,0.7926481400000001,0.8065844200000001,0.799574639 +25.850,0.8202692,0.8123251,0.7806740000000001,0.7897306000000001,0.80394975,0.7991804020000001 +25.900,0.8210042000000001,0.8108728000000001,0.7802449,0.7887262,0.8019159800000001,0.799328241 +25.950,0.8207592,0.80879113,0.7794822,0.7892045000000001,0.79951242,0.7989832800000001 +26.000,0.8205632,0.80695154,0.7793868,0.79063932,0.7978022,0.7990818430000001 +26.050,0.8199752,0.80462785,0.7791008,0.7921698500000001,0.79590709,0.7987861700000001 +26.100,0.8195831,0.80269144,0.7795299,0.79446565,0.79484398,0.798934 +26.150,0.8181131,0.800367744,0.7796252,0.7965223,0.79373464,0.79858905 +26.200,0.8171331000000001,0.79852815,0.7804356,0.7989137500000001,0.7934110900000001,0.79873689 +26.250,0.815516,0.79620446,0.7809600000000001,0.8003486230000001,0.79304131,0.7983919300000001 +26.300,0.8141440000000001,0.79421964,0.7823425,0.8016400100000001,0.7933186400000001,0.79858905 +26.350,0.811988,0.79204118,0.7834867000000001,0.80173567,0.7935035300000001,0.7982440900000001 +26.400,0.8098809100000001,0.79063728,0.7853459,0.80168784,0.79410442,0.7983919300000001 +26.450,0.8071368400000001,0.7887493000000001,0.7870621000000001,0.8004442810000001,0.7944742,0.79819481 +26.500,0.80502978,0.7873454000000001,0.7895888000000001,0.799248554,0.7953062000000001,0.7983919300000001 +26.550,0.80267772,0.7856994,0.79211546,0.7971919000000001,0.79558353,0.79814553 +26.600,0.8009626820000001,0.7846828,0.7954048800000001,0.7953744,0.7961382,0.79849049 +26.650,0.7990026330000001,0.7834242,0.7984559400000001,0.7929829500000001,0.7961382,0.79829337 +26.700,0.7972385900000001,0.7827948,0.8020314000000001,0.79102196,0.79623064,0.79868761 +26.750,0.79459252,0.7818266,0.80532082,0.7886783,0.7956759800000001,0.7985397700000001 +26.800,0.7924364700000001,0.7814877,0.80942068,0.7869565000000001,0.7953062000000001,0.7989832800000001 +26.850,0.7903294200000001,0.7808584000000001,0.8129008000000001,0.7851868000000001,0.79428931,0.7989832800000001 +26.900,0.7889084000000001,0.7807616000000001,0.8166669000000001,0.7842781000000001,0.79354975,0.799525359 +26.950,0.7868503,0.7804227,0.819575,0.7835128,0.7923479800000001,0.799525359 +27.000,0.7853803,0.7809068,0.8226260000000001,0.7834650000000001,0.7914235300000001,0.800116715 +27.050,0.7840083,0.7812457,0.8247236,0.7833693,0.7901293100000001,0.800313833 +27.100,0.7835182,0.7821655000000001,0.8266782,0.7836563000000001,0.7893435000000001,0.800855909 +27.150,0.7826852000000001,0.7827948,0.8274887000000001,0.7835128,0.7884191,0.8012008700000001 +27.200,0.7822932,0.7839567000000001,0.8280131000000001,0.7837519000000001,0.7883266,0.8019400600000001 +27.250,0.7815582000000001,0.7848765000000001,0.8274887000000001,0.7833215,0.7882342,0.8023343000000001 +27.300,0.7815582000000001,0.7865708,0.8266782,0.783178,0.7888351,0.8030734900000001 +27.350,0.7811172000000001,0.7879263000000001,0.8250097000000001,0.7825562,0.7893898,0.8035170100000001 +27.400,0.7812642000000001,0.79000795,0.8235318,0.7824127000000001,0.79059153,0.8043547600000001 +27.450,0.7813622,0.7916539,0.8214342,0.7820301000000001,0.7917008600000001,0.80484756 +27.500,0.7822932,0.79373554,0.8194796000000001,0.7823649,0.7935959800000001,0.80573459 +27.550,0.7832732,0.79538149,0.81681,0.7825562,0.79558353,0.80627667 +27.600,0.7843513000000001,0.79755995,0.8146647,0.7832258000000001,0.7982182,0.80726226 +27.650,0.7850863,0.79973841,0.8122334,0.7835606,0.80076042,0.80785361 +27.700,0.7863113,0.80240097,0.8104695000000001,0.7841346,0.8038110900000001,0.8088392000000001 +27.750,0.7874383,0.80457944,0.8083242100000001,0.7846129000000001,0.8067230900000001,0.80928272 +27.800,0.7891044,0.80729041,0.80665566,0.7860477,0.8100973,0.8101698 +27.850,0.7907214300000001,0.80951728,0.80446271,0.7885827000000001,0.8132404000000001,0.8106626 +27.900,0.7928284800000001,0.811841,0.8026511500000001,0.79317426,0.8167533,0.8113525 +27.950,0.7946905200000001,0.8135837,0.8006965620000001,0.79886592,0.819388,0.8116481000000001 +28.000,0.7968465800000001,0.8155686000000001,0.799218706,0.8059924500000001,0.8225311000000001,0.8121409 +28.050,0.79865962,0.8169241,0.7973594700000001,0.8140756,0.8252120000000001,0.8121902000000001 +28.100,0.800717676,0.8183279,0.79573859,0.8241675,0.8278466,0.8124366000000001 +28.150,0.80253072,0.8191509,0.79359332,0.8344986000000001,0.8294644,0.8122395 +28.200,0.8049317800000001,0.8200223000000001,0.79197244,0.8435861,0.8303426,0.8121902000000001 +28.250,0.8070878300000001,0.8204096000000001,0.79006553,0.8485125,0.8298804000000001,0.8117960000000001 +28.300,0.8094889000000001,0.8208937000000001,0.7884923,0.8479385,0.8284938000000001,0.811451 +28.350,0.8114489,0.8207969,0.7865854,0.8417686000000001,0.8254431000000001,0.8107118000000001 +28.400,0.813752,0.8206516,0.7850599,0.8328724000000001,0.8212831,0.8100219000000001 +28.450,0.815222,0.8198287000000001,0.783439,0.8229718,0.8157364,0.8089377600000001 +28.500,0.8165451,0.8189089,0.7822949,0.8142191000000001,0.80949642,0.80790289 +28.550,0.8176721,0.8174566000000001,0.7809600000000001,0.80761864,0.80256309,0.8064245000000001 +28.600,0.8189461,0.8161979,0.7801496,0.8047489000000001,0.79590709,0.8052417900000001 +28.650,0.8197302000000001,0.8145519,0.7789578,0.80106606,0.7892511000000001,0.80371413 +28.700,0.8207102000000001,0.8128576000000001,0.7782427000000001,0.79762237,0.7834271,0.80228502 +28.750,0.8210532,0.8108244,0.7774322,0.7965701300000001,0.7782502,0.80070807 +28.800,0.8210532,0.80898477,0.7770985,0.7971919000000001,0.7743675000000001,0.799328241 +28.850,0.8206612000000001,0.80690313,0.7767172,0.79757454,0.7715018,0.7977020100000001 +28.900,0.8204652,0.80496672,0.7768602,0.79824414,0.7700689000000001,0.79637146 +28.950,0.8193871,0.80269144,0.7771462,0.79857895,0.7696529000000001,0.7948438 +29.000,0.8184561,0.8006097950000001,0.7778613000000001,0.7994877,0.7706698000000001,0.79371037 +29.050,0.8171821,0.7982861,0.7784334,0.8004442810000001,0.7726111,0.7924291 +29.100,0.8158101000000001,0.79649492,0.7795775,0.80149652,0.7760315000000001,0.79149278 +29.150,0.814046,0.79431646,0.7807217000000001,0.80168784,0.7802840000000001,0.79050719 +29.200,0.8122820000000001,0.79257369,0.7827239,0.8014486900000001,0.7855071,0.7899158000000001 +29.250,0.8097339,0.79078251,0.7847739,0.8006355970000001,0.79086886,0.7893245000000001 +29.300,0.8079208600000001,0.7892818,0.7873482,0.799631187,0.7967853100000001,0.7890288000000001 +29.350,0.8057648000000001,0.7876358,0.7898748,0.79786151,0.8023319800000001,0.7887824 +29.400,0.8034127400000001,0.7861835,0.79297357,0.79590052,0.80755509,0.7889302 +29.450,0.8011096900000001,0.7848765000000001,0.7960723000000001,0.79379604,0.8113915,0.7890781 +29.500,0.7992476390000001,0.7837146,0.799647761,0.7920263700000001,0.8138875000000001,0.7896202 +29.550,0.79689558,0.782698,0.8028895100000001,0.7904480100000001,0.8143035000000001,0.7901622300000001 +29.600,0.79483753,0.7820203,0.8064172900000001,0.7893001000000001,0.8129169,0.79114783 +29.650,0.7927794800000001,0.7813425,0.80965904,0.7882,0.80963509,0.7919363 +29.700,0.7907214300000001,0.7809068,0.8129961,0.7871478000000001,0.80519775,0.79306973 +29.750,0.7884674,0.7805195,0.8157612000000001,0.7859999000000001,0.7996510870000001,0.7941046 +29.800,0.7869973,0.7804227,0.8185262000000001,0.785139,0.79391953,0.79558299 +29.850,0.7856253000000001,0.7805195,0.8205761,0.7843737000000001,0.7880493000000001,0.7967657 +29.900,0.7844983000000001,0.7811005000000001,0.8226260000000001,0.7842781000000001,0.7831035000000001,0.79814553 +29.950,0.7834202,0.7817298,0.8238655,0.7842781000000001,0.7788973,0.799377521 +30.000,0.7826852000000001,0.7826496000000001,0.8249620000000001,0.7843737000000001,0.7761702,0.800855909 +30.050,0.7819992,0.7836662000000001,0.8254864000000001,0.7839911,0.7746449000000001,0.8019400600000001 +30.100,0.7814602,0.7850217,0.8257248,0.7834650000000001,0.7746449000000001,0.80317205 +30.150,0.7808722,0.7863288,0.8252480000000001,0.7827954,0.7756155,0.80420692 +30.200,0.7808722,0.7878779,0.8247236,0.7825084,0.778204,0.8052417900000001 +30.250,0.7808722,0.7894754,0.8235318,0.7822692,0.7818555,0.80588243 +30.300,0.7811662,0.79136344,0.8222923000000001,0.7822692,0.7866626000000001,0.80676946 +30.350,0.7817052,0.79310621,0.8205284,0.7819823,0.7921630900000001,0.8073608200000001 +30.400,0.7826362,0.79509103,0.8186215,0.7819344,0.7981257500000001,0.80790289 +30.450,0.7835672,0.79707585,0.8164762000000001,0.7822214000000001,0.8038110900000001,0.80829713 +30.500,0.7846453,0.799447949,0.8143786000000001,0.783178,0.80889553,0.80869137 +30.550,0.7855273,0.80196528,0.8120904,0.7841824000000001,0.8128244,0.80903632 +30.600,0.7869973,0.80457944,0.8100404,0.7849477,0.8161986,0.8093812800000001 +30.650,0.7885164,0.80695154,0.8077044600000001,0.7859521,0.8182786000000001,0.8096276800000001 +30.700,0.79028041,0.80932364,0.8054161700000001,0.7885827000000001,0.8197115,0.8100219000000001 +30.750,0.7919954600000001,0.8113569,0.80298486,0.79303078,0.8201738000000001,0.810219 +30.800,0.7941515100000001,0.8132449,0.8008872530000001,0.7991528960000001,0.8207746,0.8106626 +30.850,0.7962095600000001,0.8148908,0.79850361,0.80632726,0.8212831,0.8107611 +30.900,0.79821861,0.8165852,0.79654903,0.8151756,0.8223,0.8110075 +30.950,0.799884655,0.8180375000000001,0.79440376,0.8253154,0.8234555,0.8111553 +31.000,0.8019427100000001,0.8192477,0.79259219,0.8362682,0.8251195,0.8113525 +31.050,0.80429477,0.8200707,0.79054226,0.8454993000000001,0.8265986000000001,0.8111553 +31.100,0.8067448300000001,0.8207001,0.7887307,0.8507126,0.8277542,0.8111061 +31.150,0.80875388,0.8208453,0.7868715000000001,0.8496604000000001,0.8278466,0.8107118000000001 +31.200,0.8112529,0.8208453,0.7852983,0.8431556,0.8265062000000001,0.810219 +31.250,0.813311,0.8204096000000001,0.7836774000000001,0.8333028,0.8231782000000001,0.80928272 +31.300,0.8147810000000001,0.8198287000000001,0.7824379,0.8228761,0.8183711,0.8081985700000001 +31.350,0.8160551,0.8188121,0.7812461,0.8132625,0.8120386,0.80662162 +31.400,0.8177211000000001,0.8177954000000001,0.7803880000000001,0.80670989,0.8054750900000001,0.80509396 +31.450,0.8189461,0.8162463000000001,0.7793392,0.8035053400000001,0.79854175,0.80287637 +31.500,0.8197302000000001,0.8148908,0.7787671,0.799583358,0.7923479800000001,0.80070807 +31.550,0.8204162,0.8131965000000001,0.7781950000000001,0.7964266400000001,0.7867089,0.7980962500000001 +31.600,0.8213472,0.8116957,0.7781950000000001,0.7965223,0.7824564,0.7955337100000001 +31.650,0.8214452000000001,0.80975933,0.7780043000000001,0.7973832200000001,0.7790822000000001,0.7925769300000001 +31.700,0.8210532,0.80767769,0.7780520000000001,0.79824414,0.7773720000000001,0.7899158000000001 +31.750,0.8199752,0.80530559,0.7778136,0.79853112,0.7764013000000001,0.7870576 +31.800,0.8193381000000001,0.80317554,0.7783857000000001,0.799583358,0.7768173,0.7845936 +31.850,0.8183091,0.800803436,0.7788148,0.80106606,0.778204,0.7820311000000001 +31.900,0.8175251,0.7987217900000001,0.7797682,0.80264442,0.7807000000000001,0.7801092000000001 +31.950,0.8160551,0.79659174,0.7805310000000001,0.80297922,0.7835658000000001,0.7781873 +32.000,0.814389,0.7948005600000001,0.7819135,0.8028357300000001,0.7870786000000001,0.7771031 +32.050,0.8127230000000001,0.79291256,0.7833437000000001,0.8023574400000001,0.79068397,0.7762161000000001 +32.100,0.8108119,0.79136344,0.7852506,0.80216613,0.79442798,0.7759697000000001 +32.150,0.8079208600000001,0.7896207000000001,0.7871575000000001,0.80106606,0.7978946400000001,0.7761668 +32.200,0.8055688000000001,0.7882168,0.7897795000000001,0.799344213,0.8012226400000001,0.7771524000000001 +32.250,0.8036087500000001,0.7867161,0.79249685,0.79700059,0.8038110900000001,0.7785322 +32.300,0.8018937100000001,0.7854574000000001,0.7956909200000001,0.79518308,0.80593731,0.7805527000000001 +32.350,0.7995416470000001,0.7840051,0.7988373200000001,0.7934612400000001,0.80704664,0.7829674000000001 +32.400,0.7970915900000001,0.7829885000000001,0.80236511,0.7918828800000001,0.80741642,0.7860227000000001 +32.450,0.79478853,0.7819234,0.80574988,0.79020886,0.8066306400000001,0.7892752000000001 +32.500,0.79302448,0.7812941,0.80961137,0.7888696000000001,0.8049666400000001,0.7929711700000001 +32.550,0.7910154300000001,0.7806648,0.8127578000000001,0.7876739,0.8019622000000001,0.79661786 +32.600,0.7891534,0.7804227,0.8159995,0.7869565000000001,0.79798709,0.800510952 +32.650,0.7871443,0.7803259,0.8187645,0.7861434,0.7931337500000001,0.80415764 +32.700,0.7855763,0.7806164,0.8215296000000001,0.7859999000000001,0.7882804,0.80775505 +32.750,0.7841553,0.7810036,0.8233411,0.7859999000000001,0.7833809,0.8106626 +32.800,0.7832732,0.7819234,0.8250573,0.7861912000000001,0.7796369000000001,0.8130280000000001 +32.850,0.7823422,0.7825044000000001,0.8258678,0.7857129,0.7768173,0.8145064000000001 +32.900,0.7819012000000001,0.7835694000000001,0.8265352,0.7852346,0.7755693,0.8158369000000001 +32.950,0.7813132,0.7845376000000001,0.8262492,0.7841824000000001,0.7754769,0.8161819 +33.000,0.7811172000000001,0.7858447000000001,0.8258201000000001,0.7834171000000001,0.7769560000000001,0.8161819 +33.050,0.7808232,0.7870065,0.8245806,0.7826997,0.779452,0.8151963000000001 +33.100,0.7809702000000001,0.7887493000000001,0.8234841,0.7826997,0.7834733,0.8140136 +33.150,0.7810682,0.79044364,0.8219109,0.7824127000000001,0.7881418,0.8118453 +33.200,0.7817052,0.79257369,0.8204331,0.7823171,0.7938270900000001,0.8096276800000001 +33.250,0.7820972,0.79436487,0.8183355,0.7818866,0.7994661980000001,0.8070651400000001 +33.300,0.7831262,0.7968337900000001,0.8165239000000001,0.7822214000000001,0.8049666400000001,0.8047982800000001 +33.350,0.7841063,0.799109077,0.8142833,0.7826997,0.80935776,0.8024821400000001 +33.400,0.7854293,0.80177164,0.8123287,0.7836563000000001,0.8127782,0.80070807 +33.450,0.7866543,0.80380487,0.8100881000000001,0.7841824000000001,0.8147658000000001,0.7990818430000001 +33.500,0.7883204,0.80593492,0.80803817,0.7854738,0.8161062,0.7983919300000001 +33.550,0.7899374,0.80791974,0.80560686,0.7878174,0.8166609,0.79794841 +33.600,0.7920444600000001,0.8100982,0.8036046,0.79221768,0.8173542,0.79844121 +33.650,0.79390651,0.8118894000000001,0.8013639800000001,0.79771802,0.8177702,0.799131122 +33.700,0.79596456,0.8140678,0.7995047430000001,0.80470107,0.8186946,0.800855909 +33.750,0.7977286,0.8157622,0.79740714,0.8130233,0.8198502000000001,0.80272853 +33.800,0.8000316588,0.8174082,0.79550023,0.8235935000000001,0.8218840000000001,0.80534035 +33.850,0.8018937100000001,0.8184732,0.79335496,0.8351203,0.8238715000000001,0.80790289 +33.900,0.80414776,0.819393,0.7915433900000001,0.8458819000000001,0.8259978,0.8109089 +33.950,0.8063528200000001,0.8197803,0.7893504,0.8519562,0.8274306,0.8135208 +34.000,0.80919489,0.8198771,0.7875865000000001,0.8518127000000001,0.8281702000000001,0.8160833000000001 +34.050,0.8111059,0.8194898,0.785632,0.8446862,0.8269684,0.8175617000000001 +34.100,0.812968,0.8192477,0.7842971000000001,0.8340681000000001,0.8242875000000001,0.8191386 +34.150,0.814536,0.8186668,0.7828193,0.8224935,0.8196653,0.8197793 +34.200,0.8162511,0.8181343000000001,0.7818181000000001,0.8124494,0.8140262,0.8198778 +34.250,0.8172311000000001,0.8169241,0.7806263,0.8052750200000001,0.8070928700000001,0.818843 +34.300,0.8187501,0.815859,0.7800066,0.80221396,0.8000670871000001,0.8172167 +34.350,0.8198282,0.8143583000000001,0.7791485,0.7982919700000001,0.7928102,0.8142107000000001 +34.400,0.8208572000000001,0.812906,0.7788148,0.79542223,0.7868013,0.8106133000000001 +34.450,0.8210042000000001,0.811018,0.7782427000000001,0.7952787400000001,0.781532,0.8060795500000001 +34.500,0.8212002,0.80927523,0.7782427000000001,0.7965223,0.7780191000000001,0.80144726 +34.550,0.8212982000000001,0.80714518,0.7780997000000001,0.79714408,0.7754769,0.79627291 +34.600,0.8209062,0.805354,0.7784334,0.7981963200000001,0.7745524,0.79134494 +34.650,0.8195341,0.80303031,0.7787194000000001,0.799726845,0.7746449000000001,0.7863677 +34.700,0.8187011000000001,0.80104549,0.7795775,0.8023574400000001,0.7764475000000001,0.7820311000000001 +34.750,0.8176231,0.79881862,0.7801496,0.80436626,0.7791284,0.7780394 +34.800,0.8164471000000001,0.79702744,0.7812937000000001,0.80565765,0.7830573000000001,0.7749841000000001 +34.850,0.814536,0.7949942000000001,0.7823902,0.8054185,0.7871711,0.7721752000000001 +34.900,0.812772,0.79329985,0.7842018000000001,0.80517936,0.7918395300000001,0.7704504000000001 +34.950,0.8104199000000001,0.79121821,0.7858703,0.8044140900000001,0.7961382,0.7693662 +35.000,0.80850887,0.7896207000000001,0.7880633,0.8032183700000001,0.8003906430000001,0.7697112 +35.050,0.8058138,0.7879263000000001,0.79020855,0.8004921100000001,0.80390353,0.7706968000000001 +35.100,0.80346174,0.7867645000000001,0.7931642600000001,0.79757454,0.80686175,0.7729636000000001 +35.150,0.80135469,0.7853606000000001,0.7961199800000001,0.79494394,0.80852575,0.7758219000000001 +35.200,0.799737651,0.7842471000000001,0.799600089,0.79312643,0.80949642,0.7796657 +35.250,0.79728759,0.7829400000000001,0.8026511500000001,0.7912132700000001,0.80926531,0.783953 +35.300,0.7948865300000001,0.7822139,0.8060835900000001,0.7896349,0.8082022000000001,0.7891274 +35.350,0.79248547,0.7813909,0.80927766,0.7876739,0.80575242,0.79449884 +35.400,0.79096643,0.7811973000000001,0.8128531,0.7866217000000001,0.80256309,0.800412392 +35.450,0.7889574,0.7806648,0.8159518,0.7855694000000001,0.7980795300000001,0.80598099 +35.500,0.7871443,0.7805195,0.8190983000000001,0.7849955000000001,0.7933186400000001,0.8115496000000001 +35.550,0.7852333,0.7804227,0.8214342,0.7845650000000001,0.788188,0.8163790000000001 +35.600,0.7841553,0.7809068,0.8237225,0.7847563000000001,0.7837506000000001,0.8202721000000001 +35.650,0.7829792000000001,0.7812457,0.8250097000000001,0.7847085,0.7796831000000001,0.8223911 +35.700,0.7822442000000001,0.7822139,0.8262015,0.784852,0.7770484000000001,0.8238202 +35.750,0.7813132,0.7830369,0.8264875,0.7845172,0.7754769,0.8236231 +35.800,0.7807742000000001,0.7842955,0.8267259,0.7840389000000001,0.7755693,0.8226868 +35.850,0.7801862,0.7853606000000001,0.8261062,0.7827954,0.7769098000000001,0.820075 +35.900,0.7802352,0.7868129,0.8254864000000001,0.7820779,0.7799604,0.8166747000000001 +35.950,0.7801372000000001,0.7882652,0.8241039,0.7815518,0.7839355,0.8122395 +36.000,0.7804802000000001,0.79025,0.8227691,0.7815996000000001,0.7891124,0.80770577 +36.050,0.7807252,0.792138,0.8209575,0.781217,0.79452042,0.8031227700000001 +36.100,0.7816072000000001,0.79446169,0.8192413000000001,0.7810735000000001,0.8003444200000001,0.799131122 +36.150,0.7824402,0.79634969,0.8170483000000001,0.781217,0.8054750900000001,0.79573083 +36.200,0.7835672,0.79877021,0.8149507,0.7821736,0.8100049,0.79331613 +36.250,0.7845473000000001,0.800948667,0.8124241000000001,0.782891,0.8131018000000001,0.7914435000000001 +36.300,0.7858213000000001,0.80351441,0.8103741000000001,0.7835606,0.8152742000000001,0.7911971000000001 +36.350,0.7869973,0.80554764,0.80808584,0.7842302000000001,0.8163373,0.79188702 +36.400,0.7887614000000001,0.80791974,0.80598824,0.7868130000000001,0.8171693,0.7940060400000001 +36.450,0.79081943,0.80985616,0.8035569300000001,0.79111761,0.8174466,0.7967657 +36.500,0.7931714900000001,0.811841,0.80155467,0.7971919000000001,0.8179551,0.8006587900000001 +36.550,0.79498453,0.8135837,0.799361725,0.80422278,0.8186022000000001,0.8049461200000001 +36.600,0.79694458,0.815617,0.7975978300000001,0.8130712,0.8200351,0.8100219000000001 +36.650,0.79875763,0.8170693,0.7954525600000001,0.8233544,0.8219302,0.814802 +36.700,0.8010606800000001,0.8185216000000001,0.79349797,0.8346421,0.8242875000000001,0.8194836000000001 +36.750,0.8032167400000001,0.8191993000000001,0.7913050300000001,0.8444470000000001,0.8260902,0.8230810000000001 +36.800,0.8057648000000001,0.8198287000000001,0.7895888000000001,0.8504257000000001,0.8276618,0.8253479 +36.850,0.8078718500000001,0.8200223000000001,0.7878249,0.8503778000000001,0.8281702000000001,0.8263335 +36.900,0.8101749,0.8200707,0.7863947,0.8444470000000001,0.827708,0.8261363 +36.950,0.812184,0.8197318,0.7847739,0.8342116,0.8254431000000001,0.8245594 +37.000,0.814046,0.8194414000000001,0.7835344000000001,0.8234022000000001,0.8215604000000001,0.8220461 +37.050,0.815369,0.8186668,0.7821995,0.8131668000000001,0.8160138,0.8181038 +37.100,0.8170351,0.8176502,0.7812937000000001,0.80613594,0.80958887,0.8134222000000001 +37.150,0.8184071,0.8161495000000001,0.7802449,0.80293139,0.80247064,0.80824785 +37.200,0.8195341,0.814794,0.7795299,0.7994877,0.79581464,0.8035170100000001 +37.250,0.8200242000000001,0.8131481,0.7789578,0.7965223,0.7895746,0.7987861700000001 +37.300,0.8205142000000001,0.8115989,0.7786717000000001,0.7966179600000001,0.7846289000000001,0.7945481200000001 +37.350,0.8207102000000001,0.80975933,0.7781950000000001,0.7976702,0.7806538000000001,0.7904579100000001 +37.400,0.8208082000000001,0.80801656,0.7782904,0.7987224300000001,0.7781578,0.7874026000000001 +37.450,0.8203672000000001,0.80593492,0.7783380000000001,0.7992007250000001,0.7766786,0.7849386 +37.500,0.8201222,0.80404692,0.7789101,0.8004921100000001,0.7769560000000001,0.7833616 +37.550,0.8190931,0.80186846,0.7791961000000001,0.80221396,0.7782502,0.7820311000000001 +37.600,0.8178191,0.7998836410000001,0.7801496,0.80407929,0.7809773000000001,0.7812426 +37.650,0.8161041,0.79765677,0.7810554000000001,0.80455758,0.7843053000000001,0.7805034000000001 +37.700,0.81483,0.79581718,0.7825809,0.80422278,0.7883729,0.7805527000000001 +37.750,0.8130660000000001,0.79373554,0.7841064000000001,0.8033140200000001,0.79262531,0.7807005 +37.800,0.8111059,0.79208959,0.786061,0.80288356,0.79720131,0.7812919 +37.850,0.80865587,0.79020159,0.7879203,0.80183132,0.8014537500000001,0.7816369000000001 +37.900,0.8066468200000001,0.7886041,0.7903039000000001,0.7999181613,0.80561375,0.7823268000000001 +37.950,0.8044417700000001,0.7868613000000001,0.7927828800000001,0.7973353900000001,0.8090342,0.7829674000000001 +38.000,0.80253072,0.7856510000000001,0.7958339400000001,0.7953744,0.8118075,0.7841008 +38.050,0.799884655,0.7842471000000001,0.79878965,0.7934134100000001,0.8134253,0.7852343 +38.100,0.79743459,0.7833273000000001,0.80198373,0.7918350500000001,0.8140262,0.7869591 +38.150,0.79498453,0.7822623000000001,0.80508246,0.79016103,0.8128244,0.7890288000000001 +38.200,0.79292648,0.781633,0.80865792,0.7891088,0.8101898000000001,0.7917391800000001 +38.250,0.7907214300000001,0.7810520000000001,0.8120904,0.7881044,0.80602975,0.7949916300000001 +38.300,0.7889574,0.78081,0.8154751,0.7875783000000001,0.800945309,0.798934 +38.350,0.7871933000000001,0.7804227,0.8183355,0.7867173000000001,0.7950750900000001,0.8029749300000001 +38.400,0.7857723000000001,0.7806164,0.8210052000000001,0.7862869,0.7893898,0.80731154 +38.450,0.7842043000000001,0.7807132,0.8229597000000001,0.7859999000000001,0.7839355,0.8113525 +38.500,0.7828812,0.7812457,0.8247236,0.7861434,0.7795444,0.8151963000000001 +38.550,0.7817052,0.7818266,0.8256771,0.7859521,0.7762626,0.8183995000000001 +38.600,0.7810192,0.7828916,0.8265829,0.7858086000000001,0.7746449000000001,0.8212084000000001 +38.650,0.7802352,0.7840051,0.8265829,0.7849955000000001,0.7745062,0.8225389000000001 +38.700,0.7799412,0.785409,0.8262968,0.7840867,0.7758466,0.8224897 +38.750,0.7796472000000001,0.7867645000000001,0.8251527000000001,0.7832258000000001,0.7783426,0.8203706000000001 +38.800,0.7799902000000001,0.7884588,0.8238655,0.7830345000000001,0.7825489,0.8172167 +38.850,0.7802842000000001,0.79015318,0.8221016000000001,0.7827954,0.7874022,0.8128801 +38.900,0.7809212000000001,0.792138,0.8203854,0.7826519000000001,0.79317998,0.80839569 +38.950,0.7812642000000001,0.79392918,0.8182878,0.7822214000000001,0.799050198,0.80346773 +39.000,0.7823422,0.79625287,0.8164762000000001,0.7825084,0.8049204200000001,0.798934 +39.050,0.7833222000000001,0.79857656,0.8144263,0.7830345000000001,0.8098662000000001,0.7945974 +39.100,0.7847433,0.80138436,0.8124717,0.7838476000000001,0.8138875000000001,0.7914435000000001 +39.150,0.7856743,0.80385328,0.8102311000000001,0.7841824000000001,0.8163835,0.7893738 +39.200,0.7869483,0.80637062,0.80822886,0.7850912000000001,0.8180475,0.788881 +39.250,0.7884674,0.80840385,0.8060359100000001,0.7872435000000001,0.8188333000000001,0.7895709000000001 +39.300,0.79077043,0.8105339,0.8040336600000001,0.79164373,0.819388,0.7917391800000001 +39.350,0.79302448,0.8123251,0.8018883800000001,0.79714408,0.8194804,0.79504091 +39.400,0.79532754,0.8142615000000001,0.79979078,0.8036488300000001,0.8198964000000001,0.799870317 +39.450,0.7972385900000001,0.8158106,0.7975024900000001,0.8116363000000001,0.8204511000000001,0.80529107 +39.500,0.799198638,0.8175534000000001,0.7954048800000001,0.8219195,0.8216991,0.8111553 +39.550,0.8009626820000001,0.8188605,0.7931642600000001,0.8333985,0.8230858000000001,0.8165268 +39.600,0.80316774,0.8199255000000001,0.7912573500000001,0.8444470000000001,0.8247960000000001,0.8216519 +39.650,0.80542179,0.8204096000000001,0.7893028000000001,0.8517649,0.8262751,0.8251015 +39.700,0.80801886,0.8206516,0.7877772000000001,0.8532476,0.8275231000000001,0.8269248 +39.750,0.8102239,0.8203612,0.785918,0.8480342000000001,0.8275693000000001,0.8269741 +39.800,0.812527,0.8201191000000001,0.7844878000000001,0.8385640000000001,0.8261364,0.8254464 +39.850,0.8141930000000001,0.8195866,0.7830576,0.8273242000000001,0.8227622,0.8220954 +39.900,0.8159571,0.8190541,0.7820565,0.8168496000000001,0.8180013,0.8180052 +39.950,0.8171821,0.8179891,0.7809600000000001,0.8088622000000001,0.8117613,0.8127816 +40.000,0.8185541000000001,0.8169725000000001,0.7801973000000001,0.80517936,0.8048742000000001,0.8074101 +40.050,0.8193871,0.8153749,0.7793868,0.8008269140000001,0.79775598,0.80203862 +40.100,0.8207102000000001,0.8140194000000001,0.7789578,0.79810066,0.7914235300000001,0.7975049000000001 +40.150,0.8210042000000001,0.812083,0.7783857000000001,0.7980050000000001,0.7858769000000001,0.7934146900000001 +40.200,0.8210042000000001,0.8102918,0.7781473000000001,0.7989137500000001,0.7818555,0.7901622300000001 +40.250,0.8209552,0.80821021,0.7779090000000001,0.799439871,0.7788511,0.7875997 +40.300,0.8209552,0.80646744,0.7780997000000001,0.8008269140000001,0.7774644,0.7859242000000001 +40.350,0.8203182,0.80433738,0.7783380000000001,0.80293139,0.7772333,0.7846922000000001 +40.400,0.8193871,0.8024978,0.7790531,0.80532285,0.7783889,0.7840516000000001 +40.450,0.8181621,0.800416154,0.7797205,0.8065664,0.7805613,0.7835588 +40.500,0.8170841000000001,0.79857656,0.7807693,0.80718818,0.7838893,0.7831645 +40.550,0.81532,0.79649492,0.7818658,0.8074273200000001,0.7874946,0.7826717000000001 +40.600,0.813409,0.79470374,0.783439,0.80761864,0.7916084200000001,0.7822775000000001 +40.650,0.8111059,0.79252528,0.7850599,0.8068055500000001,0.7957222,0.7815383 +40.700,0.8089988800000001,0.7907341,0.7871575000000001,0.8052271900000001,0.8002519760000001,0.7811441 +40.750,0.80669583,0.7888945,0.7894458000000001,0.80264442,0.80441198,0.7806513 +40.800,0.8047847800000001,0.7873454000000001,0.7922584800000001,0.8000138195000001,0.8082484200000001,0.7807005 +40.850,0.80272673,0.7858931,0.7951665200000001,0.79747888,0.8108831000000001,0.7807005 +40.900,0.800521671,0.7846828,0.79840827,0.7955178900000001,0.812732,0.7814890000000001 +40.950,0.7977776000000001,0.7833757,0.8015070000000001,0.79350907,0.8131480000000001,0.7827703 +41.000,0.79562155,0.7824075,0.80484409,0.79197854,0.812316,0.7852836000000001 +41.050,0.79326949,0.7814393000000001,0.80808584,0.7905914900000001,0.8099124200000001,0.7885853 +41.100,0.7911134400000001,0.7809068,0.8116136,0.7894914000000001,0.80626087,0.7931190100000001 +41.150,0.7889574,0.7802775000000001,0.8145693,0.7882479,0.80131509,0.79814553 +41.200,0.7874873,0.7802291,0.8175251,0.7875783000000001,0.79581464,0.80386197 +41.250,0.7860173,0.7800838,0.8199087,0.7872435000000001,0.7899444,0.8093812800000001 +41.300,0.7844003,0.7806164,0.8223877,0.7873869,0.7847213000000001,0.8147528 +41.350,0.7829792000000001,0.7810520000000001,0.8242946,0.7872435000000001,0.7799604,0.8191386 +41.400,0.7820482000000001,0.7819234,0.8258201000000001,0.7871956,0.7767249,0.8227361000000001 +41.450,0.7812642000000001,0.7825528,0.8263922,0.7866217000000001,0.7745986,0.8245594 +41.500,0.7808722,0.7836662000000001,0.8267259,0.7858086000000001,0.7740902000000001,0.8248551000000001 +41.550,0.7802842000000001,0.7847796,0.8262015,0.7848042000000001,0.774876,0.8228346000000001 +41.600,0.7801372000000001,0.7862319000000001,0.8255341,0.7843259,0.7773258000000001,0.8194343000000001 +41.650,0.7799412,0.7879747,0.8241516000000001,0.7837998,0.7810235000000001,0.8146542 +41.700,0.7803822,0.79015318,0.8227691,0.7833693,0.7860155000000001,0.8094798400000001 +41.750,0.7808722,0.79223482,0.8208621,0.7827475,0.79151598,0.8039605200000001 +41.800,0.7818522,0.79421964,0.8191459000000001,0.7826519000000001,0.79766353,0.799131122 +41.850,0.7826362,0.79610764,0.817239,0.7827475,0.8034875300000001,0.7949423600000001 +41.900,0.7837143000000001,0.79852815,0.8153798000000001,0.7834171000000001,0.80875686,0.79203486 +41.950,0.7847923,0.800997077,0.8131392000000001,0.7834171000000001,0.8125471000000001,0.7898173 +42.000,0.7863603,0.80370805,0.8110892000000001,0.7835128,0.8151355,0.7889302 +42.050,0.7878794,0.80603174,0.8088486100000001,0.7841824000000001,0.8162911,0.7890781 +42.100,0.7897414,0.80859749,0.80698937,0.7867173000000001,0.8167533,0.7908521500000001 +42.150,0.7912114400000001,0.8107275,0.8049394400000001,0.7905914900000001,0.8163373,0.79380892 +42.200,0.7933184900000001,0.8125671000000001,0.80298486,0.7957570300000001,0.8161986,0.7977020100000001 +42.250,0.79518054,0.813971,0.8006965620000001,0.8015921800000001,0.8162449,0.80213718 +42.300,0.7972385900000001,0.8155686000000001,0.7984559400000001,0.8091491700000001,0.8169382000000001,0.80726226 +42.350,0.7993946430000001,0.8169725000000001,0.7959769600000001,0.8185715,0.8182324000000001,0.8121902000000001 +42.400,0.8016977000000001,0.8185216000000001,0.79383168,0.8299070000000001,0.8204511000000001,0.8168718 +42.450,0.8036087500000001,0.8196834000000001,0.7916864100000001,0.841099,0.822716,0.8206663000000001 +42.500,0.8055688000000001,0.8205064000000001,0.7900178600000001,0.8497560000000001,0.8249346000000001,0.8238695 +42.550,0.8075288500000001,0.8205064000000001,0.7883016,0.8531041,0.8264138000000001,0.8255943 +42.600,0.8100769000000001,0.8202159,0.7868238000000001,0.8503778000000001,0.8271071000000001,0.8260378 +42.650,0.81238,0.8195866,0.7852983,0.8421512,0.8260902,0.8249044000000001 +42.700,0.814389,0.8191993000000001,0.7841064000000001,0.8315332000000001,0.8233631,0.8226375 +42.750,0.8158591000000001,0.8187636,0.7828669,0.820389,0.8185560000000001,0.8186458000000001 +42.800,0.8175251,0.8183279,0.7821042,0.8112536,0.8124084,0.8138657 +42.850,0.8183091,0.8172629,0.7811507000000001,0.80580114,0.8052902000000001,0.80829713 +42.900,0.8193871,0.8160043,0.7805787000000001,0.80264442,0.7980795300000001,0.80287637 +42.950,0.8204652,0.813971,0.7798159,0.7992007250000001,0.79137731,0.79760345 +43.000,0.8209552,0.8122283,0.7791961000000001,0.79838763,0.7857382,0.7930204500000001 +43.050,0.8211512000000001,0.8103403,0.7786241,0.79853112,0.7811622,0.7889302 +43.100,0.8211512000000001,0.8086459,0.7784810000000001,0.79881809,0.7780191000000001,0.7858749 +43.150,0.8204652,0.80656426,0.7785764000000001,0.799774674,0.7761702,0.7837559000000001 +43.200,0.8203672000000001,0.80457944,0.7791008,0.80197481,0.7758466,0.7827703 +43.250,0.8197792,0.80240097,0.7796252,0.80384014,0.7766786,0.7824253000000001 +43.300,0.8185541000000001,0.800367744,0.7805310000000001,0.8050358700000001,0.7788973,0.7826225 +43.350,0.8173291,0.7982861,0.7814368,0.8057533100000001,0.7820866000000001,0.7828689000000001 +43.400,0.8156631000000001,0.79634969,0.7828193,0.8066142300000001,0.7860155000000001,0.7832138000000001 +43.450,0.8135070000000001,0.79436487,0.7842971000000001,0.80699686,0.7902679800000001,0.7833616 +43.500,0.8115959,0.79257369,0.7862517,0.8065664,0.7947977500000001,0.7834602 +43.550,0.8094889000000001,0.79063728,0.7882063,0.80494021,0.799281309,0.7831153000000001 +43.600,0.80723484,0.7888461,0.79058993,0.80288356,0.8035799800000001,0.7827703 +43.650,0.80488278,0.7870549,0.7930689200000001,0.8003964520000001,0.80723153,0.7820804 +43.700,0.8028247300000001,0.7855542000000001,0.7960723000000001,0.79790934,0.8101898000000001,0.7816861 +43.750,0.800717676,0.7841987,0.7990756880000001,0.7954700600000001,0.8120386,0.7812919 +43.800,0.79875763,0.7831337,0.80226976,0.79350907,0.8128244,0.781834 +43.850,0.79596456,0.7819234,0.80532082,0.79173939,0.8121773,0.7832138000000001 +43.900,0.7933184900000001,0.7810520000000001,0.80861024,0.79035235,0.8102360000000001,0.7859242000000001 +43.950,0.79130944,0.7802291,0.8117566,0.7887740000000001,0.80686175,0.7898666000000001 +44.000,0.7897414,0.7798418,0.8149984,0.7876261,0.80251687,0.79509019 +44.050,0.7880264,0.7795513000000001,0.8180971,0.7866695,0.7973399800000001,0.80115159 +44.100,0.7862623000000001,0.7797934000000001,0.8210528,0.7860956,0.7916546400000001,0.80790289 +44.150,0.7846943000000001,0.7801807000000001,0.8236272,0.7856173000000001,0.7860155000000001,0.8143585000000001 +44.200,0.7834692000000001,0.78081,0.8258678,0.7855216,0.7810698,0.8201735 +44.250,0.7823422,0.7814393000000001,0.8275363,0.7852346,0.7770022000000001,0.8245101 +44.300,0.7812152,0.7823107,0.8287282,0.7846607000000001,0.77446,0.8268263 +44.350,0.7803822,0.7833273000000001,0.8290142,0.7837998,0.7732582,0.8267770000000001 +44.400,0.7800392,0.7844408,0.8288235,0.783178,0.7737666000000001,0.824658 +44.450,0.7799412,0.7858931,0.82787,0.7827475,0.7756618000000001,0.8201242000000001 +44.500,0.7800882,0.7877811,0.8265352,0.7822692,0.7790822000000001,0.8145556 +44.550,0.7802842000000001,0.7898627,0.8248667000000001,0.7816953000000001,0.7837044000000001,0.8083464100000001 +44.600,0.7808722,0.79194436,0.8231981,0.7811692,0.7893898,0.80228502 +44.650,0.7815582000000001,0.79383236,0.8212435,0.7807865,0.7953524200000001,0.7970121 +44.700,0.7825872,0.79586559,0.8192889000000001,0.7809779,0.80131509,0.7931190100000001 +44.750,0.7836162000000001,0.79799564,0.8171913000000001,0.7813127000000001,0.80649198,0.79026079 +44.800,0.7850863,0.800367744,0.8151414,0.7817909000000001,0.8107906,0.7885853 +44.850,0.7863113,0.80278826,0.8131392000000001,0.7819823,0.8137489,0.7881910000000001 +44.900,0.7875363000000001,0.80520877,0.8112322000000001,0.7825562,0.8154591,0.7892752000000001 +44.950,0.7890064,0.80753246,0.80918232,0.7842302000000001,0.8161062,0.79164062 +45.000,0.79096643,0.80966251,0.8069417000000001,0.7874348,0.8161986,0.79504091 +45.050,0.79326949,0.8116473000000001,0.80446271,0.7920742000000001,0.8162449,0.798934 +45.100,0.7954745400000001,0.8135353000000001,0.80193605,0.79810066,0.8163835,0.80317205 +45.150,0.7973855900000001,0.8152781,0.79945707,0.80532285,0.8169382000000001,0.80775505 +45.200,0.799198638,0.8168756,0.79716878,0.8144104000000001,0.8181400000000001,0.8125352 +45.250,0.8011096900000001,0.8182311,0.79511885,0.8253154,0.8198964000000001,0.8168225 +45.300,0.80316774,0.8192477,0.7932119400000001,0.8367465000000001,0.8218840000000001,0.8207156 +45.350,0.80522579,0.8199739,0.7913050300000001,0.8463124000000001,0.8236866,0.8239188000000001 +45.400,0.8076758500000001,0.8203612,0.7893981,0.8508083000000001,0.8250733,0.8261856000000001 +45.450,0.80992991,0.8203612,0.7875389,0.8492299000000001,0.8257204,0.8267770000000001 +45.500,0.8117920000000001,0.8199739,0.7858703,0.8420077,0.8248884000000001,0.8259392000000001 +45.550,0.813605,0.819393,0.7843925,0.8316767,0.8223462,0.8232781 +45.600,0.815467,0.8186668,0.7832483,0.8209629,0.8180475,0.8191386 +45.650,0.8167411,0.8179891,0.7822949,0.8117319000000001,0.8124546,0.8138164 +45.700,0.8179661,0.8171661,0.7816275,0.8059924500000001,0.80593731,0.8078043300000001 +45.750,0.8192401,0.8160043,0.7809124000000001,0.80297922,0.79909642,0.8018415000000001 +45.800,0.8200242000000001,0.8144551000000001,0.7801973000000001,0.79886592,0.79257909,0.79632218 +45.850,0.8204162,0.8126639,0.7794822,0.79700059,0.7868475,0.79139422 +45.900,0.8204652,0.8105823000000001,0.7787671,0.79762237,0.7821329,0.7874026000000001 +45.950,0.8203182,0.80874272,0.7784334,0.79853112,0.7787124000000001,0.7844951 +46.000,0.8204162,0.80680631,0.7784334,0.799296384,0.7767711,0.782721 +46.050,0.8199262,0.80482149,0.7789578,0.8007790840000001,0.7760778,0.7821789 +46.100,0.8190931,0.8029819,0.7796729,0.80259659,0.7766786,0.7822775000000001 +46.150,0.8181621,0.800900257,0.7803403000000001,0.8043184400000001,0.7783889,0.7828196000000001 +46.200,0.8169861,0.7987217900000001,0.7810554000000001,0.8053706700000001,0.7812084,0.7833124 +46.250,0.8157121,0.79668856,0.7819612,0.80584897,0.7848138,0.7837066 +46.300,0.814046,0.79475215,0.7833437000000001,0.8062316,0.78902,0.783953 +46.350,0.81189,0.79281574,0.7851076,0.8062316,0.7935035300000001,0.783953 +46.400,0.8094889000000001,0.79107298,0.7872052,0.8052271900000001,0.7980795300000001,0.7837066 +46.450,0.8070388300000001,0.7894754,0.7895888000000001,0.80312271,0.8023782,0.7832138000000001 +46.500,0.80507878,0.7877811,0.7921631400000001,0.8006355970000001,0.80626087,0.7825732000000001 +46.550,0.80297173,0.7861835,0.79497583,0.79838763,0.8094502,0.7820311000000001 +46.600,0.800619673,0.7846828,0.7978362,0.79618749,0.8117151,0.781834 +46.650,0.79831662,0.7833757,0.800934926,0.79384387,0.8128706,0.7824253000000001 +46.700,0.7963075700000001,0.7823591000000001,0.8041290000000001,0.79193071,0.812732,0.7839037000000001 +46.750,0.79420051,0.7817782000000001,0.8075614400000001,0.7905436600000001,0.8111142,0.786417 +46.800,0.79175045,0.7811973000000001,0.8111369,0.7892523,0.80797109,0.7899651000000001 +46.850,0.7896434,0.78081,0.8146647,0.7881044,0.8035799800000001,0.7946959600000001 +46.900,0.7875853,0.7804227,0.8178588,0.7870043,0.7982182,0.800313833 +46.950,0.7858703,0.7803743000000001,0.8207668,0.7861434,0.79248664,0.8061781100000001 +47.000,0.7844003,0.7806164,0.8232458,0.7854738,0.7867551,0.8118453 +47.050,0.7833222000000001,0.7811005000000001,0.8253911,0.7852346,0.7817631,0.8170689000000001 +47.100,0.7822932,0.7817298,0.8270596,0.785139,0.7776031,0.8214548 +47.150,0.7814602,0.7826496000000001,0.8280607000000001,0.7847563000000001,0.774876,0.8244116 +47.200,0.7807742000000001,0.783763,0.8284421000000001,0.7838476000000001,0.7733044,0.8253479 +47.250,0.7803332000000001,0.7850701,0.8282991000000001,0.7829388,0.7733506,0.8242637 +47.300,0.7798432000000001,0.7865708,0.8276794000000001,0.7821257,0.7746911000000001,0.8207649 +47.350,0.7795001,0.7883136000000001,0.8267259,0.7816475,0.7778342,0.8159355 +47.400,0.7795001,0.79020159,0.8252480000000001,0.7811213,0.7821329,0.8102683 +47.450,0.7803332000000001,0.79238005,0.8235795,0.7806430000000001,0.7876333000000001,0.80455188 +47.500,0.7815582000000001,0.79436487,0.8215772,0.7802126,0.7935035300000001,0.7992789610000001 +47.550,0.7827832,0.79649492,0.819575,0.7802604000000001,0.799604865,0.7949423600000001 +47.600,0.7838123,0.79857656,0.8173344,0.7804996000000001,0.8049666400000001,0.79149278 +47.650,0.7849393,0.800997077,0.8152844,0.78093,0.8094502,0.7891274 +47.700,0.7862623000000001,0.80327236,0.8129961,0.781217,0.812732,0.7879939 +47.750,0.7877324,0.80569287,0.8108509,0.7816953000000001,0.8148120000000001,0.7884374000000001 +47.800,0.7890554000000001,0.80791974,0.80846722,0.7826997,0.8156902,0.79021151 +47.850,0.7906724300000001,0.8102434000000001,0.80613126,0.7854738,0.8158751000000001,0.79336541 +47.900,0.79263247,0.8124219,0.8036046,0.79001754,0.8158751000000001,0.7972585000000001 +47.950,0.79498453,0.8143583000000001,0.80117329,0.7959483500000001,0.8162911,0.8019400600000001 +48.000,0.79718959,0.8161495000000001,0.7985989600000001,0.80293139,0.8170769,0.8066709000000001 +48.050,0.799198638,0.8176018,0.79645369,0.8116363000000001,0.8184635,0.8116974 +48.100,0.800913681,0.8188605,0.79435608,0.8223022000000001,0.8202200000000001,0.8161819 +48.150,0.8030207300000001,0.8200707,0.7926398600000001,0.8342116,0.8223462,0.8202721000000001 +48.200,0.80522579,0.8209421,0.79073295,0.8443514000000001,0.8243338,0.823426 +48.250,0.80738184,0.8215714000000001,0.7889214,0.8503778000000001,0.8259053000000001,0.8259392000000001 +48.300,0.8092438900000001,0.8215714000000001,0.7867761000000001,0.8497082,0.8264138000000001,0.8266784 +48.350,0.8114979,0.8213294,0.7850122,0.8432513,0.8256742,0.8262842 +48.400,0.813311,0.8206516,0.783296,0.8328246,0.8232244000000001,0.8241159 +48.450,0.815173,0.8201191000000001,0.7821519,0.8218239,0.8190644,0.8204692 +48.500,0.8168881,0.8191993000000001,0.7812461,0.8121146,0.8134715,0.8151963000000001 +48.550,0.8183091,0.8183764,0.7803880000000001,0.80594462,0.8071390900000001,0.8090856000000001 +48.600,0.8189461,0.8172145000000001,0.7793868,0.80293139,0.800205754,0.80272853 +48.650,0.8197302000000001,0.8158106,0.7785287000000001,0.7990572380000001,0.7935959800000001,0.7967657 +48.700,0.8204652,0.8138258,0.7775276000000001,0.79647447,0.7876333000000001,0.79139422 +48.750,0.8208572000000001,0.812083,0.7769555,0.7970484200000001,0.78278,0.787304 +48.800,0.8207102000000001,0.8101466,0.7766695,0.7982919700000001,0.7787124000000001,0.784298 +48.850,0.8210042000000001,0.80845226,0.7770032,0.7992007250000001,0.7764938,0.7826225 +48.900,0.8206612000000001,0.80651585,0.7774322,0.8003486230000001,0.7753844000000001,0.7818833000000001 +48.950,0.8200242000000001,0.80467626,0.7780520000000001,0.80264442,0.7758004000000001,0.7822282 +49.000,0.8189461,0.80254621,0.7783857000000001,0.80470107,0.7773258000000001,0.7829181000000001 +49.050,0.8176721,0.800464564,0.7791961000000001,0.80613594,0.7801915,0.7839037000000001 +49.100,0.8160061000000001,0.79838292,0.7801496,0.8063750900000001,0.7838431,0.7844951 +49.150,0.814879,0.79668856,0.7816275,0.8066142300000001,0.7880955000000001,0.7850372 +49.200,0.813311,0.7948005600000001,0.783296,0.8065664,0.7923942,0.7848893 +49.250,0.8109589,0.79315462,0.7853459,0.80589679,0.7969702000000001,0.7844458000000001 +49.300,0.80801886,0.79126662,0.7875389,0.80388797,0.8011302,0.7834109 +49.350,0.80591181,0.7896207000000001,0.79011321,0.80140086,0.8051515300000001,0.7824746 +49.400,0.80395176,0.7877327000000001,0.79273521,0.7986746100000001,0.8084333100000001,0.7813412000000001 +49.450,0.8017957,0.7861835,0.79588161,0.79637881,0.8111142,0.7807498 +49.500,0.799345642,0.784586,0.7988850000000001,0.7941786700000001,0.8126395000000001,0.7807498 +49.550,0.7969935800000001,0.7836178,0.80231744,0.792409,0.8131480000000001,0.7820804 +49.600,0.79459252,0.7826496000000001,0.8057022,0.79068715,0.8119924000000001,0.7844951 +49.650,0.79248547,0.7822139,0.80951603,0.7893479,0.80940398,0.7883389000000001 +49.700,0.79042742,0.7816814000000001,0.8130915000000001,0.7880087,0.8054288700000001,0.7932668500000001 +49.750,0.7887614000000001,0.7815846000000001,0.8166193,0.7870043,0.8003444200000001,0.7990818430000001 +49.800,0.7867523000000001,0.7813425,0.819575,0.7859042,0.7944742,0.8050446800000001 +49.850,0.7852823,0.7814877,0.8224830000000001,0.7853303,0.7886964000000001,0.8111061 +49.900,0.7842043000000001,0.7816814000000001,0.8247713000000001,0.7847563000000001,0.7832422,0.8165761 +49.950,0.7833712,0.7822623000000001,0.8267259,0.7846607000000001,0.7788049,0.8214548 +50.000,0.7822442000000001,0.7829400000000001,0.8276794000000001,0.7842302000000001,0.7754306000000001,0.8245594 +50.050,0.7814112000000001,0.7840535000000001,0.8281561000000001,0.7836563000000001,0.7737666000000001,0.8258407000000001 +50.100,0.7807252,0.7853606000000001,0.8277747000000001,0.7825084,0.7735355,0.8249044000000001 +50.150,0.7807252,0.7870065,0.8270596,0.7817431,0.7748298,0.8219969 +50.200,0.7807742000000001,0.7886525000000001,0.8261062,0.7811213,0.7775106,0.8170196000000001 +50.250,0.7809212000000001,0.79049205,0.8249620000000001,0.78093,0.781532,0.8114017 +50.300,0.7812642000000001,0.79228323,0.8232935,0.7804039,0.7864315000000001,0.80543891 +50.350,0.7823422,0.79426805,0.8217203000000001,0.7801648,0.79211686,0.7998210370000001 +50.400,0.7833712,0.79610764,0.8196227,0.7799734,0.79775598,0.79504091 +50.450,0.7848413000000001,0.79838292,0.8177157,0.7804039,0.8033026400000001,0.7916899000000001 +50.500,0.7860173,0.8006097950000001,0.8154274,0.7808822000000001,0.80773998,0.7892259 +50.550,0.7874383,0.80307872,0.8135205000000001,0.7814561000000001,0.8114840000000001,0.7880925000000001 +50.600,0.7887614000000001,0.80525718,0.8114706,0.7817431,0.8135178000000001,0.7879447 +50.650,0.7902314100000001,0.80762928,0.8092299900000001,0.7829867,0.8148582000000001,0.7894230000000001 +50.700,0.79155445,0.80971092,0.80689402,0.7855216,0.8151355,0.7921334200000001 +50.750,0.79326949,0.8117926000000001,0.80460573,0.7899697,0.8154129,0.79612507 +50.800,0.79518054,0.8134385000000001,0.80184071,0.7953265700000001,0.8155053000000001,0.8004616720000001 +50.850,0.7976796,0.8151329,0.799314052,0.8019269800000001,0.8162911,0.80558675 +50.900,0.7999336563,0.816682,0.79664438,0.80996226,0.817308,0.8107118000000001 +50.950,0.80218771,0.8181827,0.7944991,0.8204847000000001,0.8191569000000001,0.8159847 +51.000,0.8040007600000001,0.8192477,0.79235383,0.8320115,0.8210982,0.8204692 +51.050,0.8058138,0.8202159,0.79068528,0.8430121,0.8235018000000001,0.8243130000000001 +51.100,0.80757785,0.8207969,0.7886354000000001,0.850043,0.8253969000000001,0.8266291 +51.150,0.8100279,0.8211357,0.7869191000000001,0.8511909000000001,0.8266449,0.8277133000000001 +51.200,0.8120860000000001,0.8208453,0.7851076,0.8456906000000001,0.8263213,0.8272205 +51.250,0.8142910000000001,0.8206516,0.7839634000000001,0.8361248,0.8245186,0.8252493000000001 +51.300,0.8161531000000001,0.8199739,0.7827239,0.8246458000000001,0.8207284,0.8212084000000001 +51.350,0.8182601,0.8191993000000001,0.7816751000000001,0.8144104000000001,0.8155978,0.8159847 +51.400,0.8197302000000001,0.8180859,0.7803403000000001,0.8068055500000001,0.80926531,0.80943056 +51.450,0.8206122,0.8172629,0.7792438,0.80340968,0.80256309,0.8028270900000001 +51.500,0.8209062,0.815859,0.7781473000000001,0.7991050670000001,0.7956759800000001,0.79661786 +51.550,0.8217882000000001,0.8145035,0.7774322,0.79609184,0.7896209000000001,0.79149278 +51.600,0.8221802,0.8125187,0.7766695,0.79604401,0.7841666,0.7872055 +51.650,0.8223762,0.8107275,0.7764788,0.79723973,0.7799604,0.784298 +51.700,0.8225232,0.80874272,0.7764311,0.79786151,0.7768635,0.7824253000000001 +51.750,0.8224742,0.80699995,0.7769078,0.799248554,0.7753844000000001,0.7821789 +51.800,0.8210532,0.80482149,0.7771939,0.80149652,0.7751533,0.7824746 +51.850,0.8202202000000001,0.80303031,0.7777660000000001,0.80403146,0.77654,0.7834602 +51.900,0.8192401,0.800851846,0.7782427000000001,0.8054185,0.7788049,0.7843472 +51.950,0.8176231,0.79891544,0.7792915,0.80604028,0.7823640000000001,0.7853328 +52.000,0.8156631000000001,0.79688221,0.7806263,0.80632726,0.7863853000000001,0.7855792 +52.050,0.814389,0.7949942000000001,0.7826763,0.80670989,0.79082264,0.7856285000000001 +52.100,0.8122330000000001,0.7930093800000001,0.7847739,0.80604028,0.7951675300000001,0.7847908 +52.150,0.8094889000000001,0.79126662,0.7872052,0.8046054100000001,0.7996510870000001,0.783953 +52.200,0.8067448300000001,0.7893302000000001,0.7894458000000001,0.80230961,0.80367242,0.7826225 +52.250,0.80488278,0.7876358,0.7921631400000001,0.799822503,0.80750887,0.7815876 +52.300,0.80253072,0.7857963,0.7948328100000001,0.7971919000000001,0.8102822000000001,0.7806513 +52.350,0.800129661,0.7843924,0.79826525,0.79513526,0.8125009000000001,0.7807005 +52.400,0.7978756100000001,0.7828432000000001,0.80141165,0.7929829500000001,0.8131018000000001,0.7813412000000001 +52.450,0.79616056,0.7821171,0.80508246,0.7914045900000001,0.8125933000000001,0.7836080000000001 +52.500,0.79390651,0.7813425,0.8085625700000001,0.7897784,0.8103284000000001,0.7867619 +52.550,0.79170145,0.7810520000000001,0.8123764,0.7886305,0.8068155300000001,0.7912956600000001 +52.600,0.7892514,0.7806164,0.8158088,0.7872913,0.80186975,0.7965193 +52.650,0.7875363000000001,0.7807132,0.8193843000000001,0.7867173000000001,0.79650798,0.8026299800000001 +52.700,0.7857723000000001,0.7805679000000001,0.8221016000000001,0.7859521,0.7904990900000001,0.80844497 +52.750,0.7846943000000001,0.7809068,0.8248667000000001,0.7857129,0.7852760000000001,0.8143092000000001 +52.800,0.7834692000000001,0.7811005000000001,0.8264875,0.7853303,0.7802378000000001,0.8190894000000001 +52.850,0.7827832,0.7820203,0.8276794000000001,0.7851868000000001,0.7769098000000001,0.8229332 +52.900,0.7818522,0.7828916,0.82787,0.7840867,0.7744138,0.8246087000000001 +52.950,0.7814112000000001,0.7843924,0.8279654000000001,0.7830823,0.7738591,0.8245594 +53.000,0.7806272000000001,0.7858931,0.8272503,0.7818866,0.7744138,0.8219969 +53.050,0.7805292,0.7876843,0.8267259,0.7815040000000001,0.7768173,0.8180052 +53.100,0.7804802000000001,0.7891366000000001,0.8253434000000001,0.7808344,0.7802378000000001,0.8124859000000001 +53.150,0.7812642000000001,0.7910245600000001,0.8240562,0.7804996000000001,0.7852298000000001,0.8070651400000001 +53.200,0.7818032,0.79267051,0.8221016000000001,0.7796865000000001,0.79068397,0.8013979800000001 +53.250,0.7828322000000001,0.79475215,0.8202901,0.7795908,0.7966928600000001,0.7970121 +53.300,0.7838123,0.79659174,0.8180495000000001,0.7796865000000001,0.80219331,0.79316829 +53.350,0.7854293,0.799060667,0.8163332000000001,0.7804996000000001,0.80727775,0.79090143 +53.400,0.7866543,0.80128754,0.8142356000000001,0.7806909000000001,0.8107906,0.7891766 +53.450,0.7882224,0.80390169,0.8125194,0.7811692,0.8134253,0.7891766 +53.500,0.7894474,0.80593492,0.8100404,0.7818866,0.8145346,0.79011295 +53.550,0.79135844,0.80801656,0.80760911,0.7845650000000001,0.8153204000000001,0.7925769300000001 +53.600,0.7929754800000001,0.80966251,0.8045580600000001,0.7885827000000001,0.8152280000000001,0.7955337100000001 +53.650,0.7951315400000001,0.8116957,0.80222209,0.79427433,0.8155515000000001,0.79947608 +53.700,0.79689558,0.8132449,0.799647761,0.8004921100000001,0.8158289000000001,0.80346773 +53.750,0.79929664,0.8152781,0.79769318,0.80862305,0.8170306,0.8083464100000001 +53.800,0.80140369,0.8165368000000001,0.79526187,0.8184758000000001,0.8185098000000001,0.8127816 +53.850,0.80385375,0.8180859,0.79335496,0.8302418,0.8208671000000001,0.8174139 +53.900,0.80547079,0.8188605,0.79097132,0.8411468000000001,0.8229471,0.8209620000000001 +53.950,0.8075288500000001,0.8198771,0.7891121000000001,0.8491821,0.8251658000000001,0.8242637 +54.000,0.80934189,0.8200223000000001,0.7869191000000001,0.8512866,0.8263213,0.8259885 +54.050,0.811743,0.8202159,0.7854889,0.8474603,0.826876,0.8265799 +54.100,0.8135560000000001,0.8196834000000001,0.7837727000000001,0.8383249,0.8253969000000001,0.8249536000000001 +54.150,0.8158591000000001,0.8194414000000001,0.7826763,0.8274199000000001,0.8226698,0.8221447000000001 +54.200,0.8175741000000001,0.81857,0.7811984000000001,0.816467,0.817724,0.8170196000000001 +54.250,0.8194851000000001,0.8179407000000001,0.7801973000000001,0.8085273900000001,0.8119924000000001,0.8113032 +54.300,0.8202202000000001,0.8166336000000001,0.7788148,0.80427061,0.8049666400000001,0.804749 +54.350,0.8212002,0.8157622,0.7780043000000001,0.8005877680000001,0.7982644200000001,0.79888472 +54.400,0.8214452000000001,0.8139226,0.7770032,0.79700059,0.7914235300000001,0.7932668500000001 +54.450,0.8221802,0.8123735000000001,0.7768602,0.79695276,0.7858769000000001,0.7889795000000001 +54.500,0.8220332,0.8102434000000001,0.7765265,0.7973353900000001,0.7807462000000001,0.7853821000000001 +54.550,0.8222292,0.80874272,0.7768602,0.7980050000000001,0.7774182000000001,0.7835095000000001 +54.600,0.8216902,0.8067579,0.7767648,0.7987702600000001,0.7751071,0.7822282 +54.650,0.8213962,0.80516036,0.7774799,0.8013052,0.7747373000000001,0.7824253000000001 +54.700,0.8204162,0.80283667,0.7778613000000001,0.80340968,0.7752920000000001,0.7828196000000001 +54.750,0.8195341,0.800900257,0.7789101,0.8050358700000001,0.7776031,0.783953 +54.800,0.8174271000000001,0.7986733800000001,0.7795299,0.80546633,0.7806538000000001,0.7844458000000001 +54.850,0.8159571,0.79702744,0.7810554000000001,0.80632726,0.7849062,0.7851357 +54.900,0.813948,0.7948973800000001,0.7825332,0.8065664,0.78902,0.7848893 +54.950,0.812184,0.79325144,0.7850122,0.8063750900000001,0.79378086,0.7844951 +55.000,0.8095869,0.79121821,0.7872052,0.8046054100000001,0.79798709,0.7834602 +55.050,0.80718584,0.7899111000000001,0.7899225000000001,0.8027400800000001,0.80251687,0.7827703 +55.100,0.80429477,0.7879747,0.79221081,0.8000138195000001,0.80616842,0.7814397000000001 +55.150,0.80213871,0.7866676,0.79535721,0.79771802,0.8096813100000001,0.7807991000000001 +55.200,0.7997866530000001,0.7847312000000001,0.7983606000000001,0.79523091,0.8118538000000001,0.7802078 +55.250,0.7975816,0.7836178,0.8020314000000001,0.7934612400000001,0.8132404000000001,0.7808977 +55.300,0.7952785400000001,0.7824075,0.80522548,0.7916915600000001,0.8127782,0.7823268000000001 +55.350,0.7937105000000001,0.7819718,0.8089439500000001,0.7905436600000001,0.8112066,0.7855299 +55.400,0.7912604400000001,0.7813425,0.8121857,0.789061,0.8079248600000001,0.7893738 +55.450,0.7888594,0.7812457,0.8159042000000001,0.7881522000000001,0.8038110900000001,0.7945974 +55.500,0.7864093000000001,0.7808584000000001,0.8190029000000001,0.7871,0.7984493100000001,0.800116715 +55.550,0.7850373,0.7811489,0.8221493000000001,0.7869565000000001,0.79285642,0.80647378 +55.600,0.7835672,0.7811489,0.8242946,0.7861434,0.7869400000000001,0.8123381000000001 +55.650,0.7824892,0.7817782000000001,0.8263922,0.7861434,0.7819018000000001,0.8181531000000001 +55.700,0.7813622,0.7822623000000001,0.8273933,0.7858086000000001,0.7776031,0.8222925000000001 +55.750,0.7809212000000001,0.7834242,0.8282038,0.7852346,0.7751071,0.8249536000000001 +55.800,0.7800882,0.7844408,0.82787,0.7839433,0.7736280000000001,0.8250522 +55.850,0.7796962000000001,0.7859899,0.8275363,0.7829867,0.7740440000000001,0.8233274 +55.900,0.7793531,0.7871033000000001,0.8262968,0.7820779,0.7755231,0.8193358000000001 +55.950,0.7797942,0.7887009,0.8251527000000001,0.7817431,0.7788511,0.8145556 +56.000,0.7799902000000001,0.79025,0.8234841,0.78093,0.783196,0.80859281 +56.050,0.7806762,0.79238005,0.8219586000000001,0.7803561,0.7886964000000001,0.80292565 +56.100,0.7816562,0.79417123,0.8200040000000001,0.77983,0.79428931,0.7975541700000001 +56.150,0.7830282000000001,0.79649492,0.8183832000000001,0.7801648,0.800205754,0.79361181 +56.200,0.7844493,0.79862497,0.8164286000000001,0.7804517000000001,0.8051515300000001,0.79060575 +56.250,0.7861643,0.80119072,0.8148554,0.7808344,0.8094502,0.7890288000000001 +56.300,0.7871933000000001,0.80327236,0.8125671000000001,0.7807865,0.8122698,0.7881910000000001 +56.350,0.7887614000000001,0.80549923,0.8104218000000001,0.7817909000000001,0.8141649000000001,0.7889302 +56.400,0.79018241,0.807242,0.80760911,0.7840389000000001,0.8149044000000001,0.79080287 +56.450,0.79233847,0.80946887,0.80513013,0.7882957,0.8154129,0.7941538800000001 +56.500,0.79424951,0.8112116,0.8025558,0.7934612400000001,0.8154591,0.79794841 +56.550,0.79640557,0.8131965000000001,0.8006012170000001,0.799870332,0.8163373,0.8023835800000001 +56.600,0.79836562,0.8147456000000001,0.7983606000000001,0.8074273200000001,0.8171693,0.8069665800000001 +56.650,0.8008156780000001,0.8165368000000001,0.79631067,0.8174236,0.8189258,0.8119438 +56.700,0.8028247300000001,0.8176986,0.7940700500000001,0.8289026,0.8206822,0.8163297 +56.750,0.80468678,0.8186668,0.79202012,0.8404294,0.8229933,0.8204692 +56.800,0.80640182,0.8192477,0.7897795000000001,0.8489429000000001,0.8247498000000001,0.8236724000000001 +56.850,0.80870487,0.8198287000000001,0.7880633,0.8519562,0.8262751,0.8259885 +56.900,0.8107629000000001,0.8199255000000001,0.7862994000000001,0.8484168000000001,0.8265986000000001,0.8265799 +56.950,0.8128700000000001,0.8198771,0.7849169,0.8401424000000001,0.8259053000000001,0.8256928 +57.000,0.81483,0.8192962,0.7833437000000001,0.8290939,0.8231320000000001,0.8227361000000001 +57.050,0.8172311000000001,0.8186668,0.7820088000000001,0.8184758000000001,0.8188795000000001,0.8182023 +57.100,0.8186031,0.8176502,0.7805787000000001,0.8098187800000001,0.8129169,0.8122888 +57.150,0.8198282,0.816682,0.7794345,0.80532285,0.8062146400000001,0.80598099 +57.200,0.8207592,0.8152781,0.7782904,0.80135303,0.799050198,0.799525359 +57.250,0.8214452000000001,0.8138742,0.7776229,0.79762237,0.7923479800000001,0.79405532 +57.300,0.8219352000000001,0.8121314000000001,0.7770985,0.79695276,0.786108,0.7895216 +57.350,0.8224742,0.8105823000000001,0.7769078,0.7975267100000001,0.7809773000000001,0.7862691 +57.400,0.8224742,0.80883954,0.7767172,0.7976702,0.7771871,0.7840023 +57.450,0.8222292,0.80709677,0.7767648,0.7987702600000001,0.7749222,0.7828196000000001 +57.500,0.8216412000000001,0.80506354,0.7768602,0.8007790840000001,0.7739515,0.7824253000000001 +57.550,0.8207102000000001,0.80317554,0.7774322,0.80297922,0.7745986,0.7828689000000001 +57.600,0.8195341,0.8010939,0.7778613000000001,0.80407929,0.7763551000000001,0.7835588 +57.650,0.8183581000000001,0.799109077,0.7788148,0.80470107,0.7793595000000001,0.7843472 +57.700,0.8165451,0.79707585,0.7800542,0.80513153,0.783196,0.7847908 +57.750,0.8146340000000001,0.79538149,0.7821995,0.80565765,0.7876795000000001,0.7848893 +57.800,0.8125760000000001,0.79378395,0.7843925,0.80532285,0.79225553,0.7843472 +57.850,0.8102729000000001,0.79218641,0.7868715000000001,0.80427061,0.7970164200000001,0.7836080000000001 +57.900,0.8075288500000001,0.79044364,0.7893028000000001,0.8022617900000001,0.80149998,0.7825732000000001 +57.950,0.80512779,0.7888461,0.7922584800000001,0.8001094780000001,0.80579864,0.7816369000000001 +58.000,0.80272673,0.7870549,0.79521419,0.79771802,0.80935776,0.7807991000000001 +58.050,0.80047267,0.7855542000000001,0.79850361,0.7956613800000001,0.8120849,0.7804049000000001 +58.100,0.79792461,0.7840535000000001,0.80160235,0.7935569,0.8134253,0.7805527000000001 +58.150,0.7958665500000001,0.7830853000000001,0.8049394400000001,0.7920742000000001,0.8135178000000001,0.7821789 +58.200,0.7935635000000001,0.7823107,0.80803817,0.7907828100000001,0.8119462000000001,0.7848893 +58.250,0.79135844,0.7820687000000001,0.811709,0.7896827000000001,0.80917287,0.7889302 +58.300,0.7890064,0.7817298,0.8151891,0.788487,0.8049204200000001,0.79361181 +58.350,0.7870463000000001,0.781633,0.8188599000000001,0.7877696000000001,0.7999746426000001,0.7992789610000001 +58.400,0.7849883000000001,0.7814393000000001,0.8218633000000001,0.7869087,0.79410442,0.80509396 +58.450,0.7835182,0.7815846000000001,0.8246760000000001,0.7864782,0.788188,0.8112046 +58.500,0.7821952000000001,0.7818266,0.8266306,0.786239,0.7825026,0.8165268 +58.550,0.7814112000000001,0.7826012,0.8282514000000001,0.7861434,0.7779729000000001,0.8208142 +58.600,0.7807742000000001,0.7834726000000001,0.8287758000000001,0.7853303,0.7746449000000001,0.823426 +58.650,0.7804312,0.7848281,0.8289665,0.7844215,0.7731195000000001,0.8244608 +58.700,0.7798432000000001,0.7860383000000001,0.8281561000000001,0.7831302,0.772796,0.8229332 +58.750,0.7796962000000001,0.7874906,0.8272503,0.7823171,0.7742751,0.8200257000000001 +58.800,0.7795491,0.7888461,0.8258201000000001,0.7815040000000001,0.7771871,0.8152948000000001 +58.850,0.7799902000000001,0.79063728,0.8244853000000001,0.7809779,0.7817169,0.8101698 +58.900,0.7804802000000001,0.792138,0.8226260000000001,0.7800213,0.7870324000000001,0.80420692 +58.950,0.7814112000000001,0.79412282,0.8209098,0.7796865000000001,0.7930875300000001,0.7990818430000001 +59.000,0.7822932,0.79615605,0.8186692000000001,0.7796386000000001,0.7990039760000001,0.7947945200000001 +59.050,0.7839103000000001,0.79852815,0.81681,0.7802604000000001,0.80459687,0.79178846 +59.100,0.7853803,0.800706616,0.814474,0.7803082,0.80889553,0.7896694000000001 +59.150,0.7868503,0.80312713,0.8125194,0.7804039,0.8124084,0.788881 +59.200,0.7878304,0.80501513,0.8100404,0.7808822000000001,0.8144884,0.788881 +59.250,0.7894474,0.80714518,0.8078474800000001,0.7831302,0.8158751000000001,0.79050719 +59.300,0.7912604400000001,0.80908159,0.80532082,0.7867652,0.8161524,0.79316829 +59.350,0.7937105000000001,0.8113085000000001,0.80331856,0.79226551,0.8163835,0.79691354 +59.400,0.79576855,0.8130028,0.800982599,0.79853112,0.8166609,0.80070807 +59.450,0.7978756100000001,0.8147456000000001,0.7990280150000001,0.8063750900000001,0.8178164,0.80519251 +59.500,0.7997866530000001,0.8160527000000001,0.79664438,0.8159409000000001,0.8191569000000001,0.80972624 +59.550,0.80204071,0.817505,0.7945944500000001,0.8278025,0.8211906,0.8146049000000001 +59.600,0.80409876,0.8184248000000001,0.79221081,0.8394250000000001,0.8230395,0.8186951 +59.650,0.8063528200000001,0.8194898,0.7902562200000001,0.8487038,0.8250271,0.8223418 +59.700,0.80841087,0.8198771,0.788254,0.8521953,0.8262289,0.8247565 +59.750,0.8106159000000001,0.8203128000000001,0.7867761000000001,0.8493256,0.8267835,0.8261856000000001 +59.800,0.812331,0.8198771,0.7852029,0.8411946,0.8257204,0.8256436 +59.850,0.814585,0.8194414000000001,0.7840111000000001,0.8304331,0.8231320000000001,0.8236231 +59.900,0.8166431000000001,0.81857,0.7825332,0.8193368000000001,0.8184635,0.8193358000000001 +59.950,0.8182601,0.8178438,0.7815321000000001,0.8106319000000001,0.812732,0.8138657 +60.000,0.8195341,0.8167788,0.7804356,0.8054185,0.80584487,0.8073608200000001 +60.050,0.8212982000000001,0.8159074000000001,0.7797205,0.8014486900000001,0.79886531,0.80115159 +60.100,0.8219842,0.8143583000000001,0.7788148,0.79728756,0.7918395300000001,0.79533659 +60.150,0.8223762,0.8128576000000001,0.7782904,0.79633098,0.7858769000000001,0.79075359 +60.200,0.8225232,0.8108728000000001,0.7775276000000001,0.7965701300000001,0.7806075,0.7868605000000001 +60.250,0.8230132,0.80932364,0.7773369,0.7969049300000001,0.7770022000000001,0.784298 +60.300,0.8226702,0.80758087,0.7770985,0.7975267100000001,0.77446,0.7824746 +60.350,0.8222292,0.80603174,0.7774322,0.799631187,0.7739053,0.7819325 +60.400,0.8211022,0.80390169,0.7776706,0.80178349,0.7744138,0.7819325 +60.450,0.8199752,0.80196528,0.7784334,0.8035053400000001,0.7766324,0.7827703 +60.500,0.8185541000000001,0.79964159,0.7791485,0.8039836300000001,0.7796831000000001,0.7833616 +60.550,0.8169861,0.79775359,0.7806263,0.80470107,0.7839818000000001,0.7841008 +60.600,0.8147810000000001,0.79576877,0.7822472,0.80517936,0.7882804,0.7841994 +60.650,0.8128700000000001,0.79417123,0.7844401000000001,0.8054185,0.7930875300000001,0.7841994 +60.700,0.8102239,0.79238005,0.7864901000000001,0.8041271200000001,0.79747864,0.7834109 +60.750,0.8078228500000001,0.79092774,0.7892074,0.80250093,0.80205464,0.7826717000000001 +60.800,0.8053727900000001,0.789185,0.79173408,0.8001094780000001,0.80579864,0.7815383 +60.850,0.80331474,0.7877327000000001,0.7948328100000001,0.79795717,0.80935776,0.7808484 +60.900,0.800766677,0.7858447000000001,0.79769318,0.79556572,0.8116689,0.7801585000000001 +60.950,0.79851262,0.7845376000000001,0.80131631,0.79384387,0.8131942000000001,0.7804049000000001 +61.000,0.7962585600000001,0.7831821,0.8047010800000001,0.7921698500000001,0.8130555,0.7814397000000001 +61.050,0.7943475200000001,0.7825528,0.8085149,0.79116544,0.8117151,0.7840023 +61.100,0.7916524500000001,0.7818750000000001,0.8119473,0.7899219000000001,0.80875686,0.7875504000000001 +61.150,0.7893004,0.7817782000000001,0.8156181,0.7889175,0.80482798,0.79233054 +61.200,0.7871443,0.7813909,0.8187645,0.7879131,0.7996510870000001,0.79765273 +61.250,0.7857233,0.7815846000000001,0.8218633000000001,0.7873391000000001,0.7940582,0.8037634100000001 +61.300,0.7842043000000001,0.7815846000000001,0.8242946,0.7866217000000001,0.7880031000000001,0.8097755200000001 +61.350,0.7829792000000001,0.7821655000000001,0.8267259,0.7863825000000001,0.7825489,0.8155905 +61.400,0.7816072000000001,0.7825528,0.8282514000000001,0.7859999000000001,0.7778342,0.8201735 +61.450,0.7807742000000001,0.783521,0.8294432,0.7854260000000001,0.7747835000000001,0.8234753 +61.500,0.7799412,0.7843924,0.8293956,0.7840389000000001,0.7730271000000001,0.8245101 +61.550,0.7796962000000001,0.7858447000000001,0.8289665,0.7828432000000001,0.7730733000000001,0.8236231 +61.600,0.7794021000000001,0.7871517,0.8276794000000001,0.7816475,0.7742751,0.8204199000000001 +61.650,0.7795491,0.7889913000000001,0.8263445,0.7809779,0.7772795,0.816034 +61.700,0.7797452,0.79068569,0.8243899,0.7800691000000001,0.7814395000000001,0.8104162 +61.750,0.7805782,0.7926221,0.8227214,0.7793995,0.7868938000000001,0.8049461200000001 +61.800,0.7812642000000001,0.79436487,0.8205761,0.7788255000000001,0.79276398,0.799673198 +61.850,0.7822932,0.79649492,0.8187645,0.7791125000000001,0.7990039760000001,0.79558299 +61.900,0.7833712,0.79852815,0.8165239000000001,0.7793995,0.8044582,0.7923798200000001 +61.950,0.7849883000000001,0.800997077,0.814617,0.7800691000000001,0.80912664,0.7904579100000001 +62.000,0.7861643,0.80317554,0.812138,0.7802126,0.8122698,0.7893245000000001 +62.050,0.7876343,0.80549923,0.8100881000000001,0.7810257,0.8144422,0.7894723 +62.100,0.7891044,0.80729041,0.8075614400000001,0.7829867,0.8154129,0.7905564700000001 +62.150,0.79130944,0.80937205,0.80532082,0.7869565000000001,0.8159675000000001,0.7930204500000001 +62.200,0.7933184900000001,0.8111632000000001,0.8028895100000001,0.79193071,0.8158751000000001,0.7960757900000001 +62.250,0.7954255400000001,0.8131481,0.8008395810000001,0.7984354600000001,0.8164760000000001,0.7999688759000001 +62.300,0.7974836000000001,0.8146488000000001,0.79855129,0.8059924500000001,0.8172155000000001,0.8039605200000001 +62.350,0.7999336563,0.8161495000000001,0.7965967,0.8157018,0.8187871,0.80859281 +62.400,0.8017957,0.8173597,0.79440376,0.8271329000000001,0.8204511000000001,0.8130773 +62.450,0.8038047500000001,0.8186184000000001,0.7926398600000001,0.8391858000000001,0.8224386,0.8175617000000001 +62.500,0.8059608100000001,0.8192962,0.79058993,0.8486560000000001,0.8241489000000001,0.8212577000000001 +62.550,0.80865587,0.8199255000000001,0.7889691000000001,0.8529128,0.825628,0.8243623 +62.600,0.8108119,0.8199255000000001,0.7871575000000001,0.8503778000000001,0.8261826000000001,0.8259392000000001 +62.650,0.8129190000000001,0.8198287000000001,0.7858227000000001,0.8426295,0.825628,0.8259392000000001 +62.700,0.814879,0.8191509,0.7844878000000001,0.8316288000000001,0.8231320000000001,0.8239188000000001 +62.750,0.8169861,0.8184732,0.7836297000000001,0.8204847000000001,0.818972,0.8200257000000001 +62.800,0.8184071,0.8175534000000001,0.7824379,0.8108232000000001,0.8130555,0.8145556 +62.850,0.8199752,0.8167788,0.7814368,0.80532285,0.80649198,0.8084942500000001 +62.900,0.8209552,0.815617,0.7802449,0.80125738,0.799281309,0.80213718 +62.950,0.8223762,0.8143099,0.7793868,0.7970484200000001,0.79262531,0.7964207400000001 +63.000,0.8227682000000001,0.8126155,0.7782904,0.7955178900000001,0.7863391000000001,0.7911971000000001 +63.050,0.8229152000000001,0.811018,0.7776706,0.79618749,0.7813471000000001,0.787304 +63.100,0.8229642,0.80917841,0.7771462,0.79647447,0.7772795,0.7840516000000001 +63.150,0.8229642,0.80758087,0.7772416000000001,0.7970962500000001,0.7749222,0.7822282 +63.200,0.8223272,0.80569287,0.7771939,0.7989137500000001,0.7737666000000001,0.7812426 +63.250,0.8215922,0.80404692,0.7777183000000001,0.80178349,0.7743213,0.7812426 +63.300,0.8203672000000001,0.8020621,0.7781950000000001,0.80379232,0.7758929,0.7817354000000001 +63.350,0.8188481000000001,0.800270923,0.7792438,0.8047489000000001,0.7789898000000001,0.7826225 +63.400,0.8168391,0.79823769,0.7803403000000001,0.80498804,0.7828262,0.7832631000000001 +63.450,0.815075,0.79654333,0.7819135,0.8057054800000001,0.7874022,0.7837559000000001 +63.500,0.812821,0.79460692,0.7836297000000001,0.80580114,0.79188575,0.7836573 +63.550,0.8106649,0.7930577900000001,0.7860133,0.8050837000000001,0.7966466400000001,0.7834109 +63.600,0.80836187,0.79136344,0.78854,0.80307488,0.8009915320000001,0.782721 +63.650,0.80640182,0.7899111000000001,0.7913527,0.800970401,0.8051053100000001,0.7822282 +63.700,0.80351075,0.7882652,0.7939747,0.7987224300000001,0.8084795300000001,0.7813905000000001 +63.750,0.80101168,0.7869097,0.7971211,0.7966179600000001,0.8111604,0.7810948000000001 +63.800,0.79861062,0.7853606000000001,0.8003628530000001,0.7943699900000001,0.8126395000000001,0.7810948000000001 +63.850,0.79650357,0.7841503000000001,0.80422435,0.7927916300000001,0.8129169,0.7822775000000001 +63.900,0.7941025100000001,0.7830853000000001,0.8079905000000001,0.7914045900000001,0.8115302000000001,0.7843472 +63.950,0.79224046,0.7825044000000001,0.8118997,0.79020886,0.80889553,0.7877475 +64.000,0.79003541,0.7818750000000001,0.8153321,0.7887262,0.80482798,0.7918377400000001 +64.050,0.7881244000000001,0.7818750000000001,0.8186692000000001,0.7878652,0.7999284204,0.79715994 +64.100,0.7861153000000001,0.7818266,0.8214819000000001,0.7869565000000001,0.79419686,0.80277781 +64.150,0.7845473000000001,0.7821655000000001,0.8241992,0.7863825000000001,0.7883729,0.80878993 +64.200,0.7829792000000001,0.7821655000000001,0.8263922,0.7857129,0.78278,0.8141614 +64.250,0.7819992,0.7825528,0.8282514000000001,0.7854260000000001,0.778204,0.8191386 +64.300,0.7810682,0.7830853000000001,0.8291572,0.7846607000000001,0.7748298,0.8223911 +64.350,0.7806762,0.7842955,0.8294909,0.7836085,0.7731658,0.8243130000000001 +64.400,0.7798432000000001,0.7855058,0.8287758000000001,0.7823171,0.7728422,0.8237709000000001 +64.450,0.7795001,0.7871517,0.8277270000000001,0.7814083000000001,0.7741826,0.8212577000000001 +64.500,0.7794021000000001,0.7885556,0.8259631000000001,0.7807387,0.7768173,0.8168718 +64.550,0.7801862,0.79039523,0.8243423000000001,0.7802126,0.781116,0.8117467 +64.600,0.7808232,0.79208959,0.8223400000000001,0.7793995,0.7862929000000001,0.80603027 +64.650,0.7816072000000001,0.79407441,0.8204331,0.7790647,0.79230175,0.80070807 +64.700,0.7823422,0.79605923,0.8180971,0.7791125000000001,0.79817198,0.7960757900000001 +64.750,0.7836162000000001,0.7982861,0.8159995,0.7798778000000001,0.8038110900000001,0.7927247700000001 +64.800,0.7848413000000001,0.8005129740000001,0.8137589000000001,0.7805474,0.8083870900000001,0.79021151 +64.850,0.7862133,0.80293349,0.8118997,0.7810735000000001,0.8119462000000001,0.7890781 +64.900,0.7874873,0.80511195,0.80975439,0.7815996000000001,0.8140262,0.7887824 +64.950,0.7893984000000001,0.80729041,0.80765679,0.7836085,0.8152742000000001,0.7900143900000001 +65.000,0.79130944,0.80922682,0.80513013,0.7872913,0.815644,0.79208414 +65.050,0.7935635000000001,0.8114053,0.80284184,0.7927438,0.8160138,0.7955337100000001 +65.100,0.7955235500000001,0.8133417000000001,0.8006488900000001,0.7989137500000001,0.8163835,0.799229682 +65.150,0.7977776000000001,0.8150845,0.7988373200000001,0.8066620600000001,0.8174004,0.8035170100000001 +65.200,0.7995416470000001,0.8164399,0.7969780900000001,0.8160366,0.8187871,0.80785361 +65.250,0.8014526900000001,0.8178923,0.79526187,0.8275155000000001,0.8207746,0.8124859000000001 +65.300,0.80351075,0.8191025000000001,0.7932119400000001,0.8388988000000001,0.8226698,0.8167239000000001 +65.350,0.80591181,0.8200707,0.7914003700000001,0.8479385,0.8245186,0.8207156 +65.400,0.8079208600000001,0.8203128000000001,0.7895411,0.8515735000000001,0.8258129000000001,0.8236724000000001 +65.450,0.8101749,0.8203612,0.7880156,0.8490386000000001,0.8264600000000001,0.8257914000000001 +65.500,0.8122330000000001,0.8199255000000001,0.7863947,0.840812,0.825628,0.8260378 +65.550,0.8143400000000001,0.8195866,0.7851552,0.8301461,0.8233169,0.8245594 +65.600,0.8160061000000001,0.8188605,0.7838204000000001,0.8190976000000001,0.818972,0.8209127 +65.650,0.8180151,0.8181343000000001,0.7827239,0.8102492,0.8133791,0.8158862 +65.700,0.8193871,0.8170693,0.7814368,0.8050358700000001,0.8066306400000001,0.80972624 +65.750,0.8207102000000001,0.8161011,0.7802926,0.8015921800000001,0.7996510870000001,0.8035170100000001 +65.800,0.8217882000000001,0.8146004,0.7790531,0.7970962500000001,0.79257909,0.79740634 +65.850,0.8225232,0.8129544000000001,0.7781473000000001,0.7956613800000001,0.7864778,0.79208414 +65.900,0.8228172,0.8110664000000001,0.7772416000000001,0.7962353200000001,0.7810235000000001,0.7874519000000001 +65.950,0.8231602,0.80927523,0.7769555,0.7971919000000001,0.7770946000000001,0.7841008 +66.000,0.8228172,0.80748405,0.7766695,0.79776585,0.7744138,0.7815876 +66.050,0.8224742,0.80578969,0.7767648,0.799392042,0.7734893,0.7804542000000001 +66.100,0.8217882000000001,0.80385328,0.7770509000000001,0.8018791500000001,0.7738591,0.7800599 +66.150,0.8209062,0.80196528,0.7778136,0.80436626,0.7758929,0.7806513 +66.200,0.8192401,0.7998836410000001,0.7786241,0.8052750200000001,0.7789435,0.7813412000000001 +66.250,0.8176231,0.79814087,0.7798636000000001,0.80532285,0.7831035000000001,0.7824253000000001 +66.300,0.8156631000000001,0.79610764,0.7811984000000001,0.80513153,0.7875409,0.7830167 +66.350,0.8137030000000001,0.79436487,0.7832483,0.80517936,0.7923479800000001,0.7833616 +66.400,0.8115959,0.79267051,0.7855843,0.80436626,0.7969702000000001,0.7832631000000001 +66.450,0.8095869,0.79141185,0.788397,0.80259659,0.80149998,0.783066 +66.500,0.8071368400000001,0.79005636,0.7911620100000001,0.8001094780000001,0.80538264,0.7825239 +66.550,0.8045397700000001,0.7887009,0.7941653900000001,0.79790934,0.80889553,0.7821297 +66.600,0.8018937100000001,0.7870549,0.7970734300000001,0.7956613800000001,0.8113453,0.7815876 +66.650,0.7997866530000001,0.7855542000000001,0.800553544,0.79384387,0.8129169,0.7817354000000001 +66.700,0.7975816,0.7841987,0.80398598,0.7920263700000001,0.8129169,0.7825732000000001 +66.750,0.79537654,0.7832305,0.80779981,0.79087847,0.8116226000000001,0.7846922000000001 +66.800,0.79287748,0.7824075,0.811566,0.7897306000000001,0.80875686,0.7876983000000001 +66.850,0.7907214300000001,0.7818750000000001,0.8152844,0.7888218,0.8048742000000001,0.7918377400000001 +66.900,0.7886144,0.781633,0.8185739000000001,0.7877217000000001,0.7997897540000001,0.7965193 +66.950,0.7867523000000001,0.7817782000000001,0.8216249,0.7871,0.79428931,0.80198934 +67.000,0.7847433,0.7819234,0.8240562,0.7865260000000001,0.7882804,0.8075086600000001 +67.050,0.7835182,0.7824075,0.8262968,0.786239,0.7829186,0.8130773 +67.100,0.7821952000000001,0.7827948,0.8277747000000001,0.7858086000000001,0.7782964,0.8178081 +67.150,0.7809702000000001,0.783521,0.8287758000000001,0.7853303,0.7751071,0.8214548 +67.200,0.7797942,0.7841987,0.8288712,0.7842302000000001,0.7734431,0.8232781 +67.250,0.7793041,0.7853122,0.8284898,0.7830823,0.7733044,0.8231303000000001 +67.300,0.7793531,0.7865708,0.8273933,0.7818866,0.7745062,0.8207649 +67.350,0.7796962000000001,0.7881199000000001,0.8260108,0.781217,0.7772795,0.8169211000000001 +67.400,0.7801372000000001,0.7898143000000001,0.8242469,0.7804996000000001,0.7814395000000001,0.8117960000000001 +67.450,0.7806762,0.79175072,0.8224830000000001,0.7797821,0.7868475,0.80637523 +67.500,0.7813622,0.79359031,0.8204808,0.7791125000000001,0.79271775,0.800954468 +67.550,0.7824402,0.79572036,0.8182401,0.7791125000000001,0.7988190900000001,0.79647002 +67.600,0.7834692000000001,0.797802,0.8159042000000001,0.7797343,0.80436575,0.7930204500000001 +67.650,0.7845963,0.800174103,0.8138066,0.7805952,0.80898798,0.79075359 +67.700,0.7856743,0.80225574,0.8116613,0.7809779,0.8122235000000001,0.7893245000000001 +67.750,0.7873403,0.80443421,0.80980206,0.7816475,0.814396,0.7890288000000001 +67.800,0.7891534,0.80641903,0.8078474800000001,0.7834171000000001,0.8155053000000001,0.7899158000000001 +67.850,0.7912604400000001,0.80859749,0.8057975500000001,0.7871478000000001,0.8159213000000001,0.79203486 +67.900,0.79302448,0.8107275,0.80350926,0.7925046600000001,0.81606,0.7951394700000001 +67.950,0.79522954,0.8128576000000001,0.80141165,0.7991528960000001,0.8164298000000001,0.79888472 +68.000,0.7976796,0.814794,0.799361725,0.8069490300000001,0.817308,0.80287637 +68.050,0.8001786620000001,0.8163915,0.7975978300000001,0.8166583000000001,0.8187871,0.80731154 +68.100,0.80199171,0.8175534000000001,0.79578627,0.8279938,0.8205898,0.8118453 +68.150,0.80375575,0.8187152000000001,0.7940700500000001,0.8393771000000001,0.8225311000000001,0.8162804 +68.200,0.8057648000000001,0.8197318,0.79221081,0.8482733,0.8244262,0.8201735 +68.250,0.8080678600000001,0.820458,0.7903039000000001,0.8518127000000001,0.8259515000000001,0.8233767000000001 +68.300,0.8100769000000001,0.8205548,0.7883493,0.8488473000000001,0.8265986000000001,0.8255450000000001 +68.350,0.8122820000000001,0.8202644,0.7866808000000001,0.8406207,0.8259053000000001,0.8260871000000001 +68.400,0.8146340000000001,0.819635,0.7852029,0.8297157000000001,0.8234555,0.8247072000000001 +68.450,0.8163981,0.8189089,0.7839634000000001,0.8188106,0.8190644,0.8214055 +68.500,0.8178191,0.8179407000000001,0.7824379,0.8102014000000001,0.8131480000000001,0.8164283 +68.550,0.8192891,0.8171177000000001,0.7810077000000001,0.8053706700000001,0.80630709,0.8106133000000001 +68.600,0.8204652,0.8161979,0.7795299,0.80178349,0.79909642,0.8040098 +68.650,0.8214452000000001,0.8147456000000001,0.7783380000000001,0.7975267100000001,0.7921630900000001,0.79775129 +68.700,0.8223762,0.8130028,0.7773369,0.79599618,0.7858769000000001,0.79203486 +68.750,0.8228662,0.8111632000000001,0.7765741,0.79666579,0.7806538000000001,0.7875011000000001 +68.800,0.8226702,0.80932364,0.7761928,0.79743105,0.7766324,0.783953 +68.850,0.8224742,0.80748405,0.7761451,0.7980050000000001,0.7740902000000001,0.7815876 +68.900,0.8222782000000001,0.80554764,0.7763358,0.7991528960000001,0.7730733000000001,0.7804049000000001 +68.950,0.8212982000000001,0.80375646,0.7766218,0.80135303,0.7737666000000001,0.7801092000000001 +69.000,0.8200242000000001,0.80186846,0.7771939,0.8036966600000001,0.7758004000000001,0.7805527000000001 +69.050,0.8187501,0.79973841,0.7781950000000001,0.80513153,0.7790822000000001,0.7815383 +69.100,0.8172311000000001,0.79765677,0.7794822,0.8053706700000001,0.7833346,0.7826717000000001 +69.150,0.8155650000000001,0.79552672,0.7811507000000001,0.80513153,0.7880955000000001,0.7832631000000001 +69.200,0.813801,0.79363872,0.7831053,0.80498804,0.79299509,0.7836573 +69.250,0.8117920000000001,0.792138,0.7854413,0.80450975,0.79775598,0.7834602 +69.300,0.8095869,0.79078251,0.7882063,0.80297922,0.80219331,0.7830167 +69.350,0.80694083,0.7895723,0.79111434,0.8006355970000001,0.8062146400000001,0.7824253000000001 +69.400,0.8044417700000001,0.7884104000000001,0.7941653900000001,0.79824414,0.8094502,0.781834 +69.450,0.8020897100000001,0.7872486000000001,0.79716878,0.79599618,0.8119000000000001,0.7814397000000001 +69.500,0.800227664,0.7860867,0.8001721620000001,0.79413084,0.8131480000000001,0.7814890000000001 +69.550,0.79797361,0.7848281,0.80336624,0.79226551,0.8130555,0.7822775000000001 +69.600,0.79562155,0.7837146,0.80689402,0.7909263000000001,0.8113453,0.784298 +69.650,0.7936125,0.7827948,0.8107078000000001,0.7899697,0.8082946400000001,0.7874519000000001 +69.700,0.79155445,0.7821171,0.8143786000000001,0.7889653000000001,0.80408842,0.79178846 +69.750,0.7889574,0.7817298,0.8177634,0.7881044,0.79886531,0.7967657 +69.800,0.7867523000000001,0.7816814000000001,0.8207191,0.7874826,0.7930875300000001,0.8020879000000001 +69.850,0.7852333,0.7819234,0.8231981,0.7870521,0.7871249,0.80790289 +69.900,0.7840573000000001,0.7824075,0.8252957000000001,0.7866695,0.7817631,0.8134715 +69.950,0.7824402,0.7830369,0.8267259,0.7859999000000001,0.7771409,0.8182516000000001 +70.000,0.7810682,0.7838598,0.8277270000000001,0.7853781000000001,0.7741364000000001,0.8219476 + diff --git a/docs/Project.toml b/docs/Project.toml index 3821d09db..10845875a 100644 --- a/docs/Project.toml +++ b/docs/Project.toml @@ -3,6 +3,7 @@ Changelog = "5217a498-cd5d-4ec6-b8c2-9b85a09b6e3e" Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4" OrdinaryDiffEqTsit5 = "b1df2697-797e-41e3-8120-5422d3b24e4a" Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80" +SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf" SummationByPartsOperators = "9f78cca6-572e-554e-b819-917d2f1cf240" TrixiBase = "9a0f1c46-06d5-4909-a5a3-ce25d3fa3284" @@ -11,5 +12,6 @@ Changelog = "1" Documenter = "1" OrdinaryDiffEqTsit5 = "1" Plots = "1.9" +SparseArrays = "1" SummationByPartsOperators = "0.5.41" TrixiBase = "0.1.1" diff --git a/docs/src/dingemans.md b/docs/src/dingemans.md index d0ee0c810..68bd04865 100644 --- a/docs/src/dingemans.md +++ b/docs/src/dingemans.md @@ -23,7 +23,7 @@ Let's implement the Dingemans experiment and compare the performance of differen First, we load the necessary packages: ```@example dingemans -using DispersiveShallowWater, OrdinaryDiffEqTsit5, Plots +using DispersiveShallowWater, SummationByPartsOperators, SparseArrays, OrdinaryDiffEqTsit5, Plots ``` Next, we set up the different equation systems we want to compare: @@ -43,7 +43,7 @@ sgn = SerreGreenNaghdiEquations1D(bathymetry_type = bathymetry_variable, # Hyperbolic approximation of Serre-Green-Naghdi equations hysgn = HyperbolicSerreGreenNaghdiEquations1D(bathymetry_type = bathymetry_mild_slope, - lambda = 100.0, gravity = 9.81, eta0 = 0.8) + lambda = 200.0, gravity = 9.81, eta0 = 0.8) # for actual simulations a higher lambda (~500) is recommended # it is chosen so low to be able to see the difference between it # and the SGN equation. @@ -63,16 +63,27 @@ We create a computational domain that is large enough to contain the entire expe ```@example dingemans coordinates_min = -138.0 coordinates_max = 46.0 -N = 512 +N = 1024 mesh = Mesh1D(coordinates_min, coordinates_max, N) nothing # hide ``` -For the spatial discretization, we use fourth-order accurate [summation-by-parts operators](@ref sbp_operators): +For the spatial discretization, we use sixth-order accurate [summation-by-parts operators](@ref sbp_operators). Specifically, [upwind SBP operators](@ref upwind_sbp) are utilized for all equation systems +except the hyperbolic Serre-Green-Naghdi equations. The hyperbolic Serre-Green-Naghdi equations involve only first derivatives, making central discretization appropriate, whereas the other models include +higher-order derivatives that benefit from the upwind approach. It is important to note that when odd-order upwind operators are applied symmetrically to construct a central second-derivative operator, +the resulting operator achieves one order higher accuracy (since it is even). ```@example dingemans -accuracy_order = 4 -solver = Solver(mesh, accuracy_order) +accuracy_order = 6 +solver_central = Solver(mesh, accuracy_order) + +D1 = upwind_operators(periodic_derivative_operator; + derivative_order=1, + accuracy_order=accuracy_order - 1, + xmin=xmin(mesh), xmax=xmax(mesh), + N=nnodes(mesh)) +D2 = sparse(D1.plus) * sparse(D1.minus) +solver_upwind = Solver(D1, D2) nothing # hide ``` @@ -87,13 +98,13 @@ nothing # hide Now we create semidiscretizations for each equation system. Each semidiscretization bundles the mesh, equations, initial condition, solver, and boundary conditions: ```@example dingemans -semi_bbmbbm = Semidiscretization(mesh, bbmbbm, initial_condition, solver, +semi_bbmbbm = Semidiscretization(mesh, bbmbbm, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_sk = Semidiscretization(mesh, sk, initial_condition, solver, +semi_sk = Semidiscretization(mesh, sk, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_sgn = Semidiscretization(mesh, sgn, initial_condition, solver, +semi_sgn = Semidiscretization(mesh, sgn, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_hysgn = Semidiscretization(mesh, hysgn, initial_condition, solver, +semi_hysgn = Semidiscretization(mesh, hysgn, initial_condition, solver_central, boundary_conditions = boundary_conditions) nothing # hide ``` @@ -114,7 +125,7 @@ sol_hysgn = solve(ode_hysgn, Tsit5(); options...) nothing # hide ``` -## Visualization and Comparison +## Visualization of Temporal Evolution For proper comparison, we need to account for the fact that the BBM-BBM equations use a different reference level (``\eta_0 = 0``) compared to the other equations. We create a custom conversion function which allows us to easily shift the BBM-BBM results: @@ -147,16 +158,16 @@ for time_val in times for (i, (semi, sol, label, conversion, linestyle)) in enumerate(models) plot!(p, semi => sol, - step=step_idx, - label=label, - conversion=conversion, - plot_bathymetry=true, - legend=false, - title="Dingemans at t = $(time_val)", - suptitle="", - linewidth=[2 1], # 1 for the bathymetry - linestyles=[linestyle :solid], # :solid for the bathymetry - color=[i :black], # black for the bathymetry + step=step_idx, + label=label, + conversion=conversion, + plot_bathymetry=true, + legend=false, + title="Dingemans at t = $(time_val)", + suptitle="", + linewidth=[2 1], # 1 for the bathymetry + linestyles=[linestyle :solid], # :solid for the bathymetry + color=[i :black], # black for the bathymetry ) end @@ -169,32 +180,86 @@ legend_plot = plot(legend=:top, framestyle=:none, legendfontsize=11) for (i, (_, _, label, _, linestyle)) in enumerate(models) plot!(legend_plot, [], [], label=label, linestyles=linestyle, linewidth=2, color=i) end -plot!(legend_plot, [], [], label="Bathymetry", color=:black,) +legend_plot_bathymetry = plot(legend_plot, [], [], label="Bathymetry", color=:black,) xlims_zoom = [(-25, 15), (0, 40), (5, 45), (-100, -60)] snapshot_plots_zoom = [plot(snapshot_plots[i], xlims=xlims_zoom[i], ylims=(0.75, 0.85), title="Zoomed in at t = $(times[i])") for i in 1:4] # Combine all plots -all_plots = [snapshot_plots..., legend_plot, snapshot_plots_zoom...] +all_plots = [snapshot_plots..., legend_plot_bathymetry, snapshot_plots_zoom...] plot(all_plots..., - size=(900, 1100), - layout=@layout([a b; c d; e{0.14h}; f g; h i]), + size=(900, 1100), + layout=@layout([a b; c d; e{0.14h}; f g; h i]), ) -savefig("dingemans_comparison.png") # hide +savefig("dingemans_temporal.png") # hide nothing # hide ``` -![Dingemans comparison](dingemans_comparison.png) +![Dingemans temporal](dingemans_temporal.png) The results show how different dispersive wave models capture the wave evolution over the trapezoidal bathymetry. +## Comparison with Experimental Data + +During the experiment of Dingemans, the wave evolution was recorded at six gauges along the flume. These measurements provide a reference for validating the numerical models. +We compare the numerical results with the experimental data at the corresponding gauge locations. + +```@example dingemans +t_values, x_values, experimental_data = data_dingemans() + +tlims = [(20, 30), (25, 35), (30, 40), (35, 45), (40, 50), (45, 55)] + +snapshot_plots_time = [] +for (j, x_val) in enumerate(x_values) + p = plot(t_values, experimental_data[:, j], + xlims=tlims[j], + ylims=(0.765, 0.865), + label="experimental data", + linestyle=:dash, + color=:gray, + markershape=:circle, + markercolor=:gray, + markersize=1 + ) + for (i, (semi, sol, label, conversion, linestyle)) in enumerate(models) + plot!(p, semi => sol, x_val, + conversion=conversion, + label=label, + linestyle=linestyle, + color=i, + suptitle="", + legend=false, + title="Gauge at x = $(x_val)", + linewidth=2 + ) + + end + + push!(snapshot_plots_time, p) + +end + +legend_plot_data = plot(legend_plot, [], [], label="experimental data", linestyle=:dash, color=:gray, markershape=:circle, markercolor=:gray, markersize=1) + +all_plots2 = [snapshot_plots_time..., legend_plot_data] +plot(all_plots2..., layout=@layout([a b; c d; e f; g{0.16h}]), + size=(900, 900), suptitle="") + +savefig("dingemans_experimental.png") # hide +nothing # hide +``` + +![Dingemans experimental](dingemans_experimental.png) + +In this setup, the numerical results from the Svärd-Kalisch equations can capture the wave evolution more accurately compared to the other models. + ## [Plain program](@id overview-plain-program-dingemans) Here follows a version of the program without any comments. ```julia -using DispersiveShallowWater, OrdinaryDiffEqTsit5, Plots +using DispersiveShallowWater, SummationByPartsOperators, SparseArrays, OrdinaryDiffEqTsit5, Plots # BBM-BBM equations with variable bathymetry bbmbbm = BBMBBMEquations1D(bathymetry_type = bathymetry_variable, @@ -210,7 +275,7 @@ sgn = SerreGreenNaghdiEquations1D(bathymetry_type = bathymetry_variable, # Hyperbolic approximation of Serre-Green-Naghdi equations hysgn = HyperbolicSerreGreenNaghdiEquations1D(bathymetry_type = bathymetry_mild_slope, - lambda = 100.0, gravity = 9.81, eta0 = 0.8) + lambda = 200.0, gravity = 9.81, eta0 = 0.8) # for actual simulations a higher lambda (~500) is recommended # it is chosen so low to be able to see the difference between it # and the SGN equation. @@ -220,22 +285,30 @@ boundary_conditions = boundary_condition_periodic coordinates_min = -138.0 coordinates_max = 46.0 -N = 512 +N = 1024 mesh = Mesh1D(coordinates_min, coordinates_max, N) -accuracy_order = 4 -solver = Solver(mesh, accuracy_order) +accuracy_order = 6 +solver_central = Solver(mesh, accuracy_order) + +D1 = upwind_operators(periodic_derivative_operator; + derivative_order=1, + accuracy_order=accuracy_order - 1, + xmin=xmin(mesh), xmax=xmax(mesh), + N=nnodes(mesh)) +D2 = sparse(D1.plus) * sparse(D1.minus) +solver_upwind = Solver(D1, D2) tspan = (0.0, 70.0) saveat = range(tspan..., length = 500) -semi_bbmbbm = Semidiscretization(mesh, bbmbbm, initial_condition, solver, +semi_bbmbbm = Semidiscretization(mesh, bbmbbm, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_sk = Semidiscretization(mesh, sk, initial_condition, solver, +semi_sk = Semidiscretization(mesh, sk, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_sgn = Semidiscretization(mesh, sgn, initial_condition, solver, +semi_sgn = Semidiscretization(mesh, sgn, initial_condition, solver_upwind, boundary_conditions = boundary_conditions) -semi_hysgn = Semidiscretization(mesh, hysgn, initial_condition, solver, +semi_hysgn = Semidiscretization(mesh, hysgn, initial_condition, solver_central, boundary_conditions = boundary_conditions) ode_bbmbbm = semidiscretize(semi_bbmbbm, tspan) @@ -273,16 +346,16 @@ for time_val in times for (i, (semi, sol, label, conversion, linestyle)) in enumerate(models) plot!(p, semi => sol, - step=step_idx, - label=label, - conversion=conversion, - plot_bathymetry=true, - legend=false, - title="Dingemans at t = $(time_val)", - suptitle="", - linewidth=[2 1], # 1 for the bathymetry - linestyles=[linestyle :solid], # :solid for the bathymetry - color=[i :black], # black for the bathymetry + step=step_idx, + label=label, + conversion=conversion, + plot_bathymetry=true, + legend=false, + title="Dingemans at t = $(time_val)", + suptitle="", + linewidth=[2 1], # 1 for the bathymetry + linestyles=[linestyle :solid], # :solid for the bathymetry + color=[i :black], # black for the bathymetry ) end @@ -295,17 +368,57 @@ legend_plot = plot(legend=:top, framestyle=:none, legendfontsize=11) for (i, (_, _, label, _, linestyle)) in enumerate(models) plot!(legend_plot, [], [], label=label, linestyles=linestyle, linewidth=2, color=i) end -plot!(legend_plot, [], [], label="Bathymetry", color=:black,) +legend_plot_bathymetry = plot(legend_plot, [], [], label="Bathymetry", color=:black,) xlims_zoom = [(-25, 15), (0, 40), (5, 45), (-100, -60)] snapshot_plots_zoom = [plot(snapshot_plots[i], xlims=xlims_zoom[i], ylims=(0.75, 0.85), title="Zoomed in at t = $(times[i])") for i in 1:4] # Combine all plots -all_plots = [snapshot_plots..., legend_plot, snapshot_plots_zoom...] +all_plots = [snapshot_plots..., legend_plot_bathymetry, snapshot_plots_zoom...] plot(all_plots..., - size=(900, 1100), - layout=@layout([a b; c d; e{0.14h}; f g; h i]), + size=(900, 1100), + layout=@layout([a b; c d; e{0.14h}; f g; h i]), ) + +t_values, x_values, experimental_data = data_dingemans() + +tlims = [(20, 30), (25, 35), (30, 40), (35, 45), (40, 50), (45, 55)] + +snapshot_plots_time = [] +for (j, x_val) in enumerate(x_values) + p = plot(t_values, experimental_data[:, j], + xlims=tlims[j], + ylims=(0.765, 0.865), + label="experimental data", + linestyle=:dash, + color=:gray, + markershape=:circle, + markercolor=:gray, + markersize=1 + ) + for (i, (semi, sol, label, conversion, linestyle)) in enumerate(models) + plot!(p, semi => sol, x_val, + conversion=conversion, + label=label, + linestyle=linestyle, + color=i, + suptitle="", + legend=false, + title="Gauge at x = $(x_val)", + linewidth=2 + ) + + end + + push!(snapshot_plots_time, p) + +end + +legend_plot_data = plot(legend_plot, [], [], label="experimental data", linestyle=:dash, color=:gray, markershape=:circle, markercolor=:gray, markersize=1) + +all_plots2 = [snapshot_plots_time..., legend_plot_data] +plot(all_plots2..., layout=@layout([a b; c d; e f; g{0.16h}]), + size=(900, 900), suptitle="") ``` ### References diff --git a/docs/src/ref.md b/docs/src/ref.md index 1d6598c60..c57e03b88 100644 --- a/docs/src/ref.md +++ b/docs/src/ref.md @@ -66,6 +66,13 @@ Pages = Main.CALLBACKS_STEP_FILES ## Utilities +### Experimental data + +```@autodocs +Modules = [DispersiveShallowWater] +Pages = ["experimental_data.jl"] +``` + ### Utilities from DispersiveShallowWater.jl ```@autodocs diff --git a/src/DispersiveShallowWater.jl b/src/DispersiveShallowWater.jl index 630364f1b..f7821acd8 100644 --- a/src/DispersiveShallowWater.jl +++ b/src/DispersiveShallowWater.jl @@ -16,6 +16,7 @@ See also: [DispersiveShallowWater.jl](https://github.com/NumericalMathematics/Di module DispersiveShallowWater using BandedMatrices: BandedMatrix +using DelimitedFiles: readdlm using DiffEqBase: DiffEqBase, terminate! using FastBroadcast: @.. using ForwardDiff: ForwardDiff @@ -56,6 +57,7 @@ using TimerOutputs: TimerOutputs, print_timer, reset_timer! @reexport using TrixiBase: trixi_include using TrixiBase: TrixiBase, @trixi_timeit, timer +include("experimental_data.jl") include("boundary_conditions.jl") include("mesh.jl") include("equations/equations.jl") @@ -66,7 +68,9 @@ include("callbacks_step/callbacks_step.jl") include("visualization.jl") include("util.jl") -export examples_dir, get_examples, default_example, convergence_test +export examples_dir, get_examples, default_example, convergence_test, data_dir + +export data_dingemans export AbstractShallowWaterEquations, KdVEquation1D, diff --git a/src/experimental_data.jl b/src/experimental_data.jl new file mode 100644 index 000000000..bffc16396 --- /dev/null +++ b/src/experimental_data.jl @@ -0,0 +1,21 @@ +""" + data_dingemans() + +Load the experimental data for the Dingemans experiment. Returns the time values, x-coordinates of the six wave gauges, +and experimental data, which is a matrix of wave heights at each gauge location (columns) over time (rows). + +- Dingemans (1994) + Comparison of computations with Boussinesq-like models and laboratory measurements + [URL: https://resolver.tudelft.nl/uuid:c2091d53-f455-48af-a84b-ac86680455e9](https://resolver.tudelft.nl/uuid:c2091d53-f455-48af-a84b-ac86680455e9) +- Dingemans (1997): + Water Wave Propagation Over Uneven Bottoms (In 2 Parts). + [DOI: 10.1142/1241](https://doi.org/10.1142/1241) +""" +function data_dingemans() + path_dingemans = joinpath(data_dir(), "Dingemans.csv") + all_data, _ = readdlm(path_dingemans, ','; header = true) + t_values = all_data[:, 1] + x_values = (3.04, 9.44, 20.04, 26.04, 30.44, 37.04) + experimental_data = all_data[:, 2:end] + return t_values, x_values, experimental_data +end diff --git a/src/util.jl b/src/util.jl index da8f57523..bc97bb986 100644 --- a/src/util.jl +++ b/src/util.jl @@ -1,7 +1,7 @@ """ examples_dir() -Return the directory where the example files provided with DispersiveShallowWater.jl are located. If DispersiveShallowWater is +Return the directory where the example files provided with DispersiveShallowWater.jl are located. If DispersiveShallowWater.jl is installed as a regular package (with `]add DispersiveShallowWater`), these files are read-only and should *not* be modified. To find out which files are available, use, e.g., `readdir`. @@ -47,6 +47,19 @@ function default_example() joinpath(examples_dir(), "bbm_bbm_1d", "bbm_bbm_1d_basic.jl") end +""" + data_dir() + +Return the directory where the data files provided with DispersiveShallowWater.jl are located. If DispersiveShallowWater.jl is +installed as a regular package (with `]add DispersiveShallowWater`), these files are read-only and should *not* be +modified. To find out which files are available, use, e.g., `readdir`. + +```@example +readdir(data_dir()) +``` +""" +data_dir() = pkgdir(DispersiveShallowWater, "data")::String + function convergence_test(example::AbstractString, iterations_or_Ns; kwargs...) convergence_test(Main, example::AbstractString, iterations_or_Ns; kwargs...) end diff --git a/test/test_unit.jl b/test/test_unit.jl index 528d1e813..c4fb475f7 100644 --- a/test/test_unit.jl +++ b/test/test_unit.jl @@ -532,7 +532,11 @@ end end @testitem "util" setup=[Setup] begin + @test_nowarn examples_dir() @test_nowarn get_examples() + @test_nowarn data_dir() + t, x, experimental_data = data_dingemans() + @test size(experimental_data) == (length(t), length(x)) accuracy_orders = [2, 4, 6] for accuracy_order in accuracy_orders