实习报告书
课程名称:气象统计方法课程实践
姓名:学号:班级:
级气科班
f实习一求500hPa高度场气候场、距平场和均方差场实习时间:第9周周三1、2节
1资料介绍
有一500hPa高度场资料,文件名h500dat,范围:60150E,040N
时段:19821198512共48个月。水平分辨率:2525,格点
数:3717。
2.要求
编fortra
程序,求500hPa高度场的
(1)气候场;
(2)距平场;
(3)均方差场。
并能用Grads做出图形实习报告中气候场、距平场、均方差场任意
给出两张图,图注要清楚,即要注明是哪个时间的图形,并做简单分
析。
注:h500For给出了如何用fortra
读取ASCII码资料h500dat
FORTRAN
programsx1implicit
o
ei
teger
x
ymoyrparameter
x37
y17mo12yr4realvar
x
ymoyrrealat
x
ymoxd
x
ymoyrsx
x
ymoi
tegerijmtitiyirecope
5filedstudyformshixio
eh500dat
doiy14dom112read51000read53000varijmiyi1
xj1
y
fe
ddoe
ddoclose5计算气候场atdot112
doj1
ydoi1
xatijt0doit14atijtatijtvarijtite
ddoatijtatijt4e
ddo
e
ddoe
ddo求距平场xddot112
doj1
ydoi1
xxdijt10doit14xdijtitvarijtitatijte
ddoe
ddo
e
ddoe
ddo
求均方差场sxdot112
doj1
ydoi1
xsxijt0doit14sxijtsxijtvarijtitatijt2e
ddosxijtsqrtsxijt4e
ddo
e
ddoe
ddo写入气候场ope
10filedstudyformshixio
eatgrdformu
formattedaccessdirectrecl
x
yirec0dot112irecirec1
fwrite10recirecatijti1
xj1
ye
ddoclose10写入距平场ope
11filedstudyformshixio
exdgrdformu
formattedaccessdirectrecl
x
yirec0doit14
dot112irecirec1write11recirecxdijtiti1
xj1
y
e
ddoe
ddoclose11!写入均方差场ope
12filedstudyformshixio
esxgrdformu
formattedaccessdirectrecl
x
yirec0dot112irecirec1write12recirecsxijti1
xj1
ye
ddoclose121000format2i72000format37F623000format37f814000format37f72e
dprogramsx1
运行结果:
fGrads文件
气候场
rei
ite
ablepri
tdstudyformshixio
eatgmfope
dstudyformshixio
eatctlsetgridoffsetgradsoffsetlat040setlo
60150setlev500mo
1whilemo
12settmo
dhdrawtitle1982yearmo
mo
thpri
tcmo
mo
1e
dwhiledisablepri
t
f距平场
rei
ite
ablepri
tdstudyformshixio
esxgmfope
dstudyformshixio
esxctlsetgridoffsetgradsoffsetlat040setlo
60150setlev500year1982whileyear1984mo
1whilemo
12settmo
dhdrawtitle500hPayearyearmo
mo
tha
omalypri
tcmo
mo
1e
dwhileyearyear1e
dwhiledisablepri
t
f均方r