数值天气预报实习一

实习内容（一）

要求：对比分析以下三种空间差分格式的稳定性。

1. 根据空间前差和后差子程序，补充完整中央差格式子程序。

2. 数值试验：

1）使用时间前差格式，空间前差格式：修改相速度进行对比，绘图展示稳定和不稳定的情况，理解该格式稳定性的条件。

2）使用时间前差格式，空间后差格式：修改空间格距或时间步长进行对比，绘图展示稳定和不稳定的情况，理解该格式稳定性的条件。3）使用时间前差格式，空间中央差格式：修改空间格距或时间步长并绘图展示，理解该格式的不稳定性。

1.补充程序

添加存储地址

修改部分，调用前差，后差及中央差格式时修改

修改空间格距和相速度考察其稳定性

1. 使用时间前差格式，空间前差格式

修改相速度进行对比，绘图展示稳定和不稳定的情况，理解该格式稳定性的条件。

理解与分析：使用时间前差格式，空间前差格式，理论得当-1≤β≤0时，前差格式稳定，图中所示的β不在稳定范围内，此时不稳定，修改相速度，差分解的误差随时间增加，实践与理论相互印证。

1. 使用时间前差格式，空间后差格式

修改空间格距或时间步长进行对比，绘图展示稳定和不稳定的情况，理解该格式稳定性的条件。仅修改空间格距。

理解与分析：使用时间前差格式，空间后差格式，理论得当0≤β≤1时，前差格式稳定，图中所示的β在稳定范围内，理论上来说此时应该是稳定的，修改空间步长，由这三张图可知该格式下，随时间增加，差分解的误差逐渐减小，所以该格式稳定，与理论相互印证。

1. 使用时间前差格式，空间中央差格式

修改空间格距或时间步长并绘图展示，理解该格式的不稳定性，仅修改空间格距。

理解与分析：使用时间前差格式，空间中央差格式，理论得该格式绝对不稳定，修改空间步长，由这三张图分析得，差分解得误差随时间增加，不稳定，所以该格式下，绝对不稳定。