2: build an archive of Atmospheric models at different T, B, a  (magnetic inclination angle)
By using the matrix I we can associate at every patch of the neutron star surface the frequency dependent emissivity.                   
First compute all models spanning (so far!):
Then interpolate on a common grid and store the 6-D matrix: