module SrtmKgb27
  # KURUCZ
  #     The following values were obtained using the "low resolution"
  #     version of the Kurucz solar source function.  For unknown reasons,
  #     the total irradiance in this band differs from the corresponding
  #     total in the "high-resolution" version of the Kurucz function.
  #     Therefore, below these values are scaled by the factor SCALEKUR.
  sfluxref = [14.0526    , 11.4794    , 8.72590    , 5.56966    ,3.80927    , 1.57690    , 1.15099    , 1.10012    ,0.658212   , 5.86859e-02, 5.56186e-02, 4.68040e-02,3.64897e-02, 3.58053e-02, 1.38130e-02, 1.90193e-03]

  #     Rayleigh extinction coefficient at v = 2925 cm-1.
  rayl = [3.44534e-06,4.14480e-06,4.95069e-06,5.81204e-06,6.69748e-06,7.56488e-06,8.36344e-06,9.04135e-06,9.58324e-06,9.81542e-06,9.75119e-06,9.74533e-06,9.74139e-06,9.73525e-06,9.73577e-06,9.73618e-06]

  scalekur = 50.15/48.37

  layreffr = 32

  #     ------------------------------------------------------------------

  #     The array KA contains absorption coefs at the 16 chosen g-values
  #     for a range of pressure levels> ~100mb:and binary:temperatures
  #     species parameters (see taumol.f for definition).  The first
  #     index in the array, JS, runs from 1 to 9, and corresponds to
  #     different values of the binary species parameter.  For instance,
  #     JS=1 refers to dry air, JS = 2 corresponds to the paramter value 1/8,
  #     JS = 3 corresponds to the parameter value 2/8, etc.  The second index
  #     in the array, JT, which runs from 1 to 5, corresponds to different
  #     temperatures.  More specifically, JT = 3 means that the data are for
  #     the reference temperature TREF for this  pressure level:JT = 2 refers
  #     to TREF-15, JT = 1 is for TREF-30:and JT = 5:JT = 4 is for TREF+15
  #     is for TREF+30.  The third index:runs from 1 to 13 and refers:JP
  #     to the JPth reference pressure level (see taumol.f for these levels
  #     in mb).  The fourth index, IG, goes from 1 to 16, and indicates
  #     which g-interval the absorption coefficients are for.
  #     -----------------------------------------------------------------
  ka = zeros(5,13,16)
  ka[:, 1, 1] = [0.22907e+00,0.25625e+00,0.28779e+00,0.32376e+00,0.36426e+00]
  ka[:, 2, 1] = [0.21913e+00,0.24445e+00,0.27422e+00,0.30832e+00,0.34694e+00]
  ka[:, 3, 1] = [0.20602e+00,0.22864e+00,0.25574e+00,0.28721e+00,0.32310e+00]
  ka[:, 4, 1] = [0.19379e+00,0.21352e+00,0.23771e+00,0.26643e+00,0.29940e+00]
  ka[:, 5, 1] = [0.18369e+00,0.20061e+00,0.22200e+00,0.24786e+00,0.27816e+00]
  ka[:, 6, 1] = [0.17517e+00,0.18920e+00,0.20771e+00,0.23069e+00,0.25817e+00]
  ka[:, 7, 1] = [0.16888e+00,0.17988e+00,0.19559e+00,0.21577e+00,0.24042e+00]
  ka[:, 8, 1] = [0.16442e+00,0.17273e+00,0.18565e+00,0.20315e+00,0.22512e+00]
  ka[:, 9, 1] = [0.16159e+00,0.16738e+00,0.17752e+00,0.19242e+00,0.21179e+00]
  ka[:,10, 1] = [0.16037e+00,0.16402e+00,0.17204e+00,0.18463e+00,0.20183e+00]
  ka[:,11, 1] = [0.16036e+00,0.16399e+00,0.17198e+00,0.18455e+00,0.20173e+00]
  ka[:,12, 1] = [0.16036e+00,0.16399e+00,0.17198e+00,0.18455e+00,0.20173e+00]
  ka[:,13, 1] = [0.16036e+00,0.16399e+00,0.17198e+00,0.18455e+00,0.20173e+00]
  ka[:, 1, 2] = [0.20716e+01,0.21925e+01,0.23327e+01,0.24921e+01,0.26726e+01]
  ka[:, 2, 2] = [0.20272e+01,0.21401e+01,0.22723e+01,0.24238e+01,0.25950e+01]
  ka[:, 3, 2] = [0.19684e+01,0.20697e+01,0.21902e+01,0.23301e+01,0.24892e+01]
  ka[:, 4, 2] = [0.19133e+01,0.20021e+01,0.21101e+01,0.22375e+01,0.23842e+01]
  ka[:, 5, 2] = [0.18673e+01,0.19441e+01,0.20400e+01,0.21553e+01,0.22898e+01]
  ka[:, 6, 2] = [0.18277e+01,0.18925e+01,0.19760e+01,0.20788e+01,0.22009e+01]
  ka[:, 7, 2] = [0.17963e+01,0.18497e+01,0.19214e+01,0.20121e+01,0.21222e+01]
  ka[:, 8, 2] = [0.17749e+01,0.18156e+01,0.18763e+01,0.19555e+01,0.20540e+01]
  ka[:, 9, 2] = [0.17600e+01,0.17885e+01,0.18388e+01,0.19071e+01,0.19943e+01]
  ka[:,10, 2] = [0.17526e+01,0.17729e+01,0.18121e+01,0.18716e+01,0.19496e+01]
  ka[:,11, 2] = [0.17525e+01,0.17727e+01,0.18119e+01,0.18712e+01,0.19491e+01]
  ka[:,12, 2] = [0.17525e+01,0.17727e+01,0.18119e+01,0.18712e+01,0.19491e+01]
  ka[:,13, 2] = [0.17525e+01,0.17727e+01,0.18119e+01,0.18712e+01,0.19491e+01]
  ka[:, 1, 3] = [0.12189e+02,0.12619e+02,0.13110e+02,0.13663e+02,0.14275e+02]
  ka[:, 2, 3] = [0.12030e+02,0.12434e+02,0.12899e+02,0.13427e+02,0.14015e+02]
  ka[:, 3, 3] = [0.11816e+02,0.12183e+02,0.12611e+02,0.13101e+02,0.13653e+02]
  ka[:, 4, 3] = [0.11611e+02,0.11938e+02,0.12327e+02,0.12777e+02,0.13289e+02]
  ka[:, 5, 3] = [0.11437e+02,0.11726e+02,0.12076e+02,0.12488e+02,0.12961e+02]
  ka[:, 6, 3] = [0.11284e+02,0.11533e+02,0.11843e+02,0.12215e+02,0.12649e+02]
  ka[:, 7, 3] = [0.11160e+02,0.11370e+02,0.11641e+02,0.11975e+02,0.12370e+02]
  ka[:, 8, 3] = [0.11063e+02,0.11237e+02,0.11471e+02,0.11768e+02,0.12126e+02]
  ka[:, 9, 3] = [0.10990e+02,0.11129e+02,0.11327e+02,0.11588e+02,0.11910e+02]
  ka[:,10, 3] = [0.10945e+02,0.11053e+02,0.11223e+02,0.11454e+02,0.11746e+02]
  ka[:,11, 3] = [0.10945e+02,0.11052e+02,0.11222e+02,0.11452e+02,0.11744e+02]
  ka[:,12, 3] = [0.10945e+02,0.11052e+02,0.11222e+02,0.11452e+02,0.11744e+02]
  ka[:,13, 3] = [0.10945e+02,0.11052e+02,0.11222e+02,0.11452e+02,0.11744e+02]
  ka[:, 1, 4] = [0.55057e+02,0.56143e+02,0.57364e+02,0.58719e+02,0.60208e+02]
  ka[:, 2, 4] = [0.54647e+02,0.55677e+02,0.56842e+02,0.58142e+02,0.59575e+02]
  ka[:, 3, 4] = [0.54090e+02,0.55039e+02,0.56123e+02,0.57342e+02,0.58694e+02]
  ka[:, 4, 4] = [0.53548e+02,0.54410e+02,0.55408e+02,0.56539e+02,0.57805e+02]
  ka[:, 5, 4] = [0.53074e+02,0.53853e+02,0.54766e+02,0.55813e+02,0.56995e+02]
  ka[:, 6, 4] = [0.52644e+02,0.53336e+02,0.54162e+02,0.55123e+02,0.56218e+02]
  ka[:, 7, 4] = [0.52279e+02,0.52887e+02,0.53629e+02,0.54505e+02,0.55516e+02]
  ka[:, 8, 4] = [0.51979e+02,0.52507e+02,0.53168e+02,0.53965e+02,0.54895e+02]
  ka[:, 9, 4] = [0.51734e+02,0.52183e+02,0.52767e+02,0.53485e+02,0.54337e+02]
  ka[:,10, 4] = [0.51565e+02,0.51948e+02,0.52467e+02,0.53120e+02,0.53907e+02]
  ka[:,11, 4] = [0.51563e+02,0.51946e+02,0.52464e+02,0.53116e+02,0.53902e+02]
  ka[:,12, 4] = [0.51563e+02,0.51946e+02,0.52464e+02,0.53116e+02,0.53902e+02]
  ka[:,13, 4] = [0.51563e+02,0.51946e+02,0.52464e+02,0.53116e+02,0.53902e+02]
  ka[:, 1, 5] = [0.17794e+03,0.17973e+03,0.18164e+03,0.18366e+03,0.18581e+03]
  ka[:, 2, 5] = [0.17724e+03,0.17898e+03,0.18083e+03,0.18281e+03,0.18491e+03]
  ka[:, 3, 5] = [0.17624e+03,0.17791e+03,0.17970e+03,0.18160e+03,0.18363e+03]
  ka[:, 4, 5] = [0.17523e+03,0.17682e+03,0.17853e+03,0.18036e+03,0.18231e+03]
  ka[:, 5, 5] = [0.17429e+03,0.17581e+03,0.17744e+03,0.17920e+03,0.18107e+03]
  ka[:, 6, 5] = [0.17338e+03,0.17482e+03,0.17638e+03,0.17805e+03,0.17985e+03]
  ka[:, 7, 5] = [0.17253e+03,0.17390e+03,0.17539e+03,0.17699e+03,0.17871e+03]
  ka[:, 8, 5] = [0.17177e+03,0.17307e+03,0.17448e+03,0.17601e+03,0.17767e+03]
  ka[:, 9, 5] = [0.17107e+03,0.17230e+03,0.17364e+03,0.17511e+03,0.17669e+03]
  ka[:,10, 5] = [0.17052e+03,0.17169e+03,0.17298e+03,0.17438e+03,0.17591e+03]
  ka[:,11, 5] = [0.17051e+03,0.17168e+03,0.17297e+03,0.17437e+03,0.17590e+03]
  ka[:,12, 5] = [0.17051e+03,0.17168e+03,0.17297e+03,0.17437e+03,0.17590e+03]
  ka[:,13, 5] = [0.17051e+03,0.17168e+03,0.17297e+03,0.17437e+03,0.17590e+03]
  ka[:, 1, 6] = [0.37680e+05,0.36801e+05,0.36011e+05,0.35298e+05,0.34715e+05]
  ka[:, 2, 6] = [0.39286e+03,0.39463e+03,0.39616e+03,0.39748e+03,0.39858e+03]
  ka[:, 3, 6] = [0.39169e+03,0.39359e+03,0.39526e+03,0.39671e+03,0.39794e+03]
  ka[:, 4, 6] = [0.39034e+03,0.39239e+03,0.39421e+03,0.39580e+03,0.39717e+03]
  ka[:, 5, 6] = [0.38895e+03,0.39113e+03,0.39309e+03,0.39483e+03,0.39634e+03]
  ka[:, 6, 6] = [0.38742e+03,0.38975e+03,0.39185e+03,0.39373e+03,0.39539e+03]
  ka[:, 7, 6] = [0.38585e+03,0.38832e+03,0.39056e+03,0.39258e+03,0.39438e+03]
  ka[:, 8, 6] = [0.38426e+03,0.38686e+03,0.38924e+03,0.39140e+03,0.39333e+03]
  ka[:, 9, 6] = [0.38264e+03,0.38538e+03,0.38789e+03,0.39017e+03,0.39223e+03]
  ka[:,10, 6] = [0.38123e+03,0.38408e+03,0.38670e+03,0.38909e+03,0.39126e+03]
  ka[:,11, 6] = [0.38122e+03,0.38406e+03,0.38668e+03,0.38908e+03,0.39125e+03]
  ka[:,12, 6] = [0.38122e+03,0.38406e+03,0.38668e+03,0.38908e+03,0.39125e+03]
  ka[:,13, 6] = [0.38122e+03,0.38406e+03,0.38668e+03,0.38908e+03,0.39125e+03]
  ka[:, 1, 7] = [0.18202e+07,0.17771e+07,0.17382e+07,0.17034e+07,0.16716e+07]
  ka[:, 2, 7] = [0.63027e+03,0.63194e+03,0.63388e+03,0.63607e+03,0.63849e+03]
  ka[:, 3, 7] = [0.62938e+03,0.63090e+03,0.63268e+03,0.63472e+03,0.63701e+03]
  ka[:, 4, 7] = [0.62853e+03,0.62989e+03,0.63150e+03,0.63337e+03,0.63550e+03]
  ka[:, 5, 7] = [0.62781e+03,0.62901e+03,0.63046e+03,0.63217e+03,0.63413e+03]
  ka[:, 6, 7] = [0.62717e+03,0.62821e+03,0.62950e+03,0.63104e+03,0.63283e+03]
  ka[:, 7, 7] = [0.62666e+03,0.62753e+03,0.62867e+03,0.63005e+03,0.63168e+03]
  ka[:, 8, 7] = [0.62627e+03,0.62698e+03,0.62795e+03,0.62919e+03,0.63067e+03]
  ka[:, 9, 7] = [0.62600e+03,0.62653e+03,0.62735e+03,0.62844e+03,0.62978e+03]
  ka[:,10, 7] = [0.62583e+03,0.62624e+03,0.62692e+03,0.62788e+03,0.62910e+03]
  ka[:,11, 7] = [0.62583e+03,0.62623e+03,0.62691e+03,0.62787e+03,0.62909e+03]
  ka[:,12, 7] = [0.62583e+03,0.62623e+03,0.62691e+03,0.62787e+03,0.62909e+03]
  ka[:,13, 7] = [0.62583e+03,0.62623e+03,0.62691e+03,0.62787e+03,0.62909e+03]
  ka[:, 1, 8] = [0.55488e+07,0.54180e+07,0.53005e+07,0.51943e+07,0.50978e+07]
  ka[:, 2, 8] = [0.82143e+03,0.82232e+03,0.82327e+03,0.82429e+03,0.82541e+03]
  ka[:, 3, 8] = [0.82092e+03,0.82177e+03,0.82269e+03,0.82366e+03,0.82472e+03]
  ka[:, 4, 8] = [0.82041e+03,0.82121e+03,0.82209e+03,0.82303e+03,0.82402e+03]
  ka[:, 5, 8] = [0.81996e+03,0.82070e+03,0.82152e+03,0.82243e+03,0.82339e+03]
  ka[:, 6, 8] = [0.81952e+03,0.82021e+03,0.82098e+03,0.82184e+03,0.82277e+03]
  ka[:, 7, 8] = [0.81917e+03,0.81977e+03,0.82048e+03,0.82129e+03,0.82218e+03]
  ka[:, 8, 8] = [0.81881e+03,0.81938e+03,0.82004e+03,0.82080e+03,0.82163e+03]
  ka[:, 9, 8] = [0.81846e+03,0.81906e+03,0.81965e+03,0.82035e+03,0.82114e+03]
  ka[:,10, 8] = [0.81821e+03,0.81877e+03,0.81934e+03,0.82000e+03,0.82074e+03]
  ka[:,11, 8] = [0.81821e+03,0.81877e+03,0.81934e+03,0.81999e+03,0.82074e+03]
  ka[:,12, 8] = [0.81821e+03,0.81877e+03,0.81934e+03,0.81999e+03,0.82074e+03]
  ka[:,13, 8] = [0.81821e+03,0.81877e+03,0.81935e+03,0.81999e+03,0.82074e+03]
  ka[:, 1, 9] = [0.10204e+08,0.99477e+07,0.97209e+07,0.95103e+07,0.92960e+07]
  ka[:, 2, 9] = [0.95102e+03,0.95144e+03,0.95133e+03,0.95066e+03,0.94964e+03]
  ka[:, 3, 9] = [0.95048e+03,0.95118e+03,0.95136e+03,0.95103e+03,0.95020e+03]
  ka[:, 4, 9] = [0.94969e+03,0.95077e+03,0.95126e+03,0.95124e+03,0.95074e+03]
  ka[:, 5, 9] = [0.94872e+03,0.95015e+03,0.95101e+03,0.95129e+03,0.95105e+03]
  ka[:, 6, 9] = [0.94752e+03,0.94929e+03,0.95052e+03,0.95116e+03,0.95126e+03]
  ka[:, 7, 9] = [0.94612e+03,0.94825e+03,0.94981e+03,0.95083e+03,0.95126e+03]
  ka[:, 8, 9] = [0.94461e+03,0.94704e+03,0.94894e+03,0.95029e+03,0.95107e+03]
  ka[:, 9, 9] = [0.94291e+03,0.94567e+03,0.94791e+03,0.94957e+03,0.95069e+03]
  ka[:,10, 9] = [0.94135e+03,0.94441e+03,0.94689e+03,0.94884e+03,0.95022e+03]
  ka[:,11, 9] = [0.94133e+03,0.94440e+03,0.94688e+03,0.94883e+03,0.95021e+03]
  ka[:,12, 9] = [0.94133e+03,0.94438e+03,0.94688e+03,0.94883e+03,0.95021e+03]
  ka[:,13, 9] = [0.94133e+03,0.94438e+03,0.94687e+03,0.94883e+03,0.95021e+03]
  ka[:, 1,10] = [0.12764e+08,0.12408e+08,0.12012e+08,0.11438e+08,0.11217e+08]
  ka[:, 2,10] = [0.99589e+03,0.99493e+03,0.99263e+03,0.98974e+03,0.98742e+03]
  ka[:, 3,10] = [0.99532e+03,0.99571e+03,0.99447e+03,0.99207e+03,0.98779e+03]
  ka[:, 4,10] = [0.99446e+03,0.99559e+03,0.99536e+03,0.99386e+03,0.99088e+03]
  ka[:, 5,10] = [0.99288e+03,0.99493e+03,0.99558e+03,0.99487e+03,0.99336e+03]
  ka[:, 6,10] = [0.99069e+03,0.99366e+03,0.99519e+03,0.99534e+03,0.99428e+03]
  ka[:, 7,10] = [0.98811e+03,0.99191e+03,0.99430e+03,0.99510e+03,0.99504e+03]
  ka[:, 8,10] = [0.98531e+03,0.98996e+03,0.99300e+03,0.99479e+03,0.99513e+03]
  ka[:, 9,10] = [0.98264e+03,0.98737e+03,0.99138e+03,0.99390e+03,0.99503e+03]
  ka[:,10,10] = [0.98010e+03,0.98517e+03,0.98974e+03,0.99281e+03,0.99470e+03]
  ka[:,11,10] = [0.98010e+03,0.98509e+03,0.98985e+03,0.99280e+03,0.99468e+03]
  ka[:,12,10] = [0.98035e+03,0.98524e+03,0.98973e+03,0.99279e+03,0.99468e+03]
  ka[:,13,10] = [0.98012e+03,0.98526e+03,0.98973e+03,0.99279e+03,0.99468e+03]
  ka[:, 1,11] = [0.12266e+08,0.12000e+08,0.11724e+08,0.11828e+08,0.11523e+08]
  ka[:, 2,11] = [0.10001e+04,0.99984e+03,0.99745e+03,0.99535e+03,0.98599e+03]
  ka[:, 3,11] = [0.99963e+03,0.99927e+03,0.99848e+03,0.99596e+03,0.99390e+03]
  ka[:, 4,11] = [0.99904e+03,0.99942e+03,0.99870e+03,0.99747e+03,0.99542e+03]
  ka[:, 5,11] = [0.99826e+03,0.99927e+03,0.99919e+03,0.99821e+03,0.99627e+03]
  ka[:, 6,11] = [0.99730e+03,0.99882e+03,0.99938e+03,0.99902e+03,0.99764e+03]
  ka[:, 7,11] = [0.99551e+03,0.99810e+03,0.99919e+03,0.99961e+03,0.99863e+03]
  ka[:, 8,11] = [0.99324e+03,0.99688e+03,0.99871e+03,0.99940e+03,0.99942e+03]
  ka[:, 9,11] = [0.99058e+03,0.99489e+03,0.99773e+03,0.99912e+03,0.99954e+03]
  ka[:,10,11] = [0.98802e+03,0.99297e+03,0.99656e+03,0.99861e+03,0.99937e+03]
  ka[:,11,11] = [0.98800e+03,0.99310e+03,0.99655e+03,0.99860e+03,0.99936e+03]
  ka[:,12,11] = [0.98771e+03,0.99296e+03,0.99669e+03,0.99859e+03,0.99936e+03]
  ka[:,13,11] = [0.98800e+03,0.99296e+03,0.99655e+03,0.99859e+03,0.99936e+03]
  ka[:, 1,12] = [0.12503e+08,0.11946e+08,0.11557e+08,0.11503e+08,0.11901e+08]
  ka[:, 2,12] = [0.10015e+04,0.10004e+04,0.99932e+03,0.99373e+03,0.98591e+03]
  ka[:, 3,12] = [0.10024e+04,0.10021e+04,0.10003e+04,0.99817e+03,0.99146e+03]
  ka[:, 4,12] = [0.10022e+04,0.10027e+04,0.10021e+04,0.99956e+03,0.99610e+03]
  ka[:, 5,12] = [0.10024e+04,0.10024e+04,0.10028e+04,0.10015e+04,0.99826e+03]
  ka[:, 6,12] = [0.10017e+04,0.10025e+04,0.10025e+04,0.10023e+04,0.10008e+04]
  ka[:, 7,12] = [0.10005e+04,0.10020e+04,0.10028e+04,0.10028e+04,0.10019e+04]
  ka[:, 8,12] = [0.10000e+04,0.10010e+04,0.10024e+04,0.10029e+04,0.10023e+04]
  ka[:, 9,12] = [0.99906e+03,0.10004e+04,0.10018e+04,0.10029e+04,0.10028e+04]
  ka[:,10,12] = [0.99791e+03,0.99989e+03,0.10010e+04,0.10024e+04,0.10031e+04]
  ka[:,11,12] = [0.99769e+03,0.99967e+03,0.10008e+04,0.10024e+04,0.10029e+04]
  ka[:,12,12] = [0.99787e+03,0.99985e+03,0.10008e+04,0.10025e+04,0.10029e+04]
  ka[:,13,12] = [0.99768e+03,0.99984e+03,0.10012e+04,0.10026e+04,0.10029e+04]
  ka[:, 1,13] = [0.12359e+08,0.12469e+08,0.12868e+08,0.12776e+08,0.12590e+08]
  ka[:, 2,13] = [0.10045e+04,0.10033e+04,0.99840e+03,0.99184e+03,0.98861e+03]
  ka[:, 3,13] = [0.10050e+04,0.10041e+04,0.10024e+04,0.99701e+03,0.99083e+03]
  ka[:, 4,13] = [0.10053e+04,0.10043e+04,0.10035e+04,0.10017e+04,0.99538e+03]
  ka[:, 5,13] = [0.10048e+04,0.10056e+04,0.10038e+04,0.10026e+04,0.10008e+04]
  ka[:, 6,13] = [0.10042e+04,0.10054e+04,0.10054e+04,0.10039e+04,0.10017e+04]
  ka[:, 7,13] = [0.10029e+04,0.10046e+04,0.10055e+04,0.10046e+04,0.10033e+04]
  ka[:, 8,13] = [0.10002e+04,0.10039e+04,0.10051e+04,0.10054e+04,0.10046e+04]
  ka[:, 9,13] = [0.99825e+03,0.10019e+04,0.10043e+04,0.10051e+04,0.10052e+04]
  ka[:,10,13] = [0.99595e+03,0.99983e+03,0.10034e+04,0.10050e+04,0.10051e+04]
  ka[:,11,13] = [0.99621e+03,0.99979e+03,0.10034e+04,0.10050e+04,0.10054e+04]
  ka[:,12,13] = [0.99597e+03,0.99979e+03,0.10031e+04,0.10047e+04,0.10054e+04]
  ka[:,13,13] = [0.99623e+03,0.99978e+03,0.10031e+04,0.10047e+04,0.10054e+04]
  ka[:, 1,14] = [0.13636e+08,0.13818e+08,0.13518e+08,0.13266e+08,0.13048e+08]
  ka[:, 2,14] = [0.10049e+04,0.99946e+03,0.99774e+03,0.99774e+03,0.99771e+03]
  ka[:, 3,14] = [0.10051e+04,0.10068e+04,0.10007e+04,0.99775e+03,0.99778e+03]
  ka[:, 4,14] = [0.10071e+04,0.10061e+04,0.10050e+04,0.10005e+04,0.99776e+03]
  ka[:, 5,14] = [0.10067e+04,0.10072e+04,0.10061e+04,0.10044e+04,0.99972e+03]
  ka[:, 6,14] = [0.10049e+04,0.10070e+04,0.10074e+04,0.10057e+04,0.10037e+04]
  ka[:, 7,14] = [0.10033e+04,0.10063e+04,0.10072e+04,0.10067e+04,0.10050e+04]
  ka[:, 8,14] = [0.10025e+04,0.10048e+04,0.10068e+04,0.10071e+04,0.10062e+04]
  ka[:, 9,14] = [0.99880e+03,0.10031e+04,0.10060e+04,0.10071e+04,0.10069e+04]
  ka[:,10,14] = [0.99610e+03,0.10023e+04,0.10047e+04,0.10067e+04,0.10072e+04]
  ka[:,11,14] = [0.99601e+03,0.10023e+04,0.10047e+04,0.10067e+04,0.10072e+04]
  ka[:,12,14] = [0.99598e+03,0.10023e+04,0.10051e+04,0.10067e+04,0.10076e+04]
  ka[:,13,14] = [0.99596e+03,0.10027e+04,0.10047e+04,0.10067e+04,0.10072e+04]
  ka[:, 1,15] = [0.14285e+08,0.13947e+08,0.13644e+08,0.13317e+08,0.13015e+08]
  ka[:, 2,15] = [0.99932e+03,0.10011e+04,0.10018e+04,0.10007e+04,0.99706e+03]
  ka[:, 3,15] = [0.10097e+04,0.99743e+03,0.10015e+04,0.10018e+04,0.99920e+03]
  ka[:, 4,15] = [0.10099e+04,0.10076e+04,0.10030e+04,0.10017e+04,0.10013e+04]
  ka[:, 5,15] = [0.10075e+04,0.10079e+04,0.10073e+04,0.10035e+04,0.10018e+04]
  ka[:, 6,15] = [0.10068e+04,0.10088e+04,0.10068e+04,0.10090e+04,0.10033e+04]
  ka[:, 7,15] = [0.10057e+04,0.10073e+04,0.10090e+04,0.10087e+04,0.10062e+04]
  ka[:, 8,15] = [0.10010e+04,0.10054e+04,0.10076e+04,0.10079e+04,0.10072e+04]
  ka[:, 9,15] = [0.99853e+03,0.10054e+04,0.10071e+04,0.10079e+04,0.10078e+04]
  ka[:,10,15] = [0.99742e+03,0.10007e+04,0.10064e+04,0.10076e+04,0.10080e+04]
  ka[:,11,15] = [0.99635e+03,0.10006e+04,0.10064e+04,0.10087e+04,0.10090e+04]
  ka[:,12,15] = [0.99635e+03,0.10006e+04,0.10064e+04,0.10087e+04,0.10069e+04]
  ka[:,13,15] = [0.99635e+03,0.99960e+03,0.10064e+04,0.10076e+04,0.10080e+04]
  ka[:, 1,16] = [0.14326e+08,0.13988e+08,0.13684e+08,0.13409e+08,0.13021e+08]
  ka[:, 2,16] = [0.10008e+04,0.10027e+04,0.10036e+04,0.10037e+04,0.99673e+03]
  ka[:, 3,16] = [0.99286e+03,0.10016e+04,0.10031e+04,0.10038e+04,0.10035e+04]
  ka[:, 4,16] = [0.99278e+03,0.10078e+04,0.10023e+04,0.10035e+04,0.10038e+04]
  ka[:, 5,16] = [0.10074e+04,0.10079e+04,0.10075e+04,0.10028e+04,0.10037e+04]
  ka[:, 6,16] = [0.10069e+04,0.10002e+04,0.10079e+04,0.99212e+03,0.10032e+04]
  ka[:, 7,16] = [0.10060e+04,0.10072e+04,0.10004e+04,0.10003e+04,0.10063e+04]
  ka[:, 8,16] = [0.10009e+04,0.10066e+04,0.10075e+04,0.10080e+04,0.10074e+04]
  ka[:, 9,16] = [0.99860e+03,0.10057e+04,0.10071e+04,0.10078e+04,0.10079e+04]
  ka[:,10,16] = [0.98908e+03,0.10006e+04,0.10065e+04,0.10075e+04,0.10080e+04]
  ka[:,11,16] = [0.99644e+03,0.10006e+04,0.10065e+04,0.10000e+04,0.10005e+04]
  ka[:,12,16] = [0.99644e+03,0.10006e+04,0.10065e+04,0.10000e+04,0.10080e+04]
  ka[:,13,16] = [0.99644e+03,0.10006e+04,0.10065e+04,0.10075e+04,0.10080e+04]

  #     -----------------------------------------------------------------
  #     The array KB contains absorption coefs at the 16 chosen g-values
  #     for a range of pressure levels < ~100mb and temperatures. The first
  #     index in the array, JT, which runs from 1 to 5, corresponds to
  #     different temperatures.  More specifically, JT = 3 means that the
  #     data are for the reference temperature TREF for this pressure
  #     level, JT = 2 refers to the temperature TREF-15, JT = 1 is for
  #     TREF-30, JT = 4 is for TREF+15:and JT = 5 is for TREF+30.
  #     The second index, JP, runs from 13 to 59 and refers to the JPth
  #     reference pressure level (see taumol.f for the value of these
  #     pressure levels in mb).  The third index, IG, goes from 1 to 16,
  #     and tells us which g-interval the absorption coefficients are for.
  #     -----------------------------------------------------------------
  kb = zeros(5,47,16)
  kb[:,13-12, 1] = [0.16036e+00,0.16399e+00,0.17198e+00,0.18455e+00,0.20173e+00]
  kb[:,14-12, 1] = [0.16050e+00,0.16454e+00,0.17292e+00,0.18592e+00,0.20351e+00]
  kb[:,15-12, 1] = [0.16071e+00,0.16519e+00,0.17400e+00,0.18750e+00,0.20554e+00]
  kb[:,16-12, 1] = [0.16098e+00,0.16593e+00,0.17517e+00,0.18921e+00,0.20772e+00]
  kb[:,17-12, 1] = [0.16129e+00,0.16670e+00,0.17642e+00,0.19093e+00,0.20991e+00]
  kb[:,18-12, 1] = [0.16164e+00,0.16750e+00,0.17771e+00,0.19267e+00,0.21211e+00]
  kb[:,19-12, 1] = [0.16203e+00,0.16832e+00,0.17901e+00,0.19443e+00,0.21431e+00]
  kb[:,20-12, 1] = [0.16255e+00,0.16934e+00,0.18059e+00,0.19653e+00,0.21694e+00]
  kb[:,21-12, 1] = [0.16313e+00,0.17044e+00,0.18225e+00,0.19873e+00,0.21967e+00]
  kb[:,22-12, 1] = [0.16416e+00,0.17227e+00,0.18498e+00,0.20229e+00,0.22406e+00]
  kb[:,23-12, 1] = [0.16539e+00,0.17431e+00,0.18796e+00,0.20612e+00,0.22875e+00]
  kb[:,24-12, 1] = [0.16684e+00,0.17666e+00,0.19125e+00,0.21032e+00,0.23385e+00]
  kb[:,25-12, 1] = [0.16856e+00,0.17939e+00,0.19493e+00,0.21495e+00,0.23943e+00]
  kb[:,26-12, 1] = [0.17067e+00,0.18260e+00,0.19919e+00,0.22024e+00,0.24578e+00]
  kb[:,27-12, 1] = [0.17306e+00,0.18613e+00,0.20377e+00,0.22588e+00,0.25248e+00]
  kb[:,28-12, 1] = [0.17570e+00,0.18994e+00,0.20865e+00,0.23184e+00,0.25952e+00]
  kb[:,29-12, 1] = [0.17882e+00,0.19417e+00,0.21399e+00,0.23828e+00,0.26708e+00]
  kb[:,30-12, 1] = [0.18223e+00,0.19870e+00,0.21964e+00,0.24505e+00,0.27492e+00]
  kb[:,31-12, 1] = [0.18602e+00,0.20364e+00,0.22572e+00,0.25229e+00,0.28325e+00]
  kb[:,32-12, 1] = [0.19018e+00,0.20895e+00,0.23219e+00,0.25994e+00,0.29200e+00]
  kb[:,33-12, 1] = [0.19466e+00,0.21461e+00,0.23902e+00,0.26795e+00,0.30114e+00]
  kb[:,34-12, 1] = [0.19880e+00,0.21976e+00,0.24520e+00,0.27509e+00,0.30932e+00]
  kb[:,35-12, 1] = [0.20129e+00,0.22283e+00,0.24886e+00,0.27931e+00,0.31412e+00]
  kb[:,36-12, 1] = [0.20170e+00,0.22334e+00,0.24946e+00,0.28000e+00,0.31492e+00]
  kb[:,37-12, 1] = [0.19920e+00,0.22026e+00,0.24579e+00,0.27578e+00,0.31010e+00]
  kb[:,38-12, 1] = [0.19662e+00,0.21706e+00,0.24197e+00,0.27135e+00,0.30504e+00]
  kb[:,39-12, 1] = [0.19414e+00,0.21396e+00,0.23824e+00,0.26704e+00,0.30010e+00]
  kb[:,40-12, 1] = [0.19003e+00,0.20877e+00,0.23197e+00,0.25968e+00,0.29171e+00]
  kb[:,41-12, 1] = [0.18597e+00,0.20357e+00,0.22564e+00,0.25219e+00,0.28314e+00]
  kb[:,42-12, 1] = [0.18220e+00,0.19865e+00,0.21958e+00,0.24498e+00,0.27484e+00]
  kb[:,43-12, 1] = [0.17798e+00,0.19304e+00,0.21257e+00,0.23657e+00,0.26510e+00]
  kb[:,44-12, 1] = [0.17397e+00,0.18746e+00,0.20548e+00,0.22798e+00,0.25496e+00]
  kb[:,45-12, 1] = [0.17053e+00,0.18239e+00,0.19890e+00,0.21989e+00,0.24535e+00]
  kb[:,46-12, 1] = [0.16739e+00,0.17754e+00,0.19245e+00,0.21183e+00,0.23568e+00]
  kb[:,47-12, 1] = [0.16454e+00,0.17292e+00,0.18592e+00,0.20351e+00,0.22556e+00]
  kb[:,48-12, 1] = [0.16239e+00,0.16903e+00,0.18011e+00,0.19590e+00,0.21615e+00]
  kb[:,49-12, 1] = [0.16095e+00,0.16584e+00,0.17504e+00,0.18902e+00,0.20747e+00]
  kb[:,50-12, 1] = [0.16025e+00,0.16348e+00,0.17109e+00,0.18322e+00,0.20000e+00]
  kb[:,51-12, 1] = [0.16021e+00,0.16179e+00,0.16781e+00,0.17821e+00,0.19335e+00]
  kb[:,52-12, 1] = [0.16078e+00,0.16068e+00,0.16511e+00,0.17387e+00,0.18731e+00]
  kb[:,53-12, 1] = [0.16197e+00,0.16018e+00,0.16299e+00,0.17017e+00,0.18186e+00]
  kb[:,54-12, 1] = [0.16362e+00,0.16026e+00,0.16155e+00,0.16729e+00,0.17738e+00]
  kb[:,55-12, 1] = [0.16570e+00,0.16083e+00,0.16064e+00,0.16498e+00,0.17365e+00]
  kb[:,56-12, 1] = [0.16825e+00,0.16187e+00,0.16020e+00,0.16311e+00,0.17041e+00]
  kb[:,57-12, 1] = [0.17129e+00,0.16340e+00,0.16023e+00,0.16169e+00,0.16761e+00]
  kb[:,58-12, 1] = [0.17462e+00,0.16529e+00,0.16069e+00,0.16077e+00,0.16535e+00]
  kb[:,59-12, 1] = [0.17607e+00,0.16617e+00,0.16100e+00,0.16051e+00,0.16456e+00]
  kb[:,13-12, 2] = [0.17525e+01,0.17727e+01,0.18119e+01,0.18712e+01,0.19491e+01]
  kb[:,14-12, 2] = [0.17536e+01,0.17755e+01,0.18165e+01,0.18775e+01,0.19571e+01]
  kb[:,15-12, 2] = [0.17549e+01,0.17787e+01,0.18219e+01,0.18847e+01,0.19662e+01]
  kb[:,16-12, 2] = [0.17565e+01,0.17819e+01,0.18278e+01,0.18925e+01,0.19761e+01]
  kb[:,17-12, 2] = [0.17583e+01,0.17850e+01,0.18336e+01,0.19003e+01,0.19859e+01]
  kb[:,18-12, 2] = [0.17603e+01,0.17891e+01,0.18396e+01,0.19082e+01,0.19957e+01]
  kb[:,19-12, 2] = [0.17624e+01,0.17934e+01,0.18457e+01,0.19162e+01,0.20056e+01]
  kb[:,20-12, 2] = [0.17652e+01,0.17986e+01,0.18530e+01,0.19257e+01,0.20174e+01]
  kb[:,21-12, 2] = [0.17683e+01,0.18042e+01,0.18607e+01,0.19356e+01,0.20296e+01]
  kb[:,22-12, 2] = [0.17736e+01,0.18133e+01,0.18732e+01,0.19516e+01,0.20492e+01]
  kb[:,23-12, 2] = [0.17796e+01,0.18234e+01,0.18868e+01,0.19689e+01,0.20702e+01]
  kb[:,24-12, 2] = [0.17857e+01,0.18347e+01,0.19018e+01,0.19877e+01,0.20929e+01]
  kb[:,25-12, 2] = [0.17946e+01,0.18475e+01,0.19185e+01,0.20085e+01,0.21178e+01]
  kb[:,26-12, 2] = [0.18053e+01,0.18623e+01,0.19377e+01,0.20322e+01,0.21460e+01]
  kb[:,27-12, 2] = [0.18172e+01,0.18785e+01,0.19583e+01,0.20574e+01,0.21757e+01]
  kb[:,28-12, 2] = [0.18303e+01,0.18959e+01,0.19802e+01,0.20839e+01,0.22069e+01]
  kb[:,29-12, 2] = [0.18448e+01,0.19150e+01,0.20042e+01,0.21126e+01,0.22404e+01]
  kb[:,30-12, 2] = [0.18606e+01,0.19355e+01,0.20295e+01,0.21428e+01,0.22754e+01]
  kb[:,31-12, 2] = [0.18780e+01,0.19577e+01,0.20566e+01,0.21749e+01,0.23125e+01]
  kb[:,32-12, 2] = [0.18969e+01,0.19816e+01,0.20855e+01,0.22088e+01,0.23514e+01]
  kb[:,33-12, 2] = [0.19172e+01,0.20069e+01,0.21160e+01,0.22443e+01,0.23920e+01]
  kb[:,34-12, 2] = [0.19359e+01,0.20300e+01,0.21434e+01,0.22762e+01,0.24282e+01]
  kb[:,35-12, 2] = [0.19471e+01,0.20437e+01,0.21597e+01,0.22949e+01,0.24495e+01]
  kb[:,36-12, 2] = [0.19490e+01,0.20460e+01,0.21624e+01,0.22980e+01,0.24530e+01]
  kb[:,37-12, 2] = [0.19377e+01,0.20322e+01,0.21461e+01,0.22792e+01,0.24317e+01]
  kb[:,38-12, 2] = [0.19261e+01,0.20179e+01,0.21291e+01,0.22595e+01,0.24093e+01]
  kb[:,39-12, 2] = [0.19149e+01,0.20040e+01,0.21125e+01,0.22402e+01,0.23873e+01]
  kb[:,40-12, 2] = [0.18963e+01,0.19807e+01,0.20845e+01,0.22076e+01,0.23501e+01]
  kb[:,41-12, 2] = [0.18778e+01,0.19574e+01,0.20563e+01,0.21745e+01,0.23120e+01]
  kb[:,42-12, 2] = [0.18604e+01,0.19352e+01,0.20292e+01,0.21425e+01,0.22750e+01]
  kb[:,43-12, 2] = [0.18409e+01,0.19099e+01,0.19978e+01,0.21050e+01,0.22315e+01]
  kb[:,44-12, 2] = [0.18217e+01,0.18845e+01,0.19660e+01,0.20667e+01,0.21867e+01]
  kb[:,45-12, 2] = [0.18046e+01,0.18613e+01,0.19364e+01,0.20306e+01,0.21441e+01]
  kb[:,46-12, 2] = [0.17886e+01,0.18389e+01,0.19072e+01,0.19945e+01,0.21011e+01]
  kb[:,47-12, 2] = [0.17755e+01,0.18165e+01,0.18775e+01,0.19571e+01,0.20559e+01]
  kb[:,48-12, 2] = [0.17644e+01,0.17970e+01,0.18508e+01,0.19228e+01,0.20139e+01]
  kb[:,49-12, 2] = [0.17563e+01,0.17816e+01,0.18271e+01,0.18916e+01,0.19749e+01]
  kb[:,50-12, 2] = [0.17517e+01,0.17702e+01,0.18074e+01,0.18651e+01,0.19413e+01]
  kb[:,51-12, 2] = [0.17504e+01,0.17611e+01,0.17908e+01,0.18420e+01,0.19113e+01]
  kb[:,52-12, 2] = [0.17523e+01,0.17547e+01,0.17783e+01,0.18212e+01,0.18839e+01]
  kb[:,53-12, 2] = [0.17566e+01,0.17510e+01,0.17675e+01,0.18028e+01,0.18589e+01]
  kb[:,54-12, 2] = [0.17623e+01,0.17507e+01,0.17598e+01,0.17881e+01,0.18381e+01]
  kb[:,55-12, 2] = [0.17691e+01,0.17524e+01,0.17544e+01,0.17777e+01,0.18201e+01]
  kb[:,56-12, 2] = [0.17776e+01,0.17562e+01,0.17512e+01,0.17682e+01,0.18040e+01]
  kb[:,57-12, 2] = [0.17881e+01,0.17616e+01,0.17505e+01,0.17606e+01,0.17897e+01]
  kb[:,58-12, 2] = [0.17996e+01,0.17678e+01,0.17520e+01,0.17552e+01,0.17794e+01]
  kb[:,59-12, 2] = [0.18045e+01,0.17707e+01,0.17530e+01,0.17536e+01,0.17756e+01]
  kb[:,13-12, 3] = [0.10945e+02,0.11052e+02,0.11222e+02,0.11452e+02,0.11744e+02]
  kb[:,14-12, 3] = [0.10952e+02,0.11065e+02,0.11240e+02,0.11476e+02,0.11774e+02]
  kb[:,15-12, 3] = [0.10961e+02,0.11080e+02,0.11261e+02,0.11503e+02,0.11807e+02]
  kb[:,16-12, 3] = [0.10971e+02,0.11097e+02,0.11284e+02,0.11533e+02,0.11844e+02]
  kb[:,17-12, 3] = [0.10981e+02,0.11115e+02,0.11307e+02,0.11562e+02,0.11879e+02]
  kb[:,18-12, 3] = [0.10992e+02,0.11132e+02,0.11331e+02,0.11592e+02,0.11915e+02]
  kb[:,19-12, 3] = [0.11003e+02,0.11149e+02,0.11354e+02,0.11622e+02,0.11951e+02]
  kb[:,20-12, 3] = [0.11017e+02,0.11170e+02,0.11382e+02,0.11657e+02,0.11994e+02]
  kb[:,21-12, 3] = [0.11031e+02,0.11192e+02,0.11412e+02,0.11694e+02,0.12038e+02]
  kb[:,22-12, 3] = [0.11056e+02,0.11228e+02,0.11460e+02,0.11754e+02,0.12109e+02]
  kb[:,23-12, 3] = [0.11085e+02,0.11267e+02,0.11511e+02,0.11817e+02,0.12184e+02]
  kb[:,24-12, 3] = [0.11118e+02,0.11312e+02,0.11568e+02,0.11886e+02,0.12266e+02]
  kb[:,25-12, 3] = [0.11154e+02,0.11361e+02,0.11630e+02,0.11962e+02,0.12354e+02]
  kb[:,26-12, 3] = [0.11196e+02,0.11418e+02,0.11702e+02,0.12048e+02,0.12455e+02]
  kb[:,27-12, 3] = [0.11243e+02,0.11480e+02,0.11778e+02,0.12138e+02,0.12560e+02]
  kb[:,28-12, 3] = [0.11294e+02,0.11545e+02,0.11859e+02,0.12234e+02,0.12670e+02]
  kb[:,29-12, 3] = [0.11351e+02,0.11617e+02,0.11946e+02,0.12336e+02,0.12788e+02]
  kb[:,30-12, 3] = [0.11411e+02,0.11694e+02,0.12038e+02,0.12443e+02,0.12910e+02]
  kb[:,31-12, 3] = [0.11478e+02,0.11776e+02,0.12136e+02,0.12557e+02,0.13040e+02]
  kb[:,32-12, 3] = [0.11549e+02,0.11864e+02,0.12239e+02,0.12677e+02,0.13175e+02]
  kb[:,33-12, 3] = [0.11626e+02,0.11956e+02,0.12348e+02,0.12801e+02,0.13316e+02]
  kb[:,34-12, 3] = [0.11695e+02,0.12040e+02,0.12446e+02,0.12913e+02,0.13442e+02]
  kb[:,35-12, 3] = [0.11737e+02,0.12089e+02,0.12503e+02,0.12979e+02,0.13515e+02]
  kb[:,36-12, 3] = [0.11744e+02,0.12098e+02,0.12513e+02,0.12989e+02,0.13528e+02]
  kb[:,37-12, 3] = [0.11702e+02,0.12048e+02,0.12455e+02,0.12924e+02,0.13454e+02]
  kb[:,38-12, 3] = [0.11659e+02,0.11996e+02,0.12395e+02,0.12855e+02,0.13376e+02]
  kb[:,39-12, 3] = [0.11617e+02,0.11946e+02,0.12336e+02,0.12787e+02,0.13300e+02]
  kb[:,40-12, 3] = [0.11547e+02,0.11861e+02,0.12236e+02,0.12673e+02,0.13171e+02]
  kb[:,41-12, 3] = [0.11477e+02,0.11775e+02,0.12134e+02,0.12555e+02,0.13038e+02]
  kb[:,42-12, 3] = [0.11411e+02,0.11693e+02,0.12037e+02,0.12442e+02,0.12909e+02]
  kb[:,43-12, 3] = [0.11335e+02,0.11598e+02,0.11923e+02,0.12309e+02,0.12757e+02]
  kb[:,44-12, 3] = [0.11261e+02,0.11503e+02,0.11807e+02,0.12172e+02,0.12599e+02]
  kb[:,45-12, 3] = [0.11193e+02,0.11414e+02,0.11697e+02,0.12042e+02,0.12448e+02]
  kb[:,46-12, 3] = [0.11129e+02,0.11328e+02,0.11588e+02,0.11911e+02,0.12295e+02]
  kb[:,47-12, 3] = [0.11065e+02,0.11240e+02,0.11476e+02,0.11774e+02,0.12133e+02]
  kb[:,48-12, 3] = [0.11012e+02,0.11163e+02,0.11374e+02,0.11647e+02,0.11981e+02]
  kb[:,49-12, 3] = [0.10969e+02,0.11095e+02,0.11282e+02,0.11530e+02,0.11839e+02]
  kb[:,50-12, 3] = [0.10938e+02,0.11040e+02,0.11204e+02,0.11429e+02,0.11715e+02]
  kb[:,51-12, 3] = [0.10916e+02,0.10996e+02,0.11138e+02,0.11340e+02,0.11604e+02]
  kb[:,52-12, 3] = [0.10900e+02,0.10960e+02,0.11079e+02,0.11259e+02,0.11500e+02]
  kb[:,53-12, 3] = [0.10894e+02,0.10932e+02,0.11028e+02,0.11186e+02,0.11405e+02]
  kb[:,54-12, 3] = [0.10895e+02,0.10912e+02,0.10989e+02,0.11127e+02,0.11325e+02]
  kb[:,55-12, 3] = [0.10903e+02,0.10900e+02,0.10958e+02,0.11075e+02,0.11255e+02]
  kb[:,56-12, 3] = [0.10921e+02,0.10894e+02,0.10933e+02,0.11031e+02,0.11191e+02]
  kb[:,57-12, 3] = [0.10944e+02,0.10894e+02,0.10914e+02,0.10993e+02,0.11134e+02]
  kb[:,58-12, 3] = [0.10972e+02,0.10901e+02,0.10902e+02,0.10963e+02,0.11084e+02]
  kb[:,59-12, 3] = [0.10984e+02,0.10906e+02,0.10898e+02,0.10952e+02,0.11066e+02]
  kb[:,13-12, 4] = [0.51563e+02,0.51946e+02,0.52464e+02,0.53116e+02,0.53902e+02]
  kb[:,14-12, 4] = [0.51593e+02,0.51987e+02,0.52517e+02,0.53182e+02,0.53980e+02]
  kb[:,15-12, 4] = [0.51626e+02,0.52035e+02,0.52579e+02,0.53256e+02,0.54069e+02]
  kb[:,16-12, 4] = [0.51664e+02,0.52087e+02,0.52644e+02,0.53337e+02,0.54163e+02]
  kb[:,17-12, 4] = [0.51701e+02,0.52138e+02,0.52710e+02,0.53416e+02,0.54257e+02]
  kb[:,18-12, 4] = [0.51740e+02,0.52191e+02,0.52776e+02,0.53496e+02,0.54351e+02]
  kb[:,19-12, 4] = [0.51779e+02,0.52244e+02,0.52843e+02,0.53576e+02,0.54444e+02]
  kb[:,20-12, 4] = [0.51827e+02,0.52307e+02,0.52922e+02,0.53671e+02,0.54555e+02]
  kb[:,21-12, 4] = [0.51877e+02,0.52373e+02,0.53004e+02,0.53770e+02,0.54669e+02]
  kb[:,22-12, 4] = [0.51959e+02,0.52481e+02,0.53136e+02,0.53927e+02,0.54851e+02]
  kb[:,23-12, 4] = [0.52049e+02,0.52596e+02,0.53278e+02,0.54094e+02,0.55044e+02]
  kb[:,24-12, 4] = [0.52148e+02,0.52723e+02,0.53431e+02,0.54274e+02,0.55252e+02]
  kb[:,25-12, 4] = [0.52259e+02,0.52862e+02,0.53599e+02,0.54471e+02,0.55477e+02]
  kb[:,26-12, 4] = [0.52387e+02,0.53022e+02,0.53790e+02,0.54693e+02,0.55730e+02]
  kb[:,27-12, 4] = [0.52525e+02,0.53191e+02,0.53992e+02,0.54926e+02,0.55995e+02]
  kb[:,28-12, 4] = [0.52672e+02,0.53371e+02,0.54203e+02,0.55170e+02,0.56271e+02]
  kb[:,29-12, 4] = [0.52833e+02,0.53565e+02,0.54430e+02,0.55430e+02,0.56565e+02]
  kb[:,30-12, 4] = [0.53003e+02,0.53768e+02,0.54668e+02,0.55702e+02,0.56870e+02]
  kb[:,31-12, 4] = [0.53186e+02,0.53986e+02,0.54920e+02,0.55988e+02,0.57190e+02]
  kb[:,32-12, 4] = [0.53381e+02,0.54216e+02,0.55184e+02,0.56287e+02,0.57525e+02]
  kb[:,33-12, 4] = [0.53587e+02,0.54457e+02,0.55460e+02,0.56599e+02,0.57871e+02]
  kb[:,34-12, 4] = [0.53773e+02,0.54673e+02,0.55707e+02,0.56876e+02,0.58179e+02]
  kb[:,35-12, 4] = [0.53883e+02,0.54800e+02,0.55852e+02,0.57039e+02,0.58359e+02]
  kb[:,36-12, 4] = [0.53901e+02,0.54822e+02,0.55876e+02,0.57066e+02,0.58389e+02]
  kb[:,37-12, 4] = [0.53791e+02,0.54694e+02,0.55731e+02,0.56903e+02,0.58208e+02]
  kb[:,38-12, 4] = [0.53676e+02,0.54560e+02,0.55578e+02,0.56731e+02,0.58018e+02]
  kb[:,39-12, 4] = [0.53563e+02,0.54429e+02,0.55429e+02,0.56563e+02,0.57832e+02]
  kb[:,40-12, 4] = [0.53375e+02,0.54208e+02,0.55175e+02,0.56277e+02,0.57513e+02]
  kb[:,41-12, 4] = [0.53184e+02,0.53983e+02,0.54916e+02,0.55984e+02,0.57186e+02]
  kb[:,42-12, 4] = [0.53001e+02,0.53766e+02,0.54665e+02,0.55699e+02,0.56866e+02]
  kb[:,43-12, 4] = [0.52790e+02,0.53513e+02,0.54370e+02,0.55362e+02,0.56488e+02]
  kb[:,44-12, 4] = [0.52577e+02,0.53254e+02,0.54066e+02,0.55012e+02,0.56093e+02]
  kb[:,45-12, 4] = [0.52379e+02,0.53011e+02,0.53777e+02,0.54678e+02,0.55713e+02]
  kb[:,46-12, 4] = [0.52184e+02,0.52768e+02,0.53486e+02,0.54339e+02,0.55326e+02]
  kb[:,47-12, 4] = [0.51987e+02,0.52517e+02,0.53182e+02,0.53980e+02,0.54913e+02]
  kb[:,48-12, 4] = [0.51812e+02,0.52288e+02,0.52898e+02,0.53643e+02,0.54522e+02]
  kb[:,49-12, 4] = [0.51659e+02,0.52081e+02,0.52637e+02,0.53327e+02,0.54152e+02]
  kb[:,50-12, 4] = [0.51535e+02,0.51906e+02,0.52412e+02,0.53052e+02,0.53826e+02]
  kb[:,51-12, 4] = [0.51433e+02,0.51755e+02,0.52211e+02,0.52802e+02,0.53527e+02]
  kb[:,52-12, 4] = [0.51349e+02,0.51622e+02,0.52029e+02,0.52571e+02,0.53247e+02]
  kb[:,53-12, 4] = [0.51285e+02,0.51507e+02,0.51865e+02,0.52358e+02,0.52985e+02]
  kb[:,54-12, 4] = [0.51240e+02,0.51417e+02,0.51730e+02,0.52178e+02,0.52760e+02]
  kb[:,55-12, 4] = [0.51211e+02,0.51346e+02,0.51615e+02,0.52020e+02,0.52559e+02]
  kb[:,56-12, 4] = [0.51194e+02,0.51288e+02,0.51514e+02,0.51875e+02,0.52371e+02]
  kb[:,57-12, 4] = [0.51190e+02,0.51245e+02,0.51426e+02,0.51746e+02,0.52198e+02]
  kb[:,58-12, 4] = [0.51200e+02,0.51216e+02,0.51357e+02,0.51635e+02,0.52046e+02]
  kb[:,59-12, 4] = [0.51207e+02,0.51207e+02,0.51333e+02,0.51594e+02,0.51989e+02]
  kb[:,13-12, 5] = [0.17051e+03,0.17168e+03,0.17297e+03,0.17437e+03,0.17590e+03]
  kb[:,14-12, 5] = [0.17061e+03,0.17179e+03,0.17309e+03,0.17451e+03,0.17604e+03]
  kb[:,15-12, 5] = [0.17073e+03,0.17192e+03,0.17323e+03,0.17466e+03,0.17621e+03]
  kb[:,16-12, 5] = [0.17085e+03,0.17205e+03,0.17338e+03,0.17482e+03,0.17638e+03]
  kb[:,17-12, 5] = [0.17097e+03,0.17218e+03,0.17352e+03,0.17497e+03,0.17655e+03]
  kb[:,18-12, 5] = [0.17109e+03,0.17232e+03,0.17366e+03,0.17513e+03,0.17671e+03]
  kb[:,19-12, 5] = [0.17121e+03,0.17245e+03,0.17381e+03,0.17528e+03,0.17688e+03]
  kb[:,20-12, 5] = [0.17134e+03,0.17260e+03,0.17397e+03,0.17547e+03,0.17708e+03]
  kb[:,21-12, 5] = [0.17149e+03,0.17276e+03,0.17415e+03,0.17565e+03,0.17728e+03]
  kb[:,22-12, 5] = [0.17172e+03,0.17301e+03,0.17442e+03,0.17594e+03,0.17759e+03]
  kb[:,23-12, 5] = [0.17195e+03,0.17327e+03,0.17470e+03,0.17625e+03,0.17792e+03]
  kb[:,24-12, 5] = [0.17221e+03,0.17355e+03,0.17500e+03,0.17658e+03,0.17827e+03]
  kb[:,25-12, 5] = [0.17248e+03,0.17385e+03,0.17533e+03,0.17693e+03,0.17865e+03]
  kb[:,26-12, 5] = [0.17279e+03,0.17418e+03,0.17569e+03,0.17732e+03,0.17906e+03]
  kb[:,27-12, 5] = [0.17311e+03,0.17453e+03,0.17606e+03,0.17772e+03,0.17949e+03]
  kb[:,28-12, 5] = [0.17344e+03,0.17488e+03,0.17645e+03,0.17813e+03,0.17994e+03]
  kb[:,29-12, 5] = [0.17379e+03,0.17526e+03,0.17686e+03,0.17857e+03,0.18040e+03]
  kb[:,30-12, 5] = [0.17414e+03,0.17565e+03,0.17727e+03,0.17902e+03,0.18088e+03]
  kb[:,31-12, 5] = [0.17452e+03,0.17605e+03,0.17771e+03,0.17948e+03,0.18137e+03]
  kb[:,32-12, 5] = [0.17491e+03,0.17647e+03,0.17816e+03,0.17996e+03,0.18188e+03]
  kb[:,33-12, 5] = [0.17531e+03,0.17690e+03,0.17862e+03,0.18045e+03,0.18241e+03]
  kb[:,34-12, 5] = [0.17566e+03,0.17728e+03,0.17903e+03,0.18089e+03,0.18287e+03]
  kb[:,35-12, 5] = [0.17586e+03,0.17750e+03,0.17926e+03,0.18114e+03,0.18313e+03]
  kb[:,36-12, 5] = [0.17590e+03,0.17754e+03,0.17930e+03,0.18118e+03,0.18318e+03]
  kb[:,37-12, 5] = [0.17569e+03,0.17732e+03,0.17906e+03,0.18093e+03,0.18291e+03]
  kb[:,38-12, 5] = [0.17547e+03,0.17708e+03,0.17881e+03,0.18066e+03,0.18263e+03]
  kb[:,39-12, 5] = [0.17526e+03,0.17685e+03,0.17857e+03,0.18040e+03,0.18235e+03]
  kb[:,40-12, 5] = [0.17489e+03,0.17646e+03,0.17814e+03,0.17995e+03,0.18187e+03]
  kb[:,41-12, 5] = [0.17451e+03,0.17605e+03,0.17770e+03,0.17948e+03,0.18137e+03]
  kb[:,42-12, 5] = [0.17414e+03,0.17565e+03,0.17727e+03,0.17901e+03,0.18087e+03]
  kb[:,43-12, 5] = [0.17369e+03,0.17516e+03,0.17675e+03,0.17846e+03,0.18028e+03]
  kb[:,44-12, 5] = [0.17323e+03,0.17465e+03,0.17620e+03,0.17787e+03,0.17965e+03]
  kb[:,45-12, 5] = [0.17277e+03,0.17416e+03,0.17567e+03,0.17729e+03,0.17904e+03]
  kb[:,46-12, 5] = [0.17230e+03,0.17365e+03,0.17511e+03,0.17669e+03,0.17840e+03]
  kb[:,47-12, 5] = [0.17179e+03,0.17309e+03,0.17451e+03,0.17604e+03,0.17770e+03]
  kb[:,48-12, 5] = [0.17130e+03,0.17255e+03,0.17392e+03,0.17541e+03,0.17702e+03]
  kb[:,49-12, 5] = [0.17083e+03,0.17204e+03,0.17336e+03,0.17480e+03,0.17636e+03]
  kb[:,50-12, 5] = [0.17041e+03,0.17157e+03,0.17285e+03,0.17424e+03,0.17576e+03]
  kb[:,51-12, 5] = [0.17002e+03,0.17113e+03,0.17237e+03,0.17372e+03,0.17519e+03]
  kb[:,52-12, 5] = [0.16964e+03,0.17071e+03,0.17190e+03,0.17321e+03,0.17464e+03]
  kb[:,53-12, 5] = [0.16927e+03,0.17031e+03,0.17145e+03,0.17272e+03,0.17411e+03]
  kb[:,54-12, 5] = [0.16896e+03,0.16995e+03,0.17106e+03,0.17228e+03,0.17363e+03]
  kb[:,55-12, 5] = [0.16866e+03,0.16962e+03,0.17069e+03,0.17188e+03,0.17318e+03]
  kb[:,56-12, 5] = [0.16839e+03,0.16930e+03,0.17033e+03,0.17148e+03,0.17275e+03]
  kb[:,57-12, 5] = [0.16812e+03,0.16899e+03,0.16999e+03,0.17110e+03,0.17233e+03]
  kb[:,58-12, 5] = [0.16788e+03,0.16871e+03,0.16968e+03,0.17075e+03,0.17195e+03]
  kb[:,59-12, 5] = [0.16779e+03,0.16861e+03,0.16955e+03,0.17062e+03,0.17180e+03]
  kb[:,13-12, 6] = [0.38122e+03,0.38406e+03,0.38668e+03,0.38908e+03,0.39125e+03]
  kb[:,14-12, 6] = [0.38148e+03,0.38431e+03,0.38691e+03,0.38929e+03,0.39144e+03]
  kb[:,15-12, 6] = [0.38178e+03,0.38458e+03,0.38716e+03,0.38951e+03,0.39164e+03]
  kb[:,16-12, 6] = [0.38209e+03,0.38487e+03,0.38742e+03,0.38975e+03,0.39186e+03]
  kb[:,17-12, 6] = [0.38239e+03,0.38514e+03,0.38767e+03,0.38998e+03,0.39206e+03]
  kb[:,18-12, 6] = [0.38268e+03,0.38541e+03,0.38792e+03,0.39020e+03,0.39226e+03]
  kb[:,19-12, 6] = [0.38297e+03,0.38568e+03,0.38816e+03,0.39042e+03,0.39246e+03]
  kb[:,20-12, 6] = [0.38330e+03,0.38598e+03,0.38844e+03,0.39067e+03,0.39268e+03]
  kb[:,21-12, 6] = [0.38363e+03,0.38629e+03,0.38872e+03,0.39092e+03,0.39291e+03]
  kb[:,22-12, 6] = [0.38414e+03,0.38675e+03,0.38914e+03,0.39131e+03,0.39325e+03]
  kb[:,23-12, 6] = [0.38466e+03,0.38723e+03,0.38958e+03,0.39170e+03,0.39360e+03]
  kb[:,24-12, 6] = [0.38520e+03,0.38772e+03,0.39002e+03,0.39210e+03,0.39395e+03]
  kb[:,25-12, 6] = [0.38575e+03,0.38823e+03,0.39048e+03,0.39251e+03,0.39432e+03]
  kb[:,26-12, 6] = [0.38635e+03,0.38877e+03,0.39098e+03,0.39295e+03,0.39470e+03]
  kb[:,27-12, 6] = [0.38694e+03,0.38932e+03,0.39146e+03,0.39339e+03,0.39509e+03]
  kb[:,28-12, 6] = [0.38753e+03,0.38985e+03,0.39194e+03,0.39381e+03,0.39546e+03]
  kb[:,29-12, 6] = [0.38813e+03,0.39039e+03,0.39243e+03,0.39424e+03,0.39583e+03]
  kb[:,30-12, 6] = [0.38871e+03,0.39092e+03,0.39290e+03,0.39466e+03,0.39620e+03]
  kb[:,31-12, 6] = [0.38930e+03,0.39145e+03,0.39338e+03,0.39508e+03,0.39655e+03]
  kb[:,32-12, 6] = [0.38988e+03,0.39197e+03,0.39384e+03,0.39548e+03,0.39690e+03]
  kb[:,33-12, 6] = [0.39045e+03,0.39248e+03,0.39429e+03,0.39587e+03,0.39723e+03]
  kb[:,34-12, 6] = [0.39093e+03,0.39291e+03,0.39467e+03,0.39620e+03,0.39751e+03]
  kb[:,35-12, 6] = [0.39120e+03,0.39316e+03,0.39488e+03,0.39639e+03,0.39766e+03]
  kb[:,36-12, 6] = [0.39125e+03,0.39320e+03,0.39492e+03,0.39642e+03,0.39769e+03]
  kb[:,37-12, 6] = [0.39098e+03,0.39295e+03,0.39471e+03,0.39623e+03,0.39754e+03]
  kb[:,38-12, 6] = [0.39068e+03,0.39269e+03,0.39447e+03,0.39603e+03,0.39737e+03]
  kb[:,39-12, 6] = [0.39039e+03,0.39243e+03,0.39424e+03,0.39583e+03,0.39720e+03]
  kb[:,40-12, 6] = [0.38986e+03,0.39196e+03,0.39382e+03,0.39547e+03,0.39689e+03]
  kb[:,41-12, 6] = [0.38929e+03,0.39144e+03,0.39337e+03,0.39507e+03,0.39655e+03]
  kb[:,42-12, 6] = [0.38871e+03,0.39092e+03,0.39290e+03,0.39466e+03,0.39619e+03]
  kb[:,43-12, 6] = [0.38797e+03,0.39025e+03,0.39230e+03,0.39413e+03,0.39574e+03]
  kb[:,44-12, 6] = [0.38715e+03,0.38951e+03,0.39164e+03,0.39354e+03,0.39522e+03]
  kb[:,45-12, 6] = [0.38631e+03,0.38874e+03,0.39094e+03,0.39292e+03,0.39468e+03]
  kb[:,46-12, 6] = [0.38538e+03,0.38789e+03,0.39018e+03,0.39224e+03,0.39407e+03]
  kb[:,47-12, 6] = [0.38431e+03,0.38691e+03,0.38929e+03,0.39144e+03,0.39336e+03]
  kb[:,48-12, 6] = [0.38320e+03,0.38589e+03,0.38836e+03,0.39060e+03,0.39262e+03]
  kb[:,49-12, 6] = [0.38206e+03,0.38484e+03,0.38739e+03,0.38972e+03,0.39183e+03]
  kb[:,50-12, 6] = [0.38095e+03,0.38382e+03,0.38646e+03,0.38887e+03,0.39106e+03]
  kb[:,51-12, 6] = [0.37985e+03,0.38279e+03,0.38552e+03,0.38802e+03,0.39029e+03]
  kb[:,52-12, 6] = [0.37873e+03,0.38175e+03,0.38455e+03,0.38713e+03,0.38949e+03]
  kb[:,53-12, 6] = [0.37762e+03,0.38067e+03,0.38355e+03,0.38621e+03,0.38865e+03]
  kb[:,54-12, 6] = [0.37661e+03,0.37965e+03,0.38261e+03,0.38535e+03,0.38786e+03]
  kb[:,55-12, 6] = [0.37559e+03,0.37867e+03,0.38169e+03,0.38450e+03,0.38708e+03]
  kb[:,56-12, 6] = [0.37459e+03,0.37768e+03,0.38074e+03,0.38362e+03,0.38628e+03]
  kb[:,57-12, 6] = [0.37353e+03,0.37673e+03,0.37977e+03,0.38272e+03,0.38545e+03]
  kb[:,58-12, 6] = [0.37251e+03,0.37577e+03,0.37885e+03,0.38185e+03,0.38465e+03]
  kb[:,59-12, 6] = [0.37212e+03,0.37539e+03,0.37847e+03,0.38150e+03,0.38432e+03]
  kb[:,13-12, 7] = [0.62583e+03,0.62623e+03,0.62691e+03,0.62787e+03,0.62909e+03]
  kb[:,14-12, 7] = [0.62586e+03,0.62628e+03,0.62699e+03,0.62797e+03,0.62921e+03]
  kb[:,15-12, 7] = [0.62589e+03,0.62634e+03,0.62708e+03,0.62809e+03,0.62935e+03]
  kb[:,16-12, 7] = [0.62593e+03,0.62640e+03,0.62717e+03,0.62821e+03,0.62950e+03]
  kb[:,17-12, 7] = [0.62597e+03,0.62647e+03,0.62727e+03,0.62833e+03,0.62965e+03]
  kb[:,18-12, 7] = [0.62601e+03,0.62654e+03,0.62737e+03,0.62846e+03,0.62980e+03]
  kb[:,19-12, 7] = [0.62605e+03,0.62661e+03,0.62746e+03,0.62858e+03,0.62995e+03]
  kb[:,20-12, 7] = [0.62610e+03,0.62670e+03,0.62758e+03,0.62873e+03,0.63012e+03]
  kb[:,21-12, 7] = [0.62615e+03,0.62679e+03,0.62770e+03,0.62888e+03,0.63031e+03]
  kb[:,22-12, 7] = [0.62625e+03,0.62694e+03,0.62790e+03,0.62913e+03,0.63060e+03]
  kb[:,23-12, 7] = [0.62636e+03,0.62710e+03,0.62812e+03,0.62939e+03,0.63091e+03]
  kb[:,24-12, 7] = [0.62648e+03,0.62729e+03,0.62835e+03,0.62968e+03,0.63125e+03]
  kb[:,25-12, 7] = [0.62663e+03,0.62749e+03,0.62861e+03,0.62999e+03,0.63162e+03]
  kb[:,26-12, 7] = [0.62681e+03,0.62773e+03,0.62891e+03,0.63035e+03,0.63203e+03]
  kb[:,27-12, 7] = [0.62700e+03,0.62799e+03,0.62923e+03,0.63072e+03,0.63247e+03]
  kb[:,28-12, 7] = [0.62721e+03,0.62826e+03,0.62956e+03,0.63112e+03,0.63292e+03]
  kb[:,29-12, 7] = [0.62745e+03,0.62856e+03,0.62993e+03,0.63154e+03,0.63341e+03]
  kb[:,30-12, 7] = [0.62770e+03,0.62888e+03,0.63031e+03,0.63198e+03,0.63392e+03]
  kb[:,31-12, 7] = [0.62798e+03,0.62922e+03,0.63071e+03,0.63245e+03,0.63446e+03]
  kb[:,32-12, 7] = [0.62828e+03,0.62958e+03,0.63114e+03,0.63295e+03,0.63503e+03]
  kb[:,33-12, 7] = [0.62860e+03,0.62997e+03,0.63159e+03,0.63347e+03,0.63562e+03]
  kb[:,34-12, 7] = [0.62889e+03,0.63031e+03,0.63199e+03,0.63393e+03,0.63614e+03]
  kb[:,35-12, 7] = [0.62906e+03,0.63052e+03,0.63223e+03,0.63421e+03,0.63644e+03]
  kb[:,36-12, 7] = [0.62909e+03,0.63055e+03,0.63227e+03,0.63425e+03,0.63649e+03]
  kb[:,37-12, 7] = [0.62891e+03,0.63035e+03,0.63203e+03,0.63397e+03,0.63619e+03]
  kb[:,38-12, 7] = [0.62873e+03,0.63013e+03,0.63178e+03,0.63369e+03,0.63586e+03]
  kb[:,39-12, 7] = [0.62856e+03,0.62992e+03,0.63154e+03,0.63341e+03,0.63555e+03]
  kb[:,40-12, 7] = [0.62827e+03,0.62957e+03,0.63113e+03,0.63293e+03,0.63501e+03]
  kb[:,41-12, 7] = [0.62798e+03,0.62922e+03,0.63071e+03,0.63245e+03,0.63445e+03]
  kb[:,42-12, 7] = [0.62770e+03,0.62887e+03,0.63030e+03,0.63198e+03,0.63391e+03]
  kb[:,43-12, 7] = [0.62739e+03,0.62848e+03,0.62983e+03,0.63143e+03,0.63328e+03]
  kb[:,44-12, 7] = [0.62708e+03,0.62808e+03,0.62935e+03,0.63086e+03,0.63263e+03]
  kb[:,45-12, 7] = [0.62679e+03,0.62771e+03,0.62889e+03,0.63032e+03,0.63200e+03]
  kb[:,46-12, 7] = [0.62653e+03,0.62735e+03,0.62844e+03,0.62978e+03,0.63137e+03]
  kb[:,47-12, 7] = [0.62629e+03,0.62699e+03,0.62797e+03,0.62921e+03,0.63070e+03]
  kb[:,48-12, 7] = [0.62608e+03,0.62667e+03,0.62755e+03,0.62868e+03,0.63007e+03]
  kb[:,49-12, 7] = [0.62592e+03,0.62640e+03,0.62716e+03,0.62820e+03,0.62948e+03]
  kb[:,50-12, 7] = [0.62580e+03,0.62619e+03,0.62684e+03,0.62778e+03,0.62897e+03]
  kb[:,51-12, 7] = [0.62580e+03,0.62602e+03,0.62657e+03,0.62741e+03,0.62850e+03]
  kb[:,52-12, 7] = [0.62578e+03,0.62589e+03,0.62633e+03,0.62707e+03,0.62807e+03]
  kb[:,53-12, 7] = [0.62574e+03,0.62578e+03,0.62614e+03,0.62677e+03,0.62768e+03]
  kb[:,54-12, 7] = [0.62576e+03,0.62580e+03,0.62600e+03,0.62652e+03,0.62734e+03]
  kb[:,55-12, 7] = [0.62585e+03,0.62578e+03,0.62588e+03,0.62632e+03,0.62705e+03]
  kb[:,56-12, 7] = [0.62587e+03,0.62576e+03,0.62579e+03,0.62615e+03,0.62678e+03]
  kb[:,57-12, 7] = [0.62591e+03,0.62574e+03,0.62580e+03,0.62601e+03,0.62655e+03]
  kb[:,58-12, 7] = [0.62591e+03,0.62583e+03,0.62578e+03,0.62590e+03,0.62636e+03]
  kb[:,59-12, 7] = [0.62590e+03,0.62586e+03,0.62578e+03,0.62586e+03,0.62628e+03]
  kb[:,13-12, 8] = [0.81821e+03,0.81877e+03,0.81935e+03,0.81999e+03,0.82074e+03]
  kb[:,14-12, 8] = [0.81826e+03,0.81882e+03,0.81939e+03,0.82006e+03,0.82081e+03]
  kb[:,15-12, 8] = [0.81831e+03,0.81888e+03,0.81946e+03,0.82012e+03,0.82089e+03]
  kb[:,16-12, 8] = [0.81836e+03,0.81895e+03,0.81952e+03,0.82021e+03,0.82098e+03]
  kb[:,17-12, 8] = [0.81842e+03,0.81901e+03,0.81959e+03,0.82028e+03,0.82106e+03]
  kb[:,18-12, 8] = [0.81847e+03,0.81907e+03,0.81965e+03,0.82036e+03,0.82116e+03]
  kb[:,19-12, 8] = [0.81853e+03,0.81913e+03,0.81972e+03,0.82043e+03,0.82123e+03]
  kb[:,20-12, 8] = [0.81860e+03,0.81919e+03,0.81980e+03,0.82052e+03,0.82133e+03]
  kb[:,21-12, 8] = [0.81868e+03,0.81925e+03,0.81988e+03,0.82061e+03,0.82143e+03]
  kb[:,22-12, 8] = [0.81879e+03,0.81936e+03,0.82001e+03,0.82076e+03,0.82160e+03]
  kb[:,23-12, 8] = [0.81889e+03,0.81947e+03,0.82015e+03,0.82091e+03,0.82177e+03]
  kb[:,24-12, 8] = [0.81902e+03,0.81960e+03,0.82030e+03,0.82108e+03,0.82195e+03]
  kb[:,25-12, 8] = [0.81914e+03,0.81974e+03,0.82046e+03,0.82126e+03,0.82214e+03]
  kb[:,26-12, 8] = [0.81927e+03,0.81990e+03,0.82063e+03,0.82145e+03,0.82236e+03]
  kb[:,27-12, 8] = [0.81941e+03,0.82006e+03,0.82082e+03,0.82166e+03,0.82258e+03]
  kb[:,28-12, 8] = [0.81957e+03,0.82024e+03,0.82101e+03,0.82188e+03,0.82281e+03]
  kb[:,29-12, 8] = [0.81971e+03,0.82042e+03,0.82122e+03,0.82210e+03,0.82305e+03]
  kb[:,30-12, 8] = [0.81988e+03,0.82061e+03,0.82143e+03,0.82233e+03,0.82330e+03]
  kb[:,31-12, 8] = [0.82006e+03,0.82081e+03,0.82166e+03,0.82258e+03,0.82355e+03]
  kb[:,32-12, 8] = [0.82025e+03,0.82103e+03,0.82189e+03,0.82283e+03,0.82380e+03]
  kb[:,33-12, 8] = [0.82044e+03,0.82124e+03,0.82213e+03,0.82308e+03,0.82408e+03]
  kb[:,34-12, 8] = [0.82062e+03,0.82144e+03,0.82234e+03,0.82330e+03,0.82432e+03]
  kb[:,35-12, 8] = [0.82072e+03,0.82155e+03,0.82246e+03,0.82342e+03,0.82446e+03]
  kb[:,36-12, 8] = [0.82074e+03,0.82157e+03,0.82248e+03,0.82345e+03,0.82448e+03]
  kb[:,37-12, 8] = [0.82063e+03,0.82146e+03,0.82236e+03,0.82332e+03,0.82434e+03]
  kb[:,38-12, 8] = [0.82053e+03,0.82134e+03,0.82223e+03,0.82319e+03,0.82419e+03]
  kb[:,39-12, 8] = [0.82042e+03,0.82122e+03,0.82210e+03,0.82305e+03,0.82405e+03]
  kb[:,40-12, 8] = [0.82024e+03,0.82102e+03,0.82188e+03,0.82282e+03,0.82380e+03]
  kb[:,41-12, 8] = [0.82006e+03,0.82081e+03,0.82165e+03,0.82257e+03,0.82354e+03]
  kb[:,42-12, 8] = [0.81988e+03,0.82061e+03,0.82143e+03,0.82233e+03,0.82329e+03]
  kb[:,43-12, 8] = [0.81967e+03,0.82037e+03,0.82117e+03,0.82204e+03,0.82299e+03]
  kb[:,44-12, 8] = [0.81946e+03,0.82013e+03,0.82089e+03,0.82174e+03,0.82267e+03]
  kb[:,45-12, 8] = [0.81926e+03,0.81989e+03,0.82062e+03,0.82144e+03,0.82234e+03]
  kb[:,46-12, 8] = [0.81907e+03,0.81965e+03,0.82034e+03,0.82114e+03,0.82201e+03]
  kb[:,47-12, 8] = [0.81882e+03,0.81939e+03,0.82006e+03,0.82081e+03,0.82165e+03]
  kb[:,48-12, 8] = [0.81858e+03,0.81917e+03,0.81978e+03,0.82050e+03,0.82130e+03]
  kb[:,49-12, 8] = [0.81836e+03,0.81894e+03,0.81953e+03,0.82019e+03,0.82097e+03]
  kb[:,50-12, 8] = [0.81817e+03,0.81871e+03,0.81929e+03,0.81993e+03,0.82067e+03]
  kb[:,51-12, 8] = [0.81790e+03,0.81850e+03,0.81910e+03,0.81967e+03,0.82039e+03]
  kb[:,52-12, 8] = [0.81779e+03,0.81830e+03,0.81888e+03,0.81945e+03,0.82012e+03]
  kb[:,53-12, 8] = [0.81766e+03,0.81812e+03,0.81866e+03,0.81924e+03,0.81986e+03]
  kb[:,54-12, 8] = [0.81744e+03,0.81785e+03,0.81846e+03,0.81906e+03,0.81964e+03]
  kb[:,55-12, 8] = [0.81733e+03,0.81777e+03,0.81829e+03,0.81886e+03,0.81944e+03]
  kb[:,56-12, 8] = [0.81741e+03,0.81767e+03,0.81813e+03,0.81867e+03,0.81925e+03]
  kb[:,57-12, 8] = [0.81785e+03,0.81749e+03,0.81788e+03,0.81848e+03,0.81908e+03]
  kb[:,58-12, 8] = [0.81817e+03,0.81731e+03,0.81779e+03,0.81832e+03,0.81889e+03]
  kb[:,59-12, 8] = [0.81824e+03,0.81736e+03,0.81776e+03,0.81826e+03,0.81883e+03]
  kb[:,13-12, 9] = [0.94133e+03,0.94438e+03,0.94687e+03,0.94883e+03,0.95021e+03]
  kb[:,14-12, 9] = [0.94163e+03,0.94463e+03,0.94708e+03,0.94898e+03,0.95031e+03]
  kb[:,15-12, 9] = [0.94195e+03,0.94490e+03,0.94731e+03,0.94914e+03,0.95042e+03]
  kb[:,16-12, 9] = [0.94229e+03,0.94518e+03,0.94753e+03,0.94930e+03,0.95052e+03]
  kb[:,17-12, 9] = [0.94262e+03,0.94544e+03,0.94774e+03,0.94946e+03,0.95061e+03]
  kb[:,18-12, 9] = [0.94293e+03,0.94570e+03,0.94795e+03,0.94960e+03,0.95069e+03]
  kb[:,19-12, 9] = [0.94323e+03,0.94595e+03,0.94814e+03,0.94974e+03,0.95078e+03]
  kb[:,20-12, 9] = [0.94357e+03,0.94624e+03,0.94835e+03,0.94989e+03,0.95087e+03]
  kb[:,21-12, 9] = [0.94393e+03,0.94654e+03,0.94857e+03,0.95003e+03,0.95096e+03]
  kb[:,22-12, 9] = [0.94445e+03,0.94695e+03,0.94888e+03,0.95025e+03,0.95105e+03]
  kb[:,23-12, 9] = [0.94499e+03,0.94737e+03,0.94919e+03,0.95045e+03,0.95113e+03]
  kb[:,24-12, 9] = [0.94549e+03,0.94778e+03,0.94948e+03,0.95063e+03,0.95120e+03]
  kb[:,25-12, 9] = [0.94602e+03,0.94819e+03,0.94978e+03,0.95082e+03,0.95125e+03]
  kb[:,26-12, 9] = [0.94658e+03,0.94861e+03,0.95006e+03,0.95096e+03,0.95127e+03]
  kb[:,27-12, 9] = [0.94710e+03,0.94900e+03,0.95033e+03,0.95109e+03,0.95126e+03]
  kb[:,28-12, 9] = [0.94759e+03,0.94937e+03,0.95056e+03,0.95118e+03,0.95123e+03]
  kb[:,29-12, 9] = [0.94811e+03,0.94973e+03,0.95077e+03,0.95124e+03,0.95116e+03]
  kb[:,30-12, 9] = [0.94856e+03,0.95003e+03,0.95095e+03,0.95127e+03,0.95105e+03]
  kb[:,31-12, 9] = [0.94901e+03,0.95032e+03,0.95108e+03,0.95127e+03,0.95090e+03]
  kb[:,32-12, 9] = [0.94939e+03,0.95057e+03,0.95118e+03,0.95121e+03,0.95071e+03]
  kb[:,33-12, 9] = [0.94975e+03,0.95079e+03,0.95125e+03,0.95115e+03,0.95050e+03]
  kb[:,34-12, 9] = [0.95004e+03,0.95095e+03,0.95128e+03,0.95105e+03,0.95030e+03]
  kb[:,35-12, 9] = [0.95020e+03,0.95102e+03,0.95127e+03,0.95098e+03,0.95016e+03]
  kb[:,36-12, 9] = [0.95021e+03,0.95103e+03,0.95127e+03,0.95096e+03,0.95023e+03]
  kb[:,37-12, 9] = [0.95006e+03,0.95096e+03,0.95127e+03,0.95104e+03,0.95027e+03]
  kb[:,38-12, 9] = [0.94990e+03,0.95087e+03,0.95126e+03,0.95110e+03,0.95040e+03]
  kb[:,39-12, 9] = [0.94972e+03,0.95077e+03,0.95124e+03,0.95117e+03,0.95053e+03]
  kb[:,40-12, 9] = [0.94939e+03,0.95057e+03,0.95118e+03,0.95123e+03,0.95072e+03]
  kb[:,41-12, 9] = [0.94899e+03,0.95032e+03,0.95108e+03,0.95127e+03,0.95090e+03]
  kb[:,42-12, 9] = [0.94856e+03,0.95004e+03,0.95094e+03,0.95127e+03,0.95105e+03]
  kb[:,43-12, 9] = [0.94799e+03,0.94963e+03,0.95073e+03,0.95123e+03,0.95118e+03]
  kb[:,44-12, 9] = [0.94729e+03,0.94914e+03,0.95041e+03,0.95112e+03,0.95126e+03]
  kb[:,45-12, 9] = [0.94655e+03,0.94858e+03,0.95004e+03,0.95095e+03,0.95127e+03]
  kb[:,46-12, 9] = [0.94568e+03,0.94792e+03,0.94958e+03,0.95070e+03,0.95122e+03]
  kb[:,47-12, 9] = [0.94463e+03,0.94708e+03,0.94898e+03,0.95031e+03,0.95108e+03]
  kb[:,48-12, 9] = [0.94348e+03,0.94616e+03,0.94829e+03,0.94984e+03,0.95084e+03]
  kb[:,49-12, 9] = [0.94226e+03,0.94515e+03,0.94748e+03,0.94928e+03,0.95051e+03]
  kb[:,50-12, 9] = [0.94103e+03,0.94413e+03,0.94669e+03,0.94868e+03,0.95011e+03]
  kb[:,51-12, 9] = [0.93978e+03,0.94307e+03,0.94580e+03,0.94802e+03,0.94965e+03]
  kb[:,52-12, 9] = [0.93828e+03,0.94192e+03,0.94488e+03,0.94728e+03,0.94912e+03]
  kb[:,53-12, 9] = [0.93678e+03,0.94072e+03,0.94388e+03,0.94646e+03,0.94851e+03]
  kb[:,54-12, 9] = [0.93537e+03,0.93955e+03,0.94287e+03,0.94565e+03,0.94789e+03]
  kb[:,55-12, 9] = [0.93380e+03,0.93820e+03,0.94186e+03,0.94483e+03,0.94722e+03]
  kb[:,56-12, 9] = [0.93205e+03,0.93688e+03,0.94080e+03,0.94394e+03,0.94652e+03]
  kb[:,57-12, 9] = [0.92997e+03,0.93554e+03,0.93969e+03,0.94299e+03,0.94574e+03]
  kb[:,58-12, 9] = [0.92808e+03,0.93414e+03,0.93844e+03,0.94204e+03,0.94498e+03]
  kb[:,59-12, 9] = [0.92793e+03,0.93341e+03,0.93792e+03,0.94165e+03,0.94465e+03]
  kb[:,13-12,10] = [0.98036e+03,0.98526e+03,0.98973e+03,0.99279e+03,0.99468e+03]
  kb[:,14-12,10] = [0.98067e+03,0.98568e+03,0.99005e+03,0.99301e+03,0.99477e+03]
  kb[:,15-12,10] = [0.98125e+03,0.98614e+03,0.99040e+03,0.99324e+03,0.99485e+03]
  kb[:,16-12,10] = [0.98185e+03,0.98662e+03,0.99076e+03,0.99348e+03,0.99493e+03]
  kb[:,17-12,10] = [0.98242e+03,0.98721e+03,0.99109e+03,0.99369e+03,0.99498e+03]
  kb[:,18-12,10] = [0.98296e+03,0.98755e+03,0.99141e+03,0.99390e+03,0.99502e+03]
  kb[:,19-12,10] = [0.98350e+03,0.98800e+03,0.99172e+03,0.99408e+03,0.99505e+03]
  kb[:,20-12,10] = [0.98411e+03,0.98852e+03,0.99205e+03,0.99428e+03,0.99507e+03]
  kb[:,21-12,10] = [0.98464e+03,0.98894e+03,0.99238e+03,0.99447e+03,0.99496e+03]
  kb[:,22-12,10] = [0.98547e+03,0.98983e+03,0.99286e+03,0.99470e+03,0.99509e+03]
  kb[:,23-12,10] = [0.98617e+03,0.99062e+03,0.99331e+03,0.99487e+03,0.99508e+03]
  kb[:,24-12,10] = [0.98718e+03,0.99115e+03,0.99373e+03,0.99499e+03,0.99503e+03]
  kb[:,25-12,10] = [0.98813e+03,0.99180e+03,0.99413e+03,0.99494e+03,0.99492e+03]
  kb[:,26-12,10] = [0.98918e+03,0.99244e+03,0.99450e+03,0.99508e+03,0.99475e+03]
  kb[:,27-12,10] = [0.99010e+03,0.99304e+03,0.99478e+03,0.99509e+03,0.99499e+03]
  kb[:,28-12,10] = [0.99090e+03,0.99357e+03,0.99495e+03,0.99505e+03,0.99433e+03]
  kb[:,29-12,10] = [0.99167e+03,0.99393e+03,0.99505e+03,0.99495e+03,0.99385e+03]
  kb[:,30-12,10] = [0.99249e+03,0.99447e+03,0.99508e+03,0.99477e+03,0.99348e+03]
  kb[:,31-12,10] = [0.99290e+03,0.99477e+03,0.99509e+03,0.99452e+03,0.99326e+03]
  kb[:,32-12,10] = [0.99360e+03,0.99496e+03,0.99505e+03,0.99419e+03,0.99287e+03]
  kb[:,33-12,10] = [0.99410e+03,0.99506e+03,0.99493e+03,0.99381e+03,0.99210e+03]
  kb[:,34-12,10] = [0.99447e+03,0.99508e+03,0.99465e+03,0.99347e+03,0.99121e+03]
  kb[:,35-12,10] = [0.99453e+03,0.99509e+03,0.99465e+03,0.99330e+03,0.99083e+03]
  kb[:,36-12,10] = [0.99467e+03,0.99509e+03,0.99463e+03,0.99327e+03,0.98992e+03]
  kb[:,37-12,10] = [0.99451e+03,0.99508e+03,0.99475e+03,0.99344e+03,0.99124e+03]
  kb[:,38-12,10] = [0.99429e+03,0.99507e+03,0.99486e+03,0.99365e+03,0.99174e+03]
  kb[:,39-12,10] = [0.99405e+03,0.99505e+03,0.99495e+03,0.99373e+03,0.99219e+03]
  kb[:,40-12,10] = [0.99347e+03,0.99496e+03,0.99505e+03,0.99420e+03,0.99289e+03]
  kb[:,41-12,10] = [0.99302e+03,0.99477e+03,0.99509e+03,0.99452e+03,0.99315e+03]
  kb[:,42-12,10] = [0.99237e+03,0.99435e+03,0.99508e+03,0.99477e+03,0.99348e+03]
  kb[:,43-12,10] = [0.99148e+03,0.99394e+03,0.99491e+03,0.99498e+03,0.99395e+03]
  kb[:,44-12,10] = [0.99040e+03,0.99324e+03,0.99485e+03,0.99508e+03,0.99453e+03]
  kb[:,45-12,10] = [0.98910e+03,0.99241e+03,0.99449e+03,0.99508e+03,0.99476e+03]
  kb[:,46-12,10] = [0.98735e+03,0.99138e+03,0.99388e+03,0.99490e+03,0.99500e+03]
  kb[:,47-12,10] = [0.98565e+03,0.99006e+03,0.99301e+03,0.99477e+03,0.99509e+03]
  kb[:,48-12,10] = [0.98383e+03,0.98834e+03,0.99196e+03,0.99423e+03,0.99507e+03]
  kb[:,49-12,10] = [0.98174e+03,0.98651e+03,0.99073e+03,0.99346e+03,0.99492e+03]
  kb[:,50-12,10] = [0.97955e+03,0.98480e+03,0.98924e+03,0.99257e+03,0.99470e+03]
  kb[:,51-12,10] = [0.97755e+03,0.98298e+03,0.98767e+03,0.99154e+03,0.99399e+03]
  kb[:,52-12,10] = [0.97554e+03,0.98111e+03,0.98597e+03,0.99036e+03,0.99323e+03]
  kb[:,53-12,10] = [0.97320e+03,0.97892e+03,0.98413e+03,0.98888e+03,0.99232e+03]
  kb[:,54-12,10] = [0.97086e+03,0.97720e+03,0.98263e+03,0.98733e+03,0.99135e+03]
  kb[:,55-12,10] = [0.96881e+03,0.97539e+03,0.98096e+03,0.98583e+03,0.99029e+03]
  kb[:,56-12,10] = [0.96735e+03,0.97336e+03,0.97906e+03,0.98435e+03,0.98900e+03]
  kb[:,57-12,10] = [0.96407e+03,0.97115e+03,0.97740e+03,0.98282e+03,0.98752e+03]
  kb[:,58-12,10] = [0.96154e+03,0.96893e+03,0.97574e+03,0.98129e+03,0.98611e+03]
  kb[:,59-12,10] = [0.95571e+03,0.96861e+03,0.97503e+03,0.98063e+03,0.98557e+03]
  kb[:,13-12,11] = [0.98771e+03,0.99296e+03,0.99655e+03,0.99859e+03,0.99936e+03]
  kb[:,14-12,11] = [0.98850e+03,0.99335e+03,0.99679e+03,0.99871e+03,0.99939e+03]
  kb[:,15-12,11] = [0.98905e+03,0.99377e+03,0.99705e+03,0.99882e+03,0.99942e+03]
  kb[:,16-12,11] = [0.98962e+03,0.99419e+03,0.99730e+03,0.99893e+03,0.99945e+03]
  kb[:,17-12,11] = [0.99016e+03,0.99444e+03,0.99754e+03,0.99903e+03,0.99947e+03]
  kb[:,18-12,11] = [0.99067e+03,0.99497e+03,0.99775e+03,0.99911e+03,0.99948e+03]
  kb[:,19-12,11] = [0.99117e+03,0.99532e+03,0.99795e+03,0.99918e+03,0.99949e+03]
  kb[:,20-12,11] = [0.99173e+03,0.99572e+03,0.99816e+03,0.99924e+03,0.99949e+03]
  kb[:,21-12,11] = [0.99228e+03,0.99610e+03,0.99835e+03,0.99929e+03,0.99946e+03]
  kb[:,22-12,11] = [0.99309e+03,0.99663e+03,0.99862e+03,0.99936e+03,0.99953e+03]
  kb[:,23-12,11] = [0.99388e+03,0.99698e+03,0.99900e+03,0.99942e+03,0.99925e+03]
  kb[:,24-12,11] = [0.99466e+03,0.99758e+03,0.99904e+03,0.99947e+03,0.99905e+03]
  kb[:,25-12,11] = [0.99542e+03,0.99800e+03,0.99919e+03,0.99963e+03,0.99878e+03]
  kb[:,26-12,11] = [0.99618e+03,0.99839e+03,0.99930e+03,0.99945e+03,0.99847e+03]
  kb[:,27-12,11] = [0.99683e+03,0.99872e+03,0.99939e+03,0.99934e+03,0.99751e+03]
  kb[:,28-12,11] = [0.99741e+03,0.99897e+03,0.99946e+03,0.99913e+03,0.99749e+03]
  kb[:,29-12,11] = [0.99792e+03,0.99917e+03,0.99949e+03,0.99884e+03,0.99707e+03]
  kb[:,30-12,11] = [0.99821e+03,0.99929e+03,0.99946e+03,0.99850e+03,0.99637e+03]
  kb[:,31-12,11] = [0.99871e+03,0.99938e+03,0.99949e+03,0.99810e+03,0.99532e+03]
  kb[:,32-12,11] = [0.99899e+03,0.99961e+03,0.99926e+03,0.99761e+03,0.99460e+03]
  kb[:,33-12,11] = [0.99918e+03,0.99949e+03,0.99880e+03,0.99700e+03,0.99413e+03]
  kb[:,34-12,11] = [0.99929e+03,0.99946e+03,0.99850e+03,0.99635e+03,0.99329e+03]
  kb[:,35-12,11] = [0.99934e+03,0.99941e+03,0.99830e+03,0.99590e+03,0.99292e+03]
  kb[:,36-12,11] = [0.99935e+03,0.99940e+03,0.99826e+03,0.99583e+03,0.99287e+03]
  kb[:,37-12,11] = [0.99930e+03,0.99946e+03,0.99846e+03,0.99628e+03,0.99327e+03]
  kb[:,38-12,11] = [0.99924e+03,0.99949e+03,0.99866e+03,0.99670e+03,0.99350e+03]
  kb[:,39-12,11] = [0.99917e+03,0.99949e+03,0.99899e+03,0.99707e+03,0.99451e+03]
  kb[:,40-12,11] = [0.99898e+03,0.99946e+03,0.99913e+03,0.99777e+03,0.99505e+03]
  kb[:,41-12,11] = [0.99871e+03,0.99939e+03,0.99935e+03,0.99810e+03,0.99546e+03]
  kb[:,42-12,11] = [0.99835e+03,0.99929e+03,0.99947e+03,0.99851e+03,0.99637e+03]
  kb[:,43-12,11] = [0.99780e+03,0.99913e+03,0.99949e+03,0.99893e+03,0.99722e+03]
  kb[:,44-12,11] = [0.99704e+03,0.99882e+03,0.99942e+03,0.99928e+03,0.99778e+03]
  kb[:,45-12,11] = [0.99613e+03,0.99837e+03,0.99930e+03,0.99947e+03,0.99849e+03]
  kb[:,46-12,11] = [0.99492e+03,0.99787e+03,0.99911e+03,0.99949e+03,0.99897e+03]
  kb[:,47-12,11] = [0.99335e+03,0.99679e+03,0.99871e+03,0.99939e+03,0.99936e+03]
  kb[:,48-12,11] = [0.99156e+03,0.99560e+03,0.99810e+03,0.99923e+03,0.99951e+03]
  kb[:,49-12,11] = [0.98955e+03,0.99414e+03,0.99727e+03,0.99894e+03,0.99946e+03]
  kb[:,50-12,11] = [0.98749e+03,0.99256e+03,0.99630e+03,0.99847e+03,0.99918e+03]
  kb[:,51-12,11] = [0.98496e+03,0.99085e+03,0.99509e+03,0.99784e+03,0.99916e+03]
  kb[:,52-12,11] = [0.98304e+03,0.98897e+03,0.99370e+03,0.99701e+03,0.99883e+03]
  kb[:,53-12,11] = [0.98104e+03,0.98693e+03,0.99213e+03,0.99600e+03,0.99833e+03]
  kb[:,54-12,11] = [0.97906e+03,0.98440e+03,0.99053e+03,0.99485e+03,0.99771e+03]
  kb[:,55-12,11] = [0.97789e+03,0.98286e+03,0.98886e+03,0.99361e+03,0.99696e+03]
  kb[:,56-12,11] = [0.97514e+03,0.98117e+03,0.98708e+03,0.99224e+03,0.99622e+03]
  kb[:,57-12,11] = [0.97242e+03,0.97945e+03,0.98472e+03,0.99072e+03,0.99500e+03]
  kb[:,58-12,11] = [0.96969e+03,0.97814e+03,0.98316e+03,0.98916e+03,0.99384e+03]
  kb[:,59-12,11] = [0.96859e+03,0.97725e+03,0.98261e+03,0.98851e+03,0.99335e+03]
  kb[:,13-12,12] = [0.99768e+03,0.99984e+03,0.10010e+04,0.10024e+04,0.10029e+04]
  kb[:,14-12,12] = [0.99795e+03,0.99992e+03,0.10011e+04,0.10025e+04,0.10029e+04]
  kb[:,15-12,12] = [0.99880e+03,0.99996e+03,0.10013e+04,0.10026e+04,0.10031e+04]
  kb[:,16-12,12] = [0.99852e+03,0.10003e+04,0.10015e+04,0.10028e+04,0.10029e+04]
  kb[:,17-12,12] = [0.99876e+03,0.10000e+04,0.10016e+04,0.10027e+04,0.10029e+04]
  kb[:,18-12,12] = [0.99898e+03,0.10001e+04,0.10018e+04,0.10028e+04,0.10028e+04]
  kb[:,19-12,12] = [0.99918e+03,0.10003e+04,0.10021e+04,0.10028e+04,0.10028e+04]
  kb[:,20-12,12] = [0.99940e+03,0.10005e+04,0.10026e+04,0.10029e+04,0.10031e+04]
  kb[:,21-12,12] = [0.99959e+03,0.10007e+04,0.10022e+04,0.10029e+04,0.10026e+04]
  kb[:,22-12,12] = [0.99984e+03,0.10010e+04,0.10024e+04,0.10030e+04,0.10023e+04]
  kb[:,23-12,12] = [0.99996e+03,0.10013e+04,0.10024e+04,0.10029e+04,0.10025e+04]
  kb[:,24-12,12] = [0.10000e+04,0.10016e+04,0.10029e+04,0.10029e+04,0.10023e+04]
  kb[:,25-12,12] = [0.10003e+04,0.10023e+04,0.10029e+04,0.10026e+04,0.10019e+04]
  kb[:,26-12,12] = [0.10007e+04,0.10022e+04,0.10029e+04,0.10026e+04,0.10015e+04]
  kb[:,27-12,12] = [0.10012e+04,0.10025e+04,0.10031e+04,0.10024e+04,0.10010e+04]
  kb[:,28-12,12] = [0.10015e+04,0.10027e+04,0.10029e+04,0.10022e+04,0.10005e+04]
  kb[:,29-12,12] = [0.10019e+04,0.10028e+04,0.10028e+04,0.10019e+04,0.99992e+03]
  kb[:,30-12,12] = [0.10022e+04,0.10029e+04,0.10026e+04,0.10017e+04,0.99965e+03]
  kb[:,31-12,12] = [0.10025e+04,0.10030e+04,0.10022e+04,0.10011e+04,0.99843e+03]
  kb[:,32-12,12] = [0.10027e+04,0.10029e+04,0.10020e+04,0.10007e+04,0.99748e+03]
  kb[:,33-12,12] = [0.10029e+04,0.10028e+04,0.10021e+04,0.10000e+04,0.99658e+03]
  kb[:,34-12,12] = [0.10031e+04,0.10026e+04,0.10017e+04,0.99945e+03,0.99563e+03]
  kb[:,35-12,12] = [0.10030e+04,0.10025e+04,0.10031e+04,0.99888e+03,0.99496e+03]
  kb[:,36-12,12] = [0.10030e+04,0.10025e+04,0.10020e+04,0.99880e+03,0.99482e+03]
  kb[:,37-12,12] = [0.10029e+04,0.10026e+04,0.10015e+04,0.99939e+03,0.99547e+03]
  kb[:,38-12,12] = [0.10029e+04,0.10027e+04,0.10017e+04,0.99958e+03,0.99634e+03]
  kb[:,39-12,12] = [0.10028e+04,0.10028e+04,0.10017e+04,0.99993e+03,0.99613e+03]
  kb[:,40-12,12] = [0.10027e+04,0.10029e+04,0.10022e+04,0.10003e+04,0.99696e+03]
  kb[:,41-12,12] = [0.10027e+04,0.10029e+04,0.10024e+04,0.10013e+04,0.99857e+03]
  kb[:,42-12,12] = [0.10022e+04,0.10029e+04,0.10028e+04,0.10021e+04,0.99929e+03]
  kb[:,43-12,12] = [0.10018e+04,0.10028e+04,0.10028e+04,0.10020e+04,0.10001e+04]
  kb[:,44-12,12] = [0.10013e+04,0.10026e+04,0.10029e+04,0.10023e+04,0.10016e+04]
  kb[:,45-12,12] = [0.10007e+04,0.10022e+04,0.10029e+04,0.10026e+04,0.10015e+04]
  kb[:,46-12,12] = [0.10014e+04,0.10016e+04,0.10028e+04,0.10028e+04,0.10020e+04]
  kb[:,47-12,12] = [0.99993e+03,0.10012e+04,0.10025e+04,0.10029e+04,0.10024e+04]
  kb[:,48-12,12] = [0.99937e+03,0.10005e+04,0.10022e+04,0.10029e+04,0.10027e+04]
  kb[:,49-12,12] = [0.99853e+03,0.10002e+04,0.10015e+04,0.10026e+04,0.10029e+04]
  kb[:,50-12,12] = [0.99741e+03,0.99993e+03,0.10008e+04,0.10023e+04,0.10029e+04]
  kb[:,51-12,12] = [0.99665e+03,0.99912e+03,0.10002e+04,0.10018e+04,0.10028e+04]
  kb[:,52-12,12] = [0.99385e+03,0.99826e+03,0.10000e+04,0.10013e+04,0.10025e+04]
  kb[:,53-12,12] = [0.99099e+03,0.99714e+03,0.99965e+03,0.10007e+04,0.10022e+04]
  kb[:,54-12,12] = [0.98826e+03,0.99612e+03,0.99903e+03,0.10002e+04,0.10018e+04]
  kb[:,55-12,12] = [0.98397e+03,0.99379e+03,0.99822e+03,0.10000e+04,0.10013e+04]
  kb[:,56-12,12] = [0.98006e+03,0.99120e+03,0.99722e+03,0.99969e+03,0.10006e+04]
  kb[:,57-12,12] = [0.97631e+03,0.98840e+03,0.99624e+03,0.99911e+03,0.10002e+04]
  kb[:,58-12,12] = [0.97279e+03,0.98508e+03,0.99422e+03,0.99838e+03,0.10001e+04]
  kb[:,59-12,12] = [0.97171e+03,0.98327e+03,0.99315e+03,0.99818e+03,0.99997e+03]
  kb[:,13-12,13] = [0.99623e+03,0.99978e+03,0.10034e+04,0.10050e+04,0.10054e+04]
  kb[:,14-12,13] = [0.99674e+03,0.10002e+04,0.10034e+04,0.10050e+04,0.10053e+04]
  kb[:,15-12,13] = [0.99652e+03,0.10007e+04,0.10036e+04,0.10051e+04,0.10051e+04]
  kb[:,16-12,13] = [0.99785e+03,0.10007e+04,0.10039e+04,0.10049e+04,0.10053e+04]
  kb[:,17-12,13] = [0.99816e+03,0.10017e+04,0.10043e+04,0.10052e+04,0.10052e+04]
  kb[:,18-12,13] = [0.99792e+03,0.10021e+04,0.10042e+04,0.10053e+04,0.10051e+04]
  kb[:,19-12,13] = [0.99826e+03,0.10023e+04,0.10042e+04,0.10053e+04,0.10051e+04]
  kb[:,20-12,13] = [0.99869e+03,0.10028e+04,0.10038e+04,0.10053e+04,0.10044e+04]
  kb[:,21-12,13] = [0.99914e+03,0.10028e+04,0.10048e+04,0.10054e+04,0.10048e+04]
  kb[:,22-12,13] = [0.99989e+03,0.10035e+04,0.10050e+04,0.10053e+04,0.10046e+04]
  kb[:,23-12,13] = [0.10011e+04,0.10037e+04,0.10051e+04,0.10053e+04,0.10041e+04]
  kb[:,24-12,13] = [0.10018e+04,0.10041e+04,0.10050e+04,0.10052e+04,0.10037e+04]
  kb[:,25-12,13] = [0.10026e+04,0.10040e+04,0.10056e+04,0.10050e+04,0.10036e+04]
  kb[:,26-12,13] = [0.10031e+04,0.10048e+04,0.10054e+04,0.10050e+04,0.10030e+04]
  kb[:,27-12,13] = [0.10047e+04,0.10050e+04,0.10051e+04,0.10045e+04,0.10025e+04]
  kb[:,28-12,13] = [0.10039e+04,0.10052e+04,0.10052e+04,0.10041e+04,0.10018e+04]
  kb[:,29-12,13] = [0.10044e+04,0.10053e+04,0.10051e+04,0.10037e+04,0.10011e+04]
  kb[:,30-12,13] = [0.10048e+04,0.10054e+04,0.10048e+04,0.10029e+04,0.99977e+03]
  kb[:,31-12,13] = [0.10050e+04,0.10053e+04,0.10045e+04,0.10025e+04,0.99968e+03]
  kb[:,32-12,13] = [0.10052e+04,0.10050e+04,0.10041e+04,0.10015e+04,0.99880e+03]
  kb[:,33-12,13] = [0.10053e+04,0.10050e+04,0.10033e+04,0.10008e+04,0.99675e+03]
  kb[:,34-12,13] = [0.10051e+04,0.10048e+04,0.10028e+04,0.10000e+04,0.99529e+03]
  kb[:,35-12,13] = [0.10054e+04,0.10047e+04,0.10003e+04,0.99992e+03,0.99481e+03]
  kb[:,36-12,13] = [0.10054e+04,0.10046e+04,0.10017e+04,0.99987e+03,0.99445e+03]
  kb[:,37-12,13] = [0.10054e+04,0.10048e+04,0.10030e+04,0.99995e+03,0.99517e+03]
  kb[:,38-12,13] = [0.10054e+04,0.10049e+04,0.10034e+04,0.10006e+04,0.99603e+03]
  kb[:,39-12,13] = [0.10053e+04,0.10051e+04,0.10037e+04,0.10011e+04,0.99693e+03]
  kb[:,40-12,13] = [0.10052e+04,0.10052e+04,0.10041e+04,0.10018e+04,0.99860e+03]
  kb[:,41-12,13] = [0.10048e+04,0.10053e+04,0.10045e+04,0.10023e+04,0.99951e+03]
  kb[:,42-12,13] = [0.10048e+04,0.10054e+04,0.10046e+04,0.10024e+04,0.10003e+04]
  kb[:,43-12,13] = [0.10043e+04,0.10053e+04,0.10051e+04,0.10038e+04,0.10013e+04]
  kb[:,44-12,13] = [0.10036e+04,0.10051e+04,0.10053e+04,0.10044e+04,0.10012e+04]
  kb[:,45-12,13] = [0.10028e+04,0.10048e+04,0.10054e+04,0.10048e+04,0.10031e+04]
  kb[:,46-12,13] = [0.10008e+04,0.10042e+04,0.10053e+04,0.10052e+04,0.10038e+04]
  kb[:,47-12,13] = [0.10002e+04,0.10039e+04,0.10050e+04,0.10053e+04,0.10045e+04]
  kb[:,48-12,13] = [0.99859e+03,0.10025e+04,0.10043e+04,0.10054e+04,0.10050e+04]
  kb[:,49-12,13] = [0.99776e+03,0.10009e+04,0.10041e+04,0.10052e+04,0.10053e+04]
  kb[:,50-12,13] = [0.99565e+03,0.99944e+03,0.10034e+04,0.10049e+04,0.10054e+04]
  kb[:,51-12,13] = [0.99271e+03,0.99811e+03,0.10022e+04,0.10043e+04,0.10053e+04]
  kb[:,52-12,13] = [0.99056e+03,0.99717e+03,0.10006e+04,0.10036e+04,0.10051e+04]
  kb[:,53-12,13] = [0.98830e+03,0.99502e+03,0.99909e+03,0.10028e+04,0.10048e+04]
  kb[:,54-12,13] = [0.98610e+03,0.99275e+03,0.99821e+03,0.10021e+04,0.10042e+04]
  kb[:,55-12,13] = [0.98386e+03,0.99045e+03,0.99729e+03,0.10008e+04,0.10038e+04]
  kb[:,56-12,13] = [0.98150e+03,0.98846e+03,0.99517e+03,0.99919e+03,0.10028e+04]
  kb[:,57-12,13] = [0.97902e+03,0.98636e+03,0.99303e+03,0.99832e+03,0.10022e+04]
  kb[:,58-12,13] = [0.97656e+03,0.98426e+03,0.99079e+03,0.99735e+03,0.10008e+04]
  kb[:,59-12,13] = [0.97503e+03,0.98337e+03,0.99004e+03,0.99647e+03,0.10002e+04]
  kb[:,13-12,14] = [0.99596e+03,0.10023e+04,0.10047e+04,0.10067e+04,0.10072e+04]
  kb[:,14-12,14] = [0.99625e+03,0.10024e+04,0.10052e+04,0.10068e+04,0.10071e+04]
  kb[:,15-12,14] = [0.99662e+03,0.10025e+04,0.10054e+04,0.10069e+04,0.10071e+04]
  kb[:,16-12,14] = [0.99703e+03,0.10026e+04,0.10056e+04,0.10070e+04,0.10070e+04]
  kb[:,17-12,14] = [0.99784e+03,0.10032e+04,0.10054e+04,0.10071e+04,0.10070e+04]
  kb[:,18-12,14] = [0.99947e+03,0.10032e+04,0.10060e+04,0.10071e+04,0.10069e+04]
  kb[:,19-12,14] = [0.10001e+04,0.10036e+04,0.10062e+04,0.10071e+04,0.10068e+04]
  kb[:,20-12,14] = [0.10009e+04,0.10037e+04,0.10064e+04,0.10076e+04,0.10067e+04]
  kb[:,21-12,14] = [0.10015e+04,0.10046e+04,0.10065e+04,0.10072e+04,0.10065e+04]
  kb[:,22-12,14] = [0.10023e+04,0.10047e+04,0.10068e+04,0.10071e+04,0.10063e+04]
  kb[:,23-12,14] = [0.10021e+04,0.10055e+04,0.10069e+04,0.10071e+04,0.10060e+04]
  kb[:,24-12,14] = [0.10028e+04,0.10059e+04,0.10071e+04,0.10069e+04,0.10057e+04]
  kb[:,25-12,14] = [0.10033e+04,0.10062e+04,0.10068e+04,0.10068e+04,0.10052e+04]
  kb[:,26-12,14] = [0.10043e+04,0.10066e+04,0.10072e+04,0.10061e+04,0.10047e+04]
  kb[:,27-12,14] = [0.10033e+04,0.10068e+04,0.10071e+04,0.10062e+04,0.10042e+04]
  kb[:,28-12,14] = [0.10057e+04,0.10070e+04,0.10070e+04,0.10058e+04,0.10036e+04]
  kb[:,29-12,14] = [0.10062e+04,0.10075e+04,0.10068e+04,0.10053e+04,0.10030e+04]
  kb[:,30-12,14] = [0.10065e+04,0.10072e+04,0.10069e+04,0.10048e+04,0.10022e+04]
  kb[:,31-12,14] = [0.10068e+04,0.10071e+04,0.10062e+04,0.10042e+04,0.10014e+04]
  kb[:,32-12,14] = [0.10070e+04,0.10070e+04,0.10058e+04,0.10035e+04,0.99925e+03]
  kb[:,33-12,14] = [0.10071e+04,0.10068e+04,0.10053e+04,0.10029e+04,0.99808e+03]
  kb[:,34-12,14] = [0.10072e+04,0.10065e+04,0.10048e+04,0.10021e+04,0.99788e+03]
  kb[:,35-12,14] = [0.10072e+04,0.10064e+04,0.10045e+04,0.10016e+04,0.99700e+03]
  kb[:,36-12,14] = [0.10072e+04,0.10063e+04,0.10044e+04,0.10015e+04,0.99731e+03]
  kb[:,37-12,14] = [0.10072e+04,0.10065e+04,0.10047e+04,0.10021e+04,0.99778e+03]
  kb[:,38-12,14] = [0.10072e+04,0.10067e+04,0.10050e+04,0.10026e+04,0.99820e+03]
  kb[:,39-12,14] = [0.10071e+04,0.10068e+04,0.10053e+04,0.10030e+04,0.99855e+03]
  kb[:,40-12,14] = [0.10070e+04,0.10070e+04,0.10058e+04,0.10036e+04,0.99965e+03]
  kb[:,41-12,14] = [0.10068e+04,0.10071e+04,0.10062e+04,0.10042e+04,0.10010e+04]
  kb[:,42-12,14] = [0.10065e+04,0.10072e+04,0.10065e+04,0.10048e+04,0.10022e+04]
  kb[:,43-12,14] = [0.10061e+04,0.10071e+04,0.10069e+04,0.10055e+04,0.10031e+04]
  kb[:,44-12,14] = [0.10054e+04,0.10069e+04,0.10071e+04,0.10061e+04,0.10039e+04]
  kb[:,45-12,14] = [0.10046e+04,0.10066e+04,0.10072e+04,0.10065e+04,0.10048e+04]
  kb[:,46-12,14] = [0.10023e+04,0.10060e+04,0.10075e+04,0.10069e+04,0.10055e+04]
  kb[:,47-12,14] = [0.10024e+04,0.10044e+04,0.10068e+04,0.10071e+04,0.10062e+04]
  kb[:,48-12,14] = [0.10006e+04,0.10039e+04,0.10063e+04,0.10072e+04,0.10067e+04]
  kb[:,49-12,14] = [0.99705e+03,0.10026e+04,0.10052e+04,0.10070e+04,0.10070e+04]
  kb[:,50-12,14] = [0.99574e+03,0.10014e+04,0.10040e+04,0.10066e+04,0.10072e+04]
  kb[:,51-12,14] = [0.99447e+03,0.99966e+03,0.10033e+04,0.10065e+04,0.10071e+04]
  kb[:,52-12,14] = [0.99299e+03,0.99670e+03,0.10025e+04,0.10054e+04,0.10069e+04]
  kb[:,53-12,14] = [0.99138e+03,0.99552e+03,0.10013e+04,0.10044e+04,0.10065e+04]
  kb[:,54-12,14] = [0.98979e+03,0.99462e+03,0.99874e+03,0.10030e+04,0.10060e+04]
  kb[:,55-12,14] = [0.98814e+03,0.99291e+03,0.99630e+03,0.10029e+04,0.10050e+04]
  kb[:,56-12,14] = [0.98639e+03,0.99150e+03,0.99562e+03,0.10014e+04,0.10045e+04]
  kb[:,57-12,14] = [0.98454e+03,0.98998e+03,0.99438e+03,0.99904e+03,0.10032e+04]
  kb[:,58-12,14] = [0.98268e+03,0.98844e+03,0.99315e+03,0.99688e+03,0.10025e+04]
  kb[:,59-12,14] = [0.98190e+03,0.98779e+03,0.99262e+03,0.99638e+03,0.10024e+04]
  kb[:,13-12,15] = [0.99739e+03,0.10017e+04,0.10064e+04,0.10076e+04,0.10080e+04]
  kb[:,14-12,15] = [0.99781e+03,0.10021e+04,0.10065e+04,0.10077e+04,0.10090e+04]
  kb[:,15-12,15] = [0.99827e+03,0.10029e+04,0.10067e+04,0.10077e+04,0.10079e+04]
  kb[:,16-12,15] = [0.99875e+03,0.10040e+04,0.10068e+04,0.10078e+04,0.10089e+04]
  kb[:,17-12,15] = [0.99816e+03,0.10041e+04,0.10080e+04,0.10079e+04,0.10078e+04]
  kb[:,18-12,15] = [0.99861e+03,0.10065e+04,0.10071e+04,0.10090e+04,0.10078e+04]
  kb[:,19-12,15] = [0.10001e+04,0.10056e+04,0.10072e+04,0.10079e+04,0.10077e+04]
  kb[:,20-12,15] = [0.99952e+03,0.10059e+04,0.10074e+04,0.10090e+04,0.10087e+04]
  kb[:,21-12,15] = [0.10000e+04,0.10061e+04,0.10085e+04,0.10080e+04,0.10085e+04]
  kb[:,22-12,15] = [0.10012e+04,0.10064e+04,0.10076e+04,0.10090e+04,0.10073e+04]
  kb[:,23-12,15] = [0.10032e+04,0.10067e+04,0.10088e+04,0.10079e+04,0.10070e+04]
  kb[:,24-12,15] = [0.10052e+04,0.10070e+04,0.10079e+04,0.10089e+04,0.10067e+04]
  kb[:,25-12,15] = [0.10057e+04,0.10083e+04,0.10079e+04,0.10098e+04,0.10062e+04]
  kb[:,26-12,15] = [0.10061e+04,0.10085e+04,0.10080e+04,0.10075e+04,0.10068e+04]
  kb[:,27-12,15] = [0.10066e+04,0.10087e+04,0.10111e+04,0.10072e+04,0.10051e+04]
  kb[:,28-12,15] = [0.10069e+04,0.10089e+04,0.10079e+04,0.10068e+04,0.10054e+04]
  kb[:,29-12,15] = [0.10072e+04,0.10069e+04,0.10077e+04,0.10063e+04,0.10034e+04]
  kb[:,30-12,15] = [0.10075e+04,0.10111e+04,0.10075e+04,0.10068e+04,0.10088e+04]
  kb[:,31-12,15] = [0.10077e+04,0.10080e+04,0.10072e+04,0.10061e+04,0.10007e+04]
  kb[:,32-12,15] = [0.10078e+04,0.10089e+04,0.10068e+04,0.10043e+04,0.10013e+04]
  kb[:,33-12,15] = [0.10079e+04,0.10077e+04,0.10063e+04,0.10044e+04,0.10009e+04]
  kb[:,34-12,15] = [0.10090e+04,0.10075e+04,0.10058e+04,0.10025e+04,0.10002e+04]
  kb[:,35-12,15] = [0.10080e+04,0.10073e+04,0.10054e+04,0.10021e+04,0.99972e+03]
  kb[:,36-12,15] = [0.10090e+04,0.10073e+04,0.10054e+04,0.10020e+04,0.99964e+03]
  kb[:,37-12,15] = [0.10080e+04,0.10075e+04,0.10057e+04,0.10024e+04,0.10001e+04]
  kb[:,38-12,15] = [0.10090e+04,0.10076e+04,0.10060e+04,0.10040e+04,0.10006e+04]
  kb[:,39-12,15] = [0.10079e+04,0.10077e+04,0.10063e+04,0.10034e+04,0.10010e+04]
  kb[:,40-12,15] = [0.10089e+04,0.10089e+04,0.10068e+04,0.10054e+04,0.10014e+04]
  kb[:,41-12,15] = [0.10087e+04,0.10090e+04,0.10082e+04,0.10051e+04,0.10018e+04]
  kb[:,42-12,15] = [0.10085e+04,0.10090e+04,0.10075e+04,0.10058e+04,0.10025e+04]
  kb[:,43-12,15] = [0.10071e+04,0.10079e+04,0.10078e+04,0.10086e+04,0.10037e+04]
  kb[:,44-12,15] = [0.10067e+04,0.10088e+04,0.10079e+04,0.10071e+04,0.10048e+04]
  kb[:,45-12,15] = [0.10061e+04,0.10075e+04,0.10090e+04,0.10096e+04,0.10057e+04]
  kb[:,46-12,15] = [0.10054e+04,0.10071e+04,0.10068e+04,0.10078e+04,0.10065e+04]
  kb[:,47-12,15] = [0.10020e+04,0.10065e+04,0.10087e+04,0.10079e+04,0.10093e+04]
  kb[:,48-12,15] = [0.99938e+03,0.10058e+04,0.10073e+04,0.10079e+04,0.10087e+04]
  kb[:,49-12,15] = [0.99869e+03,0.10038e+04,0.10079e+04,0.10099e+04,0.10089e+04]
  kb[:,50-12,15] = [0.99593e+03,0.10003e+04,0.10062e+04,0.10075e+04,0.10090e+04]
  kb[:,51-12,15] = [0.99417e+03,0.99877e+03,0.10055e+04,0.10061e+04,0.10090e+04]
  kb[:,52-12,15] = [0.99334e+03,0.99717e+03,0.10030e+04,0.10077e+04,0.10077e+04]
  kb[:,53-12,15] = [0.99034e+03,0.99548e+03,0.99989e+03,0.10071e+04,0.10085e+04]
  kb[:,54-12,15] = [0.98844e+03,0.99385e+03,0.99849e+03,0.10064e+04,0.10081e+04]
  kb[:,55-12,15] = [0.98654e+03,0.99323e+03,0.99708e+03,0.10005e+04,0.10066e+04]
  kb[:,56-12,15] = [0.98558e+03,0.99048e+03,0.99560e+03,0.99999e+03,0.10071e+04]
  kb[:,57-12,15] = [0.98349e+03,0.98867e+03,0.99404e+03,0.99866e+03,0.10054e+04]
  kb[:,58-12,15] = [0.98039e+03,0.98688e+03,0.99250e+03,0.99837e+03,0.10028e+04]
  kb[:,59-12,15] = [0.97954e+03,0.98614e+03,0.99186e+03,0.99679e+03,0.10019e+04]
  kb[:,13-12,16] = [0.98905e+03,0.99319e+03,0.10065e+04,0.10075e+04,0.10080e+04]
  kb[:,14-12,16] = [0.98947e+03,0.10009e+04,0.10067e+04,0.10076e+04,0.10005e+04]
  kb[:,15-12,16] = [0.98993e+03,0.10052e+04,0.10068e+04,0.10076e+04,0.10080e+04]
  kb[:,16-12,16] = [0.99040e+03,0.10054e+04,0.10069e+04,0.10077e+04,0.10005e+04]
  kb[:,17-12,16] = [0.99823e+03,0.10056e+04,0.99955e+03,0.10078e+04,0.10080e+04]
  kb[:,18-12,16] = [0.99867e+03,0.99832e+03,0.10071e+04,0.10003e+04,0.10079e+04]
  kb[:,19-12,16] = [0.99168e+03,0.10059e+04,0.10072e+04,0.10078e+04,0.10079e+04]
  kb[:,20-12,16] = [0.99955e+03,0.10061e+04,0.10073e+04,0.99297e+03,0.10003e+04]
  kb[:,21-12,16] = [0.10000e+04,0.10063e+04,0.99994e+03,0.10079e+04,0.10002e+04]
  kb[:,22-12,16] = [0.10007e+04,0.10066e+04,0.10075e+04,0.10005e+04,0.10075e+04]
  kb[:,23-12,16] = [0.10052e+04,0.10068e+04,0.10002e+04,0.10080e+04,0.10072e+04]
  kb[:,24-12,16] = [0.10056e+04,0.10070e+04,0.10078e+04,0.10005e+04,0.10068e+04]
  kb[:,25-12,16] = [0.10060e+04,0.99977e+03,0.10079e+04,0.99294e+03,0.10064e+04]
  kb[:,26-12,16] = [0.10064e+04,0.99996e+03,0.10080e+04,0.10077e+04,0.99839e+03]
  kb[:,27-12,16] = [0.10067e+04,0.10001e+04,0.98561e+03,0.10074e+04,0.10052e+04]
  kb[:,28-12,16] = [0.10069e+04,0.10003e+04,0.10080e+04,0.10070e+04,0.99700e+03]
  kb[:,29-12,16] = [0.10072e+04,0.10078e+04,0.10079e+04,0.10065e+04,0.10036e+04]
  kb[:,30-12,16] = [0.10074e+04,0.98555e+03,0.10002e+04,0.99845e+03,0.95900e+03]
  kb[:,31-12,16] = [0.10076e+04,0.10080e+04,0.10074e+04,0.99776e+03,0.10037e+04]
  kb[:,32-12,16] = [0.10077e+04,0.10005e+04,0.10070e+04,0.10044e+04,0.10038e+04]
  kb[:,33-12,16] = [0.10079e+04,0.10079e+04,0.10064e+04,0.99602e+03,0.10038e+04]
  kb[:,34-12,16] = [0.10005e+04,0.10077e+04,0.10059e+04,0.10036e+04,0.10037e+04]
  kb[:,35-12,16] = [0.10080e+04,0.10075e+04,0.10055e+04,0.10037e+04,0.10036e+04]
  kb[:,36-12,16] = [0.10005e+04,0.10075e+04,0.10055e+04,0.10037e+04,0.10036e+04]
  kb[:,37-12,16] = [0.10080e+04,0.10077e+04,0.10058e+04,0.10037e+04,0.10037e+04]
  kb[:,38-12,16] = [0.10004e+04,0.10078e+04,0.10062e+04,0.99569e+03,0.10037e+04]
  kb[:,39-12,16] = [0.10078e+04,0.10079e+04,0.10065e+04,0.10036e+04,0.10038e+04]
  kb[:,40-12,16] = [0.10003e+04,0.10005e+04,0.10070e+04,0.99698e+03,0.10038e+04]
  kb[:,41-12,16] = [0.10001e+04,0.10005e+04,0.99992e+03,0.10052e+04,0.10037e+04]
  kb[:,42-12,16] = [0.99994e+03,0.10005e+04,0.10077e+04,0.10059e+04,0.10036e+04]
  kb[:,43-12,16] = [0.10071e+04,0.10078e+04,0.10079e+04,0.99172e+03,0.10038e+04]
  kb[:,44-12,16] = [0.10068e+04,0.10002e+04,0.10080e+04,0.10072e+04,0.10049e+04]
  kb[:,45-12,16] = [0.10063e+04,0.10074e+04,0.10005e+04,0.99276e+03,0.10059e+04]
  kb[:,46-12,16] = [0.10057e+04,0.10071e+04,0.10078e+04,0.10079e+04,0.10067e+04]
  kb[:,47-12,16] = [0.10009e+04,0.10067e+04,0.10001e+04,0.10080e+04,0.99247e+03]
  kb[:,48-12,16] = [0.99941e+03,0.10061e+04,0.10073e+04,0.10079e+04,0.10004e+04]
  kb[:,49-12,16] = [0.99034e+03,0.10054e+04,0.99943e+03,0.99276e+03,0.10005e+04]
  kb[:,50-12,16] = [0.99601e+03,0.10003e+04,0.10064e+04,0.10074e+04,0.10005e+04]
  kb[:,51-12,16] = [0.99422e+03,0.99882e+03,0.10058e+04,0.10071e+04,0.10003e+04]
  kb[:,52-12,16] = [0.98496e+03,0.99726e+03,0.10013e+04,0.99930e+03,0.10076e+04]
  kb[:,53-12,16] = [0.99030e+03,0.99555e+03,0.99990e+03,0.99882e+03,0.99991e+03]
  kb[:,54-12,16] = [0.98836e+03,0.98653e+03,0.99856e+03,0.99827e+03,0.99962e+03]
  kb[:,55-12,16] = [0.98643e+03,0.98485e+03,0.99716e+03,0.10012e+04,0.10067e+04]
  kb[:,56-12,16] = [0.97714e+03,0.99044e+03,0.99567e+03,0.99999e+03,0.99885e+03]
  kb[:,57-12,16] = [0.97511e+03,0.98860e+03,0.99409e+03,0.99872e+03,0.10058e+04]
  kb[:,58-12,16] = [0.98036e+03,0.98678e+03,0.99251e+03,0.99003e+03,0.10052e+04]
  kb[:,59-12,16] = [0.97953e+03,0.98603e+03,0.99185e+03,0.99687e+03,0.10009e+04]
