AIMOS - Atmospheric Ionization Module OSnabrueck - version 1.6


                 README
                 ======
1. Files included (when downloading)
2. Data format
3. Handle with care-important information
4. Publication/Citation
Appendix A. Changes and information about obsolete versions


1.) README.txt is the file you are reading at the moment.
    grid_*.txt is the grid you chose. The shows the borders of every cell; the columns are: 
                lower geographic latitude border
                upper geographic latitude border
                lower geographic longitude border
                upper geographic longitude border
                The cells are separated by a NEWLINE.
    pressure_*.txt includes the half-level pressure of the cells. 
                This represents the mean pressure in the cell which is NOT the same 
                as the one in the (vertical) middle of the cell!
                The original ionization rates are calculated for a 67-layer resolution, as given from HAMMONIA.
                    original pressure at half level in [Pa] ={
                    100000., 98644.1, 95846.7, 91271.7, 84957.5, 77188.8, 68393.3, 59058.8, 49786.9, 41211.,
                    33758.4, 27490.3, 22246.9, 17886.2, 14282.2, 11324.1, 8914.13, 6966.1, 5403.64, 4159.87,
                    3177.39, 2407.46, 1809.01, 1347.74, 995.278, 728.35, 528.058, 379.176, 269.585, 189.722,
                    132.121,  91.0181, 62.0069, 41.7605, 27.7944, 18.3654, 12.135, 8.01833, 5.29817, 3.50081,
                    2.31319,  1.52846, 1.00994, 0.667326, 0.440941, 0.291355, 0.192515, 0.127206, 0.0840524, 0.0555383,
                    0.0366973, 0.0242481, 0.0160221, 0.0105867, 6.99526e-03, 4.62218e-03, 3.05414e-03, 2.01805e-03, 1.33344e-03, 8.81081e-04,
                    5.82181e-04, 3.84681e-04, 2.52932e-04, 1.62338e-04, 9.80145e-05,  5.26505e-05, 1.70941e-05
                    };   // 67-layers
               If you choose another vertical resolution it will be interpolated from these values.
    *_AIMOS_version.txt contains just the used AIMOS version number.
    *_ionisationrates_year*_doy*.dat these are the ionization rates. For further description read section 2.


2.) Data format
Year and doy (day of year) are given in the filename.

The file itself has 6+"amount of your pressure cells" columns.
Note: If you have one 7+"amount of your pressure cells" columns your pressure file is ending with a newline and therefore not correctly defined.
      You can easily identify the problem by the last two columns in the ionization file which will be identical.
      Solution: Make a new pressure definition without a newline in the end. (Or use the previous ionization files skipping the last column.)

The first 6 columns are indices. In detail:
column 1: interval
          The time is divided into intervals.
          The number of intervals depends on what you chose.
          If you choose 2h, there will be 12 intervals and the first interval
          represents 0h-2h (UTC).

column 2: particle species
          p: proton
          e: electron
          a: alpha

column 3-6: spatial information:
             3: lower geographic latitude border
             4: upper geographic latitude border
             5: lower geographic longitude border
             6: upper geographic longitude border
             Unit: degree. minutes and seconds etc are converted to decimal.


The following columns represent the ionization rates (ion pair production)
for every pressure level. Units: [ion pair production/(m^3 s)]

The pressure levels used are the ones you chose in the pressure definition. The half-level
 pressure for this data set is stored in the pressure_*.txt.

The altitude is not explicitly given here. It varies depending on latitude, season and solar influence.
The according altitude, which was given by Jens Kieser from the HAMMONIA Model (MPI Hamburg), was
used for calculation purpose.

In the case of total column calculation the unis is [ion pair production/m^2 s].

Small note: along with the update to version 1.6 the number digits have been reduced in order to minimize the file size.
None of the erased digits has an impact on data quality.


