气象统计方法实习BD - 360文档中心

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

实习一：气候场、距平场、均方差场

编程如下:

parameter(ii=37,jj=17,mon=12,year=4)

real var(ii,jj,mon,year),ave(ii,jj,mon),jp(ii,jj,mon,year)

real s(ii,jj,mon)

integer i,j,iy,m

open(5,file='d:\ex1\h500.dat')

open(6,file='d:\ex1\ave.grd',form='binary')

open(7,file='d:\ex1\jp.grd',form='binary')

open(8,file='d:\ex1\s.grd',form='binary')

open(12,file='d:\ex1\outall.grd',form='binary'

open(9,file='d:\ex1\ave.txt')

open(10,file='d:\ex1\jp.txt')

open(11,file='d:\ex1\s.txt')

！读数据

DO iy=1,4

do m=1,12

!ccc read h500

read(5,1000)

read(5,2000) ((var(i,j,m,iy),i=1,ii),j=1,jj)

enddo

！计算气候场

do j=1,jj

do i=1,ii

do m=1,12

ave(i,j,m)=var(i,j,m,1)+var(i,j,m,2)+var(i,j,m,3)+var(i,j,m,4)

ave(i,j,m)=ave(i,j,m)/4.0

enddo

！计算距平场

do iy=1,4

do m=1,12

do j=1,jj

do i=1,ii

jp(i,j,m,iy)=var(i,j,m,iy)-ave(i,j,m)

enddo

！计算均方差场

do j=1,jj

do i=1,ii

do m=1,12

s(i,j,m)=jp(i,j,m,1)*jp(i,j,m,1)+jp(i,j,m,2)*jp(i,j,m,2)+jp(i,j /,m,3)*jp(i,j,m,3)+jp(i,j,m,4)*jp(i,j,m,4)

s(i,j,m)=s(i,j,m)/4.0

s(i,j,m)=sqrt(s(i,j,m))

enddo

do iy=1,4

do m=1,12

write(6)((ave(i,j,m),i=1,ii),j=1,jj)

write(7)((jp(i,j,m,iy),i=1,ii),j=1,jj)

write(8)((s(i,j,m),i=1,ii),j=1,jj)

write(9,2000)((ave(i,j,m),i=1,ii),j=1,jj)

write(10,2000)((jp(i,j,m,iy),i=1,ii),j=1,jj)

write(11,2000)((s(i,j,m),i=1,ii),j=1,jj)

write(12)((ave(i,j,m),i=1,ii),j=1,jj)

write(12)((jp(i,j,m,iy),i=1,ii),j=1,jj)

write(12)((s(i,j,m),i=1,ii),j=1,jj)

enddo

1000 format(2i7) 2000 format(37f8.1) close(5)

close(6)

close(7)

close(8)

close(9)

close(10)

close(11)

close(12)

end

给ave配的ctl文件：

dset ^d:\ex1\ave.grd

undef -9.99E+33

title NCEP/NCAR REANALYSIS PROJECT xdef 37 linear 60.000 2.500 ydef 17 linear 0.000 2.500

zdef 1 levels 500

tdef 12 linear JAN1982 12mo vars 1

ave 1 99 H500

endvars

给ave配的gs文件：

'reinit'

'open d:\ex1\ave.ctl'

'enable print d:\ex1\ave.gmf'

mon=1

while(mon<=12)

'set t 'mon''

'd ave'

'draw title qihouchang of 'mon' '

'print'

'c'

mon=mon+1

endwhile

'disable print'

;

气候场图：

一月份高度的气候场呈现南高北低的状态，陆地上的高度场比较稀疏，而在西太平洋上高度场比较密集。

八月份高度的气候场呈现东高西低的状态，在我国东北部以北以及印度东北部出现低压中心，而在赤道西太平洋地区出现高压中心。35°N以北高度分布很密集，而35°N以南比较稀疏。

给jp配的ctl文件：

dset ^d:\ex1\jp.grd

undef -9.99E+33

title NCEP/NCAR REANALYSIS PROJECT

xdef 37 linear 60.000 2.500

ydef 17 linear 0.000 2.500

zdef 1 levels 500

tdef 48 linear JAN1982 1mo

vars 1

jp 1 99 H500

endvars

给jp配的gs文件：

'reinit'

'open d:\ex1\jp.ctl'

'enable print d:\ex1\jp.gmf'

year=1982

while(year<=1985)