!
      subroutine protons
      implicit none
      integer,parameter,dimension(1) :: ntimes = (/192,
      real,dimension(192,1) :: fe,eo,hour
      integer,dimension(192,1) :: iday
!
!  Proton flux parameters for January 15-22, 2005, Solar Eruptions
      fe(1:192,1)=(/
     | 2.097e+01, 2.655e+01, 2.634e+01, 2.256e+01, 2.514e+01, 2.765e+01,
     | 2.350e+01, 3.386e+01, 1.357e+02, 4.459e+02, 7.942e+02, 9.935e+02,
     | 1.070e+03, 1.022e+03, 9.725e+02, 8.582e+02, 7.748e+02, 7.445e+02,
     | 4.541e+02, 5.568e+02, 5.102e+02, 3.508e+02, 3.269e+02, 3.009e+02,
     | 3.499e+02, 7.030e+02, 1.772e+03, 3.068e+03, 4.388e+03, 4.609e+03,
     | 5.968e+03, 7.123e+03, 7.821e+03, 9.291e+03, 1.239e+04, 1.402e+04,
     | 1.714e+04, 1.818e+04, 2.360e+04, 3.035e+04, 3.113e+04, 3.732e+04,
     | 4.162e+04, 3.948e+04, 3.627e+04, 3.219e+04, 3.602e+04, 3.666e+04,
     | 3.808e+04, 3.645e+04, 4.259e+04, 4.074e+04, 3.164e+04, 2.711e+04,
     | 2.589e+04, 2.739e+04, 2.795e+04, 2.605e+04, 3.593e+04, 3.596e+04,
     | 2.341e+04, 5.159e+04, 1.405e+05, 1.585e+05, 3.038e+05, 3.580e+05,
     | 3.201e+05, 2.513e+05, 2.717e+05, 2.827e+05, 2.761e+05, 2.051e+05,
     | 1.762e+05, 1.914e+05, 2.106e+05, 1.726e+05, 1.788e+05, 1.655e+05,
     | 1.682e+05, 1.624e+05, 1.410e+05, 1.393e+05, 1.348e+05, 1.418e+05,
     | 1.107e+05, 8.431e+04, 8.006e+04, 7.573e+04, 7.053e+04, 6.352e+04,
     | 5.235e+04, 5.451e+04, 4.917e+04, 4.221e+04, 3.412e+04, 2.217e+04,
     | 1.779e+04, 1.835e+04, 1.094e+04, 1.149e+04, 1.176e+04, 1.130e+04,
     | 1.121e+04, 1.075e+04, 9.043e+03, 1.036e+04, 1.056e+04, 1.119e+04,
     | 9.786e+03, 9.211e+03, 8.681e+03, 8.133e+03, 8.490e+03, 7.496e+03,
     | 5.283e+03, 4.902e+03, 5.190e+03, 4.790e+03, 5.197e+03, 4.019e+03,
     | 3.642e+03, 3.301e+03, 3.415e+03, 3.147e+03, 2.509e+03, 2.597e+03,
     | 8.029e+03, 6.580e+05, 4.506e+05, 2.462e+05, 1.679e+05, 1.631e+05,
     | 1.318e+05, 1.155e+05, 1.152e+05, 1.052e+05, 9.456e+04, 7.608e+04,
     | 6.256e+04, 4.637e+04, 3.788e+04, 3.690e+04, 2.879e+04, 2.465e+04,
     | 2.938e+04, 2.082e+04, 1.979e+04, 1.914e+04, 1.959e+04, 2.461e+04,
     | 2.510e+04, 2.601e+04, 2.575e+04, 2.513e+04, 2.054e+04, 2.117e+04,
     | 2.497e+04, 2.636e+04, 2.908e+04, 3.271e+04, 3.078e+04, 5.363e+04,
     | 1.503e+04, 7.755e+03, 6.702e+03, 5.545e+03, 5.097e+03, 4.413e+03,
     | 4.235e+03, 3.919e+03, 3.737e+03, 3.360e+03, 2.783e+03, 2.919e+03,
     | 2.609e+03, 2.448e+03, 2.297e+03, 2.085e+03, 1.899e+03, 1.771e+03,
     | 1.792e+03, 1.744e+03, 1.535e+03, 1.544e+03, 1.353e+03, 1.083e+03,
     | 9.641e+02, 8.945e+02, 7.245e+02, 6.453e+02, 6.738e+02, 6.157e+02,
     | /)
      eo(1:192,1)=(/
     | 1.245e+01, 1.349e+01, 1.324e+01, 1.541e+01, 1.561e+01, 1.345e+01,
     | 1.288e+01, 1.461e+01, 1.225e+01, 9.241e+00, 7.598e+00, 6.615e+00,
     | 6.222e+00, 5.969e+00, 5.877e+00, 5.688e+00, 5.427e+00, 5.366e+00,
     | 6.098e+00, 5.475e+00, 5.714e+00, 6.091e+00, 6.375e+00, 7.210e+00,
     | 6.966e+00, 8.032e+00, 7.727e+00, 7.397e+00, 7.031e+00, 7.069e+00,
     | 7.122e+00, 7.003e+00, 7.145e+00, 6.946e+00, 6.769e+00, 6.543e+00,
     | 6.559e+00, 6.873e+00, 6.921e+00, 6.765e+00, 6.660e+00, 6.185e+00,
     | 5.922e+00, 5.729e+00, 5.687e+00, 5.641e+00, 5.335e+00, 5.247e+00,
     | 5.310e+00, 5.183e+00, 4.971e+00, 4.869e+00, 4.885e+00, 4.916e+00,
     | 4.877e+00, 4.767e+00, 4.707e+00, 4.679e+00, 4.633e+00, 4.609e+00,
     | 5.515e+00, 7.115e+00, 9.217e+00, 9.785e+00, 9.850e+00, 9.718e+00,
     | 9.690e+00, 9.391e+00, 8.897e+00, 8.840e+00, 8.738e+00, 9.041e+00,
     | 9.203e+00, 8.874e+00, 8.580e+00, 8.590e+00, 7.831e+00, 8.002e+00,
     | 7.678e+00, 7.357e+00, 7.259e+00, 7.181e+00, 7.159e+00, 6.797e+00,
     | 6.960e+00, 7.358e+00, 7.610e+00, 8.228e+00, 8.090e+00, 8.016e+00,
     | 8.252e+00, 8.180e+00, 8.371e+00, 8.507e+00, 8.268e+00, 7.871e+00,
     | 7.389e+00, 7.255e+00, 7.024e+00, 6.707e+00, 6.513e+00, 6.478e+00,
     | 6.531e+00, 6.646e+00, 7.120e+00, 7.055e+00, 7.037e+00, 7.191e+00,
     | 7.052e+00, 7.028e+00, 7.084e+00, 7.098e+00, 7.545e+00, 7.836e+00,
     | 8.474e+00, 8.706e+00, 8.657e+00, 9.188e+00, 8.828e+00, 9.536e+00,
     | 9.620e+00, 9.734e+00, 9.564e+00, 9.438e+00, 9.897e+00, 9.405e+00,
     | 1.882e+01, 8.846e+01, 4.535e+01, 2.630e+01, 1.951e+01, 1.568e+01,
     | 1.371e+01, 1.269e+01, 1.067e+01, 9.978e+00, 1.001e+01, 1.006e+01,
     | 1.004e+01, 1.071e+01, 1.129e+01, 1.073e+01, 1.120e+01, 1.177e+01,
     | 9.737e+00, 1.117e+01, 1.122e+01, 1.068e+01, 9.776e+00, 8.290e+00,
     | 7.778e+00, 7.382e+00, 7.182e+00, 7.055e+00, 7.257e+00, 7.110e+00,
     | 6.597e+00, 6.227e+00, 5.934e+00, 5.617e+00, 5.581e+00, 4.665e+00,
     | 5.545e+00, 6.315e+00, 6.341e+00, 6.425e+00, 6.541e+00, 6.562e+00,
     | 6.651e+00, 6.782e+00, 6.760e+00, 6.805e+00, 7.002e+00, 6.795e+00,
     | 6.769e+00, 6.707e+00, 6.716e+00, 6.656e+00, 7.046e+00, 7.051e+00,
     | 7.134e+00, 7.153e+00, 7.579e+00, 7.297e+00, 7.567e+00, 7.747e+00,
     | 8.074e+00, 8.019e+00, 8.446e+00, 8.700e+00, 8.497e+00, 8.480e+00,
     | /)
      iday(1:192,1)=(/
     |   15,  15,  15,  15,  15,  15,  15,  15,  15,  15,  15,  15,
     |   15,  15,  15,  15,  15,  15,  15,  15,  15,  15,  15,  15,
     |   16,  16,  16,  16,  16,  16,  16,  16,  16,  16,  16,  16,
     |   16,  16,  16,  16,  16,  16,  16,  16,  16,  16,  16,  16,
     |   17,  17,  17,  17,  17,  17,  17,  17,  17,  17,  17,  17,
     |   17,  17,  17,  17,  17,  17,  17,  17,  17,  17,  17,  17,
     |   18,  18,  18,  18,  18,  18,  18,  18,  18,  18,  18,  18,
     |   18,  18,  18,  18,  18,  18,  18,  18,  18,  18,  18,  18,
     |   19,  19,  19,  19,  19,  19,  19,  19,  19,  19,  19,  19,
     |   19,  19,  19,  19,  19,  19,  19,  19,  19,  19,  19,  19,
     |   20,  20,  20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
     |   20,  20,  20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
     |   21,  21,  21,  21,  21,  21,  21,  21,  21,  21,  21,  21,
     |   21,  21,  21,  21,  21,  21,  21,  21,  21,  21,  21,  21,
     |   22,  22,  22,  22,  22,  22,  22,  22,  22,  22,  22,  22,
     |   22,  22,  22,  22,  22,  22,  22,  22,  22,  22,  22,  22,
     | /)
      hour(1:192,1)=(/
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     |  0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,10.5,11.5,
     | 12.5,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,
     | /)
      end subroutine protons