end


using MAT
using MatrixProductStates
BLAS.set_num_threads(1)

#Spin model parameters
const na = 70
const od = 33
const ad = 2.5*60/na #atomic distance in micro m
const cloud_sig = 23.6344/ad
const cloud_coup = exp.(-((1:na)-(na+1)/2).^2/2/cloud_sig^2)
const gam_1d = cloud_coup/sum(cloud_coup)/2*od
const gam_eg = 1.0
const del_p = 0.0
const k_wg = 0.5*pi
const rj = collect(1.0:na) #sort!(rand(na)*na)#

#Pulse shape parameters
const k_in = k_wg
const gam_exp = 2*pi*6.065e6
const f_amp_exp = sqrt.([0.0589,0.0721,0.0995,0.1637,0.3081,0.6395,1.2121,1.8898,
    2.6039,4.2050,7.1291,10.4823,17.0816,27.1592,42.8378,71.7874]/gam_exp/1e-6)
const f_amp = f_amp_exp[1]

#Rydberg specific parameters
const om = 10.0/6.065/2.0
const k_c = 0.0
const del_s = 0.0
const gam_es = 0.0

# rmax = 20 na = 70
const uu = [524.328 + 89.9671im, 524.328 - 89.9671im, 649.737 + 219.429im, 649.737 - 219.429im, 488.928]
const lam = [-0.228929 + 0.189648im, -0.228929 - 0.189648im, 0.146279 + 0.297638im, 0.146279 - 0.297638im, 0.480032]

