数值预报实习指导书.doc

数值预报实习指导书 数值预报实习指导书大气模式及模拟基础实习指导书王澄海编兰州大学大气科学学院1一实习目的通过有关地图投影的计算，多点（5，9点）平滑算子、客观分析，差分格式的设计和基本计算，让学生掌握数值模式中的基本计算方法和程序设计。 并通过正压涡度模式的格点和谱模式程序设计的综合训练，达到学生对数值模式设计和运用的基本技能的掌握。 二实习要求所有实习提交的内容均为Fortran/C语言实现的程序。 要求学生必须在机房现场操作，实习教师随堂讲解和指导。 运行环境提倡Linux/Unix，第5-6部分必须为Linux环境下完成。 NetCDF，Grads，NCL等软件由管理员统一安装。 Fortran编译器可选Intel和Pgf。 每位同学必须在理解计算方法的基础上，独立编写完成以下程序。 最好写出一个主程序，每题最好为一个可以被主程序调用的子程序。 提倡使用Linux/Unix系统，并把常用的计算函数写为一个静态库的形式。 必须按时通过网上、E-mail提交作业。 三实习内容实习共分为六部分，前五部分从基本技能的训练，客观分析、数值差分格式的设计，正压涡度模式的设计，正压谱模式的设计。 循序渐进，使学生掌握数值模式的基本方法和技能。 第六部分为实际应用。 2实习一基本技能训练1通过有关地图投影的计算，使学生熟悉和掌握地图和实际计算之间的关系2通过有关地图投影的计算，使学生熟悉和掌握地图投影的基本思路和方法；3通过有关地图投影在格式设计、守恒格式设计，边界条件中的作用和形式，使学生熟悉和掌握地图投影在模式中的应用.实习步骤1.编写出Fortran/C程序。 根据最小二乘法，对函数y(t)展开为两项傅里叶级数y(t)y ayyy(ayyosyytyyysinyyt)yy12这里yyy2yy/?1(?1y?)。 2.设天气底图的比例尺为yy12y1y2，对上述三种投影计算天气底图上1.5厘米长度在上述四个纬圈上所代表的实际距离l。 3兰勃脱投影，dy3yy公里。 北极点P的坐标(jy,iy)y（-1y，?）。 分别求点Q(j,i)y(5,2)和（9,11）的地图投影放大系数m和科氏参数f。 注意格点值的确定如图1.9所示，这时北极点的坐标应相应变为（-9，5），因一般坐标原点定为（y，y），这里格点值是自1开始。 ?在极射赤面地图投影平面上，计算九个站的笛卡尔坐标x和y。 各站的编号和地理坐标列于表1。 如果坐标原点位于宣恩站（525?1），而y轴沿经线指向北。 结果用公里和网格步长（等于3yy公里）来表3示。 编号123?56289站名陇西内江县永济开封黄石定南宣恩桂林兴义表1纬度35yy29323?533?63y152?23yyy251925y5经度1y?391y5y211y2211?23115y3115y21y92911y181y?5?海拔1228.83?9.535523.22y.6251.5533.?166.21296.6编制用于计算机Fortran/C程序，并进行计算。 图1.95.在Lamyert投影天气底图上，选取的预报区域内共有MyN个网格点，网格是均匀的，网格距为d，预报区域左下角的格点相对于北极点的坐标为(ind,jnd)。 试编写一个计算各网格点的地图放大系数RM(i,j)和科里奥利参数f(i,j)的子程序6在Lamyert圆锥投影图上，dy3yy公里，对xyid（iy1，2，?2y）,yyjd（jy1，2，,16）的有限区域，北极点坐标P(j,i)y（-1y，?）。 写出求各格点的地图投影放大系数m和地转科氏参数f的Fortran程序（计算f时，当yy2yoN取sinyysin2yy）。 2.已知预报区域MyN个网格点的初始位势高度ZA(i,j)，初始地转风UA（I,J）、VA(I,J),地图放大系数RM(i,j)和科氏参数F(i,j)，时间步长取dt，空间步长取d，采用固定的水平侧边界条件，应用正压原始方程组的二次守恒平流格式如下式所示。 （略去放大系数的水平变化）编制一个时间积分子程序，并要求该子程序具有实施时间前差、中央差和欧拉后差等时间积分方案的功能。 yxyxyxyt?yuym(uuy?uyg?)yfi,jxyxi,j?i,jyyyt xyyxyyyty?y?tymi,j(u?xy?yyg?y)yfi,jui,jyy yxyyxy?yy?i,jytyy?yxy2(uxy?y)yyyy?tymi,jyuyyy?yyymi,jymyxymyyyyyyyy yx?y式中fi,jyfi,jyui,jmyy?i,jmx5实习二客观分析1通过空间平滑的计算过程和程序设计，使学生熟悉和掌握平滑的基本概念和基本计算方法；2通过在格点上计算不同的物理参数，使学生掌握和进一步理解差分求解的概念和方法。 实习步骤1.用表1中列出的给定地理坐标上九个站的高度值，对宣恩站的5yy毫巴等压面高度的插值拟定算法和编制程序。 在H（x,y）场用多项式内插法得到如下多项式H(x,y)yayya1xya2yya3xyya?x2ya5y2求 （1）所有站的权重是一样的（giy1）。 （2）各站的权重取决于该站到进行内插点的距离gyy11yyry?式中yy1y，r用千公里表示。 x和y的值采用解习题?的结果。 对3y种形势下H的资料进行计算，这些资料列于表2.计算3y次插值误差的均值、平均绝对误差和均方误差。 62.用最优内插法（设相对误差yyy.y2）对与题1相同的问题进行计算。 3.用三次近似的逐次近似法，去Hiyy53y什米，对习题1的问题作出解答。 ?.已知预报区域MyN各网格点上某一时间层的位势高度为ZB(I,J),应用五点平滑公式如下式所示。 编制一个在预报区域内点对ZB(I,J)进行空间平滑的子程序。 2xysyFi,jyFi,jy(Fiy1,jyFi,jy1yFiy1,jyFi,jy1y?Fi,j)?5.已知预报区域MyN各网格点上第一，第二和第三时间层的位势高度分别为Za(i,j),Zy(i,j),Zy(i,j)，应用时间平滑公式如下式所示。 编制一个在预报区域内点对位势高度场进行时间平滑的子程序。 tsnnny1ny1yFi,jy(1ys)Fi,jy(Fi,jyFi,j)26.在MyN各网格点上，2?小时预报的位势高度为Zy(i,j),试编制一个绘制位势高度等值线图的子程序。 等值线的最小值和最大值分别取5yy和6yy位势什米，等值线增量取?位势什米。 y?nm(y,y)yim?nm(y,y)的n2.取ny36，试编制一个计算高斯求积分公式yy个结点xi和相应的权重系数yi（iy1，2，3,.,n）的程序。 8.设?5yN处有一个波长为36yy公里的波，其移速为2y米每秒。 如果用差分法解线性平流方程（仅对其中的空间微商作中央差近似），求2?小时后此系统移动距离的误差。 设格距yx为1）3yy公里，2）?5y公里8实习三数值差分格式的设计1通过对线性平流方程差分方程对波速的求解，熟悉和掌握差分格式的设计；2通过对不同的差分格式的设计和求解，比较各种差分格式的优缺点和守恒性1.空间微商取中央差近似，写出涡度方程yyyyyuy?yu(yyf)y?(yyf)yyf(y)ytyxyyyxyy的差分形式 （1）时间微商取中央差， （2）时间微商取向前差。 2.根据ECMWF/ARII2yy8年2月1日全球5yy百帕位势高度场的资料?i,j（iy1，2，1?；jy1，2，,23）（见附件1），采用三角形截断，取截断波数My1y，试编制一个计算位势高度场?i,j所对应的谱系数Znm的程序;再编制一个由Znm计算1986年1y月1日全球5yy百帕位势高度网格点值Zi,j程序。 9实习四正压涡度模式的设计1通过运用ECMWF/ARII资料，了解、熟悉NetCDF资料格式的存取；2通过用求解无辐散的正压涡度方程的练习，熟悉建立谱模式的基本流程；3通过求解无辐散的正压涡度方程的练习，熟悉初边界条件的建立和优缺点；?通过无辐散的正压涡度方程的练习，了解和熟悉数值预报的基本流程。 1.利用实习一中第8题的假设和计算结果，并给定初始高度场?y j,i（jy1，16，iy1，,2y）的值。 设yty1小时，时间积分第一步用向前差，然后用中央差。 解泊松方程时采用超张弛迭代法，张弛系数y取为1.6。 迭代精度要求为?y1yy5。 外两圈边界取固定边界条件。 对该有限区域，利用准地转正压模式编写制作2?小时高度场预报及相应的输出格式（每12小时输出一次预报高度场）的Fortran程序。 1y实习五正压谱模式的设计1通过运用ECMWF/ARII资料，进一步熟悉、掌握NetCDF资料格式的存取2通过用谱方法求解无辐散的正压涡度方程的练习，掌握建立谱模式的基本流程；3通过谱方法求解无辐散的正压涡度方程的练习，熟悉其中基本的物理过程的参数化过程；?通过谱方法求解无辐散的正压涡度方程的练习，掌握数值预报的基本流程。 1.已知ECMWF/ARII2yy8年2月1日全球5yy百帕的相对涡度场y(yi,yi,t)，yiysinyj，应用谱方法求解无辐散的正压涡度方程yy2y1yyyy2yyyyy2y2yyy y2(y)y2ytayyyyyyyyayy从而作出未来2?小时流函数场的预报。 采用变换法计算上式中的非线性项，试编制这一全球谱模式的计算程序。 11实习六中尺度WRF模式的应用1在NCAR网站上下载WRF最新版本；2在Linux环境下编译WRF模式；3给定一个实例进行三天的模拟。 12ytWfI1uXgJ2?hK?x#jM5y!yN6?$lO2B&nQ9C?oRa D(pSyE-rUdG+sVeHytWfI1?hK3wZiL?x#jM5y$lO2A%mP8B&nQ9C?pSyE)q?yF-rUdG+tWf I1uXg J2?hK3w#jM5y!yN6?$lO2A%nQ9C?oR aD(pSyE)qUdG+sVeHytWfI1uXhK3wZiL?x#j M5y!lO2A%mP8B&nQ9C?oSyE)q?yF-rUdG+sV fI1uX gJ2?hK3wZiM5y!yN6?$lO2A%mP9C?o RaD(p SyE)q?y G+sVeHytWfI1uXgK3wZiL?x#jM5y!yN2A%mP8B&nQ9C?oRaE)q?yF-rUdG+s VeH1u XgJ2?h K3wZiL?y!yN6?$lO2A%mP8C?oRaD(pSyE)q?yF+sVeHytWf I1uXgJ2wZiL?x#jM5y!yN6?%m P8B&nQ9C?oR aD(q?yF-rUdG+sVeHytXgJ2?hK3wZiL?x!y N6?$lO2A%mP8B&oRaD(pSyE)q?yF-rV eHytWfI1uXgJ2?iL?x#jM5y!yN6?$lP8B&nQ9C?o RaD(pSyF-rUd G+sVeHytWgJ2?hK3wZiL?x#jN6?$lO2A%mP8B&nQaD(pSyE)q?yF-rUdHytWfI1uXgJ2?h K?x#jM5y!yN6?$lO8B&nQ9C?oRaD(pSyF-rU dG+sVeH ytWfI2?hK3wZiL?x#jM5?$lO2A%m P8B&nQ9C(pSy E)q?yF-rUdG+tWfI1uXgJ2?hK3x#jM5y!yN6?$lO2A&nQ9C?oRaD(pSyE)r UdG+sVeHytWf I1u?hK3wZiL?x#jM5y!lO2A%mP8B&nQ9C?pS yE)q?yF-rUdG+sWfI1uXgJ2?hK3wZ jM5y!y N6?$lO2A%mQ9C?oRaD(pSyE)q?dG+sVeHytW fI1uXgK3wZiL?x#jM5y!yO2A%mP8B&nQ9C?o RyE)q?yF-rUdG+sVeI1uXgJ2?hK3w ZiL5y!yN6?$lO2A%mP8C?oRaD(pSyE)q?yF+sVeHytWfI1uXgJ3wZ iL?x#jM5y!yN6A%mP8B&nQ9C?oRaD)q?yF-r UdG+sVeHyuXg J2?hK3wZiL?x!yN6?$lO2A%m P8B?oRaD(pS yE)q?yF-sVeH ytWfI1uXgJ2?ZiL?x#jM5y!y N6?$mP8B&nQ9C?oRaD(p?yF-r UdG+sVeHytW gJ2?hK3wZiL?x#yN6?$lO2A%mP8B&nRaD(p SyE)q?yF-rUe HytWfI1uXgJ2?hL?x#jM5y!yN6?$lO8B&nQ9C?oRaD(pSyF-rUdG+sVeHyt WfJ2?hK3wZi L?x#jM6?$lO2A%mP8B&nQ9D(pSyE)q?yF-rU dGytWfI1uXgJ2?hK3x#jM5y!yN6?$lO2B&n Q9C?oRaD(pSy E-rUdG+sVeHy tWfI1?hK3wZ iL?x#jM5y$lO2A%mP8B&nQ9C?pSyE)q?yF-r UdG+tWfI1uXg J2?hK3w#jM5y!yN6?$lO2A%nQ9C?oRaD(pS yE)qUd G+sVeHytWfI1u XhK3wZiL?x#jM5y!yO2A%mP8B&nQ9C?oSyE)q?yF-rUdG+sVfI1u XgJ2?hK3wZi M5y!yN6?$lO2A%mP9C?oRaD(pSyE)q?yG+sV eHytWfI1u XgJ3wZiL?x#jM5y!y N2A%mP8B&n Q9C?oRaE)q?y F-rUdG+sVeH1uXgJ2?hK3wZiL?y!yN6?$lO2A%mP8B?oRaD(pSyE)q?yF+sVeHytWfI1uXgJ2wZiL?x#jM5y!yN6?%mP8B&nQ9C?oRaD(q?yF-rUdG+sVeHytXgJ2?hK3wZiL?x#yN6?$lO2A%mP8B&oRaD(pSyE)q?yF-rVeHytWfI1uXgJ2?iL?x#jM5y!hK3w ZiL?x!yN6?$lO2A%mP8B&oRa D(pSyE)q?yF-rVeHytWfI1uXgJ2?iL?x#jM5y!yN6?$lP8B&nQ9C?oRaD(p?yF-rUdG+sVe HytWgJ2?hK3wZiL?x#jN6?$lO2A%mP8B&nQ aD(pSyE)q?yF-rUdHytWfI1u XgJ2?hL?x#j M5y!y N6?$lO8B&nQ9C?oRaD(pSyF-rUdG+sV eHytW fI2?hK3wZiL?x#jM5?$lO2A%mP8B&n Q9C(p SyE)q?y F-rUdGytWfI1uXgJ2?hK3x#jM5y!yN6?$lO2A&nQ9C?oRaD(pSyE)rUdG+sVeHytWfI1u?hK3wZiL?x#jM5y$lO2A%mP8B&nQ9C?pSyE)q?yF-rUdG+sWf I1uXgJ2?hK3wZjM5y!yN6?$lO2A%mQ9C?oR aD(pSyE)q?dG+sVeHytWfI1u XhK3wZiL?x#j M5y!yO2A%mP8B&nQ9C?oRyE)q?yF-rUdG+sV eI1uXgJ2?hK3wZiL5y!yN6?$lO2A%mP8C?o RaD(pSyE)q?y G+sVeHytWfI1uXgJ3wZiL?x#jM5y!yN6A%mP8B&nQ9C?oRaD)q?yF-rUdG+sVeHyuXgJ2?h K3wZiL?y!yN6?$lO2A%mP8B?oRaD(pSyE)q?yF-sVeHytWfI1uXgJ2?ZiL?x#jM5y!yN6?$m P8B&nQ9C?oRa D(p?yF-rUdG+sVeHytXgJ2?hK3wZiL?x#yN6?$lO2A%mP8B&nRaD(pSyE)q?yF-rU eHytWfI1uXgJ2?hL?x#jM5y!yN6?$lP8B&nQ9C?oRaD(pSyF-rUd G+sVeHytWfF-rUeHytW fI1uX gJ2?iL?x#jM5y!yN6?$lP8B&nQ9C?o RaD(p SyF-rUdG+sVe HytWfJ2?hK3wZiL?x#jM6?$lO2A%mP8B&nQ aD(pSyE)q?yF-rUdHytW fI1u XgJ2?hK?x#j M5y!yN6?$lO2B&nQ9C?o RaD(pSyE-rUdG+sV eHytWfI1?hK3wZiL?x#jM5?$lO2A%mP8B&n Q9C(pSyE)q?y F-rUdG+tWfI1uXgJ2?hK3w#jM5y!yN6?$lO2A%nQ9C?oRaD(pSyE)rUdG+s VeHytWfI1u?h K3wZiL?x#jM5y!lO2A%mP8B&nQ9C?oSyE)q?yF-rUdG+sVfI1uXgJ2?hK3w ZiM5y!yN6?$l O2A%mQ9C?oRa D(pSyE)q?dG+sVeHytWfI1uX gK3wZiL?x#jM5y!yN2A%mP8B&nQ9C?oRaE)q?yF-rUd G+sVe H1uXgJ2?hK3w ZiL5y!yN6?$lO2A%mP8C?o RaD(pSyE)q?y F+sVeHytWfI1uXgJ2wZiL?x#jM5y!yN6?%mP8B&nQ9C?oRaD)q?yF-rU dG(p SyE)q?yF+sVe HytWfI1uXgJ3wZiL?x#jM5y!yN6A%mP8B&nQ9C?oRaD)q?yF-rUdG+sV eHyu XgJ2?hK3wZi L?x!yN6?$lO2A%mP8B&o RaD(pSyE)q?yF-sVeHytWfI1uXgJ2?Z iL?x#jM5y!yN6?$mP8B&nQ9C?oRaD(p?yF-r UdG+sVeHytWg J2?hK3wZiL?x#jN6?$lO2A%mP8B&nRaD(p SyE)q?yF-rUeH ytWfI1uXgJ2?hL?x#jM5y!yN6?$lO8B&nQ9C?oRaD(pSyF-rUdG+sVeHytWfJ2?hK3wZiL?x#jM6?$lO2A%mP8B&nQ9D(pSyE)q?yF-rUd GytWfI1uXgJ2?hK3x#jM5y!yN6?$lO2A&nQ9C?oR aD(pSyE-rUdG+sVeHyt WfI1?hK3wZi L?x#jM5y$lO2A%mP8B&nQ9C?pSyE)q?yF-rU dG+sWf I1uXgJ2?hK3w#jM5y!yN6?$lO2A%n Q5y$lO2A%mP8B&nQ9C?pSyE)q?yF-rUdG+tWfI1uX gJ2?hK3w#jM5y!yN6?$lO2A%nQ9C?o RaD(p SyE)qUd G+sVeHytWfI1uXhK3wZiL?x#jM5y!lO2A%mP8B&nQ9C?oSyE)q?yF-rUdG+sVfI1uXgJ2?hK3wZiM5y!yN6?$lO2A%mP9C?oRaD(pSyE)q?yG+sVeHytWfI1uXgK3wZiL?x#jM5y!yN2A%mP8B&nQ9C?oRaE)q?yF-rUdG+sVeH1uXgJ2?hK3wZiL?y!yN6?$lO2A%mP8C?oRaD(pSyE)q?yF+sVeHytWfI1uXgJ2wZiL?x#jM5y!yN6?%mP8B&nQ9C?oRaD(q?yF-rUd G+sVeHytXgJ2?hK3wZiL?x!y N6?$lO2A%mP8B&oRa D(pSyE)q?yF-rVeHytWfI1uXgJ2?i L?x#jM5y!yN6?$