天气学诊断实习程序

实习资料:该过程的观测资料和再分析资料，实习要求用 micaps 资料进行计算和分析。1、观测资料：计算范围：东经 30-160，北纬 10-80，格距 2.5*2.5。所给要素为高度场，风场，温度以及温度露点差，共11 层。时次是 2013 年 5 月 21 日 08时到 2013 年 5 月 31 日 18 时，每 12 小时一次，共 22 个时次。2、再分析资料：东经30-160，北纬 10-80，格距 2.5*2.5。所给要素为高度场，风场，温度以及比湿（单位为：kg/kg），共 17 层，时次是 2013 年 1 月 1 日 00时到 2013 年 5 月 31 日 20 时，每 6 小时一次，共 604 个时次。实习内容：1. 利用所给资料和绘图软件绘制等高线、等压线、等温线，了解当时环流形势（主要是 500 和 850hPa）。2. 计算涡度和散度（主要是 300，500，850hPa），500hPa 实测风涡度平流，计算 850hPa 温度平流并绘图，利用高空和地面形势预报方程，从热力和动力因子分析南方气旋生成、发展的原因，以及散度高低层的配置。3. 计算 850，700，500hPa 水汽通量并绘图，分析水汽通量输送较大的层次及水汽通量散度，指出水汽输送通道大致范围，并说明水汽辐合区域的位置，什么位置有利于降水的发生，原因是什么。比较水汽通量中水汽平流和风的散度项的大小4. 编制计算垂直速度程序，并绘制 500hPa 垂直速度。用第二种修正方案，其中大气层顶的垂直速度可以直接采用 0，也可以用绝热法。5. 以上任务完成后，有兴趣的同学可以编写计算流函数和势函数的程序。也可以对以上计算结果模仿文献绘制一些垂直剖面图等，从不同角度对该次过程进行分析.程序 1 读取保存数据program duqushujureal h(53,29,11,8),t(53,29,11,8),u(53,29,11,8),v(53,29,11,8),ttd(53,29,11,8)integer,parameter::nz=11,nt=8,nx=53,ny=29character z(11)*4,sj(8)*12sj(1)='13052508.000'sj(2)='13052520.000'sj(3)='13052608.000'sj(4)='13052620.000'sj(5)='13052708.000'sj(6)='13052720.000'sj(7)='13052808.000'sj(8)='13052820.000'z(1)='1000'z(2)='925'z(3)='850'z(4)='700'z(5)='500'z(6)='400'z(7)='300'z(8)='250'z(9)='200'z(10)='150'z(11)='100'do iz=1,nzdo it=1,ntopen(1,file='j:\tianzhen\micaps\height\'//trim(z(iz))//'/'//sj(it))do i=1,4read(1,*)enddodo j=ny,1,-1read(1,*) (h(i,j,iz,it),i=1,nx)enddoopen(2,file='j:\tianzhen\micaps\temper\'//trim(z(i...