# rmax = 30 na = 60
#uu = [1055.32, 1156.39 + 264.599im, 1156.39 - 264.599im, 781.661]
#lam = [-0.220397, 0.0354409 + 0.244048im, 0.0354409 - 0.244048im, 0.396262]
#uu = [15158.7, 11895.1 - 1929.93im, 11895.1 + 1929.93im, 2399.43, 16.4076]
#lam = [-0.103573, 0.015448 + 0.126927im, 0.015448 - 0.126927im, 0.30542, 0.606657]
#uu = [90866.8, 59173.0 - 18223.1im, 59173.0 + 18223.1im, 5257.17, 107.47, 1.26558]
#lam = [-0.0606795, 0.0112804 + 0.0784944im, 0.0112804 - 0.0784944im, 0.246432, 0.48336, 0.720283]

# rmax = 20 na = 60
#const uu = [0.704528 + 135.691im, 0.704528 - 135.691im, 225.917]
#const lam = [0.226209 + 0.32374im, 0.226209 - 0.32374im, 0.475997]

#uu = [ 2255.07, 686.082 + 91.2384im, 686.082 - 91.2384im, 477.07]
#lam = [-0.0998312, 0.0646318 + 0.25377im, 0.0646318-0.25377im, 0.424938]
#uu = [-508.38 - 1018.61im, -508.38 + 1018.61im, -1546.1 - 1942.74im, -1546.1 + 1942.74im,
#        28.5203]
#lam = [ -0.0374226 + 0.193987im, -0.0374226 - 0.193987im, 0.250214 + 0.0887921im,
#       0.250214 - 0.0887921im, 0.58104]
#uu = [-658.963 + 91.3732im, -658.963 - 91.3732im, 309.325 + 1085.49im, 309.325 - 1085.49im,
#       709.204, 4.33617]
#lam = [-0.0921073 + 0.171392im, -0.0921073 - 0.171392im, 0.124872 + 0.194013im,
#       0.124872 - 0.194013im, 0.382177, 0.670268]

