2023年气象统计实习报告

发布 2022-03-19 19:00:28 阅读 4257

实习报告书。

课程名称:气象统计方法课程实践。

姓名:学号:

班级: 级气科班。

\*实习一求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年,面对错综复杂的国内外形势,市委市 深入贯彻落实党的十八大会议精神,坚持 稳中求进 转中求好 能快则...

统计学专业教学实习报告

本文部分内容来自网络,本司不为其真实性负责,如有异议请及时联系,本司将予以删除。一。实习目的。实习是统计学专业教学计划的重要组成部分,是对学生进行实际统计工作能力初步训练的基本形式,是培养学生职业技能与能力的重要环节,是全面检验和提高我校教育教学质量的必要措施。实习的目的是使学生巩固和运用所学的基础...