目录实习一求500hPa高度场气候场、距平场和均方差场 -------------------------------- 31、资料介绍 ------------------------------------------------------------------------------------------------- 32．要求-------------------------------------------------------------------------------------------------------- 33、实习结果 ------------------------------------------------------------------------------------------------- 31）、FORTRAN源程序 ---------------------------------------------------------------------------------- 3(2)、grads文件 --------------------------------------------------------------- 5（3）、实习结果 ------------------------------------------------------------------------------------- 5实习二计算给定数据资料的简单相关系数和自相关系数---------------------------- 71、资料介绍 ------------------------------------------------------------------------------------------------- 72、要求-------------------------------------------------------------------------------------------------------- 73、实习结果 ------------------------------------------------------------------------------------------------- 7（1）、Fortran源程序 ------------------------------------------------------------------------------ 7（2）、程序运行结果：---------------------------------------------------------------------------- 9实习三分析中国夏季降水线性趋势的分布特征---------------------------------------- 91.资料介绍及要求： --------------------------------------------------------------------------------------- 92.实习结果 -------------------------------------------------------------------------------------------------- 9（1）.matlab程序------------------------------------------------------------------------------------- 9（2）.程序运行结果 -------------------------------------------------------------------------------- 10实习四求给定数据的一元线性回归方程------------------------------------------------ 101、资料介绍及要求 --------------------------------------------------------------------------------------- 102、实习结果 ------------------------------------------------------------------------------------------------ 11（1）、MATLAB程序----------------------------------------------------------------------------- 11 （2）、程序运行结果 ------------------------------------------------------------------------------ 11 （3）、结果分析 ------------------------------------------------------------------------------------ 12 实习五对给定的海温数据进行EOF分析----------------------------------------------- 121、资料介绍 ------------------------------------------------------------------------------------------------ 122、要求------------------------------------------------------------------------------------------------------- 133、实习结果： --------------------------------------------------------------------------------------------- 13（1）、FORTRAN源程序 ------------------------------------------------------------------------ 13 （2）空间场和时间序列的ctl文件------------------------------------------------------------- 14 （3）运行结果 --------------------------------------------------------------------------------------- 15 （4）分析---------------------------------------------------------------------------------------------- 16 实习三（附加）计算给定数据的11年滑动平均和累积距平 ---------------------- 161、资料介绍 ------------------------------------------------------------------------------------------------ 162、要求------------------------------------------------------------------------------------------------------- 163、实习结果 ------------------------------------------------------------------------------------------------ 17 实习四（附加）求给定数据的多元线性回归方程------------------------------------ 181、说明------------------------------------------------------------------------------------------------------- 182、要求------------------------------------------------------------------------------------------------------- 183、实习结果： --------------------------------------------------------------------------------------------- 18（1）Matlab源程序--------------------------------------------------------------------------------- 18 （2）运行结果 --------------------------------------------------------------------------------------- 22实习一求500hPa高度场气候场、距平场和均方差场1、资料介绍有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.时段：1982.1~1985.12共48个月。

水平分辨率：2.5*2.5，格点数：37*17。

2．要求编fortran程序，求500hPa高度场的（1）气候场；（2）距平场；（3）均方差场。

并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。

注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.3、实习结果1）、FORTRAN源程序program ex_gradsimplicit noneinteger,parameter::nx=37,ny=17,nz=4,nt=12integer i,j,iz,itreal var(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt)! Opening fileopen(10,file='g:\gradsdata\h500.dat')do iz=1,nzdo it=1,ntread(10,1000)read(10,3000) ((var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo1000 format(2i7)2000 format(37f6.2)3000 format(37f8.1)4000 format(37f7.2)close(10)!Outputopen(16,file='g:\gradsdata\h500.grd',form='binary')do iz=1,nzdo it=1,ntwrite(16) ((var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!Calculating the Climatological Fielddo it=1,ntdo i=1,nxdo j=1,nysum=0do iz=1,nzsum=sum+var(i,j,iz,it)enddocl(i,j,it)=sum/4enddoenddoenddo! Output climate-fileopen(12,file='g:\gradsdata\climate.grd',form='binary') do it=1,ntwrite(12) ((cl(i,j,it),i=1,nx),j=1,ny)enddo!Calculating the Anomalydo iz=1,nzdo it=1,ntdo i=1,nxdo j=1,nyjp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it)enddoenddoenddoenddoopen(13,file='g:\gradsdata\anomaly.grd',form='binary') !Output anomaly-filedo iz=1,nzdo it=1,ntwrite(13) ((jp(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!Calculating the Mean-square Deviationdo it=1,ntdo i=1,nxdo j=1,nysum=0do iz=1,nzsum=sum+(jp(i,j,iz,it))**2enddojfc(i,j,it)=sqrt(sum/4)enddoenddoenddo!Output mean-square deviation-fileopen(14,file='g:\gradsdata\deviation.grd',form='binary')do it=1,ntwrite(14) ((jfc(i,j,it),i=1,nx),j=1,ny)enddoend(2)、grads文件'open g:\gradsdata\***.ctl' (***为所求场对应的ctl文件名)'set lat 0 40''set lon 60 150''set lev 500''enable print g:\gradsdata\***.gmf' （***为所求场名称）i=1while(i<=48(或 12))'set t 'i'd h''print''c'i=i+1endwhile'disable print';（3）、实习结果①、原始场1982年1月1982年7月结果分析：冬季（此处以1月为代表）等高线分布整体平缓，表明高度场分布相对均匀，且北部接近极地位势高度低，赤道地区位势高度高，这与太阳直射点在1月在南半球，使北半球整体的辐射吸收随纬度增加而减小有关。