Exshall公司

用Turkel-Zwas显式大时间步长FORTRAN程序求解球坐标系中的浅水方程组。本文介绍并编制了一个FORTRAN计算机程序,将Turkel–Zwas显式大时间步长格式应用于半球正压模式,并对浅水方程的积分不变量进行了约束恢复。接着,我们详细介绍了本文中所包含的算法,特别是与T-Z格式的效率和稳定性有关的算法,以及基于变分方法的二次约束恢复方法。我们特别提供了代码中使用的高纬度滤波、夏皮罗滤波和罗伯特滤波算法的详细信息。详细说明了EXSHALL代码中的各种子程序,重点介绍了代码中实现的算法,并给出了一些主要子程序的流程图。最后,我们提供了一个可视化的例子,用真实的初始数据演示了一个4天的运行,以及一个打印出来的样本和高度场和速度场的等值线等值线。