Following r-2886-5, debugging some issues ( conservation of H2O, water_reservoir not initialized, info_pem which was not working)
Cleaning of some redundant routines (e.g., stopping criterion of the PEM) and variables renamed (e.g., alpha_criterion now transofrmed in ice_criterion)
Ice table subroutine is now in a module and recalled "compute_ice_table_equilibrium" (v.s. dynamic ice table, efforts ongoing)
Abort_pem introduced to avoid just stopping the pem with raw "stop" in the codes
conf_pem has been improved so that values are not hard coded (e.g., geothermal flux, mass of water_reservoir)
Regolith thermal properties updated in a cleaner way