| 查看: 1359 | 回复: 4 | |||
lijunjiexm铁虫 (小有名气)
|
[求助]
求助一个编译出错的简单程序，请多多指教，谢谢
|
|
在建立的一个project的Source files的其中命名为STRESSAVG的f文件如下： C C********************************************************************** C********************************************************************** C********************************************************************** C SUBROUTINE STRESSAVG(NUMEL, STORE, I1, I2, STR, NKT, > ORDWHERE, FIRST, LAST, > COUNTPLUS, COUNTMINS, COUNTFLAG) IMPLICIT NONE INCLUDE 'PARAMETERS.INC' INCLUDE 'NUMBERS.INC' INTEGER*4 ORDWHERE(NDMB) INTEGER*4 I INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) REAL*8 STORE (NDMB) REAL*8 SUM1, SUM2, NKT REAL*8 STR(2, 3, NUMEL) INTEGER*4 II, J, I1, I2, INDEX, NUMEL INTEGER*4 COUNTPLUS(NUMEL), COUNTMINS(NUMEL) LOGICAL*1 COUNTFLAG C C---------------------------------------------------------------------- C STR (i , j , k) C C ^ ^ ^ C | | |_________ element number C | |_____________ components in order xx, yx, yy C |_________________ triangle number (1 or 2) C C---------------------------------------------------------------------- C IF (COUNTFLAG) THEN DO I = 1 , NUMEL COUNTPLUS(I) = 0 COUNTMINS(I) = 0 ENDDO ENDIF IF ((I1 .EQ. 1) .AND. (I2 .EQ. 1)) THEN INDEX = 1 ELSEIF ((I1 .EQ. 1) .AND. (I2 .EQ. 2)) THEN INDEX = 2 ELSEIF ((I1 .EQ. 2) .AND. (I2 .EQ. 2)) THEN INDEX = 3 ELSE PRINT *,'Wrong indices for stress tensor' ENDIF DO J = 1 , NUMEL SUM1 = ZERO SUM2 = ZERO *vdir nodep CDIR$ IVDEP DO II = FIRST(J) , LAST(J) I = ORDWHERE(II) IF (I .GT. 0) THEN SUM1 = SUM1 + STORE(I) IF (COUNTFLAG) THEN COUNTPLUS(J) = COUNTPLUS(J) + 1 ENDIF ELSE SUM2 = SUM2 + STORE(-I) IF (COUNTFLAG) THEN COUNTMINS(J) = COUNTMINS(J) + 1 ENDIF ENDIF ENDDO STR(1 , INDEX , J) = SUM1 STR(2 , INDEX , J) = SUM2 ENDDO DO J = 1 , NUMEL STR(1 , INDEX , J) = STR(1 , INDEX , J) * NKT / > MAX(1 , COUNTPLUS(J)) STR(2 , INDEX , J) = STR(2 , INDEX , J) * NKT / > MAX(1 , COUNTMINS(J)) ENDDO IF ((INDEX .EQ. 1) .OR. (INDEX .EQ. 3)) THEN DO J = 1 , NUMEL STR(1 , INDEX , J) = STR(1 , INDEX , J) - NKT STR(2 , INDEX , J) = STR(2 , INDEX , J) - NKT ENDDO ENDIF RETURN END 而PARAMETERS.INC和NUMBERS.INC同时保存在同一个project下的Extenal Dependenciec文件夹中： PARAMETERS.INC文件为： INTEGER*4 NSD, NEN, NDOF INTEGER*4 NNPMAX, NELMAX, NEQMAX INTEGER*4 NSBDIV, NBWMAX INTEGER*4 NBUF , MAXTRACE INTEGER*4 NDMB C PARAMETER ( NSD = 2, > NEN = 4, > NDOF = 2, > NNPMAX = 3000, > NEQMAX = NDOF * NNPMAX, > NSBDIV = 110, > NBWMAX = NDOF * NSBDIV, > NDMB = 524288, C > NDMB = 65536, C > NDMB = 2048, > NBUF = 1024, > NELMAX = 3000, > MAXTRACE = 10000 ) NUMBERS.INC 文件为： REAL*8 ZERO, ONE, TWO, THREE, FOUR, FIVE, SIX, > HALF, THIRD, SIXTH, > FOURTH, TWSVTH, NINTH, TWFOUR, > PI, TWOPI, FTPI, > RSMALL PARAMETER (ZERO = 0.0 , > ONE = 1.0 , > TWO = 2.0 , > THREE = 3.0 , > FOUR = 4.0 , > FIVE = 5.0 , > SIX = 6.0 , > HALF = 1.0 / 2.0 , > THIRD = 1.0 / 3.0 , > SIXTH = 1.0 / 6.0 , > FOURTH = 1.0 / 4.0 , > TWSVTH = 1.0 / 27.0 , > NINTH = 1.0 / 9.0 , > TWFOUR = 24. , > PI = 3.14159265358979323846 , > TWOPI = 2. * 3.14159265358979323846 , > FTPI = 4. * 3.14159265358979323846 / 3. , > RSMALL = 1.0E-30 ) |
回复此楼» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
极限拉扯，如何抉择
已经有7人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助一简单matlab程序的出错问题，谢谢
已经有3人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】VASP5.2 编译出错make: *** [vasp] 错误 1,请大家帮忙分析分析,谢谢
已经有7人回复
- 【求助】f90子程序怎么编译错误？请大家帮忙，谢谢
已经有3人回复
- 【求助】请教编译过程中的两个问题【已完结】
已经有7人回复
- 【求助】为什么我编译了一个书本上的无错误的C++程序，却有这么多的错误？
已经有8人回复
- 【求助】fortran出错 编译通过 运行溢出
已经有10人回复
- 【求助】VASP编译出错
已经有5人回复
lijunjiexm
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 22.8
- 红花: 1
- 帖子: 182
- 在线: 87.1小时
- 虫号: 1198843
- 注册: 2011-02-03
- 专业: 常微分方程与动力系统
★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 13:31:08
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 13:31:08
|
在windows下Compaq Visual Fotran编译出现如下错误信息（请多多指教，谢谢）： -------------------Configuration: sucontr - Win32 Debug-------------------- Compiling Fortran... G:\Program Files\sucontr\src\stressavg.f G:\Program Files\sucontr\src\stressavg.f(14) : Error: This name does not have a type, and must have an explicit type. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------^ G:\Program Files\sucontr\src\stressavg.f(14) : Error: A specification expression is invalid. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------^ G:\Program Files\sucontr\src\stressavg.f(14) : Error: A specification expression is invalid. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------------------------^ G:\Program Files\sucontr\src\stressavg.f(17) : Error: A specification expression is invalid. [NUMEL] REAL*8 STR(2, 3, NUMEL) ---------------------------^ Error executing df.exe. stressavg.obj - 4 error(s), 0 warning(s) |
2楼2011-07-17 21:58:30
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
3楼2011-07-17 22:12:46
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2011-07-17 23:32:05
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
5楼2011-07-17 23:35:19
