MATLAB集成C代码-冰盖模型(DEISM): 深度综合冰盖模型

简介：
本项目介绍如何在MATLAB中集成C语言编写的代码，以构建和运行深度综合冰盖模型DEISM，实现跨编程语言的有效建模。 Matlab集成的C代码深度综合冰盖模型（DEISM）包含我在博士学位论文中编写的冰盖模型的代码。该模型用Matlab编写，并实现了在文献[1]中导出的冰流方程，基于文献[2]进行了数值模拟。此模型使用自动微分执行基础阻力的反演，在ISMIP-HOM实验中经过测试以验证其正确性。有关模型的具体信息可参见文献[3]或[4]。 该模型解决了文献[1]提出的深度综合冰流方程，通常称为“混合方程”。此名称源于它结合了浅冰近似（SIA）和浅层近似（SSA），这两种方法用于描述冰盖的流动特性。通过这种公式化方式，可以求解二维网格上的冰盖速度，并在此之后重建垂直速度剖面。 模型在Arakawa-C网格上使用有限差分法进行方程的数值求解。该模型可根据表面速度和地形观测反演基础阻力，且与已实现的所有三个滑动定律（线性、Budd及Schoof）兼容。文献[4]中描述了一种新程序，此程序利用Budd和Schoof定律将冰下水文学纳入到基础阻力的反演过程中。 此外，在文献[4]和[5]中，还介绍了如何将冰流模型与一个独立的冰下水文模型整合在一起。关于获取该水文模型副本的信息，请直接联系相关作者。