text:
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


text:
include("modulos/datasets.jl")
include("modulos/attributes_from_dataset.jl")
include("modulos/models_cross_validation.jl")

#=
hay tres extracciones
    · :A1 (aproximación 1)
    · :A21 (aproximación 2 extracción 1)
    · :A22 (aproximación 2 extracción 2)
=#
extraction = :A33;
dataset_name="datasets/faces.data";
change = false;
if extraction==:A1
    x = 6;
    y = 7;
elseif extraction==:A21
    x = 42;
    y = 43;
elseif (extraction==:A22 || extraction==:A23)
    x = 36;
    y = 37;
elseif extraction==:A31
    x = 48;
    y = 49;
elseif extraction==:A32
    x = 66;
    y = 67;
elseif extraction==:A33
    x = 72;
    y = 73;
elseif extraction==:A34
    x = 84;
    y = 85;
end;
if (!isfile(dataset_name) || change)
    createDataset(dataset_name,extraction);
end;

dataset = readdlm(dataset_name,',');
inputs = convert(Array{Float64,2}, dataset[:,1:x]);
targets = convert(Array{Any,1},dataset[:,y]);
seed!(1);
numFolds = 10;

# Parametros principales de la RNA y del proceso de entrenamiento
topology = [3];
learningRate = 0.01; # Tasa de aprendizaje
numMaxEpochs = 1000; # Numero maximo de ciclos de entrenamiento
validationRatio = 0.2; # Porcentaje de patrones que se usaran para validacion. Puede ser 0, para no usar validacion
maxEpochsVal = 30; # Numero de ciclos en los que si no se mejora el loss en el conjunto de validacion, se para el entrenamiento
numRepetitionsAANTraining = 50; # Numero de veces que se va a entrenar la RNA para cada fold por el hecho de ser no determinístico el entrenamiento
# Parametros del SVM
kernel = "rbf";
kernelDegree = 2;
kernelGamma = 1;
C=1000000;
# Parametros del arbol de decision
maxDepth = 15;
# Parapetros de kNN
numNeighbors = 1;
# Normalizamos las entradas, a pesar de que algunas se vayan a utilizar para test
#normalizeMinMax!(inputs);
# Entrenamos las RR.NN.AA.

modelHyperparameters = Dict();
modelHyperparameters["topology"] = topology;
modelHyperparameters["learningRate"] = learningRate;
modelHyperparameters["validationRatio"] = validationRatio;
modelHyperparameters["numExecutions"] = numRepetitionsAANTraining;
modelHyperparameters["maxEpochs"] = numMaxEpochs;
modelHyperparameters["maxEpochsVal"] = maxEpochsVal;
modelCrossValidation(:ANN, modelHyperparameters, inputs, targets, numFolds);

# Entrenamos las SVM
modelHyperparameters = Dict();
modelHyperparameters["kernel"] = kernel;
modelHyperparameters["kernelDegree"] = kernelDegree;
modelHyperparameters["kernelGamma"] = kernelGamma;
modelHyperparameters["C"] = C;
modelCrossValidation(:SVM, modelHyperparameters, inputs, targets, numFolds);

# Entrenamos los arboles de decision
modelCrossValidation(:DecisionTree, Dict("maxDepth" => maxDepth), inputs, targets, numFolds);

# Entrenamos los kNN
modelCrossValidation(:kNN, Dict("numNeighbors" => numNeighbors), inputs, targets, numFolds);

ext: jl


text:
a = VecSE2()
@test typeof(a) <: VecSE
@test typeof(a) <: AbstractVec
@test a.x == 0.0
@test a.y == 0.0
@test a.θ == 0.0

b = VecSE2(0.0,0.0,0.0)
@test b.x == 0.0
@test b.y == 0.0
@test b.θ == 0.0

c = VecSE2(VecE2(0.0,0.0),0.0)
@test c.x == 0.0
@test c.y == 0.0
@test c.θ == 0.0

@test length(a) == 3
@test a == b
@test b == c
@test a == c
@test isequal(b, a)
@test copy(a) == a
@test convert(Vector{Float64}, a) == [0.0,0.0,0.0]
@test convert(VecSE2, [1.0,2.0,3.0]) == VecSE2(1.0,2.0,3.0)
@test convert(VecE2, VecSE2(1.0,2.0,3.0)) == VecE2(1.0,2.0)
@test convert(VecE3, VecSE2(1.0,2.0,3.0)) == VecE3(1.0,2.0,3.0)

