(完整版)数值天气预报考试题答案.docVIP

兰州大学大气科学学院 大气科学 / 应用气象 专业（班） 年级 课程 数值天气预报 200 — 200 学年 第 2 学期（ A ）卷 期末考试题（ 120 分钟） 考试时间 姓 名 学 号 题号 得 分 分 数 主考教师： 隆霄 阅卷教师： ．对线性平流方程： u u 0 tx Crank-Nicolson 的差分格式为： u nj 1 unj a u nj 11 u jn 11 u jn 1 u nj 1 0 t ( 2 x 2 ) 2 x 分析其稳定性条件和截断误差。 (20 分) 解 : 令 a t , 则差分格式可以展开为 : u nj 1 u nj(u nj 11 unj 11 u nj 1 u nj 1 ) 4 x 利用冯 - 纽曼方法 ,设 u nj GeIkxi , I 1, 代入上式 , 化简可以得到增幅因子表达式为 : G 1 2 I sin k x 1 2 I sin k x 由于 | G | | 1 2 I sin k x | 1 1 2 I sin k x 所以格式绝对稳定 . 利用泰勒展开可以得到 : u nj 1 u nj unj t 2u nj ( t) 2 /(2! ) ....... t t 2 所以 u nj u nj 1 u nj 2u nj t /(2!) ....... , 时间截断误差为 O( t ) t t t 2 类似的方法把 unj 11 , unj 11 , u nj 11 , u nj 11在空间域展开 , 可以得到截断误差为 O ( x2 ) ．对线性平流方程： u c u d u 0 t x y 其中 c=2.5, d=0.2, x y 0.1, t 0.004 1 初始条件为： ui0, j 1.5 sin(2 i x) sin(2 j y) 请写出时间中央差、空间中央差分的差分格式，并用 Fortran 语言写出计算程序。 （启动时步采用时间前差格式，东西向边界取周期条件，南北向边界取固定边界条件） 。 要求 :程序中要有明确的变量说明，所写程序符合 Fortran 语言基本标准。 (20 分 ) 答：微分方程的差分格式为 启动时刻的差分格式为  : :  uin, j 1 uin, j 1 uin 1,j uin 1, j uin, j 1 uin, j 1 0 2 t c 2 x d 2 x ui1, j ui0, j c ui0 1, j ui0 1, j d ui0, j 1 ui0, j 1 0 2 t 2 x 2 x 程序略 . 3． 试分析如下差分形式的地转适应方程的解的频散性质。 （20 分） u x x fv g z x 0 其中x t v 0 fu t z x Hu x 0 t 解：设方程有如下形式的波动解  1 ( i 1 i 1 ) 2 x ^ u u eI (kid t ) ^ v veI ( kid t ) ^ zeI ( kid t ) 将此解代入差分方程组，并进行消元得 I u fv gzI sin kd / d 0 I v fu 0 I z HuI sin kd / d 0 进行消元得频率方程 2 2 2 ( ) 1 sin 2kd / 2d 相应的群速度公式为 f 2 2 sin 2kd cg 2d 由此可知当 kd=π/2 时，群速度为 0，波能量不能频散。当 kd=π时，重力惯性波变为纯惯性振荡。 2 4 ．已知 P 坐标系中的连续方程为： ( u v p ) p0 ，其中 x y 坐标系中的连续方程。 (20 分 ) 解：对任意变量 F ( x, y, p,t ) F [ x, y, ( x, y, p,t ), t ] 有： ( F ) P ( F ) F ( ) p (1) s s s 其中 s 表示 x ， y ，或者 t  P Pt P Pt Ps Pt  ，请推导出 又 ( ) p 1 ( p pt [ ( ppt ) ] p ( p pt ) ( s ) s ( 1 ) 把（ 2）代入（ 1）则有： ( F ) P ( F ) F s s s 利用（ 3）把 P 坐标系中的连续方程 ( u v ) p 0 展开有： x y p u v ? Vh 0 （ 4） ( ) x y p 又 dp d ( pt ) （ 5） dt dt 所以： p 1 1 [ ] 1 [ (Vh ? t 把（ 6）代入（ 4），整理后则有： t ? ( Vh ) 0  ) p pt (2) 2 s s s 3） ) ] （6） 7） （ 7）即是 坐标系中的连续方程。 5 ．原始方程模式与过滤模式有哪些不同之处？ (10 分 ) 答 :