气象统计分析 实习六 气象要素线性趋势分析方法

实习六 气象要素线性趋势分析方法

分析中国夏季降水线性趋势的分布特征

利用数据160zhan-rainfall-summer.txt，编写求1982-2006年中国160站各站夏季降水线性倾向率，给出分布图，并进行简单分析。read-rain.for给出了阅读资料的fortran程序，数据在文件夹中单独给出。

Fortran代码：

program SX6
    implicit none 
    integer,parameter :: m=160,n=25 
    integer :: sta(m)
    real :: jd(m),wd(m),f(m,n),t(n) 
    real :: a(m),b(m),b_(m),syear,sf0,r 
    integer :: i,j

    open(22,file='D:\data\160zhan-rainfall-summer.txt')
    read(22,*) !空读一行
    do i=1,m
        read(22,*)sta(i),jd(i),wd(i),(f(i,j),j=1,n)
    enddo
    close(22)
    do i=1,n
        t(i)=i  !将1982-2006改成时间序列1,2,3，...，n
    enddo
    do i=1,m
        call xgxs(n,t(1:n),f(i,1:n),r,syear,sf0,b_(i)) 
    end do
    call stntogrd(jd,wd,m,b_)  !调用输出函