cryotools / cosipy

Coupled snowpack and ice surface energy and mass balance model in Python
GNU General Public License v3.0
52 stars 30 forks source link

Modification to the heat equation module #52

Open MarcusGastaldello opened 1 year ago

MarcusGastaldello commented 1 year ago

Modification of the heat_equation.py module in order to implement a Neumann condition (prescribed heat flux) on the base of the computational domain (see Issue #24 'Geothermal heat flux' raised by @apfelbuzza). Modified matrix form to allow a different equation for the basal node.

Updated on 02.02.2024 to enable the user to either have a prescribed and fixed basal temperature or a basal heat flux. This is specified by the user in the 'constants.py' file.

(Code modifications contained within the 'Modified-Heat-Equation' branch)