Changed the max chemical evolution time

to 1000 times the orbital time.

Now the binding eenergy of surface species will depends on the H2 coverage on
the surface, following Garrod 2011.
Bug corrected in column making

Should not use id==0 to decide whether the cell should be indexed.
Should use using==.true.