# rmax = 20 na = 50
#uu = [1740.29, 2277.16 + 454.662im, 2277.16 - 454.662im, 945.959, 6.06471]
#lam = [-0.126487, 0.0183872 + 0.140738im, 0.0183872 - 0.140738im, 0.296793, 0.601841]
#uu = [ 57.1125, 170.91 + 211.667im, 170.91 - 211.667im, 336.281]
#lam = [ -0.298989, 0.0548962 + 0.260436im, 0.0548962 - 0.260436im, 0.382366]

# rmax = 20 na = 40
#uu = [19014.9 - 3038.05im, 19014.9 + 3038.05im, 2635.02, 47.4863, 0.531659]
#lam = [-0.0167917 + 0.0377784im, -0.0167917 - 0.0377784im, 0.155726, 0.386972, 0.653397]

# rmax = 20 na = 30
#uu = [-2652.81-1864.93im, -2652.81 + 1864.93im, 274.344, 5.91125, 0.0704647]
#lam = [0.0313486 + 0.0360902im, 0.0313486 - 0.0360902im, 0.179561, 0.407673, 0.666632]

# rmax = 20 na = 20
#uu = [1350.25, 707.234, 15.7581, 0.426333, 0.0053846]
#lam = [-0.0138335, 0.0500349, 0.19712, 0.418572, 0.672633]