3.) The ionization data is generated from satellite measurements from NOAA POES N15/N16 for 2002--2005 
    and N17/18 from 2006 onward. The higher energies are taken from NOAA GOES-10, used for 2002 and 2003,
    GOES-11 for 2004-2010, GOES-15 in 2011 and GOES-13 in 2012.
    
    There are some assumptions used in the model and there are also known problems in the satellite data, 
    so please read this carefully.

    low latitudes:
    Particle precipitation and therefore the ionization rate in low latitudes is set to zero, even though
    satellites show count rates here. However, these particles relect a trapped population and will not
    cause ionization - at least not at these latitudes. 

    SAA:
    Problems with the POES satellites are, that they move through the South Atlantic Anomaly and they
    show high count rates there. These count rates are generated by the high particle flux in the ring current 
    which is near earth due to the low magnetic field in the south atlantic. We assume that it is
    only a dip in the field and that most of the particles do not enter the atmosphere there. In fact 
    we do not know how much particles actually enter the atmosphere there, but we do know, that it is only a
    a part of the counts we see in this region. That's why we cut this area out and set the counts to zero.

    get GOES spatially:
    Another problem is that the GOES satellite has a geostationary orbit. Therefore we do not get spatial 
    information from this satellite. We need this satellite because it has high energetic protons and
    alphas. In order to implement the data from this satellite we looked at the high energetic end of the
    proton spectrum from POES, which has spatial information and we saw that these particles occur 
    homogeneously within the polar cap. So we had a look at the mean size of the polar cap of the highest
    POES proton channel and took this as the size of the polar cap for the GOES channels. This is
    an underestimate of the polar cap because the cap grows with increasing particle energy, but it seems
    to be the best to implement the GOES data.

    POES TED:
    According to David Evans from the NOAA the POES TED (the very low energy channels) have a bad shielding
    against MeV electrons. That's why we can see a second polar oval where there should not be one.
    In the case of TED electrons this effect is not severe. The count rates of the TED electrons are
    high, so the crosstalk could not change much. But in case of the TED protons we can see a second oval 
    (due to crosstalk) exactly where we assume the MeV electrons to enter the atmosphere and it shows a comparable
    intensity as the original oval. Since the only channel on POES measuring MeV electrons is contaminated with protons
    and therefore not usable for correction, we had to decide whether to wipe out all TED proton particles below
    a certain latitude or to use them as is. It's more or less the same problem as it is with the SAA. 
    Following we decided to work with the original values. So please remember that the proton ionization rates in 
    very high altitude may be increased in the area of the MeV-electron polar oval.

    energy range:
    the energy range is:
    154eV-500MeV for protons
    154eV-300keV for electrons (154eV-2.5MeV for version 1.2 and 1.7, but data quality suffers significantly, 
                               so the energy range has been reduced in versions 1.6 and 1.8)
                               Thus the lower altitude range of direct electron ionization is at about 70~km. 
    4MeV-500MeV for alphas (per particle and NOT per nuc)
    The upper energy range indicates that galactic cosmic rays are not included and so the ionization rate 
    at ground level is zero almost all the time. The electron ionization in the lower atmosphere is a result
    of bremsstrahlung.

    KNOWN PROBLEMS:
     With the transition from GOES-12 (and before) to GOES 13/15 NOAA noticed that the energy thresholds of the 
     proton channels that were published in their data books have been affected by a transcription fault.
     Instead of publishing data for an isotropic particle distribution, the published values accounting for a 
     distribution that enters the detector by the aperture only. The new energy thresholds (that should also 
     be valid for the previous GOES satellites) are quite different from the old ones for high energetic channel 
     like P7: old (165.--500.MeV), new (110.--900.MeV). 
     Given that GOES protons are widely used for 20 years now (with these old energy thresholds) we would have 
     expected a small note on their website that something "changed", but you cannot find this information online
     yet. Contacting NOAA, it looks like that GOES high energetic proton channel thresholds are subject of 
     _severe_ uncertainties. 
     J. Rodriguez (NOAA, responsible for the particle instrument on GOES) proposes to read:     
     Sandberg et al. (2014, http://onlinelibrary.wiley.com/doi/10.1002/2014GL060469/full) who suggest that the 
     highest proton channel P7 (nominal isotropic energy range 110-900MeV) on GOES-11 has an effective energy 
     range of 120-186 MeV only. As this investigation is based on IMP-8, a satellite that was degradating in 
     orbit for 27 years before the comparison with GOES-11 was possible, I am more than skeptical about exact 
     numbers, but anyway, the energy thresholds of the higher channels are uncertain - and NOAA suggests that 
     they may be significantly higher _OR_ lower than the old values.
     Assuming a typical Power-Law Spectrum with gamma=-4 the new energy thresholds would boost the high energetic
     particle population by almost one (0.9) order of magnitude. The Sandberg energy range would decline this 
     population by 1.5 orders of magnitude, and the NOAA new thresholds and Sandberg would be almost 2.5 oders of 
     magnitude apart.
     Given this uncertainty we did not update the thresholds yet but we are looking for a replacement for these 
     GOES channels.
     To keep a continous data set, years 2011 and 2012 have been processed with the old thresholds as well.
     
    
    Data gaps of Version 1.6:
    The time coverage by GOES and POES satellites is relatively good, but not absolutely without gaps.
     POES:
     In most cases small gaps within one of the POES-satellites can be accumulated by data from the other one.
     Otherwise: In order to get easily usable data we decided to take a mean value for missing 
     POES data (lower energies). This mean value is based on a 8 years (for version 1.2 and before 4 years) average for the same geomagnetic 
     disturbance and the same geographic position, so the mean value is assumed to be a good choice. But data 
     within the polar cap is NOT (can not be) emulated! Ionization in the polar cap will be treated as zero.
     GOES:
     The high energetic protons and alphas measured at GOES have gaps too, but these can not be replaced by mean values.
     When gaps are identified they will be announced here.
     
    Data gaps:
    year: 2001  doy: 60-90 GOES alpha missing
    year: 2002  doy: 29-31 MEPED electron missing
    year: 2002  doy: 90 MEPED electron missing
    year: 2002  doy: 151 MEPED electron missing
    year: 2002  doy: 212 MEPED electron missing
    year: 2002  doy: 243 MEPED electron missing
    year: 2002  doy: 304 MEPED electron missing
    year: 2003  doy: 29-31 MEPED electron missing
    year: 2003  doy: 62 interval: 11-12 (20-24h) TED/MEPED missing - approximated with Kp maps 
    year: 2003  doy: 90 MEPED electron missing
    year: 2003  doy: 117 interval: 9-12 (16-24h) GOES missing
    year: 2003  doy: 118 interval: 1 alpha channel a_a1 (0-2h) missing
    year: 2003  doy: 148-150 GOES Proton missing
    year: 2003  doy: 151 GOES Proton missing and MEPED electron missing
    year: 2003  doy: 212 MEPED electron missing
    year: 2003  doy: 243 MEPED electron missing
    year: 2003  doy: 251 interval: 4 (6-8h) MEPED proton channel mep0P4 missing
    year: 2003  doy: 302 interval: 12 (22-24h) MEPED electron missing
    year: 2003  doy: 304 MEPED electron missing
    year: 2003  doy: 344 interval: 10-12 (18-24h) GOES missing
    year: 2004  doy: 5 GOES missing
    year: 2004  doy: 30-31 MEPED electron missing
    year: 2004  doy: 38 interval: 1-11 (0-22h) TED/MEPED missing - approximated with Kp maps 
    year: 2004  doy: 38 interval: 4 (6-8h) GOES missing
    year: 2004  doy: 91 MEPED electron missing
    year: 2004  doy: 129 interval: 3-6 (4-12h) GOES missing
    year: 2004  doy: 152 MEPED electron missing
    year: 2004  doy: 173 interval: 6-8 (10-16h) GOES missing
    year: 2004  doy: 213 MEPED electron missing
    year: 2004  doy: 231 interval: 10 (18-20h) TED/MEPED missing - approximated with Kp maps 
    year: 2004  doy: 244 MEPED electron missing
    year: 2004  doy: 305 MEPED electron missing
    year: 2005  doy: 29-31 MEPED electron missing
    year: 2005  doy: 90 MEPED electron missing
    year: 2005  doy: 151 MEPED electron missing
    year: 2005  doy: 194 interval: 12 (22-24h) TED/MEPED missing - approximated with Kp maps 
    year: 2005  doy: 212 MEPED electron missing
    year: 2005  doy: 243 MEPED electron missing
    year: 2005  doy: 304 MEPED electron missing
    year: 2005  doy: 319 interval: 12 (22-24h) TED/MEPED missing - approximated with Kp maps 
    year: 2006  doy: 29-31 MEPED electron missing
    year: 2006  doy: 90 MEPED electron missing
    year: 2006  doy: 151 MEPED electron missing
    year: 2006  doy: 212 MEPED electron missing
    year: 2006  doy: 243 MEPED electron missing
    year: 2006  doy: 304 MEPED electron missing
    year: 2007  doy: 29-31 MEPED electron missing
    year: 2007  doy: 33 interval: 2-7 (2-14h) GOES missing
    year: 2007  doy: 90 MEPED electron missing
    year: 2007  doy: 145 interval: 2-12 (2-24h) TED/MEPED missing - approximated with Kp maps 
    year: 2007  doy: 151 MEPED electron missing
    year: 2007  doy: 212 MEPED electron missing
    year: 2007  doy: 243 MEPED electron missing
    year: 2007  doy: 304 MEPED electron missing
    year: 2008  doy: 30-31 MEPED electron missing
    year: 2008  doy: 91 MEPED electron missing
    year: 2008  doy: 152 MEPED electron missing
    year: 2008  doy: 183 TED/MEPED missing - full day approximated with Kp maps 
    year: 2008  doy: 220 TED/MEPED missing - full day approximated with Kp maps 
    year: 2008  doy: 244 MEPED electron missing
    year: 2008  doy: 213 MEPED electron missing
    year: 2008  doy: 276 interval: 1-6 (0-12h) TED/MEPED missing - approximated with Kp maps 
    year: 2008  doy: 305 MEPED electron missing
    year: 2009  doy: 29-31 MEPED electron missing
    year: 2009  doy: 90 MEPED electron missing
    year: 2009  doy: 151 MEPED electron missing
    year: 2009  doy: 212 MEPED electron missing
    year: 2009  doy: 243 MEPED electron missing
    year: 2009  doy: 304 MEPED electron missing
    year: 2010  doy: 91-120 GOES proton missing (no SPE during this period)
    year: 2010  doy: 130 interval: 11-12 (20-24h) GOES proton and alpha missing
    year: 2010  doy: 274-304 GOES proton missing (no SPE during this period)
    year: 2011  doy: 226 interval: 6 (10-12h) TED/MEPED missing - approximated with Kp maps
    year: 2012  doy: 60-67 TED/MEPED missing - full day approximated with Kp maps 
    
   MEPED electron missing:
     Electron flux is corrected based on Janet Green. 
     In case that the correction results in negative flux 
     or in case that high energy proton channel P7 shows counts >2, 
     MEPED electrons are neglected.
     Results in a limitation of the covered altitude range (low altitude electron ionization missing).

   TED/MEPED missing:
     All electron and low/mid proton data is missing for the whole day (or for indicated intervals).
     If possible approximated values from mean Kp maps will be used instead (and indicated here).
     Treat with care.
                   
     e.g. interval: 4 (6-8h) mep0P4 missing: 
        single channel missing
        no approximation with Kp-maps. 
        However, as long as the channel is not the lowest or highest that should not 
        be a major problem as the spectrum will be derived from the remaining channels.

   GOES missing: 
     all alpha channels and high energetic protons (6.9-500MeV) are missing. 
     Values cannot be approximated, therefore the energy range will be reduced.

   GOES proton missing: 
     high energetic protons (6.9-500MeV) are missing. 
     Values cannot be approximated, therefore the energy range will be reduced.

     
4. Publication/Citation:
   The AIMOS data are freely available for everyone. However, it should be cited properly by using the model description (of the first version 1.0):
     Wissing, J.M., and M.-B. Kallenrode, 
     Atmospheric Ionization Module OSnabrueck (AIMOS) : A 3-D model to determine atmospheric ionization by energetic charged particles from different populations 
     J. Geophys. Res., 114, A06104, doi:10.1029/2008JA013884
   Additionally a link on http://aimos.physik.uos.de is highly appreciated.
   Please send a copy of the published paper to jan.maik.wissing@uni-osnabrueck.de to let us know about the impact of the AIMOS project. 

   It might also be useful to contact the AIMOS-team before publication in case that you have questions.
   In case that intensive correspondance with the AIMOS-team is required that might be awarded with co-authorship as well.


   
Appendix A (changes in version):
   -version 1.0: First version. Documentation as printed in JGR:
      Wissing, J.M., and M.-B. Kallenrode: Atmospheric Ionization Modeule OSnabrück (AIMOS) 1: 
      A 3D model to determine atmospheric ionization by energetic charged particles from different populations, 
      J. Geophys. Res. doi:10.1029/2008JA013884

================
   -version 1.1 (April 9th 2009): 
     first: recalculation of the iosotropic pitch angle resolution with much higher resolution.
       expected benefit: Spectre with a strong population of high energetic particles (gamma=-1) were underestimated in version 1.0 by 10% in the 
       lower atmosphere and overestimated by 20% in the very high atmosphere. This will correct it.
     second: Mean precipitation maps in version 1.0 were leveled by measurement from the complete N-15 and N-16 orbits. This forces a strong 
       impact of the equatorial region on the polar oval population. Version 1.1 now cuts out the equatorial region so that only the polar oval 
       region measurements will be used for the comparison and adaption to the ionization maps. This is a MAJOR(!) improvement to the ionization 
       data! Heavily affected is high energetic proton ionization at geomagnetically quiet time.

================
   -version 1.2 (Feb. 17th 2011, outdated but still available as it has been used in multi-model intercomparisons)
     first: The subdivision of the Kp-index has been revised from 19 into 17 subdivisions. The whole data set has been recalculated, but only 
       minor effects are expected. The spatial distribution of ionization at Kp>=6 (e.g. solar events) will benefit from the revision. No changes 
       (less that 0.5%) for the other Kp-values. 
     second: The data base has been extended to almost present (2002-01-01 -- Nov. 2010). This also allows modeling of the recent solar minimum. 
     third: Starting with 2006-01-01 we now use the NOAA POES satellites N17 and N18. To sum it up we use N15/N16 for 2002--2005 and N17/18 for 
       2006--2010. Apart from that GOES-10 is used in 2001--2003 while GOES-11 is used from 2004 to present. The change of satellites was necessary 
       as old detectors show degradation.

     KNOWN PROBLEMS:
     The electron ionization rate between 50 and 70km altitude (and during SPEs also above) is overestimated!
     The reasons are: 
     (a) that the upper energy threshold of the POES MEPED mep0e3 electron detector (2.5MeV) -which 
     has been given in the instrument description- turned out to be mistaken. According to 
     Yando, et al.: "A Monte Carlo simulation of the NOAA POES Medium Energy Proton and Electron Detector instrument", JGR, 2011,
     the upper energy threshold is at least above 10MeV. Consequently the detector will count more particles than expected.
     In order to get rid of this problem I suggest not to use the electron ioniaztion rates below 70km (just cut).
     (b) Especially during periods of high energetic proton flux (e.g. SPEs) these protons cause crosstalk in the electron channels.
     This affects not only the 50-70km altitude range but also higher altitudes.

    Missing data:
     The time coverage by GOES and POES satellites is relatively good, but not absolutely without gaps.
     POES:
     In most cases small gaps within one of the POES-satellites can be accumulated by data from the other one.
     Otherwise: In order to get easily usable data we decided to take a mean value for missing 
     POES data (lower energies). This mean value is based on a 4 years average for the same geomagnetic 
     disturbance and the same geographic position, so the mean value is assumed to be a good choice. But data 
     within the polar cap is NOT (can not be) emulated! Ionization in the polar cap will be treated as zero.
     GOES:
     The high energetic protons and alphas measured at GOES have gaps too, but these can not be replaced by mean values.
     When gaps are identified they will be announced here.
    Known data gaps:
    year   doy     UTC    satellite         gap in channel                                        information
    2001    23     0-24   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2001    46     0-24   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2001    75     0-24   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2001   223     0-24   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2003    62    20-24   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2003   117    16-24   GOES-10           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2003   118     0- 2   GOES-10           low energetic alpha (a_a1)                            calculated without this data
    2003   148     0-24   GOES-10           high energetic protons (6.9-500MeV)                   calculated without this data
    2003   149     0-24   GOES-10           high energetic protons (6.9-500MeV)                   calculated without this data
    2003   150     0-24   GOES-10           high energetic protons (6.9-500MeV)                   calculated without this data
    2003   151     0-24   GOES-10           high energetic protons (6.9-500MeV)                   calculated without this data
    2003   344    18-24   GOES-10           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2004     5     0-24   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2004*   38     0-22   POES (TED/MEPED)  all electron & low/mid energetic protons              values were substituted by assumed values using the kp-index as proxy
    2004*   38     6- 8   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2004   129     4-12   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2004   173    10-16   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2004   231    18-20   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2005    64    10-12   POES (MEPED)      high energetic electrons & medium protons             values were substituted by assumed values using the kp-index as proxy
    2005   194    22-24   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2005   319    22-24   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2006   240    18-20   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2007    33     2-14   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2007*  145     2-24   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2008*  183     0-24   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy         
    2008*  220     0-24   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2008   276     0-14   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2010   100    20-24   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2010   130    20-24   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data
    2010   222    16-20   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2010   229    18-20   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2010   230     0- 2   POES (TED/MEPED)  all electron & low energetic protons                  values were substituted by assumed values using the kp-index as proxy
    2010   274     2- 4   GOES-11           high energetic protons (6.9-500MeV) and all alphas    calculated without this data

*try not to use:
- year 2001 doy 23  (Jan. 23rd, NOAA-16 data not existing) 
- year 2001 doy 46  (Feb. 15th, NOAA-16 data not existing)
- year 2001 doy 75  (Mar. 16th, NOAA-16 data not existing)
- year 2001 doy 223 (Aug. 11th)
- year 2004 doy 38  (Feb. 7th)
- year 2007 doy 145 (Apr. 25th)
- year 2008 doy 183 (July 1st)
- year 2008 doy 220 (Aug. 7th)
If you are looking for the effect of particle events on the lower atmosphere you should also exclude the periods without GOES data - or at least 
keep in mind that some ionization is missing at these periods.

=================
   -versions 1.3--1.5 have not been published

=================
   -version 1.6 (Oct. 5th 2013, default version, released for the ISSI model intercomparison)
     local time sector restructuring: now using MLT instead of LT, resulting in a more structured diurnal variation in the magnetospheric ionization
     mean precipitation maps: to reduce the impact of spikes/outliers the upper and lower 25% of the count rates of each channel 
             are neglected. This superseeds the former smoothing of the  and therefore preserves the extreme steepness (e.g. 4 orders of 
             magnitute when entering the auroral oval) of the spatial distribution.
             Based on 8 years insted of 4.
     spatial resolution: low latitude regions that are not affected by significant particle precipitation are left out in order to get a 
            higher spatial resolution in high latitudes.
     reduce spikes: scaling of the mean precipitation matrices every day is now neglects the upper and lower 5% of the count rates of each channel,
            which reduces measuring faults as well as the impact of e.g. short-term substorms that simply cannot be modeled with AIMOS - but would
            haven an unpredictable affect on the modeled ionization rates. 
     polar cap size: high energetic particles' cap size has been adjusted to the more accurate spatial resolution in high latitudes. 
            In fact this resulted in a slight increase of the polar cap size.
     MEPED electron channels: using the original MEPED channels as published by NOAA, but:
            Since the MEPED electron detector is sensitive to high energetic protons, the electron channels will be set to an error value in case that
            the omnidirectional proton channel P7 shows more than 2 counts (no unit available). This also effectively cuts out the SAA. 

=================
   -version 1.7 (Oct. 4th 2013, unpublished development version, still has issues in the corrected mep0e3 electron channel.)
     local time sector restructuring: now using MLT instead of LT, resulting in a more structured diurnal variation in the magnetospheric ionization
     mean precipitation maps: to reduce the impact of spikes/outliers the upper and lower 25% of the count rates of each channel 
             are neglected. This superseeds the former smoothing of the  and therefore preserves the extreme steepness (e.g. 4 orders of 
             magnitude when entering the auroral oval) of the spatial distribution.
             Based on 8 years insted of 4.
     spatial resolution: low latitude regions that are not affected by significant particle precipitation are left out in order to get a 
            higher spatial resolution in high latitudes.
     reduce spikes: scaling of the mean precipitation matrices every day is now neglects the upper and lower 5% of the count rates of each channel,
            which reduces measuring faults as well as the impact of e.g. short-term substorms that simply cannot be modeled with AIMOS - but would
            haven an unpredictable affect on the modeled ionization rates. 
     polar cap size: high energetic particles' cap size has been adjusted to the more accurate spatial resolution in high latitudes. 
            In fact this resulted in a slight increase of the polar cap size.
     MEPED electron channels: using Janet Green's correction of the MEPED electron channels
            Since Janet Green's corrected may create negative counts, these will be set to an error value.
            Since the MEPED electron detector is sensitive to high energetic protons, the electron channels will be set to an error value in case that
            the omnidirectional proton channel P7 shows more than 2 counts (no unit available). This also effectively cuts out the SAA.

    KNOWN PROBLEMS:
    50-70km electron ionization might be unaccurate (noticed at 69 degrees N), and even be anti-correlated with Kp. 
    This is due to problematic electron measurements in the mep0e3 channel on POES and due to missing information how the 
    electron spectrum will continue at the high energetic end (no measurements available). 
    There is no way to further improve the electron ionization without better instrumentation on new satellites.