资源描述：
《2013气象统计实习报告.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实习一求500hPa高度场气候场、距平场和均方差场一、实习内容有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.时段：1982.1~1985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17。编fortran程序，求500hPa高度场的（1）气候场；（2）距平场；（3）均方差场。并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。二、fortran程序programmainimplicitnoneintegerii,jj,mo,yrparameter(ii=3
2、7,jj=17,mo=12,yr=4)realvar(ii,jj,mo,yr)realclimate(ii,jj,mo),juping(ii,jj,mo,yr),junfang(ii,jj,mo)integeri,j,m,iyopen(5,file='d:1h500.dat')open(1,file='d:1h500aver.grd',form='binary')open(2,file='d:1h500juping.grd',form='binary')open(3,file='d:1h500junfang.grd',form='binary')doiy=1,4!读数
3、据资料h500dom=1,12read(5,1000)read(5,3000)((var(i,j,m,iy),i=1,ii),j=1,jj)EnddoEnddo1000format(2i7)2000FORMAT(37F6.2)3000format(37f8.1)4000format(37f7.2)dom=1,mo!气候场doj=1,jjdoi=1,iiclimate(i,j,m)=0doiy=1,yrclimate(i,j,m)=climate(i,j,m)+var(i,j,m,iy)enddoclimate(i,j,m)=climate(i,j,m)/4enddoenddoendd
4、owrite(1)(((climate(i,j,m),i=1,ii),j=1,jj),m=1,mo)doi=1,ii!距平场doj=1,jjdom=1,modoiy=1,yrjuping(i,j,m,iy)=var(i,j,m,iy)-climate(i,j,m)enddoenddoenddoenddowrite(2)((((juping(i,j,m,iy),i=1,ii),j=1,jj),m=1,mo),iy=1,yr)dom=1,mo!均方差场doi=1,iidoj=1,jjjunfang(i,j,m)=0doiy=1,yrjunfang(i,j,m)=junfang(i,j,m
5、)+juping(i,j,m,iy)**2enddojunfang(i,j,m)=sqrt(junfang(i,j,m)/yr)enddoenddoenddowrite(3)(((junfang(i,j,m),i=1,ii),j=1,jj),m=1,mo)End一、图形显示1、气候场：t=7分析：7月份多年平均场在东部海面有高压中心，西部有低压中心。2、距平场：T=1982.7分析：对于7月份而言，1982年的高度场偏低，大部分地区为负距平。实习二计算给定数据资料的简单相关系数和自相关系数一、实习内容根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单
6、相关系数。2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）二、fortran程序programmainimplicitnoneintegeri,j,maxireal::k=0,sumx=0,sumy=0,xa,ya,sx,sy,r1,y(20),x(20),rx(10),ry(10)real::data1(20)=(/3.4,3.3,3.2,2.9,3.4,2.8,3.6,3.0,2.8,3.0,3.1,3.0,2.9,2.7,3.5,3.2,3.1,2.8,2.9,2.9/)real::data2(20)=(/3.24,3.14,3.26,2.3
7、8,3.32,2.71,2.84,3.94,2.75,1.83,2.8,2.81,2.63,3.2,3.6,3.4,3.07,1.87,2.63,2.47/)doi=1,20sumx=sumx+data1(i)sumy=sumy+data2(i)enddoxa=sumx/20;ya=sumy/20sumx=0;sumy=0doi=1,20sumx=sumx+(data1(i)-xa)**2；sumy=sumy+(data2(i)-ya)**2enddosx
