- 1、本文档共30页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
2014年气象统计实习报告
?
?
实?习?报?告?书
课程名称:气象统计方法课程实习一 求500hPa高度场气候场、距平场和均方差场
实习时间:第周周、节For给出了如何用fortran读取ASCII码资料h500.dat.
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\h500.dat)
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)
enddo
enddo
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)
enddo
at(i,j,t)=at(i,j,t)/4
enddo
enddo
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
enddo
enddo
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
enddo
sx(i,j,t)=sqrt(sx(i,j,t)/4)
enddo
enddo
enddo
!写入气候场
open(10,file=d:\study\form\shixione\at.grd,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)
enddo
close(10)
!写入距平场
open(11,file=d:\study\form\shixione\xd.grd,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)
enddo
enddo
close(11)
!写入均方差场
open(12,file=d:\study\form\shixione\sx.grd,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)
enddo
close(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\at.gmf
open d:\study\form\shixione\at.ctl
set grid off
set grads off
set lat 0 40
set lon 60 150
set lev 500
mon=1
while(mon=12)
set t mon
d h
draw title 1982yearmonmonth
print
c
mon=mon+1
endwhile
disable print
;
距平场
reinit
enable print d:\study\form\shixione\sx.gmf
open d:\study
您可能关注的文档
- 15届二模秦淮化学试题.doc
- 15材料周建强4150112005.docx
- 15秋地大《高等数学(专下)》在线作业一答案.doc
- 15秋西南交大《混凝土与砌体结构B》在线作业二答案.doc
- 160临床之最.doc
- 15级广西中医药大学成教院《内科护理学》试卷.doc
- 161-关于大悬臂盖梁钢筋间距调整的函.doc
- 168小组工作.doc
- 16号工字钢代替ZLP630吊篮悬臂安装方案(含计算书).docx
- 16春东财《基金管理》在线作业二.doc
- 2022内科主治考试真题及答案5篇.pdf
- 2022年上海造价工程师考试真题卷二.pdf
- 2022小学学生寒假双减作业设计清单方案(详细版).pdf
- 2022年中考道德与法治二轮复习:守望精神家园 专项测试卷(部编版,含答案).pdf
- 2022年襄阳职业技术学院公共课《思想道德基础与法律修养》科目期末试卷A(有答案).pdf
- 2022年新疆师范大学计算机科学与技术专业《数据结构与算法》科目期末试卷A(有答案).pdf
- 2022年湖南冶金职业技术学院公共课《马克思主义基本原理概论》期末试卷A(有答案).pdf
- 2022信息科技课程标准解读心得体会(通用5篇).pdf
- PEP人教小学英语六年级上册单元检测试题附答案(全册).pdf
- 2022年安全生产法及相关法律知识.pdf
文档评论(0)