数值天气预报试卷 Lisa 2009.4.11 1 大气科学专业数值天气预报试卷（A） 1、 请写出 坐标的定义？并给出水平气压梯度力在 坐标下的形式？（8 分） 答： （1） 坐标的定义：（4 分） （2） 根据， 气压梯度力在 坐标下的形式： （2 分） 进一步利用状态方程有：（2 分） 2、 大气数值模式可以分为哪几类？（8 分） 答：大气模式的可以分为：斜压原始方程模式正压原始方程模式原始方程模式平衡模式线性平衡模式非地转模式准地转斜压模式准地转正压模式准地转模式过滤模式数值天气预报模式3、 已知地图投影直角坐标系中x和y方向的地图投影放大系数分别为m和n，请写出散度V、标量平流AV 的表达式。 （8 分） 答：Vmn nU X+ mV Y （4 分）；AV mUXA +nVYA （4 分） 4、 请写出A2的差分形式。 答： 2222 2 yA xAA ，2, 1, 1 22)(2sAAAxAjijiji； 21,1, 22)(2sAAAyAjijiji最终有：2, 1, 11,1, ,2 )(4)(sAAAAAAjijijijiji ji（8 分） 5、 对于一维平流方程：0xuctu （其中c为常数） ，时间、空间均采用向前差分格式， 写出对应的差分方程，并讨论其稳定性。 （15 分） 答： （1）差分方程：tuun in i 1 cxuun in i 1=0； （5 分） （2）稳定性讨论：设iIkxnn ieAu ，将其代入差分方程，可以求得增幅因子 G： xksinIxkcosG1 01111122 xkcosG数值天气预报试卷 Lisa 2009.4.11 2 当01 时，该格式为稳定的。 （10 分） 6、 请详细说明用数值方法求解准地转正压模式的具体计算步骤。 答：求解准地转正压模式的具体计算步骤如下： （1）输入初始位势高度场zi j ,0 ； （2）计算放大系数mi j ,和科氏参数fij,； （3）计算绝对涡度i j ,； （4）计算绝对涡度平流Fi j ,； （5）求解预报方程tzji ,。 （8 分） 7、 请写出方程0vFFuFy yxxxt t对应的差分格式。 （10 分） 答：根据差分算子的定义：jijix xjijijixxjijijixxjijixjijixAAdAAAAdAAAAAAAdAAAA, 1, 1, 1, 12, 1, 1,21,21,21,212121241121 jijix xjijijixxjijijixxjijixjijixAAdAAAAdAAAAAAAdAAAA, 1, 1, 1, 12, 1, 1,21,21,21,212121241121 jijix xjijijixxjijijixxjijixjijixAAdAAAAdAAAAAAAdAAAA, 1, 1, 1, 12, 1, 1,21,21,21,212121241121 )(211 ,1 ,n jin jit tFFtF)()(41 , 1, 1, 1, 1n jin jin jin jin jin jin jin jixxxFFuuFFuuxFun jin jiyyvFvFyvF1,1,)()(21 最终有差分格式： )(211 ,1 ,n jin jit tFFtF)()(41 , 1, 1, 1, 1n jin jin jin jin jin jin jin jixxxFFuuFFuuxFun jin jiyyvFvFyvF1,1,)()(21 0 （10 分） 8、 设一维变量Fi可表示为如下谐波形式：iikxcosFCF，其中C为一常数，F、k 和 分别为谐波的振幅、波数和位相，空间三点平滑算子定义为1ii1iix iFF2F2SFF ， （1）请导出空间三点平滑响应函数的表达式； （2）试证明：当 S=1/2 时，三点平滑可以完全滤除原谐波中波长为L=2x 的波动。 （15 分） 答： （1）将变量Fi的谐波形式代入空间平滑公式，整理后可得： ix ikxxkSFCFcoscos11R定义为平滑后的谐波振幅与原来谐波振幅之比，不难得到：LxSLSR2 3sin21,（8 分） （2）当 S=1/2 时，三点平滑的响应函数为：LxcosL,R2321（3 分） 对于原谐波中波长为L=2x，02cos2 ,2123 xxxR ， 因此，此时，三点平滑可以完全滤除原谐波中波长为L=2x的波动。 （4 分） 9、 简述模式初始化的概念模式？简单说明模式初始化有哪些常用的方法？（10 分） 答：（1） 用原始方程模式作预报时， 由于观测或分析资料的误差导致风场和气压场的不平衡，数值天气预报试卷 Lisa 2009.4.11 3 初始资料和数值模式之间的不平衡，直接采用未经任何处理的观测值或分 析值作为模式的初始场，容易导致高频振荡，产生计算不稳定。因此，需 要对模式的初值进行处理，即称为模式的“初始化” （5 分） ； （2）模式 常用的初始化方法包括：静力初始化、动力初始化和变化初始化三种（5 分） 。 10、什么是数值模式物理过程参数化？模式的物理过程参数化包括哪些主要的内容？（10 分） 答： （1）为了提高预报的确率，往往要在数值天气预报模式中引入一些影响大气运动的重要 物理过程，而这些过程基本上属于次网格过程，很难用显式的方法在模式 中表示出来。通常把上述次网格的物理过程，通过大尺度的物理量来描述 它们的统计效应， 并作为某些物理量的源或汇包含在大尺度运动方程组中， 这种方总统称为“参数化”的方法。 （5 分） （2）模式的物理过程参数化的主要内容包括：辐射传输、水汽凝结、边界层的处理（湍流对热量、动量和水汽的输送等） 、地形作用等等（5 分） 。 大气科学专业数值天气预报试卷（B） 1、地图投影放大系数是如何定义的？请列举出并简单说明三种常用的地图投影方式及其基 本特点。 （10 分） 答： （1） 地图放大系数或地图放大因子， 又称为映像比例尺：m=地球表面上相应的距离映像平面上的距离； （4 分） （2）三种常用的地图投影方式及其基本特点： （a）极射赤面投影，在极地和高纬度地区产 生的变形较小，这种投影方式通常用于制作极地天气图和北半球天气底图； （b）兰勃托 投影， 在中纬度地区产生的变形较小， 这种投影方式通常用于制作中纬度地区的天气图， 如亚欧天气底图； （c）而麦卡托投影，在低纬或热带地区产生的变形较小，这种投影方 式通常用于制作低纬或热带地区的天气底图。 （6 分） 2、已知地图投影直角坐标系中x和y方向的地图投影放大系数分别为m和n，请导出 P 坐标系下的连续方程0 pyv xu在地图投影坐标系中的形式。 （10 分） 答： mn nU X+ mV Y+p=0（10 分） 3、请写出)B,A(J的差分形式。 （8 分） 答： jijijijijijijijicBBAABBAA dBAJ, 1, 11,1,1,1, 1, 1241, （8 分） 数值天气预报试卷 Lisa 2009.4.11 4 4、现有如下形式的方程：0yubxuatu，时间微商用向前差分格式、空间微商用中央 差格式来构造对应的差分方程。 （10 分） 2、 答： tuun jin ji ,1 , xuuan jin ji 2, 1, 1yuubn jin ji 21,1,0（10 分） 5、根据准地转涡度方程：0pffVtp ，导出准地转正压模式的预报方程。 （12 分） 答： （1）将准地转涡度方程：0pffVtp ，应用于无辐散层上： tVfp0（4 分） （2）利用准地转近似的散度方程：f 22 gVkfV1；gf21（4 分） （3）旋转风V和相对涡度可以分别由计算的地转风Vg和地转风涡度g来代替 0fVtgg 02yxz xyzzt ，即为预报方程 （4 分） 6、五点平滑和九点平滑是如何定义的？并分别写出二者的计算表达式？（15 分） 答： （1）五点平滑的定义：对变量Fi,j同时沿x方向和y方向进行一维三点平滑运算，然后再取其平均值。其计算表达式为：jijijijijijiyx jiFFFFFSFF,1, 11, 1, ,44（7 分） （2）九点平滑的定义：对变量Fi,j先进行x方向的三点平滑，然后再对其进行y方向的三 点平滑；其计算表达式为： （8 分） jijijijijijijijijijijiyx jiFFFFFSFFFFFSSFF,1, 11, 11, 11, 12,1, 11, 1, ,444127、正压原始方程模式在地图投影坐标系中的预报方程组形式如下： xzmgvfyuvxuumtu* yzmgufyvvxvumtv* yv xu mz mz yvmz xumtz2请采用二次守恒的平流格式，写出其对应的差分方程组。 （15 分） 答：正压原始方程模式采用二次守恒的平流格式所对应的差分方程组为： 数值天气预报试卷 Lisa 2009.4.11 5 000, ,2,* ,* ,y yx xjijiyyyxxxjit tjijiy yyyyxxx jit tjijix xyyyxxx jit tvumzmzumzumzufzgvuvumvvfzguuuumu（15 分） 8、准地转正压模式、正压原始方程模式的时间步长大致是多少，并说明二者存在差别的原 因？（10 分） 答：准地转正压模式的时间步长大致为 1 小时，而正压原始方程模式的时间步长大致为 10 分钟； （5 分） 原因：原始方程模式保留了快速移动的重力惯性波，时间积分的时间步长必须取得很小。 （5 分） 9、简述模式初始化和物理过程参