资源描述:
《短期气候预测基础实习四》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、一、实习目的掌握短期气候预测因子的分析和选择,加深对夏季降水分布、坏流界常在短期气候预测屮物理机制的认识。二、实习要求>熟悉资料和方法程序(提供部分子程序);>使用GRADS对结果进行绘图输出和文字分析,完成实习报告。>三、实习内容>计算1951-2001年夏季三类年雨型合成图;>计算各类雨型的前期冬季高度场距平合成图,指出可能出现的遥和关型。四、资料介绍A、人气环流场资料:NCEP/NCAR再分析资料1948-2008年(61年)的500百帕刀平均高度场资料;资料范围为(90°S-90°N,0°-360°E);网格距为2.50X2.50;纬向格点数为144
2、,经向格点数为73。资料为GRD格式,资料从南到北、白西向东排列,每月为一个记录,按年逐月排放。B、站点资料:国家气候中心整编的6、7、8刀降水虽资料(时间段:1951~2010年,资料的格式参见readme.txt文件)C、1951〜2005年雨型分类表:ddi100一类雨型010二类雨型001三类雨型五、方法介绍1、降水距平百分率:RPj=乩二尺x100%「RRi某年夏季降水量;R为1971-2000年夏季降水多年平均值2、总体均值的t统计量:t=l“°JTn为样本量丘、S分别代表样本吻值和标准差“°为总体均值一、二、三类雨型的样本数分别为2()、17和
3、18o我们认为简单认为它们通过0.01显著性水平的t值都为2.88,通过0.05显著性水平的t值都为2.10o八、实习步羽>编程计算并绘制1951-2005年夏季三类雨型年合成图;>编程计算并绘制前期冬季500百帕高度场距平合成图,指出可能存在的遥相关型环流特征。A・源程序:1•三类雨型及前期高度场合成;programmainimplicitnonereal,dimension(160)::sumr,aver,sumrp1,sumrp2,sumrp3,ip1,rp2,rp3,sunup,averp,sumrp11,sumrp22,sumrp33,jfangc
4、ha1,jfangcha2Jfangcha3,tI,t2,t3integer::rsty(55,3),i,j,it,imoreal,dimension(160,55)::r6,r7,r8,summer,rprcal,dimcnsion(l44,73,55)::winter,jphreal::a(144,73,12,61)real,dimension(144,73)::aveh,sumh,sumjph1,sumjph2,sumjph3jph1,jph2,jph3,sumjph1l,sumjph22,sumjph33,jfangcha11Jfangcha22,j
5、fangcha33,th1,th2,th3,sumjph,avejphopen(1,file=*g:r1606*)open(2,file='g:rl607*)open(3,file=*g:r1608')open(4,file='g:ddi')open(5,file='g:rp1.txt')open(6,file='g:rp2.txt*)open(7,file='g:rp3.txf)open(8,file='g:tl.txt1)open(9,file='g:t2.txt')open(10,file='g:t3.txt*)open(11,fi
6、le='g:h5004808.grd',form='binary')opcn(l2,filc='g:th1.grd',form='binary')open(13,file='g:th2.grdform='binary')open(14,file='g:th3.grd',form=,binary')open(15,file='g:jph1,grd',form='binary,)open(l6,file='g:jph2.grd',form='binary,)opcn(17,filc='g:jph3.grd,,form='binary,)read(
7、4,*)((rsty(i,j),j=1,3)3=1,55)read(l,*)((r6(i,j),i=l,160),j=l,55)read(2,*)((r7(i,j),i=1,160),j=1,55)read(3,*)((r8(ij),i=l,160),j=l,55)doit=l,61doimo=1,12doj=1,73doi=l,144read(11)a(i,j,imo,it)cnddo;cnddo;cnddo;cnddowrite(*,*)'okl'doit=l,55doj=1,73doi=l,144winter(ijat)=(a(ijJ2,it4-2)+
8、a(i,j,l,it+2)+a(ij,2,it+2)