# pure dephasing
const gam_gg = 0.0
const gam_ss = 40e3/6.065e6
const gam_ee = 0.0

# simulation parameters
const setrand = 1 # random number generator seed
const dt = 0.01
const t_fin = 10000.0
const d = 3
const d_max = 10
const measure_int = 5
const path_data = "/home/jdouglas/data/" # save directory
const base_filename = string(path_data, "Ryd_Jumps_N", na, "_D", d_max,
    "_Tf", t_fin, "_f", round(f_amp, 3), "_dt", dt, "_traj", setrand)

# envelope of input pulse
function f(t)

   f_amp

end

# local spin operators
const hgg = [1 0 0; 0 0 0; 0 0 0]
const hee = [0 0 0; 0 1 0; 0 0 0]
const hss = [0 0 0; 0 0 0; 0 0 1]
const hge = [0 1 0; 0 0 0; 0 0 0]
const heg = hge'
const hes = [0 0 0; 0 0 1; 0 0 0]
const hse = hes'
const id = eye(3)

# types used in mps representation
const TN = Complex{Float64}
const TA = Array

function time_evolve()

    # create jump and measurement operators
    jumpleft = makempo(TN, TA, [jj->id jj->im*sqrt(gam_1d[jj]/2)*exp(im*k_wg*rj[jj])*hge;
                                jj->0 jj->id], na, d)
    jumpright = makempo(TN, TA, [jj->id jj->im*sqrt(gam_1d[jj]/2)*exp(-im*k_wg*rj[jj])*hge;
                                 jj->0 jj->id], na, d)
    jumpright[1][1, :, 2, :] = f(0.0)*id + im*sqrt(gam_1d[1]/2)*exp(-im*k_wg*rj[1])*hge
    ir_mpo = applyMPOtoMPO(jumpright, conj_mpo(jumpright))
    ir2_mpo = applyMPOtoMPO(applyMPOtoMPO(jumpright, ir_mpo), conj_mpo(jumpright))

    # step times and measurement times
    t = 0.0:dt:t_fin
    t_m = t[1:measure_int:end]
    tstep = length(t) - 1

    # initialize output measures
    p_coh = zeros(tstep)
    p_sum = zeros(tstep)
    p_acc = zeros(tstep, 2na + 4)
    e_pop = zeros(na, length(t_m))
    s_pop = zeros(na, length(t_m))
    epop = zeros(na)
    spop = zeros(na)
    I_r = zeros(TN, length(t_m))
    I2_r = zeros(TN, length(t_m))
    times = zeros(tstep, 4)
    t_r = Float64[]
    t_l = Float64[]
    t_eg = Float64[]
    pos_eg = Float64[]
    t_ss = Float64[]
    pos_ss = Float64[]
    cont = zeros(Int64, 4)

    #Construct left-canonical MPS for the initial state and temp arrays for time evolution
    A, dims =  mpsgroundstate(TN, TA, na, d_max, d)
    A_coh = similar(A)
    A_r = similar(A)
    A_l = similar(A)
    A1 = similar(A)
    A2 = similar(A)
    A3 = similar(A)
    A_coh .= copy.(A)
    A_r .= copy.(A)
    A_l .= copy.(A)
    A1 .= copy.(A)
    A2 .= copy.(A)
    A3 .= copy.(A)
    mpo_size = 4 + length(lam)
    envop = build_env(TN, TA, dims, dims, ones(na + 1)*mpo_size)
    envop_jump = build_env(TN, TA, dims, dims, ones(na + 1)*2)
    env1 = build_env(TN, TA, dims, dims)
    env2 = build_env(TN, TA, dims, dims)
    env3 = build_env(TN, TA, dims, dims)

    # initial measurements
    I_r[1] = scal_op_prod(A, ir_mpo, A)
    I2_r[1] = scal_op_prod(A, ir2_mpo, A)
    measure_excitations!((@view e_pop[:, 1]), (@view s_pop[:, 1]), A)

    # random jump variables
    srand(setrand)
    r1 = rand(tstep)
    r2 = rand(tstep)

    # create evolution operators for Runge-Kutta algorithm
    H1 = rydberg_hamiltonian(TN, TA, 1, dt/2, t[1])
    H2 = rydberg_hamiltonian(TN, TA, 0, 1, t[1] + dt/2)
    H3 = rydberg_hamiltonian(TN, TA, 1, dt/2, t[2])

    H = [H1, H2, H3]

    println("Initialisation complete")

    #EVOLUTION
    for i = 1:tstep
        time_in = time()

        # update time evolution operators
        update_rydberg_hamiltonian!(H1, 1, dt/2, t[i])
        update_rydberg_hamiltonian!(H2, 0, 1, t[i] + dt/2)
        update_rydberg_hamiltonian!(H3, 1, dt/2, t[i+1])

        # hoherent time evolution probability
        RK4stp_apply_H_half!(dt, H, A_coh, A, A1, A2, A3, env1, env2, env3, envop)
        p_coh[i] = normalize_lo!(A_coh)
        times[i, 1] = time() - time_in
        if r1[i] < p_coh[i]
            A .= copy.(A_coh)
            times[i,2] = 0
        else
            #Update waveguide jump operator
            jumpright[1][1,:,2,:] = f(t[i])*id + im*sqrt(gam_1d[1]/2)*exp(-im*k_wg*rj[1])*hge

            #Jump probabilities: decay into the waveguide
            compress_var_apply_H!(A_r, envop_jump, A, A, jumpright, 1)
            p_r = normalize_lo!(A_r)
            compress_var_apply_H!(A_l, envop_jump, A, A, jumpleft, 1)
            p_l = normalize_lo!(A_l)

            #Jump probabilities: decay outside the waveguide
            measure_excitations!(epop, spop, A)
            p_acc[i, 1] = 0.0
            for j = 1:na
                p_acc[i, j + 1] = p_acc[i, j] + dt*gam_eg*epop[j]
            end
            for j = 1:na
                p_acc[i, na + j + 1] = p_acc[i, na + j] + dt*gam_ss*spop[j]
            end
            p_acc[i, 2na + 2] = p_acc[i, 2na + 1] + dt*p_r
            p_acc[i, 2na + 3] = p_acc[i, 2na + 2] + dt*p_l
            p_sum[i] = p_acc[i, 2na + 3]
            println(p_sum[i] + p_coh[i])
            println(t[i])
            p_acc[i, :] = p_acc[i, :]/p_acc[i, 2na + 3]

            #Do jumps according to above probabilities
            if r2[i] < p_acc[i, na + 1]
                for j = 1:na
                    if r2[i] < p_acc[i, j + 1]
                        apply_site_operator!(A[j], hge)
                        leftorthonormalizationQR!(A)
                        push!(t_eg, t[i])
                        push!(pos_eg, j)
                        cont[1] += 1
                        break
                    end
                end
            elseif r2[i] < p_acc[i, 2na + 1]
                for j = 1:na
                    if r2[i] < p_acc[i, na + j + 1]
                        apply_site_operator!(A[j], hss)
                        leftorthonormalizationQR!(A)
                        push!(t_ss, t[i])
                        cont[2] += 1
                        break
                    end
                end
            elseif r2[i] < p_acc[i, 2na + 2]
                A .= copy.(A_r)
                push!(t_r, t[i])
                cont[3] += 1
            else
                A .= copy.(A_l)
                push!(t_l, t[i])
                cont[4] += 1
            end
            times[i, 2] = time() - times[i, 1] - time_in

        end

        # measurements (every measure_int time steps)
        if rem(i, measure_int) == 0

            mind = div(i, measure_int)+1
            # atomic state populations
            measure_excitations!((@view e_pop[:, mind]), (@view s_pop[:, mind]), A)
            times[i, 3] = time() - times[i, 2]  - times[i, 1] - time_in

            # output right field update
            jumpright[1][1, :, 2, :] = f(t[i + 1])*id +
                    im*sqrt(gam_1d[1]/2)*exp(-im*k_wg*rj[1])*hge

            ir_mpo[1] =  apply_site_MPOtoMPO(jumpright[1], conj_site_mpo(jumpright[1]))
            ir2_mpo[1] =  apply_site_MPOtoMPO(
                apply_site_MPOtoMPO(jumpright[1], ir_mpo[1]),
                conj_site_mpo(jumpright[1]))

            I_r[mind] = scal_op_prod(A, ir_mpo, A)
            I2_r[mind] = scal_op_prod(A, ir2_mpo, A)

            times[i,4] = time() - times[i, 3] - times[i, 2] - times[i, 1] - time_in

        end

        # saving temp file
        if rem(i, measure_int*2000) == 0

            write_data_file(string(base_filename,"_temp.mat"), t, t_m, e_pop,
                            s_pop, p_coh, p_acc, p_sum, I_r, I2_r, cont, t_r,
                            t_l, t_eg, t_ss, pos_eg, A, times)

        end

    end

    # saving output file
    write_data_file(string(base_filename,".mat"), t, t_m, e_pop,
                    s_pop, p_coh, p_acc, p_sum, I_r, I2_r, cont, t_r,
                    t_l, t_eg, t_ss, pos_eg, A, times)