@test isapprox(polar( 1.0, 0.0, 0.5), VecSE2(1.0,0.0,0.5))
@test isapprox(polar( 2.0, π/2,-0.5), VecSE2(0.0,2.0,-0.5))
@test isapprox(polar( 3.0,-π/2,   1), VecSE2(0.0,-3.0,1.0))
@test isapprox(polar(-0.5,  1π, -1), VecSE2(0.5,0.0,-1))

a = VecSE2(0.0,1.0,π/2)
b = VecSE2(0.5,2.0,-π)

@test a != b

@test a + b == VecSE2(0.5, 3.0, -π/2)
@test a + 1 == VecSE2(1.0, 2.0,  π/2)
@test a + 0.5 == VecSE2(0.5, 1.5, π/2)
@test 1 + a == VecSE2(1.0, 2.0,  π/2)
@test 0.5 + a == VecSE2(0.5, 1.5, π/2)

@test a - b == VecSE2(-0.5, -1.0, 3π/2)
@test a - 1 == VecSE2(-1.0,  0.0,  π/2)
@test a - 0.5 == VecSE2(-0.5, 0.5, π/2)
@test 1 - a == VecSE2(1.0,  0.0,  π/2)
@test 0.5 - a == VecSE2(0.5, -0.5, π/2)

@test a * 2 == VecSE2(0.0, 2.0, π/2)
@test a * 0.5 == VecSE2(0.0, 0.5, π/2)
@test 2a == VecSE2(0.0, 2.0, π/2)
@test 0.5a == VecSE2(0.0, 0.5, π/2)

@test a / 2 == VecSE2(0.0, 0.5, π/2)
@test a / 0.5 == VecSE2(0.0, 2.0, π/2)

@test b^2 == VecSE2(0.25, 4.0, -π)
@test b^0.5 == VecSE2(0.5^0.5, 2.0^0.5, -π)

@test a % 2.0 == VecSE2(0.0, 1.0, π/2)

@test clamp(VecSE2(1.0, 10.0, 0.5), 0.0, 5.0) == VecSE2(1.0, 5.0, 0.5)
@test clamp(VecSE2(-1.0, 4.0, 5.5), 0.0, 5.0) == VecSE2(0.0, 4.0, 5.5)

@test isfinite(a)
@test !isfinite(VecSE2(Inf,0))
@test !isfinite(VecSE2(0,Inf))
@test !isfinite(VecSE2(0,-Inf))
@test !isfinite(VecSE2(Inf,Inf))

@test !isinf(a)
@test isinf(VecSE2(Inf,0))
@test isinf(VecSE2(0,Inf))
@test isinf(VecSE2(0,-Inf))
@test isinf(VecSE2(Inf,Inf))

@test !isnan(a)
@test isnan(VecSE2(NaN,0))
@test isnan(VecSE2(0,NaN))
@test isnan(VecSE2(NaN,NaN))

# @test round(VecSE2(0.25,1.75)) == VecSE2(0.0,2.0)
# @test round(VecSE2(-0.25,-1.75)) == VecSE2(-0.0,-2.0)

# @test floor(VecSE2(0.25,1.75)) == VecSE2(0.0,1.0)
# @test floor(VecSE2(-0.25,-1.75)) == VecSE2(-1.0,-2.0)

# @test ceil(VecSE2(0.25,1.75)) == VecSE2(1.0,2.0)
# @test ceil(VecSE2(-0.25,-1.75)) == VecSE2(-0.0,-1.0)

# @test trunc(VecSE2(0.25,1.75)) == VecSE2(0.0,1.0)
# @test trunc(VecSE2(-0.25,-1.75)) == VecSE2(-0.0,-1.0)

# c = VecSE2(3.0,4.0)

# @test isapprox(abs(a), 1.0)
# @test isapprox(abs(c), 5.0)
# @test isapprox(hypot(a), 1.0)
# @test isapprox(hypot(c), 5.0)
# @test isapprox(abs2(a), 1.0)
# @test isapprox(abs2(c), 25.0)
# @test isapprox(norm(a), VecSE2(0.0,1.0))
# @test isapprox(norm(b), VecSE2(0.5/sqrt(4.25),2.0/sqrt(4.25)))
# @test isapprox(norm(c), VecSE2(3.0/5,4.0/5))

# @test isapprox(dist(a,a), 0.0)
# @test isapprox(dist(b,b), 0.0)
# @test isapprox(dist(c,c), 0.0)
# @test isapprox(dist(a,b), hypot(0.5, 1.0))
# @test isapprox(dist(a,c), hypot(3.0, 3.0))

