Rewrote integration to be more transparent with regards to changes. Preparing for completely adaptive mass step.
Removed superflous comment.
Too tiny dm causes rho to go to zero after 1 step.

Currently in process of fixing dynamic mass step. Problem is not in decision of new step but in that smaller dm than -1e24 makes rho go to zero after one step.

* Fixed calculation of rho from EOS - is now explicit variable in integration scheme.
Changed dm, moved main loop i area, added printout of params when something goes to zero.
Moved dm. Now use P_new and T_new to calc rho.
In process of getting adaptive dm to work.

* Fixed loop, is now while so that will always reach center of star
  unless something is zero.
* saveEvery is now deprecated - now use new var massPile to count how
  much mass has been passed and every time a change of 0.0001 of the
  starting mass is reached, a point is written to the datafile.
  This will be conventient for huge simulations.
* Added an upper dm limit, so that if the proposed dm is bigger than
  this, a limit is set. Problem now is that this always happens.

TODO: There is obviously something wrong in the calculation of the new
dm, so that it is set to be way huge. Check this, if the adaptive dm
starts working. The program will be better.

TODO: Remove all traces of the saveEvery var, it is no longer necessary.
Found reason for why calculated dm was wrong.

Used the wrong values in calculation. The values should not be multiplied by dm, this is the change, not the value deciding the step size. Should be correct now. Added a set of vars.
Adaptive dm is working.

Found that what was screwing up everything was error in the formula, and
the fact that the chosen change deciding the dm was max and not min of
all the changes. Should be min so that dm always is the smallest and in
that covering well enough all the vars.

Noticed that if the percentage is lowered from 30% to 20% or 10% the
integration starts hanging at about 14% and dm just keeps dropping. Can
maybe include lower limit to dm to avoid this.

TODO: Remove the saveEvery and N commandline arguments. Not necessary
anymore.
Save ten times more points for more smoothness when plotting.
Cleaning: Away with N and saveEvery.
Removed saveEvery and N from initexperimenter.
Removed saveEvery and N, fixed saving.
Saving to different folder works.
Some newlines in printing.
Extended spectrum with 10% in each direction.
Reading from opacityfile is correct.
Removed print and exit that shouldnt be there.
Algorithm for finding kappa should be correct.
Fixed lots of bugs. Integration now works with new kappa finding method. More effective.
Implemented argparse.