数值天气预报(5-3-5)--第3节斜压原始方程模式的数值解法.pptVIP

* * * 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 §5.3 斜压原始方程模式的数值解法 一、 基本方程组 二、 差分方程 三、 时间积分 四、 水平边界的处理 六、 地形及初值处理 七、 模式的物理过程及其参数化 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 引言 国家气象中心 1983 LFM (Limited area Fine mesh Model） 有限区域细网格模式，主要用于短期降水预报。 1990 LAFS (Limited area Analysis and Forecast System） 1996 HLAFS (High resolution LAFS) 本节将介绍LFM模式： σ坐标系的五层斜压原始方程模式； 水平格距为190.5公里，积分区域为45×37个格点，覆盖中国 及其邻近地区； 每天以世界时00时的资料为初值进行预报，向各台站提供24 及36小时的降水及诊断预报结果。 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 一、基本方程组 LFM采用 σ = p / ps 坐标及兰勃托地图投影下的通量方程组： 摩擦力和水平动量扩散 温度扩散 热源 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 水汽源 （5.151-157) 垂直边界条件 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 地面气压倾向方程： 垂直速度的计算方程： 此外还用到连续方程的两个变形方程： 从 到 垂直积分 从 到任意高度 垂直积分 （5.158) （5.159) 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 二、差分方程 模式在垂直方向上按σ等分为五层。 模式变量定义如图所示。 另外，因高空水汽含量很少，对降水影响不大，只在对流层下半部考虑凝结现象，故只在 k＝6，8，10上考虑水汽方程。 交界面 中间层 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 模式水平取规则的正方形网格，所有变量定义在同一格点上。 差分格式是按照质量守恒、动量守恒和能量守恒原则设计的。 其中： 非线性平流项采用利利(Lilly)格式； 线性项采用中央差格式； 气压梯度力项采用科比(Corby)格式。 为了书写方便，对任一气象要素A，引用以下符号： 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 模式的差分方程组： （5.160) （5.161) （5.162) （5.163) 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 （5.164) （5.165) （5.166) 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 三、时间积分 为了阻尼高频波，保证计算稳定性。时间积分采用欧拉后差和中央差交替使用的格式，时间步长为4.5分钟。 在每6小时内先做5步欧拉后差，其余时间均用中央差方案。 中央差阶段需要进行时间滤波，其滤波公式为 ： 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 LFM模式虽然采用固定边界，但每小时对第二圈格点上的变量进行一次空间平滑，保证了预报效果及计算的稳定性。 为了减小边界上波的反射所造成的影响，后来采用戴维斯(Davies)方法，修改边界附近四圈格点上的计算倾向。 设A为任一气象要素： K为给定的空间函数，由边界向内减小，至第五减小为零。 四、水平边界的处理 修正后的值 预报值 初始时刻的值 第五章 斜压原始方程模式 第3节 斜压原始方程模式的数值解法 大气科学学院闵锦忠 五、地形及初值的处理 LFM模式考虑了地形的作用（平滑地形）。 LFM模式需要的初始资料是1000、700、500、300、100百帕五个等五面的位势高度和850、700、500百帕三个等压面上的温度露点差值。 地形高度上的ps值根据等压面高度值插值求得的； 等σ面上的位势高度及温度露点差用二次拉格朗日插值计算； 温度的初值是由σ面上的位势高度通过静力方程计算得到； 比湿q的初值是由计算出的温度值及温度露点差求得； 风的初值由σ面上的位势高度通过求解平衡方程而得平衡风。 游性适等