# @test isapprox(dist2(a,a), 0.0)
# @test isapprox(dist2(b,b), 0.0)
# @test isapprox(dist2(c,c), 0.0)
# @test isapprox(dist2(a,b), 0.5^2 + 1^2)
# @test isapprox(dist2(a,c), 3^2 + 3^2)

# @test isapprox(dot(a, b), 2.0)
# @test isapprox(dot(b, c), 1.5 + 8.0)
# @test isapprox(dot(c, c), abs2(c))

# @test isapprox(proj(b, a, Float64), 2.0)
# @test isapprox(proj(c, a, Float64), 4.0)

# @test isapprox(proj(b, a, VecSE2), VecSE2(0.0, 2.0))
# @test isapprox(proj(c, a, VecSE2), VecSE2(0.0, 4.0))

# @test isapprox(lerp(VecSE2(0,0), VecSE2(2,-2), 0.25), VecSE2(0.5,-0.5))
# @test isapprox(lerp(VecSE2(2,-2), VecSE2(3,-3), 0.5), VecSE2(2.5,-2.5))

# @test isapprox(rot180(b), VecSE2(-0.5,-2.0))
# @test isapprox(rotr90(b), VecSE2( 2.0,-0.5))
# @test isapprox(rotl90(b), VecSE2(-2.0, 0.5))

# @test isapprox(rot(b, 1.0*pi), VecSE2(-0.5,-2.0))
# @test isapprox(rot(b, -π/2), VecSE2( 2.0,-0.5))
# @test isapprox(rot(b,  π/2), VecSE2(-2.0, 0.5))
# @test isapprox(rot(b, 2π), b)


text:
##interlace block operators
function isboundaryrow(A,k)
    for j=1:size(A,2)
        if isafunctional(A[k,j])
            return true
        end
    end

    return false
end

domainscompatible{T<:Operator}(A::Matrix{T})=domainscompatible(map(domainspace,A))

function spacescompatible{T<:Operator}(A::Matrix{T})
    for k=1:size(A,1)
        if !spacescompatible(map(rangespace,vec(A[k,:])))
            return false
        end
    end
    for k=1:size(A,2)
        if !spacescompatible(map(domainspace,A[:,k]))
            return false
        end
    end
    true
end

spacescompatible{T<:Operator}(A::Vector{T}) = spacescompatible(map(domainspace,A))

function domainspace{T<:Operator}(A::Matrix{T})
    @assert spacescompatible(A)

    spl=map(domainspace,vec(A[1,:]))
    if spacescompatible(spl)
        ArraySpace(first(spl),length(spl))
    elseif domainscompatible(spl)
        TupleSpace(spl)
    else
        PiecewiseSpace(spl)
    end
end

function rangespace{T<:Operator}(A::Vector{T})
    @assert spacescompatible(A)

    spl=map(rangespace,A)
    if spacescompatible(spl)
        ArraySpace(first(spl),length(spl))
    else
        TupleSpace(spl)
    end
end

function promotespaces{T<:Operator}(A::Matrix{T})
    isempty(A) && return A
    A=copy(A)#TODO: promote might have different Array type
    for j=1:size(A,2)
        A[:,j]=promotedomainspace(A[:,j])
    end
    for k=1:size(A,1)
        A[k,:]=promoterangespace(vec(A[k,:]))
    end

    # do a second loop as spaces might have been inferred
    # during range space
    for j=1:size(A,2)
        A[:,j]=promotedomainspace(A[:,j])
    end
    A
end

## Interlace operator

immutable InterlaceOperator{T,p,DS,RS,DI,RI,BI} <: Operator{T}
    ops::Array{Operator{T},p}
    domainspace::DS
    rangespace::RS
    domaininterlacer::DI
    rangeinterlacer::RI
    bandinds::BI
end

InterlaceOperator{T,p}(ops::Array{T,p},ds,rs,di,ri,bi) =
    InterlaceOperator{T,p,typeof(ds),typeof(rs),
                        typeof(di),typeof(ri),typeof(bi)}(ops,ds,rs,di,ri,bi)

