fortran编写气象要素重现期

自由行

最近在网上查阅计算气象要素重现期的帖子，很少，很想与群里与大牛们进行交流，想问问大家，我最近在做降水重现期的计算，比如降水量的重现期，自己有编程，但是不停的提示错误，想请大家一起来交流，有相关程序可以借鉴吗，谢谢。反查程序发现没有定义类型的错误


program main
parameter(n=10,pi=3.1415926)
parameter(mx1=4.43,mx2=3.33,mx3=4.93,mx4=3.43,mx5=2.84,mx6=3.95,mx7=4.49,mx8=5.32,mx9=3.93,mx10=4.23)
real avex,sx,aa,bb,mx,pp,dd(n),xx,ss,ppp

OPEN(100,FILE='jizhi.txt')
open (1,file='jizhigengdie3.txt')
do i=1,n
READ(100,*) dd(i)                           
end do
   
      
xx=0.0
ss=0.0
      
do i=1,n
      xx=xx+dd(i)
end do
avex=xx/n
do i=1,n
ss=ss+(dd(i)-avex)*(dd(i)-avex)
end do
      sx=sqrt(ss/n)
  
aa=real(1.2825/sx)
bb=avex-0.5772*(sqrt(6.0)/pi)*sx
      
    mx=mx1
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,4.43
      mx=mx2
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,3.33
mx=mx3
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,4.93
mx=mx4
call gbe(pp,aa,bb,mx)
      write(1,*)pp,aa,bb,3.43
mx=mx5
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,2.84
mx=mx6
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,3.95
    mx=mx7
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,4.49
    mx=mx8
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,5.32
     mx=mx9
call gbe(ppp,aa,bb,mx)
      write(1,*)ppp,aa,bb,3.93
     mx=mx10
call gbe(ppp,aa,bb,mx)                                
      write(1,*)ppp,aa,bb,4.23

      end
subroutine gbe(ppp,aa,bb,mx)  !输出结果
      real ppp,aa,bb,mx  
   ppp=1-1.0/exp(exp(1.0/aa(mx-bb)))  
end
                                                                              

但是程序总显示警告：Warning: In the call to GBE, actual argument #2 does not match the type and。。。

王磊

建议楼主加上implicit none试试

自由行

王磊 发表于 2017-2-9 18:36
建议楼主加上implicit none试试

好的，我马上试试

qiulin_first

大神，这个是哪种算法？

nunu18

楼主你的程序调试出来了没？