自定义模型通不过编译

xingzaiyuanfang · 发表于 2014-6-11 21:41:34

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

请教一下各位大神，我的AC8B自定义模型程序总是通不过编译报错END STATEMENT IS MISSING,Syntax error found identifier strtin when
expecting one of the end of statement, 到底错在哪里啊，求助一下！！！！
程序如下
SUBROUTINE AC8B(MC,ISLOT)
INTEGER MC,ISLOT
C
C MC =MACHINE ARRAY INDEX
C ISLOT=ARRAY ALLOCATION TQABLE INDEX
C J =STRTIN(1,ISLOT) [USES CON(J) THROUGH CON(J+20) ]
C K =STRTIN(2,ISLOT) [USES STATE(K) AND STATE(K+4) ]
C
$INSERT COMON4
C
INTRINSIC MAX,  EXP,  LOG
EXTERNAL BADMID
C
INTEGER  IB,  J,  K
REAL VERROR,  VPID,  SE,  IN, FEX,  VEMAX, A, B
REAL VE0,  VR0, VPID0
C
C DO NOTHING IN ACTIVITIES 'DOCU' AND 'DYDA' FOR NOW
C
IF (MODE.GT.4) RETURN
C
C BUS SEQUENCS NUMBER NEGATIVE IF MACHINE
C IS OFF LINE,SVS OR INDUCTION MACHINE
C
IB=NUMTRM(MC)
IF (IB.LE.0) RETURN
C
IF (MIDTRM) GO TO 900
C
C GET STARTING 'CON' AND 'STATE' INDICES
C
J=STRTIN(1,ISLOT)
K=STRTIN(2,ISLOT)
B=LOG(CON(J+20)/CON(J+18))/(CON(19)-CON(17))
A=CON(J+18)/EXP(B*CON(J+17))
C
GO TO  (100,200,300,400),MODE
C
C MODE = 1 – INITIALIZE
C
100 STATE(K)=ECOMP(MC)
VREF(MC)=ECOMP(MC)-VUEL(MC)-VOEL(MC)
STATE(K+2)=0
VE0=EFD(MC)+0.577*CON(J+11)*XADIFD(MC)
VR0=(A*EXP(B*VE0)+CON(J+13))*VE0+CON(J+12)*XADIFD(MC)
VPID0=VR0/CON(J+7)
STATE(K+4)=VE0
STATE(K+3)=VR0
STATE(K+2)=VPID0
VAR(L)=VPID0
VAR(L+1)=1-0.577*XADIFD(MC)/VE0
VAR(L+2)=(CON(J+15)-CON(J+12)*XADIFD(MC))/(CON(J+13)+A*EXP(B*VE0))
IF (VR0.LT.CON(J+10).OR.VR0.GT.CON(J+9)) THEN
   WRITE(LPDEV,307) NUMBUS(IB),MACHID(MC)
ENDIF
IF (VE0.LT.CON(J+16).OR.VE0.GT.VAR(L+2)) THEN
   WRITE(LPDEV,307) NUMBUS(IB),MACHID(MC)
ENDIF
RETURN
C
C MODE = 2 – CALCULATE DERIVATIVES
C
200 DSTATE(K)=(ECOMP(MC)-STATE(K))/CON(J)
VERROR=VREF(MC)+VOTHSG(MC)+VUEL(MC)+VOEL(MC)-STATE(K)
DSTATE(K+1)=CON(J+2)*VERRROR
DSTATE(K+2)=(CON(J+3)/CON(J+4)*VERRROR-STATE(K+2))/CON(J+4)
VPID=STATE(K+1)-STATE(K+2)+(CON(J+1)+CON(J+3)/CON(J+4))*VERROR
IF(VPID.GT.CON(J+5)) VPID=CON(J+5)
IF(VPID.LT.CON(J+6)) VPID=CON(J+6)
DSTATE(K+3)=(CON(J+7)*VPID-STATE(K+3))/CON(J+8)
IF(STATE(K+3).GT.CON(J+9)) STATE(K+3)=CON(J+9)
IF(STATE(K+3).LT.CON(J+10)) STATE(K+3)=CON(J+10)
SE=A*EXP(B*STATE(K+4))
DSTATE(K+4)=(STATE(K+3)-(SE+CON(13))*STATE(K+4)
* -CON(J+12)*XADIFD(MC))/CON(J+14)
VEMAX=(CON(J+15)-CON(J+12)*XADIFD(MC))/(CON(J+13)+SE)
IF(STATE(K+4).GT.CON(J+9)) STATE(K+3)=CON(J+9)
IF(STATE(K+4).LT.CON(J+16)) STATE(K+3)=CON(J+16)
RETURN
C
C MODE 3 – SET EFD
C
300 IN=CON(J+11)*XADIFD(MC)/STATE(K+4)
IF(IN.LE.0) THEN
   FEX=1
ELSEIF(IN.LE.0.433) THEN
   FEX=1-0.577*IN
ELSEIF(IN.LT.0.75) THEN
   FEX=SQRT(0.75-IN*IN)
ELSEIF(IN.LE.1) THEN
   FEX=1.732*(1-IN)
ELSE
   FEX=0
ENDIF
EFD(I)=FEX*STATE(K+4)
RETURN
C
C MODE = 4 – SET NINTEG
C
400 NINTEG=MAX(NINTEG,K+1)
RETURN
C
C MODEL NOT IMPLEMENTED FOR MSTR/MRUN
C
307 FORMAT('AC8B AT BUS',16,'MACHINE',A,
*  , 'INITIALIZED OUT OF LIMITS')
C
900 CALL BADMID(MC,IB,'AC8B')
RETURN
C
END

apsal · 发表于 2014-8-6 08:59:14

猜你是用的for格式写的Fortran代码。那你应该每行空出前6列，你现在只空了4列吧？

大肚腩 · 发表于 2015-1-14 00:34:06

做的自定义励磁？

kinghcyy · 发表于 2015-3-8 07:00:26

请问你怎么自建模型的？谢谢谢谢！

君弃剑 · 发表于 2016-7-24 16:34:28

楼主好人，请问你怎么自建模型的？谢

小宇1994 · 发表于 2017-10-18 11:20:28

楼主您好，我想问问您是怎样进行编译的？万分感谢！

snowlyy · 发表于 2020-5-14 00:41:35

最后结尾得 end sub routine AC8B这样？

账号		自动登录	找回密码
密码			立即加入

[讨论] 自定义模型通不过编译

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子