function InterlaceOperator{T}(ops::Matrix{Operator{T}},ds::Space,rs::Space)
    # calculate bandinds TODO: generalize
    p=size(ops,1)
    if size(ops,2) == p && all(isbanded,ops)
        l,u = 0,0
        for k=1:p,j=1:p
            l=min(l,p*bandinds(ops[k,j],1)+j-k)
        end
        for k=1:p,j=1:p
            u=max(u,p*bandinds(ops[k,j],2)+j-k)
        end
    elseif p == 1 && size(ops,2) == 2 && size(ops[1],2) == 1
        # special case for example
        l,u = min(bandinds(ops[1],1),bandinds(ops[2],1)+1),bandinds(ops[2],2)+1
    else
        l,u = (1-dimension(rs),dimension(ds)-1)  # not banded
    end

    InterlaceOperator(ops,ds,rs,
                        cache(interlacer(ds)),
                        cache(interlacer(rs)),
                        (l,u))
end

function InterlaceOperator{T}(ops::Vector{Operator{T}},ds::Space,rs::Space)
    # calculate bandinds
    p=size(ops,1)
    if all(isbanded,ops)
        l,u = 0,0
        for k=1:p
            l=min(l,p*bandinds(ops[k],1)+1-k)
        end
        for k=1:p
            u=max(u,p*bandinds(ops[k],2)+1-k)
        end
    else
        l,u = (1-dimension(rs),dimension(ds)-1)  # not banded
    end

    InterlaceOperator(ops,ds,rs,
                        InterlaceIterator(tuple(dimension(ds))),
                        cache(interlacer(rs)),
                        (l,u))
end

function InterlaceOperator{T}(opsin::Matrix{Operator{T}})
    isempty(opsin) && throw(ArgumentError("Cannot create InterlaceOperator from empty Matrix"))

    ops=promotespaces(opsin)
    # TODO: make consistent
    # if its a row vector, we assume scalar
    if size(ops,1) == 1
        InterlaceOperator(ops,domainspace(ops),rangespace(ops[1]))
    else
        InterlaceOperator(ops,domainspace(ops),rangespace(ops[:,1]))
    end
end

function InterlaceOperator{T}(opsin::Vector{Operator{T}})
    ops=promotedomainspace(opsin)
    InterlaceOperator(ops,domainspace(first(ops)),rangespace(ops))
end

InterlaceOperator{T,p}(ops::Array{T,p}) =
    InterlaceOperator(Array{Operator{mapreduce(eltype,promote_type,ops)},p}(ops))

function Base.convert{T}(::Type{Operator{T}},S::InterlaceOperator)
    if T == eltype(S)
        S
    else
        ops=Array(Operator{T},size(S.ops)...)
        for j=1:size(S.ops,2),k=1:size(S.ops,1)
            ops[k,j]=S.ops[k,j]
        end
        InterlaceOperator(ops,domainspace(S),rangespace(S),
                            S.domaininterlacer,S.rangeinterlacer,S.bandinds)
    end
end

#TODO: More efficient to save bandinds
bandinds(M::InterlaceOperator) = M.bandinds

# function colstop{T}(M::InterlaceOperator{T,2},j::Integer)
#     l = M.bandinds[1]
#
#     if isinf(l)
#         (J,ξ) = M.domaininterlacer[j]
#         ret = j
#
#         for K=1:size(M.ops,1)
#             cs=colstop(M.ops[K,J],ξ)
#
#             for k=j:size(M,1)
#                 (KK,κ)=M.rangeinterlacer[k]
#                 if K == KK && κ ≥ cs
#                     ret = max(ret,k)
#                     break
#                 end
#             end
#         end
#         return ret
#     else
#         return  min(size(M,1),j-l)
#     end
# end

# function colstop{T}(M::InterlaceOperator{T,1},j::Integer)
#     l = M.bandinds[1]
#
#     if isinf(l)
#         ret = j
#
#         for K=1:length(M.ops)
#             cs=colstop(M.ops[K],j)
#
#             for k=j:size(M,1)
#                 (KK,κ)=M.rangeinterlacer[k]
#                 if K == KK && κ ≥ cs
#                     ret = max(ret,k)
#                     break
#                 end
#             end
#         end
#         return ret
#     else
#         return  min(size(M,1),j-l)
#     end
# end

israggedbelow(M::InterlaceOperator) = all(isbandedbelow,M.ops)

function getindex{T}(op::InterlaceOperator{T,2},k::Integer,j::Integer)
    M,J = op.domaininterlacer[j]
    N,K = op.rangeinterlacer[k]
    op.ops[N,M][K,J]::T
end

