计算CAPE指数程序里的一段程序，里面用红字标的语句在读取的时候出错。

zyq

这是计算CAPE指数程序里的一段程序，里面用红字标的语句在读取的时候出错，不知道是何原因？哪位高手帮我看下，先谢谢啦！
!-----------------计算LCL处的气压，温度，假相当位温 (也就是抬升凝结高度)
l=p(1)
m=t(1)
n=td(1)
tc=m-0.976*(m-n)/(0.976-0.000833*((237.3+n)**2)/(273.16+n))  !LCL处的温度
pc=l*((273.16+tc)/(273.16+m))**3.50
ec=6.11*10**((a*tc)/(b+tc)) !水汽压
qc=0.622*ec/(pc-0.378*ec) ! q=0.622*E/(p-0.378E)   
thse=(tc+273.16)*exp(0.28586*LOG(1000/pc)+2500*(qc/(338.52-0.24*(tc+273.16)+1.24*tc)))
write(*,*)"pc is:",pc,"tc is:",tc,"ec is:",ec
write(*,*)"qc is:",qc,"thse is:",thse

mofangbao

没有明白楼主的意思，编译能通过吗？出错的话是什么错误提示呢？还有啊 变量的上下文什么的，其实，说实话，调错这东西，别人能帮看出来的也只有语法错误了，不然就是从头开始分析程序，然后那着全部的程序开始调试，所以，光看这句话是没有办法看出错误的，既然叫调试，就得逐步的调

只能说，木有发现语法错误

阿华田

什么叫读取错误？能把你的错误信息发出来了。

topmad

把错误代码发上来吧，要不，很难看的

言深深

有可能是因为1000/pc=0，在log函数里面是超出定义域的，这样的错误显示代码6201，改成LOG(1000.0/pc)试试。如果不是就请如楼上所说的，最好给出错误显示

阿华田

言深深 发表于 2011-9-27 14:00

                               
登录/注册后可看大图

有可能是因为1000/pc=0，在log函数里面是超出定义域的，这样的错误显示代码6201，改成LOG(1000.0/pc)试试。 ...

我估计就是深深说的问题。我开始也认为log函数里面的出问题，还是深深看出了问题所在，我赌你赢！

言深深

阿华田 发表于 2011-9-27 17:48

                               
登录/注册后可看大图

我估计就是深深说的问题。我开始也认为log函数里面的出问题，还是深深看出了问题所在，我赌你赢！

因为不知道出错原因，还不好说，呵呵，不过谢谢支持

wazky

改成1000.0了，可是错误的现实代码还是6201

海蓝蓝

是不是数据有问题

river

wazky 发表于 2011-9-28 08:48

                               
登录/注册后可看大图

改成1000.0了，可是错误的现实代码还是6201

你这个问题深深一下就判断出来了，你还不回复给深深帮你解决啊。直接这么发上看来深深又看不到