实习报告书。
课程名称:气象统计方法课程实践。
姓名:学号:
班级: 级气科班。
\*实习一求500hpa高度场气候场、距平场和均方差场
实习时间:第9周周三节。
1. 资料介绍。
有一500hpa高度场资料,文件名范围:60~150e,0~40n.
时段:1982.1~1985.12共48个月。水平分辨率:2.5*2.5,格点数:37*17。
2.要求。编fortran程序,求500hpa高度场的。
1)气候场;
2)距平场;
3)均方差场。
并能用grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图,图注要清楚,即要注明是哪个时间的图形,并做简单分析。
注:给出了如何用fortran读取ascii码资料。
fortran
program sx1
implicit none
integer nx,ny,mo,yr
parameter(nx=37,ny=17,mo=12,yr=4)
real var(nx,ny,mo,yr)
real at(nx,ny,mo),xd(nx,ny,mo,yr),sx(nx,ny,mo)
integer i,j,m,t,it,iy,irec
open(5,file='d:\study\form\shixione\')
do iy=1,4
do m=1,12
read(5,1000)
read(5,3000)((var(i,j,m,iy),i=1,nx),j=1,ny)
enddoenddo
close(5)
计算气候场at
do t=1,12
do j=1,ny
do i=1,nx
at(i,j,t)=0
do it=1,4
at(i,j,t)=at(i,j,t)+var(i,j,t,it)
enddoat(i,j,t)=at(i,j,t)/4
enddoenddo
enddo求距平场xd
do t=1,12
do j=1,ny
do i=1,nx
xd(i,j,t,1)=0
do it=1,4
xd(i,j,t,it)=var(i,j,t,it)-at(i,j,t)
enddo
enddoenddo
enddo!求均方差场sx
do t=1,12
do j=1,ny
do i=1,nx
sx(i,j,t)=0
do it=1,4
sx(i,j,t)=sx(i,j,t)+(var(i,j,t,it)-at(i,j,t))*2
enddosx(i,j,t)=sqrt(sx(i,j,t)/4)
enddoenddo
enddo写入气候场。
open(10,file='d:\study\form\shixione\',form='unformatted',access='direct',recl=nx*ny)
irec=0
do t=1,12
irec=irec+1
write(10,rec=irec)((at(i,j,t),i=1,nx),j=1,ny)
enddoclose(10)
写入距平场。
open(11,file='d:\study\form\shixione\',form='unformatted',access='direct',recl=nx*ny)
irec=0
do it=1,4
do t=1,12
irec=irec+1
write(11,rec=irec)((xd(i,j,t,it),i=1,nx),j=1,ny)
enddoenddo
close(11)
写入均方差场。
open(12,file='d:\study\form\shixione\',form='unformatted',access='direct',recl=nx*ny)
irec=0
do t=1,12
irec=irec+1
write(12,rec=irec)((sx(i,j,t),i=1,nx),j=1,ny)
enddoclose(12)
1000 format(2i7)
2000 format(37f6.2)
3000 format(37f8.1)
4000 format(37f7.2)
end program sx1
运行结果:grads文件。
气候场。reinit'
enable print d:\study\form\shixione\'
open d:\study\form\shixione\'
set grid off'
set grads off'
set lat 0 40'
set lon 60 150'
set lev 500'
mon=1while(mon<=12)
set t 'mon''
d h'draw title 1982year'mon'month'
print'
c'mon=mon+1
endwhile
disable print'
距平场。reinit'
enable print d:\study\form\shixione\'
open d:\study\form\shixione\'
set grid off'
set grads off'
set lat 0 40'
set lon 60 150'
set lev 500'
year=1982
while(year<=1984)
mon=1while(mon<=12)
set t 'mon''
d h'draw title 500hpa 'year'year'mon'month anomaly'
print'
c'mon=mon+1
endwhile
year=year+1
endwhile
disable print';
均方差。reinit'
2023年气象统计期末试卷
09海科。一 列举题。1.列举三个海洋或大气资料集。2.宽平稳数据时间序列的2个特点。3.多元线性回归的三种形式的矩阵表达式。4.三种滤波方法。5.三种衡量变量间线性相关程度的相关系数。6.m个变量,n次观测的样本的两种不同的矩阵表示形式。7.求非齐次线性方程的解和逆矩阵的两个变换方法。8.度量相似...
2023年统计报告
2013年 经济运行稳中较快。提示 本文原版含图表word版全文 地址附后 正式会员会看到 地址 这里只复制粘贴部分内容或目录 下面显示的字数不代表全文字数 有任何不清楚的烦请咨询本站客服。2013年,面对错综复杂的国内外形势,市委市 深入贯彻落实党的十八大会议精神,坚持 稳中求进 转中求好 能快则...
统计学专业教学实习报告
本文部分内容来自网络,本司不为其真实性负责,如有异议请及时联系,本司将予以删除。一。实习目的。实习是统计学专业教学计划的重要组成部分,是对学生进行实际统计工作能力初步训练的基本形式,是培养学生职业技能与能力的重要环节,是全面检验和提高我校教育教学质量的必要措施。实习的目的是使学生巩固和运用所学的基础...