topic Error in code? in Intel® Fortran Compiler

Error in code?

mechprog — Fri, 23 Sep 2011 10:04:10 GMT

Hello all,

I can not understand the error in the following code:

CODE

IF(JRCD.NE.0)THEN

1 WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',

2 NOEL,'INTEGRATION POINT NUMBER ',NPT

I obtain the following error messages, but do not fully under it:

ERROR MESSAGE
cyclic_fatigue_B_Reduced.for(65) : Error: Syntax error, found IDENTIFIER 'WRITE' when expecting one of: ;
1 WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',
--------^
cyclic_fatigue_B_Reduced.for(65) : Error: Syntax error, found CHARACTER_CONSTANT 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ' when expecting one of: => = . ( %
1 WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',
-------------------^
cyclic_fatigue_B_Reduced.for(66) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: :: , : )
2 NOEL,'INTEGRATION POINT NUMBER ',NPT
------------------------------------------------^

Any input would be much appreciated. Thanks.

Error in code?

Arjen_Markus — Fri, 23 Sep 2011 10:07:57 GMT

The "1" in the second statement you have shown is probably in the sixth column (given
yesterday's discussion) so itcontinues the first line. But that is already complete.

Remove that "1" and it will probably work.

Regards,

Arjen

Error in code?

mechprog — Fri, 23 Sep 2011 10:47:44 GMT

Thanks Arjen,

I already had the column thing sorted in my code (I just copy and pasted the code on the OP which lost the formatting).

Problem solved!

Code changed to

IF(JRCD.NE.0)THEN

WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',

1 NOEL,'INTEGRATION POINT NUMBER ',NPT

Error in code?

mechprog — Fri, 23 Sep 2011 13:45:19 GMT

I am curious, looking up the following code I found that DIMENSION e.g.

DIMENSION

ARRAY(15)

"DIMENSION" defines an array, in this case ARRAY(15), why is it not just called "ARRAY" and not "DIMENSION"?

Thanks

Error in code?

Arjen_Markus — Fri, 23 Sep 2011 13:54:15 GMT

DIMENSION is the keyword, ARRAY is the name of the array variable. There is no keyword ARRAY in Fortran

I do not know the historical background for using the word "DIMENSION", but to me it
seems quite apt.

Regards,

Arjen

Error in code?

Steven_L_Intel1 — Fri, 23 Sep 2011 13:56:53 GMT

DIMENSION is an attribute. As with most attributes, it can be used as a separate statement, such as:

DIMENSION ARRAY(15)

This might need a separate type declaration statement if implicit typing was not to be used.

or, in more modern Fortran, in a comma list of attributes:

REAL, DIMENSION(15) :: ARRAY

One can also use:

REAL ARRAY(15)

DIMENSION is one of those attributes where you don't have to use the keyword. In more modern code I'd expect to find it in a comma list of attributes for a single declaration. The DIMENSION statement is a holdover from very early Fortran.

Error in code?

mechprog — Fri, 23 Sep 2011 15:02:35 GMT

Thanks for the really interesting input from everyone on this..

I have created an empty text file called STATEV.txt and am trying to write the data stored in STATEV(1) (produced during Abaqus analysis)into the text file. The text fileremains blank. Can anyone help? Thanks.

SUBROUTINE
USDFLD (FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,TIME,DTIME,CMNAME,ORNAME,
NFIELD,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,NDI,NSHR,
COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
c
INCLUDE 'ABA_PARAM.INC'
c
CHARACTER*80
CMNAME,ORNAME
CHARACTER*3FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),T(3,3),TIME
DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
c
CALL GETVRM('E',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
MATLAYO,LACCFLA)
EPS=ABS(ARRAY(1))
c
FIELD(1)=EPS
c
STATEV(1)=FIELD(1)
c
OPEN (unit=1, file='STATEV.txt')
WRITE (1),STATEV(1)
c
RETURN
END

Error in code?

jimdempseyatthecove — Fri, 23 Sep 2011 15:22:22 GMT

Place a break point on the WRITE statement and see what is in STATEV(1).
What is in STATEV? What type is STATEV (INCLUDE 'ABA_PARAM.INC' may have declared STATEV as character)
Are you looking in the correct directory for STATEV.txt?

Jim Dempsey

Error in code?

mechprog — Fri, 23 Sep 2011 15:49:17 GMT

Hi Jim,

ABA_PARAM.INC is an internal file which includes one of the following in the user subroutine based on the fact that the solver is SINGLE or DOUBLE PRECISION:

parameter (nprecd=1)
or
implicit real*8(a-h,o-z)
parameter (nprecd=2)

I am using this code as a subroutine during and SIMULIA abaqus analysis, so i am not sure how I can use a break point to see what STATEV(1) actually looks like during the analysis?

STATEV(1) is defined as "An array containing the solution-dependant state variable" produced during an analysis run.

In terms of the directory of STATEV(1) I have include the path C:\......\STATEV(1) and still recieve a blank text file at the end of an analysis.

Thanks again for all your help.

Error in code?

jimdempseyatthecove — Fri, 23 Sep 2011 20:18:55 GMT

Then this means that STATEV is either a REAL*4 or REAL*8 array.
The code you provided should have

a) created an output file of the name specified within the current directory at the time of the OPEN
b) written a floating point value, presumably non-blank
c) closed the file

