短期实习考核.doc

南京信息工程大学 实验（实习）报告 实验（实习）名称 短期考核 日期 12.14 得分 指导教师 邓伟涛 系 大气科学 专业 大气科学（长望实验班） 年级 2012 班次 1 姓名 汤钰 学号 20121373001 一．实习内容 1．给定了某指数的标准化序列，从中挑选了该指数正位相年（指数 ），负位相年（指数 ），绘制该指数正、负位相年500hPa高度场合成差值图（以及t检验），确定该指数并分析它与500hPa高度场的关系； 2．按照《短期气候预测基础》教材当中该指数的定义，将该指数延长至2014年（1960~2014年），绘制该指数与地面气温的相关系数图，并对其进行分析。 二．实习步骤 1. 资料： NCEP/NCAR 1948年1月－2014年8月（800个月）的500hPa位势高度场、地面气温场月平均资料 资料范围为（900S-900N,00-3600E) 网格距为2.50×2.50,纬向格点数为144，经向格点数为73 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 2. 实习方法： 不同样本之间的均值t统计量： 分别代表第i类样本均值和标准差 ， 为第i类样本的样本量 2)不同样本之间的均值t统计量： 正位相年和负位相年的样本数分别为8和10。任意两类之间通过0.01 显著性水平的t值分别为2.92，通过0.05显著性水平的t值都为2.12。 3) 样本数为55，通过0.01显著性水平的相关系数为0.34，通过0.05显著性水平的t值都为0.26。 3. 编写程序 1) PROGRAM kaohe PARAMETER(N=144*73) dimension hh(144,73,800),h(144,73,67),r(144,73) dimension ave1(144,73),ave2(144,73),cz(144,73),jpc(144,73) dimension a1(144,73,8),a2(144,73,10),zave(144,73) dimension s1(144,73),s2(144,73),t(144,73),jp1(144,73),jp2(144,73) dimension air(144,73,800),aave(144,73,67),kave(144,73) real eap(55),a(55),b(55),c(55) real sum,d,ave,s,p,x,y,z open(10,file=f:\kaohe\data\hgt500.grb,form=binary,status=old) open(20,file=f:\kaohe\data\air.grb,form=binary,status=old) cccccccccccccccccc读数据 do k=1,800 do j=1,73 do i=1,144 read(10)hh(i,j,k) read(20)air(i,j,k) enddo;enddo;enddo write(*,*)read data ok ccccccccccccccccc计算年平均 ccccccc一年平均 do it=1,67 do j=1,73 do i=1,144 sum=0 do k=1,3 sum=sum+hh(i,j,(it-1)*12+5+k) aave(i,j,it)=aave(i,j,it)+air(i,j,(it-1)*12+5+k) end do h(i,j,it)=sum/3 aave(i,j,it)=aave(i,j,it)/3 end do;end do;end do ccccccccccc67年平均 do j=1,73 do i=1,144