实习一

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

姓名:xxx

学号:xxx 班级:xxx

实习一：气候场、距平场、均方差场 编程如下:

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 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 enddo 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 enddo enddo 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 enddo 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 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) mon=1

while(mon<=12) 'set t 'mon'' 'd jp'

'draw title jupingchang of 'year'.'mon'' 'print' 'c'

mon=mon+1

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库实习一在线全文阅读。