The important factor is "what was the current directory" at the time of the OPEN?

The current directory is not necessarily the directory where the execuitable is located. In fact the file (as specified in your program) is opened: within the current directory .AND. within the current drive
Each are seperately settable.

It would appear the the file you are looking at is not the file written.
IOW this may have been an older empty file generated at some other time in the location you are looking at (but not the location where the current file is written).

How are you launching the program? If from IDE you can check the properties (default is project folder .NOT. execuitable folder).

Jim Dempsey

Error in code?

mechprog — Sat, 24 Sep 2011 07:36:24 GMT

Hi Jim,

I am not sure what you mean by ""what was the current directory" at the time of the OPEN?"? I have created a STATEV.txt file and placed this into the working directory of the analysis with ABAQUS.

The program is being launch by ABAQUS which calls the subrotuine I have written above during an analysis run. I use the following command in the ABAQUS Command Window to launch the analysis input file (.inp) and subroutine (.for):

abaqus job=model_1 user=Cyclic_Fatigue_B_Reduced

Thanks for all you help.

Error in code?

jimdempseyatthecove — Sat, 24 Sep 2011 12:28:10 GMT

Your subroutine is not closing the file STATEV.txt. The WRITE is buffered. Should you examine STATEV.txt prior to close then it may be empty. If STATEV.txt is required to remain open (more state variables to be written) consider adding FLUSH(1).

Consider changing the I/O unit number for STATEV.txt from 1 to some number greater than 9. Lower unit numbers are typically reserved for legacy I/O devices. Ignore this if ABAQUS requires STATEV.txt to be located on unit 1. Also, do not use literal unit numbers. Use parameters from shared header or module.

integer, parameter :: uSTATEV = 10
...
OPEN(UNIT=uSTATEV, FILE='STATEV.txt', STATUS='REPLACE')

** or STATUS='OLD' if that is the case

Jim Dempsey

Error in code?

mechprog — Wed, 28 Sep 2011 10:52:51 GMT

Hi everyone,

Still having problems with this, I have updated this in line with the advice from the forum and reference to the following text:

Structured Fortran 77 for engineers and scientists (5th Ed.) Delores M.Etter

INTEGER, PARAMETER :: uSTATEV = 105

101 FORMAT (F12.6)
OPEN (UNIT=105,
1 FILE='C:\.....\STATEV.txt',
2 STATUS='NEW')
WRITE(105,101) STATEV(1)
CLOSE (UNIT=105)

Please note. I have taken into account the .for column formatting in the code.

I would really appreciate any further input or advice. Thanks.

Error in code?

TimP — Wed, 28 Sep 2011 11:53:19 GMT

Did you try changing that FILE path name to an actual path, or allowing it to use your current working directory?

Error in code?

mechprog — Wed, 28 Sep 2011 12:08:16 GMT

Thanks for you response, yes i have tried both and still no success. in the code below:

FILE='C:\.....\STATEV.txt',

I have just shorterned the text on the forum post to make it easier to read. In the real code I have tried both including the full FILE path as well as trying just:

FILE='STATEV.txt',

Still not sure whereI am going wrong. Thanks.

Error in code?

mecej4 — Wed, 28 Sep 2011 12:13:07 GMT

The WRITE statement should not have a continuation character in COL-6.

The syntax is

IF ( CONDITION ) THEN
Statements to execute if CONDITION is .TRUE.
ENDIF

Note that the first and last lines are complete statements by themselves. Likewise, the statements in the IF...ENDIF block are complete statements.

Error in code?

mechprog — Wed, 28 Sep 2011 13:06:30 GMT

Thanks mecej4, my WRITE STATEMENT does not have a continuation character in COL-6.

Sorry, did I misunderstand your post?

Error in code?

mecej4 — Wed, 28 Sep 2011 13:18:24 GMT

The initiating post of this thread contains:

IF(JRCD.NE.0)THEN

1 WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',

2 NOEL,'INTEGRATION POINT NUMBER ',NPT

If these are the statements in your program, the '1' at the beginning of the second line should not be there. The error message from the compiler is consistent with such a continuation character having been scanned.

Error in code?

mechprog — Wed, 28 Sep 2011 13:40:10 GMT

Thanks mecej4, I corrected this before. This still did no solve the problem I am having with not being able to write to the STATEV.txt file. Perhaps the problem is that there is nothing to write to the file, but I just do know about it. I will try write a simple programme to write to a file to test out the code.

[bash] 101    FORMAT (F12.6)
        OPEN (UNIT=uSTATEV,
      1 FILE='STATEV.txt',
      2 STATUS='NEW') 
        WRITE (uSTATEV,101) STATEV(1)
        CLOSE (UNIT=uSTATEV)





[/bash]

Error in code?

mecej4 — Wed, 28 Sep 2011 13:56:19 GMT

You will find this forum more helpful if you post complete code instead of fragments, and if you show the actual error messages from the compiler and runtime, instead of simply saying that something "did not work".

When you post code, please use the Syntax Highlighter Tool (the yellow pencil/pen icon).