Two changes:

1. A factor 2 is multiplied to the gas-dust collision cooling rate.
2. Use the above cell for initial condition if using the generic initial
   condition failed.
Two major changes

1. The analysing subroutine is improved.
2. A new option to use a different way of implementing H2O self-shielding.