请问运行fortran程序时 出现math—exp：overflow error 怎么办

心语

程序是求各等压面上的比湿：
do it=1,nt
  do iz=1,nz
   do ij=1,ny
            do ii=1,nx
         if(temper(ii,ij,iz,it).ge.-15)then
            AA=17.2693882
                        BB=35.86
                        else if(temper(ii,ij,iz,it).le.-40)then
                        AA=21.8745584
                        BB=7.66
                        else
                        AA=21.8745584+(17.2693882-21.8745584)/(-15-(-40.))*(temper(ii,ij,iz,it)-(-40))
                        BB=7.66+(35.86-7.66)/(-15-(-40.))*(temper(ii,ij,iz,it)-(-40))
                        endif
                        EE=6.1078*exp(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))-BB)
                        q(ii,ij,iz,it)=622*EE/(p(iz)-0.378*EE)
                enddo
          enddo
        enddo
enddo

其中temper数组中数据也就-60到30之间 ，ttd数组中数据在0到30之间，P数组的数据为1000,925,850,700,500,400,300,250,200,150,100
请问要咋改啊

river

提示说的比较清楚了，数学函数exp出错了，溢出错误。估计是你的数据超出了exp的定义域了。是不是没处理缺测值啊

心语

river 发表于 2013-5-25 14:18

                               
登录/注册后可看大图

提示说的比较清楚了，数学函数exp出错了，溢出错误。估计是你的数据超出了exp的定义域了。是不是没处理缺测 ...

应该没有缺测值  数据是老师给的，没提有缺测的

river

心语 发表于 2013-5-25 14:25

                               
登录/注册后可看大图

应该没有缺测值  数据是老师给的，没提有缺测的

自己输出到屏幕上看一下，老师不可能什么都交代给你

言深深

建议将 EE=6.1078*exp(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))-BB)
                        q(ii,ij,iz,it)=622*EE/(p(iz)-0.378*EE)
这两个公式里面的变量都print出来看看，可能出现了异常

心语

言深深 发表于 2013-5-25 14:29

                               
登录/注册后可看大图

建议将 EE=6.1078*exp(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-ttd(ii,ij ...

其他变量没错  就是EE输出前半部分是0.0000E+00,后半部分是无穷

心语

river 发表于 2013-5-25 14:29

                               
登录/注册后可看大图

自己输出到屏幕上看一下，老师不可能什么都交代给你

其他变量没错  就是EE输出前半部分是0.0000E+00,后半部分是无穷

心语

print *,(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))-BB)的结果

心语

我懂啦  这儿少了个括号exp（(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))-BB)）

river

心语 发表于 2013-5-25 16:02

                               
登录/注册后可看大图

我懂啦  这儿少了个括号exp（(AA*temper(ii,ij,iz,it)-ttd(ii,ij,iz,it))/(273.16+(temper(ii,ij,iz,it)-tt ...

自问自答了，不错