课程名称: 气象统计方法课程实践
姓名: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