end

function write_data_file(filename, t, t_m, e_pop, s_pop, p_coh, p_acc,
                         p_sum, I_r, I2_r, cont, t_r, t_l, t_eg, t_ss,
                         pos_eg, A, times)

    file = matopen(string(filename),"w")
    write(file, "TN", string(TN))
    write(file, "TA", string(TA))
    write(file, "na", na)
    write(file, "gam_1d", gam_1d)
    write(file, "gam_eg", gam_eg)
    write(file, "del_p", del_p)
    write(file, "k_wg", k_wg)
    write(file, "rj", rj)
    write(file, "k_in", k_in)
    write(file, "f_amp", f_amp)
    write(file, "f", f.(t_m))
    write(file, "om", om)
    write(file, "k_c", k_c)
    write(file, "del_s", del_s)
    write(file, "gam_es", gam_es)
    write(file, "uu", uu)
    write(file, "lam", lam)
    write(file, "gam_gg", gam_gg)
    write(file, "gam_ss", gam_ss)
    write(file, "gam_ee", gam_ee)
    write(file, "dt", dt)
    write(file, "t_fin", t_fin)
    write(file, "d_max", d_max)
    write(file, "measure_int", measure_int)

    write(file, "t", collect(t))
    write(file, "t_m", collect(t_m))
    write(file, "e_pop", e_pop)
    write(file, "s_pop", s_pop)
    write(file, "p_coh", p_coh)
