24小时热门版块排行榜    

查看: 965  |  回复: 4
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 lx5599 的 10 个金币

lx5599

银虫 (初入文坛)

[交流] 【求助】Fortran编译链接Aspen的问题

各位哥哥姐姐们,东西有点长,不过请您看下去,您的坚持说不定对我有很大的影响,谢谢了啊!请各位Fortran高手指点:下面是我编的一个动力学程序,链接到ASPEN中的,出现的问题,(本人是外行,还望各位帮忙),再鞠躬!
源程序这样的:
Program DLX
CALL ASPEN(N,       NCOMP,   NR,      NRL,     NRV,     T,      TLIQ,   TVAP,    P,       PHFRAC,&
          F,       X,       Y,       IDX,     NBOPST,  KDIAG,  STOIC,  IHLBAS,  HLDLIQ,  TIMLIQ,&
                  IHVBAS,  HLDVAP,  TIMVAP,  NINT,    INT,     NREAL,  REAL,   RATES,   RATEV,   NINTB,&
                  INTB,    NREALB,  REALB,   NIWORK,  IWORK,   NWORK,  WORK)
STOP
END
!----------------------------------------------------------------------
SUBROUTINE ASPEN(N,       NCOMP,   NR,      NRL,     NRV,     T,      TLIQ,   TVAP,    P,       PHFRAC,&
                F,       X,       Y,       IDX,     NBOPST,  KDIAG,  STOIC,  IHLBAS,  HLDLIQ,  TIMLIQ,&
                        IHVBAS,  HLDVAP,  TIMVAP,  NINT,    INT,     NREAL,  REAL,   RATES,   RATEV,   NINTB,&
                        INTB,    NREALB,  REALB,   NIWORK,  IWORK,   NWORK,  WORK)
!----------------------------------------------------------------------

IMPLICIT NONE

!DECLARE VARIABLES USED IN DIMENSIONING
!----------------------------------------------------------------------

INTEGER NCOMP, NR, NRV, NINT, NINTB, NREALB, NIWORK, NWORK, USER_IUMISS
REAL    USER_RUMISS
!#include "ppexec_user.cmn"
    EQUIVALENCE(RMISS,USER_RUMISS)
    EQUIVALENCE(IMISS,USER_IUMISS)

!DECLARE ARGUMENTS
!----------------------------------------------------------------------

INTEGER NRL(3),        IDX(NCOMP), NBOPST(6), INT(NINT), INTB(NINTB),&
       IWORK(NIWORK), N,          KDIAG,     IHLBAS,    IHVBAS, NREAL
REAL*8 PHFRAC(1),  X(NCOMP,3),     Y(NCOMP),     STOIC(NCOMP,NR),  RATES(NCOMP),  RATEL(1),&
      RATEV(NRV), REALB(NREALB),  WORK(NWORK),  T,                TLIQ,          TVAP,&
          P,          F,              HLDLIQ,       TIMLIQ,           HLDVAP,        TIMVAP

!----------------------------------------------------------------------
!DECLARE LOCAL VARIABLES

INTEGER IMISS
REAL*8 REAL(NREAL), RMISS

!BEGIN EXECUTABLE CODE
!-----------------------------------------------------------------------
REAL*8, DIMENSION (1 : 5) ::  K, Ea, a, b,REATS
K  (1:5) = (/ 0.000683, 0.001584, 0.0000010948, 0.006148, 0.91509 /)
Ea (1:5) = (/ 15578,    54112,    49646,        17404,    27131   /)
a  (1:5) = (/ 0.6785,   1.7544,   0.55373,      0.30319,  1.6426  /)
b  (1:5) = (/ 0.4550,   2.2449,   1.9672,       0.11062,  0.66696 /)
!------------------------------------------------------------------------
REATS(1)= -K(1) * EXP( -Ea(1)/8.314/(T-273.15) ) * ( F*PHFRAC(1)*Y(1) )**a(1) * ( Y(7)*P )**b(1)
REATS(2)= -K(2) * EXP( -Ea(2)/8.314/(T-273.15) ) * ( F*PHFRAC(1)*Y(2) )**a(2) * ( Y(7)*P )**b(2)
REATS(3)= -K(3) * EXP( -Ea(3)/8.314/(T-273.15) ) * ( F*PHFRAC(1)*Y(3) )**a(3) * ( Y(7)*P )**b(3)
REATS(4)= -K(4) * EXP( -Ea(4)/8.314/(T-273.15) ) * ( F*PHFRAC(1)*Y(4) )**a(4) * ( Y(7)*P )**b(4)
REATS(5)= -K(5) * EXP( -Ea(5)/8.314/(T-273.15) ) * ( F*PHFRAC(1)*Y(5) )**a(5) * ( Y(7)*P )**b(5)
!-------------------------------------------------------------------------
RETURN
END
编译时,出现警告:
Compiling Fortran...
D:\exam\exam1\exam1.f90
D:\exam\exam1\exam1.f90(2) : Warning: In the call to ASPEN, actual argument #6 does not match the type and kind of the corresponding dummy argument.
CALL ASPEN(N,       NCOMP,   NR,      NRL,     NRV,     T,      TLIQ,   TVAP,    P,       PHFRAC,&


在我用Aspen中的 aspcomp进行编译时,出现这些错误:
第一类错误:
f90: Error: Illegal character in statement label field  [P]
f90: Error: Illegal character in statement label field  [r]
f90: Error: Illegal character in statement label field  [o]

第二类:
f90: Error: First statement in file must not be continued
exam1.f90(1) : Error: Syntax error, found END-OF-STATEMENT when expecting one of
: ( : % . = =>

第三类错误:
f90: Error: First statement in file must not be continued
exam1.f90(1) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: ( : % . = =>
这是什么原因啊?怎么解决?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianlangxingaa

铁杆木虫 (著名写手)

★ ★
余泽成(金币+2,VIP+0):谢谢参与应助! 1-4 21:32
很抱歉,我对你这个程序进行了测试,编译和运行时没有发现你发现的错误

但是你这个程序没有输出,所以不太清楚是不是真的有错误,还是编译器的问题

希望在检查下

cheers
2楼2010-01-04 19:25:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★
余泽成(金币+2,VIP+0):专家辛苦了! 1-4 21:33
我编译了一下,好像也没问题,
不过看你的第一类错误,好像是需要把第一行往里挪挪的意思,但是你这是90,我不知道。

[ Last edited by maomao1210 on 2010-1-4 at 20:45 ]
3楼2010-01-04 20:41:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

许虎_lenovo

铁虫 (初入文坛)


jjdg(金币+1,VIP+0):欢迎新虫 1-19 00:43
我做过类似的问题,要仔细排查,需要你个人努力,当初我就是这么弄出来的!!
4楼2010-01-18 21:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

f41898

至尊木虫 (文坛精英)

nono2009(金币+0,VIP+0):请勿在专业区灌水。在这个版面再提醒1次! 1-21 09:53
????
5楼2010-01-20 19:35:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lx5599 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见