计算气象要素重现期

自由行

通过编程根据10个要素的值，分别计算这些要素的重现期，但是程序一直显示有10个warning，自己看了很久不知道具体怎么修改想请教下各位大神，比较急，希望大神门交流。
！计算气象要素重现期
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
!cccc  ccccccccccccccccccccccccccccccccccccccccccccc   
OPEN(100,FILE='jizhi.txt')
open (1,file='jizhigengdie3.txt')
do i=1,n
READ(100,*) dd(i)                           
end do
!c----------------------------------------------------      
      
xx=0.0
ss=0.0
      
do i=1,n
      xx=xx+dd(i)
end do
avex=xx/n
ss=ss+(dd(i)-avex)*(dd(i)-avex)
do i=1,n
end do
      sx=sqrt(ss/n)
  
aa=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(-aa(mx-bb)))  
end

！错误提示：Warning: In the call to GBE, actual argument #2 does not match the type and kind of the correspondin

言深深

错误提示的意思是gbe子程序中aa的定义和你传递下去的aa的定义不同，但是你这边似乎是一样的···
不过我想问一下
ss=ss+(dd(i)-avex)*(dd(i)-avex)
do i=1,n
end do
这是什么鬼？

自由行

言深深 发表于 2017-2-9 09:21
错误提示的意思是gbe子程序中aa的定义和你传递下去的aa的定义不同，但是你这边似乎是一样的···
不过我 ...

这个粘贴时弄错了，do i=1,n是应该写在第一行的，但是改过之后还是显示这个warnning，请教大神这个该怎么办喃