#    write(file, "p_acc", p_acc)
    write(file, "p_sum", p_sum)
    write(file, "I_r", I_r)
    write(file, "I2_r", I2_r)
    write(file, "times", sum(times,1))
    write(file, "cont", cont)
    write(file, "t_r", t_r)
    write(file, "t_l", t_l)
    write(file, "t_eg", t_eg)
    write(file, "t_ss", t_ss)
    write(file, "pos_eg", pos_eg)
    write(file, "A", collect.(A))

    close(file)

end

function rydberg_hamiltonian(::Type{TN}, ::Type{TA}, con, delt, time) where {TN, TA}

    H = Array{TA{TN, 4}, 1}(na)

    no_exp = length(lam)
    max_ind = 4 + no_exp
    drj = diff(rj)
    ph = exp.(im*k_wg*drj)
    cp = sqrt.(delt*gam_1d/2)

    H[1] = zeros(TN, 1, 3, max_ind, 3)
    H[1][1, :, 1, :] = id
    H[1][1, :, 2, :] = -cp[1]*ph[1]*heg
    H[1][1, :, 3, :] = -cp[1]*ph[1]*hge
    for exp_ind = 1:no_exp
        H[1][1, :, 3 + exp_ind, :] = sqrt(-im*delt*uu[exp_ind]*lam[exp_ind])*hss
    end
    H[1][1, :, max_ind, :] = local_h(1, con, delt, time)

    for jj = 2:(na - 1)
        H[jj] = zeros(TN, max_ind, 3, max_ind, 3)
        H[jj][1, :, 1, :] = id
        H[jj][1, :, 2, :] = -cp[jj]*ph[jj]*heg
        H[jj][1, :, 3, :] = -cp[jj]*ph[jj]*hge
        H[jj][1, :, max_ind, :] = local_h(jj, con, delt, time)
        H[jj][2, :, max_ind, :] = cp[jj]*hge
        H[jj][3, :, max_ind, :] = cp[jj]*heg
        H[jj][2, :, 2, :] = ph[jj]*id
        H[jj][3, :, 3, :] = ph[jj]*id
        for exp_ind = 1:no_exp
            H[jj][1, :, 3 + exp_ind, :] = sqrt(-im*delt*uu[exp_ind]*lam[exp_ind])*hss
            H[jj][3 + exp_ind, :, max_ind, :] = sqrt(-im*delt*uu[exp_ind]*lam[exp_ind])*hss
            H[jj][3 + exp_ind, :, 3 + exp_ind, :] = lam[exp_ind]*id
        end
        H[jj][max_ind, :, max_ind, :] = id
    end

    H[na] = zeros(TN, max_ind, 3, 1, 3)
    H[na][1, :, 1, :] = local_h(na, con, delt, time)
    H[na][2, :, 1, :] = cp[na]*hge
    H[na][3, :, 1, :] = cp[na]*heg
    for exp_ind = 1:no_exp
        H[na][3 + exp_ind, :, 1, :] = sqrt(-im*delt*uu[exp_ind]*lam[exp_ind])*hss
    end
    H[na][max_ind, :, 1, :] = id

    return H

