用户模型接口的问题!请高手看下我的错在哪了!? 下面是我的程序: C$ #2 BY: ANAVI DATE: 15-NOV-1994 FIXED TYPO INI(NINT)-->INT(NINT) C$ #1 BY: ANAVI DATE: 1-JUL-1994 NEW FOR USER MODELS C C User Unit Operation Model (or Report) Subroutine for USER C SUBROUTINE USRUSR (NSIN, NINFI, SIN1, SIN2, SIN3, 2 SIN4, SINFI, NSOUT, NINFO, SOUT1, 3 SOUT2, SOUT3, SOUT4, SINFO, NSUBS, 4 IDXSUB, ITYPE, NINT, INT, NREAL, 5 REAL, IDS, NPO, NBOPST, NIWORK, 6 IWORK, NWORK, WORK, NSIZE, SIZE, 7 INTSIZ, LD) C IMPLICIT NONE C C DECLARE VARIABLES USED IN DIMENSIONING C INTEGER NSUBS, NINT, NPO, NIWORK,NWORK, + NSIZE C #include "ppexec_user.cmn" EQUIVALENCE (RMISS, USER_RUMISS) EQUIVALENCE (IMISS, USER_IUMISS) C C C #include "dms_ncomp.cmn" C C THIS SUBROUTINE WILL SPLIT AN INLET TO TWO OUTLETS C C RETRIEVE THE SLPIT FACTOR C C C DECLARE ARGUMENTS C INTEGER IDXSUB(NSUBS),ITYPE(NSUBS), INT(NINT), + IDS(2,13), NBOPST(6,NPO), + IWORK(NIWORK),INTSIZ(NSIZE),NSIN, NINFI, + NSOUT, NINFO, NREAL, LD, I REAL*8 SIN1(1), SIN2(1), SIN3(1), + SIN4(1), SOUT1(1), SOUT2(1), + SOUT3(1), SOUT4(1), WORK(NWORK), + SIZE(NSIZE), SINFI, SINFO C C DECLARE LOCAL VARIABLES C INTEGER IMISS REAL*8 REAL(NREAL), RMISS, SPLIT C C BEGIN EXECUTABLE CODE REAL*8 POWER,PRES,T,PH2,PO2,VP,H2,CKH2,YDH2,O2,CKO2,YDO2 REAL*8 IC,VT,VANODE,VCATHODE,UF,AREA,V,H2O,CURRENT DATA PH2,PO2,UF,AREA/0.2,0.2,0.65,96.1/ POWER=REAL(1) PRES=REAL(2) T=REAL(3) H2=REAL(4) O2=REAL(5) IC=42 VP=76*LOG10(PRES/1) VT=0.008*(T-1000)*IC VANODE=172*LOG10(PH2/0.15) VCATHODE=92*LOG10(PO2/0.164) V=1.38+(VP+VT+VANODE+VCATHODE)/1000 CURRENT=POWER*1000/V YDH2=0.018655*CURRENT/UF CKH2=H2-YDH2 YDO2=0.5*YDH2 CKO2=O2-YDO2 H2O=YDH2 SOUT1(1)=CKO2 SOUT2(1)=CKH2 SOUT3(1)=H2O C SPLIT = 0D0 IF ( NREAL .GT. 0 ) SPLIT = REAL(1) IF ( SPLIT .LT. 0D0 ) SPLIT = 0D0 IF ( SPLIT .GE. RMISS ) SPLIT = 0D0 C C SPLIT FIRST INLET TO FIRST AND SECOND OUTLETS C DO 100 I = 1, NCOMP_NCC+1 SOUT1(1) = SIN1(1) * SPLIT SOUT2(1) = SIN1(1) - SOUT1(1) 100 CONTINUE C C COPY ALL INTENSIVE VARIABLES C DO 200 I = NCOMP_NCC+2, NCOMP_NCC+3 SOUT1(1) = SIN1(1) SOUT2(1) = SIN1(1) 200 CONTINUE C RETURN END 我的模型本意就是模拟H2和O2的反应,进口为H2和O2,出口是剩余的气体和生成的水。但是运行后说是我的质量不平衡,而且出口无气流,都是0。所以我怀疑是我的模型接口出了问题。看了半天用户模型手册也没看懂。我希望的结果是显示进出口的H2、O2、水的摩尔数。 下面是我模拟的截图:查看更多4个回答 . 3人已关注