topic This looks like it should be in Intel® Fortran Compiler

Fortran error

Le__Dai — Tue, 28 Aug 2018 12:42:25 GMT

Hello, i am getting the following error when trying to complied my program.

Here is my code:

*USER SUBROUTINE

SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,

1 JLTYP,TEMP,PRESS,SNAME)

include 'ABA_PARAM.INC'

DIMENSION FLUX(2),TIME(2),COORDS(3)

CHARACTER*80 SNAME

*code

! REAL SOL,TEMP,PRESS

! INTEGER JSTEP,JINC,NOEL,NPT,JLTYP

! REAL t,dt,x,y,z

! REAL FLUX1,FLUX2,FLUX3

INTEGER flag,kstep

REAL t,x,y,z

Real Q,a,b,cf,cr,pi,m,n

*times

t=time(1)

a= 6

b = 2.5

cf = 3

cr= 1

pi = 3.14159265359

*start point parameters

FLUX(1) = 0

FLUX(2) = 0

x = coords(1)

y = coords(2)

z = coords(3)

C FORM HEAT FLUX FROM SAMPLING POINTS

if(kstep.EQ.3) then

Q = 0.5*10*90

x1=x

y1=y

z1=z - 5*t

m=exp(-3*((x1)**2/(a)**2+(y1)**2/(b)**2+(z1)**2/(cf)** 2))

n=exp(-3*((x1)**2/(a)**2+(y1)**2/(b)**2+(z1)**2/(cr)** 2))

FLUX(1)=(6*1.732*Q)*((1.5/cf)*m+(.5/cr)*n)/(pi*sqrt(pi)*(a*b))

end if

return

end

this is the error I am receiving:

Welding.for(35): error #5082: Syntax error, found '=' when expecting one of: .EQV. .NEQV. .XOR. .OR. .AND. .LT. < .LE. <= .EQ. == .NE. /= .GT. > ...

y=y

------^

Welding.for(34): error #6385: The highest data type rank permitted is INTEGER(KIND=8). [X1]

x1=x1

-----------^

Welding.for(35): error #6385: The highest data type rank permitted is INTEGER(KIND=8).

y=y

-------^

compilation aborted for Welding.for (code 1)

Abaqus Error: Problem during compilation - C:\temp\Han\Welding.for

Begin MFS->SFS and SIM cleanup

8/28/2018 9:25:28 PM

Please let me know your thoughts on these few errors and how I could approach fixing them. Thanks!

The error messages you report

Steve_Lionel — Tue, 28 Aug 2018 19:11:16 GMT

The error messages you report don't match what's in the code you show. Please attach a complete source rather than typing it in.

As a side comment, your value of PI is wrong (last digit should be 6 if correctly rounded), not that it matters much since PI is implicitly single precision (as is your constant), so the extra digits won't be used.

Is this an Abaqus user

Greg_T_ — Tue, 28 Aug 2018 19:35:50 GMT

Is this an Abaqus user subroutine? The routine name and later message indicate this is a user subroutine for an Abaqus finite element analysis. Have you also tried compiling the routine outside of Abaqus, perhaps from the Intel Fortran initialized command line? It can also be worthwhile to create a small main program that can call the subroutine and pass in values to test that it is working as expected before including in the Abaqus FEA.

Regards, Greg

This looks like it should be

Lorri_M_Intel — Wed, 29 Aug 2018 13:24:13 GMT

This looks like it should be a fixed-form file (the extension ".for", the "C" in column 1). Please make sure that you are following all the fixed-form rules and that your assignment statements start at column 7 or later.

Also, follow Greg's suggestion above, so that you can confirm that your subroutine works as expected.

--Lorri