end

function local_h(jj, con, delt, time)

    con*id/na + delt*((im*del_p -
        (gam_eg + gam_ee + gam_1d[jj])/2)*hee +
        im*f(time)*sqrt(gam_1d[jj]/2)*exp(im*k_in*rj[jj])*heg -
        1/(2na)*abs2(f(time))*id + (im*del_s - gam_ss/2)*hss +
        im*om*(exp(im*k_c*rj[jj])*hes + exp(-im*k_c*rj[jj])*hse))

end

function update_rydberg_hamiltonian!(H, con, delt, time)

    no_exp = length(lam)
    max_ind = 4 + no_exp

    H[1][1, :, max_ind, :] = local_h(1, con, delt, time)

    for jj = 2:(na - 1)
        H[jj][1, :, max_ind, :] = local_h(jj, con, delt, time)
    end

    H[na][1, :, 1, :] = local_h(na, con, delt, time)

end

# calculate local expectation values, assuming mps is already left normalized
function measure_excitations!(e_pop, s_pop, A)

    # initialize right environment
    env = ones(1, 1)

    for n = na:(-1):1
        Aenv = prod_LR(A[n], env)
        e_pop[n] = real(local_exp(Aenv, A[n], hee))
        s_pop[n] = real(local_exp(Aenv, A[n], hss))
        env = update_renv(A[n], Aenv)
    end

end

time_evolve()


test_LP() = CoreModel(
    zeros(4, 3),
    zeros(4),
    ones(3),
    ones(3),
    ones(3),
    ["r$x" for x = 1:3],
    ["m$x" for x = 1:4],
)

test_simpleLP() = CoreModel(
    [
        1.0 1.0
        -1.0 1.0
    ],
    [3.0, 1.0],
    [-0.25, 1.0],
    -ones(2),
    2.0 * ones(2),
    ["r$x" for x = 1:2],
    ["m$x" for x = 1:2],
)

test_simpleLP2() = CoreModel(
    zeros(2, 2),
    [0.0, 0.0],
    [-0.25, 1.0],
    -ones(2),
    2.0 * ones(2),
    ["r$x" for x = 1:2],
    ["m$x" for x = 1:2],
)

test_sparseLP() = CoreModel(
    sprand(4000, 3000, 0.5),
    sprand(4000, 0.5),
    sprand(3000, 0.5),
    sprand(3000, 0.5),
    sprand(3000, 0.5),
    ["r$x" for x = 1:3000],
    ["m$x" for x = 1:4000],
)

test_coupledLP() = CoreModelCoupled(
    CoreModel(
        sprand(4000, 3000, 0.5),
        sprand(4000, 0.5),
        sprand(3000, 0.5),
        sprand(3000, 0.5),
        sprand(3000, 0.5),
        ["r$x" for x = 1:3000],
        ["m$x" for x = 1:4000],
    ),
    sprand(2000, 3000, 0.5),
    sprand(2000, 0.5),
    sprand(2000, 0.5),
)

test_toyModel() = CoreModel(
    [
        -1.0 1.0 0.0 0.0 0.0 0.0 0.0
        -2.0 0.0 1.0 0.0 0.0 0.0 0.0
        1.0 0.0 0.0 0.0 0.0 0.0 -1.0
        0.0 -1.0 0.0 -1.0 0.0 0.0 0.0
        0.0 0.0 -1.0 0.0 -1.0 0.0 0.0
        0.0 0.0 0.0 0.0 0.0 -1.0 1.0
    ],
    zeros(6),
    [0, 0, 0, 0, 0, 0, 1.0],
    fill(-1000.0, 7),
    fill(1000.0, 7),
    ["r1", "m1t", "m3t", "EX_m1(e)", "EX_m3(e)", "EX_biomass(e)", "biomass1"],
    ["m1[c]", "m3[c]", "m2[c]", "m1[e]", "m3[e]", "biomass[c]"],
)

const reaction_standard_gibbs_free_energies = Dict(
    #=
    ΔᵣG⁰ data from Equilibrator using the E. coli core model's reactions
    To generate this data manually, go to https://equilibrator.weizmann.ac.il/ and
    enter each reaction into the search bar, fill in the ΔG⁰ below from there. To generate
    automatically, use the eQuilibrator.jl package.
    =#
    "ACALD" => -21.268198981756314,
    "PTAr" => 8.651025243027988,
    "ALCD2x" => 17.47260052762408,
    "PDH" => -34.246780702043225,
    "PYK" => -24.48733600711958,
    "CO2t" => 0.0,
    "MALt2_2" => -6.839209921974724,
    "CS" => -39.330940148207475,
    "PGM" => -4.470553692565886,
    "TKT1" => -1.4976299544215408,
    "ACONTa" => 8.46962176350985,
    "GLNS" => -15.771242654033706,
    "ICL" => 9.53738025507404,
    "FBA" => 23.376920310319235,
    "SUCCt3" => -43.97082007726285,
    "FORt2" => -3.4211767267069124,
    "G6PDH2r" => -7.3971867270165035,
    "AKGDH" => -28.235442362320782,
    "TKT2" => -10.318436817276165,
    "FRD7" => 73.61589524884772,
    "SUCOAS" => -1.1586968559555615,
    "FBP" => -11.606887851480572,
    "ICDHyr" => 5.398885342794983,
    "AKGt2r" => 10.085299238039797,
    "GLUSy" => -47.21690395283849,
    "TPI" => 5.621932460512994,
    "FORt" => 13.509690780237293,
    "ACONTb" => -1.622946931741609,
    "GLNabc" => -30.194559792842753,
    "RPE" => -3.3881719029747615,
    "ACKr" => 14.027197131450492,
    "THD2" => -33.84686533309243,
    "PFL" => -19.814421615735,
    "RPI" => 4.477649590945703,
    "D_LACt2" => -3.4223903975852004,
    "TALA" => -0.949571985515206,
    "PPCK" => 10.659841402564751,
    "ACt2r" => -3.4196348363397995,
    "NH4t" => -13.606633927039097,
    "PGL" => -25.94931748161696,
    "NADTRHD" => -0.014869680795754903,
    "PGK" => 19.57192102020454,
    "LDH_D" => 20.04059765689044,
    "ME1" => 12.084968268076864,
    "PIt2r" => 10.415108493818785,
    "ATPS4r" => -37.570267233299816,
    "PYRt2" => -3.422891289768689,
    "GLCpts" => -45.42430981510088,
    "GLUDy" => 32.834943812395665,
    "CYTBD" => -59.700410815493775,
    "FUMt2_2" => -6.845105500839001,
    "FRUpts2" => -42.67529760694199,
    "GAPD" => 0.5307809794271634,
    "H2Ot" => -1.5987211554602254e-14,
    "PPC" => -40.81304419704113,
    "NADH16" => -80.37770501380615,
    "PFK" => -18.546314942995934,
    "MDH" => 25.912462872631522,
    "PGI" => 2.6307087407442395,
    "O2t" => 0.0,
    "ME2" => 12.099837948872533,
    "GND" => 10.312275879236381,
    "SUCCt2_2" => -6.82178244356977,
    "GLUN" => -14.381960140443113,
    "ETOHt2r" => -16.930867506944217,
    "ADK1" => 0.3893321583896068,
    "ACALDt" => -3.197442310920451e-14,
    "SUCDi" => -73.61589524884772,
    "ENO" => -3.8108376097261782,
    "MALS" => -39.22150045995042,
    "GLUt2r" => -3.499043558772904,
    "PPS" => -6.0551989457468665,
    "FUM" => -3.424133018702122,
)