在matlab中集成C代码实现的深度综合冰盖模型（DEISM）

Matlab集成的c代码深度综合冰盖模型（DEISM） 该存储库包含我在博士学位论文中编写的冰盖模型的代码。 该模型用Matlab编写，实现了在[1]中导出的冰流方程，并基于[2]进行了数值模拟。 该模型使用自动微分执行基础阻力的反演。 该模型已在ISMIP-HOM实验中进行了测试，以验证其实现。 有关模型的完整详细信息，请参见[3]或[4]。 执行 该模型解决了[1]中提出的深度综合冰流方程，通常称为“混合方程”。 该名称源于以下事实：它结合了冰盖方程的浅冰近似（SIA）和浅层近似（SSA）。 通过此公式，您可以求解2D网格上的冰盖速度，然后在此之后重建垂直速度剖面。 在Arakawa-C网格上使用有限差分法求解方程。 基于表面速度和地形的观察，该模型可用于反演基础阻力。 反转过程可与已实现的所有三个滑动定律（线性，Budd和Schoof）兼容。 在[4]中描述了一种新的程序，该程序使用Budd和Schoof定律将基础水文学纳入基础阻力反演中。 在[4]和[5]中，冰流模型与[6]中描述的冰下水文模型整合在一起。 如果您需要水文模型的副本，请与[6]的作者联系。 我很高兴分享将两者结