实用标准文案
气象统计实习报告
专业:大气科学 班级:xxxx级x班
学号:2012130xxxx 姓名:
*实习一 求500hPa高度场气候场、距平场和均方差场
一、实习结果
文档大全
实用标准文案
1981.1距平场
1981.1 500hpa高度场在欧亚大陆为正距平,在印度洋和太平洋为负距平。
1981.1气候场
1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。
文档大全
实用标准文案
1981.1均方差场
1981.1的均方差场在欧亚大陆的南侧有最小值,在10N~20N的南侧岁纬度而降低,在其北侧随纬度升高。
二、相关的fortran程序,gs文件
(1)Fortran程序:
program main
parameter(nx=37,ny=17,mo=12,yr=4) real var(nx,ny,mo,4)!数据
real vars(nx,ny,mo)!4年气候态 real jp(nx,ny, mo,4)!距平 real fc(nx,ny,mo)!方差 real jfc(nx,ny,mo)!均方差 integer i,j,m,y,irec real::summ=0.0
!*********************求均方差********************* do m=1,12 do i=1,37 do j=1,17 do y=1,4
summ=summ+jp(i,j,m,y)**2 enddo
文档大全
实用标准文案
fc(i,j,m)=summ/4!方差
jfc(i,j,m)=sqrt(fc(i,j,m))!均方差 summ =0.0 enddo enddo enddo
!******************************************
open(11,file='d:\\shixi\\vars.grd',form='binary' ) open(12,file='d:\\shixi\\jp.grd',form='binary' ) open(13,file='d:\\shixi\\jfc.grd',form='binary' ) do m=1,12
write(11)((vars(i,j,m),i=1,nx),j=1,ny) enddo
do m=1,12
write(13)((jfc(i,j,m),i=1,nx),j=1,ny) enddo do y=1,4
do m=1,12
write(12)((jp(i,j,m,y),i=1,nx),j=1,ny) enddo enddo
close(11) close(12) close(13) End
(2)Gs程序
1.距平场
'open e:\\anomaly.ctl' 'set gxout shaded' 'd xiaobao'
'set gxout contour' 'd xiaobao'
'draw title 1982年1月'
'enable print e:\\anomaly.gmf' 'print'
'disable print' 2.气候场
'open e:\\climate.ctl' 'set gxout shaded' 'd xiaobao'
'set gxout contour' 'd xiaobao'
文档大全
实用标准文案
'draw title 1982.1'
'enable print e:\\climate.gmf' 'print'
'disable print' 3.均方差场
'open e:\\deviation.ctl' 'set gxout shaded' 'd xiaolu'
'set gxout contour' 'd xiaolu'
'draw title 1982.1 deviation' 'enable print e:\\deviation.gmf' 'print'
'disable print'
*实习二 计算给定数据资料的简单相关系数和自相关系数
根据下表中年平均气温和冬季平均气温的等级数据进行下列计算: 1)计算两个气温之间的简单相关系数。
2)分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。(滞后长度τ最大取10)
一、实习结果
文档大全
气象统计实习报告材料(2)



