fortran计算假相当位温然后计算对流性稳定度指数，运行问题，不懂哇

cycycycycy

!tlcl,theta,eqt应该分别代表着抬升凝结高度处的温度，位温，假相当位温
tlcl=1.e+36
theta=1.e+36
eqt=1.e+36
do k=1,MZ
  do t=1,MT
    do j=1,MY
          do i=1,MX         
          tlcl(i,j,k,t)=55.0+2840.0/(3.5*log(temp(i,j,k,t))-log(e(i,j,k,t))-4.805)
          tt=0.2854*(1.0-0.28*q(i,j,k,t))
          theta(i,j,k,t)=temp(i,j,k,t)*((1000.0/p(k))**tt)
          enddo
    enddo
  enddo
enddo
do k=1,MZ
  do t=1,MT
    do j=1,MY
          do i=1,MX         
      eqt(i,j,k,t)=theta(i,j,k,t)*exp(((3376./tlcl(i,j,k,t))-2.54)*q(i,j,k,t)*(1.0+0.81*q(i,j,k,t)))
      enddo
    enddo
  enddo
enddo
!对流性稳定指数Convective instability index
  do t=1,MT
    do j=1,MY
          do i=1,MX         
      Iconve(i,j,t)=eqt(i,j,5,t)-eqt(i,j,3,t)   
      enddo
    enddo
  enddo
根据我查找到的gs编的fortran，可是有问题，不懂啊
http://bbs.06climate.com/forum.php?mod=viewthread&tid=22546

勿忘草

这是全部的运行程序还是其中的一个子函数？

王磊

MATH 。。log。。可能是你的log函数里面的值出现了负数，比如log（-1），建议版主在tlcl(i,j,k,t)=55.0+2840.0/(3.5*log(temp(i,j,k,t))-log(e(i,j,k,t))-4.805)前输出temp(i,j,k,t))-log(e(i,j,k,t))-4.805 的值看看书不是有这种情况。

cycycycycy

王磊 发表于 2015-12-18 11:50
MATH 。。log。。可能是你的log函数里面的值出现了负数，比如log（-1），建议版主在tlcl(i,j,k,t)=55.0+284 ...

对！我也发现这个问题了！感谢！