数值预报试验指导书.doc

《大气模式及模拟基础》 实习指导书 王澄海 编 兰州大学大气科学学院  一．实习目的： 通过有关地图投影的计算，多点（5，9点）平滑算子、客观分析，差分格式的设计和基本计算，让学生掌握数值模式中的基本计算方法和程序设计。并通过正压涡度模式的格点和谱模式程序设计的综合训练，达到学生对数值模式设计和运用的基本技能的掌握。 二． 实习要求 所有实习提交的内容均为Fortran/C语言实现的程序。要求学生必须在机房现场操作，实习教师随堂讲解和指导。运行环境提倡Linux/Unix，第5-6部分必须为Linux环境下完成。 NetCDF，Grads，NCL等软件由管理员统一安装。Fortran编译器可选Intel和Pgf。 每位同学必须在理解计算方法的基础上，独立编写完成以下程序。最好写出一个主程序，每题最好为一个可以被主程序调用的子程序。提倡使用Linux/Unix系统，并把常用的计算函数写为一个静态库的形式。 必须按时通过网上、E-mail提交作业。 三． 实习内容 实习共分为六部分，前五部分从基本技能的训练，客观分析、数值差分格式的设计，正压涡度模式的设计，正压谱模式的设计。循序渐进，使学生掌握数值模式的基本方法和技能。第六部分为实际应用。 实习一：基本技能训练： 通过有关地图投影的计算，使学生熟悉和掌握地图和实际计算之间的关系 通过有关地图投影的计算，使学生熟悉和掌握地图投影的基本思路和方法； 通过有关地图投影在格式设计、守恒格式设计，边界条件中的作用和形式，使学生熟悉和掌握地图投影在模式中的应用. 实习步骤： 1. 编写出Fortran/C程序。根据最小二乘法，对函数y(t)展开为两项傅里叶级数 y(t)= 这里=2k/(=4)。 2. 设天气底图的比例尺为=1：2，对上述三种投影计算天气底图上1.5厘米长度在上述四个纬圈上所代表的实际距离l。 3． 兰勃脱投影，d=300公里。北极点P的坐标=（-10，4）。分别求点Q(j,i)=(5,7)和（9,11）的地图投影放大系数m和科氏参数f。注意格点值的确定如图1.9所示，这时北极点的坐标应相应变为（-9，5），因一般坐标原点定为（0，0），这里格点值是自1开始。 4． 在极射赤面地图投影平面上，计算九个站的笛卡尔坐标x和y。各站的编号和地理坐标列于表1。如果坐标原点位于宣恩站（57541），而y轴沿经线指向北。结果用公里和网格步长（等于300公里）来表示。 表1 编号 站名 纬度 经度 海拔 1 陇西 3500 10439 1728.8 2 内江县 2937 10507 349.5 3 永济 3453 11027 355 4 开封 3446 11423 73.7 5 黄石 3015 11503 20.6 6 定南 2447 11502 251.5 7 宣恩 3000 10929 533.4 8 桂林 2519 11018 166.2 9 兴义 2505 10454 1296.6 编制用于计算机Fortran/C程序，并进行计算。 图1.9 5. 在Lambert投影天气底图上，选取的预报区域内共有MN个网格点，网格是均匀的，网格距为d，预报区域左下角的格点相对于北极点的坐标为。试编写一个计算各网格点的地图放大系数RM(i,j)和科里奥利参数f(i,j)的子程序 6． 在Lambert 圆锥投影图上，d=300公里，对x=id（i=1，2，…20）,y=jd（j=1，2，…,16）的有限区域，北极点坐标P(j,i)=（-10，4）。写出求各格点的地图投影放大系数m和地转科氏参数f的Fortran程序（计算f时，当N取）。 7. 已知预报区域MN个网格点的初始位势高度ZA(i,j)，初始地转风UA（I,J）、VA(I,J),地图放大系数RM(i,j)和科氏参数F(i,j)，时间步长取dt，空间步长取d，采用固定的水平侧边界条件，应用正压原始方程组的二次守恒平流格式如下式所示。（略去放大系数的水平变化）编制一个时间积分子程序，并要求该子程序具有实施时间前差、中央差和欧拉——后差等时间积分方案的功能。 式中：  实习二：客观分析 通过空间平滑的计算过程和程序设计，使学生熟悉和掌握平滑的基本概念和基本计算方法； 通过在格点上计算不同的物理参数，使学生掌握和进一步理解差分求解的概念和方法。 实习步骤： 1. 用表1中列出的给定地理坐标上九个站的高度值，对宣恩站的500 毫巴等压面高度的插值拟定算法和编制程序。在H（x,y）场用多项式内插法得到如下多项式： 求：（1）所有站的权重是一样的（=1）。 （2）各站