# the domain is not interlaced
function getindex{T}(op::InterlaceOperator{T,1},k::Integer,j::Integer)
    N,K = op.rangeinterlacer[k]
    op.ops[N][K,j]::T
end

function getindex{T}(op::InterlaceOperator{T},k::Integer)
    if size(op,1) == 1
        op[1,k]
    elseif size(op,2) == 1
        op[k,1]
    else
        error("Only implemented for row/column operators.")
    end
end

#####
# optimized copy routine for when there is a single domainspace
# and no interlacing of the columns is necessary
# this is especially important for \
######

function Base.convert{SS,PS,DI,RI,BI,T}(::Type{Matrix},
                            S::SubOperator{T,InterlaceOperator{T,1,SS,PS,DI,RI,BI}})
    kr,jr=parentindexes(S)
    P=parent(S)
    ret=Array(eltype(S),size(S,1),size(S,2))
    for k in eachindex(kr)
        K,κ=P.rangeinterlacer[kr[k]]
        @inbounds ret[k,:]=P.ops[K][κ,jr]
    end
    ret
end

function Base.convert{SS,PS,DI,RI,BI,T}(::Type{BandedMatrix},
                            S::SubOperator{T,InterlaceOperator{T,1,SS,PS,DI,RI,BI}})
    kr,jr=parentindexes(S)
    P=parent(S)
    ret=BandedMatrix(eltype(S),size(S,1),size(S,2),bandwidth(S,1),bandwidth(S,2))
    for j=1:size(ret,2),k=colrange(ret,j)
        K,κ=P.rangeinterlacer[kr[k]]
        @inbounds ret[k,j]=P.ops[K][κ,jr[j]]
    end
    ret
end

domainspace(IO::InterlaceOperator) = IO.domainspace
rangespace(IO::InterlaceOperator) = IO.rangespace

#tests whether an operator can be made into a column
iscolop(op) = isconstop(op)
iscolop(::Multiplication) = true

promotedomainspace{T}(A::InterlaceOperator{T,1},sp::Space) =
    InterlaceOperator(map(op->promotedomainspace(op,sp),A.ops))

interlace{T<:Operator}(A::Array{T}) = length(A)==1?A[1]:InterlaceOperator(A)

immutable DiagonalInterlaceOperator{OPS,DS,RS,T<:Number} <: Operator{T}
    ops::OPS
    domainspace::DS
    rangespace::RS
end

function DiagonalInterlaceOperator(v::Tuple,ds::Space,rs::Space)
    T=mapreduce(eltype,promote_type,v)
    w=map(Operator{T},v)
    DiagonalInterlaceOperator{typeof(w),typeof(ds),typeof(rs),T}(w,ds,rs)
end
DiagonalInterlaceOperator{ST<:Space}(v::Tuple,::Type{ST}) =
    DiagonalInterlaceOperator(v,ST(map(domainspace,v)),ST(map(rangespace,v)))
DiagonalInterlaceOperator(v::Vector,k...) = DiagonalInterlaceOperator(tuple(v...),k...)

Base.convert{T}(::Type{Operator{T}},op::DiagonalInterlaceOperator) =
        DiagonalInterlaceOperator(map(Operator{T},op.ops),op.domainspace,op.rangespace)

function bandinds(S::DiagonalInterlaceOperator)
    binds=map(bandinds,S.ops)
    bra=mapreduce(first,min,binds)
    brb=mapreduce(last,max,binds)
    n=length(S.ops)
    n*bra,n*brb
end

function getindex(D::DiagonalInterlaceOperator,k::Integer,j::Integer)
    n=length(D.ops)
    mk = n+mod(k,-n)
    T=eltype(D)
    if mk == n+mod(j,-n)  # same block
        k=(k-1)÷n+1  # map k and j to block coordinates
        j=(j-1)÷n+1
        D.ops[mk][k,j]::T
    else
        zero(T)
    end
end

domainspace(D::DiagonalInterlaceOperator) = D.domainspace
rangespace(D::DiagonalInterlaceOperator) = D.rangespace

## Convert Matrix operator to operators

Base.convert{OO<:Operator}(::Type{Operator},M::Array{OO}) = InterlaceOperator(M)

function choosedomainspace{T}(A::InterlaceOperator{T,1},sp::UnsetSpace)
    # this ensures correct dispatch for unino
    sps = Vector{Space}(
        filter(x->!isambiguous(x),map(choosedomainspace,A.ops)))
    if isempty(sps)
        UnsetSpace()
    else
        union(sps...)
    end
end