Fortran runtime error bad integer for item 1 in list input

To expand on francescalus’ and Vladimir F’s comments:

Catching bugs

The simplest way of debugging a program is adding a bunch of print statements. If you add these to your code, like

    implicit none
    
    !declaring the single precision for real numbers, array, and variables.
    integer :: i, firstnumber
    real*4 :: average, sum_, standarddeviation
    
    real*4, allocatable :: arrayone(:)!making sure that the array is single precision real
    
    open(unit=1,file="C:cygwin64data sets for labs and assignmentsNumbers.txt")!opening file
    
    do 
        read(1,*)firstnumber!reading the first number
        print *, 'firstnumber = ', firstnumber
    end do

    allocate(arrayone(firstnumber))!allocating the array with the first number variable obtained from previous statement
    
    close(1)
    
    do i = 2, firstnumber
    
    open(i,file="C:cygwin64data sets for labs and assignmentsNumbers.txt")
        read(i,*)
        
        print *, 'arrayone = ', arrayone
        
        sum_ = Sum(arrayone) !sums all of the numbers from arrayone
        
        print *, 'sum_= ', sum_
        
        average = Sum_/firstnumber 
        
        print *, 'average= ', average
        
        standarddeviation = sqrt(sum_ - average) * (sum_ - average)!calculates the standard deviation of the data set
        
        print *, 'standarddeviation= ', standarddeviation
        
        close(i)
    end do
    print*,"Average = "
    print*,"Standard Deviation = "

then you can see what your program is actually doing, and work out what’s going wrong.

Your current error

Your first problem is that the lines

do 
  read(1,*)firstnumber!reading the first number
end do

are performing the read operation multiple times. It will literally try to read integers from the file forever. Obviously the file is finite, so eventually it runs out of valid things to read, and throws the error you are seeing.

To fix this, you should simply get rid of the infinite loop, by removing the lines do and end do.

The empty read

Your next problem is the statement

read(i,*)

This reads from the file you just opened, but doesn’t store the results anywhere. It looks like you’re trying to populate the arrayone array, which you would do like

read(i,*) arrayone(j)

where j is the index of arrayone corresponding to the element which you want to fill.

In a loop

The other problem you have is that you’re doing a lot of things inside a loop which you should be doing outside the loop. As it stands, the code

do i = 2, firstnumber
  open(i,file="C:cygwin64data sets for labs and assignmentsNumbers.txt")
  read(i,*)
  sum_ = Sum(arrayone)
  average = Sum_/firstnumber 
  standarddeviation = sqrt(sum_ - average) * (sum_ - average)
  close(i)
end do

will:

• open the file
• read the first number in the file, but not store it anywhere
• calculate the sum, average and standard deviation of the entire array arrayone, even though most of the elements of arrayone are not yet filled.
• close the file

and it will do all of this over and over again, firstnumber times.

I suspect that what you want to do is read the firstnumber elements of arrayone from lines 2 through firstnumber+1 of the file, and only then (after the loop) calculate the sum, average, and standard deviation of arrayone.

In order to read multiple rows from your file, you should only open it once (at the statement open(unit=1...), and only close it once (using close(1), but at the end of your program). Each read statement should then read from the same file (so each should be read(1,...).

I put «Real» MATC before I even opened my project and restart two twice, but still the same error message:
ELMER SOLVER (v 7.0) STARTED AT: 2014/08/11 10:31:01
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 — , CSC — IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Library version: 7.0 (Rev: exported)
MAIN: =============================================================
MAIN:
MAIN:
MAIN: ————————————-
MAIN: Reading Model: case.sif

Loading user function library: [HeatSolve]…[HeatSolver_Init0]
LoadMesh: Scaling coordinates: 1.000E-03 1.000E-03 1.000E-03

MAIN: ————————————-

Loading user function library: [HeatSolve]…[HeatSolver_Init]
Loading user function library: [HeatSolve]…[HeatSolver]
OptimizeBandwidth: ———————————————————
OptimizeBandwidth: Computing matrix structure for: heat equation…done.
OptimizeBandwidth: Half bandwidth without optimization: 8575
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization …done.
OptimizeBandwidth: Half bandwidth after optimization: 1478
OptimizeBandwidth: ———————————————————

MAIN:
MAIN: ————————————-
MAIN: Time: 1/14 86400.000000000000
MAIN: ————————————-
MAIN:
HeatSolve:
HeatSolve:
HeatSolve: ————————————-
HeatSolve: TEMPERATURE ITERATION 1
HeatSolve: ————————————-
HeatSolve:
HeatSolve: Starting Assembly…
HeatSolve: Assembly:

: …..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 120074
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.02

1 0.3895E-02
2 0.1909E-03
3 0.2327E-04
4 0.1216E-04
5 0.5040E-05
6 0.5204E-04
7 0.1487E-05
8 0.1616E-05
9 0.2218E-06
10 0.1219E-06
11 0.2561E-06
12 0.3507E-07
13 0.7409E-07
14 0.7608E-08
15 0.8945E-08
16 0.1223E-08
17 0.1234E-08
18 0.1603E-09
19 0.2121E-09
20 0.3031E-10
20 0.3031E-10
ComputeChange: NS (ITER=1) (NRM,RELC): ( 24.900009 2.0000000 ) :: heat equation
HeatSolve: iter: 1 Assembly: (s) 5.53 5.53
HeatSolve: iter: 1 Solve: (s) 0.16 0.16
HeatSolve: Result Norm : 24.900008502275778
HeatSolve: Relative Change : 2.0000000000000000
HeatSolve:
HeatSolve:
HeatSolve: ————————————-
HeatSolve: TEMPERATURE ITERATION 2
HeatSolve: ————————————-
HeatSolve:
HeatSolve: Starting Assembly…
HeatSolve: Assembly:

At line 2841 of file Lists.f90
Fortran runtime error: Bad real number in item 1 of list input

Hello courtois,

well … quite frankly, I didn’t really know what to choose from all these options. I don’t have Fortran experience. However, in my Makefile.inc, there seems to be no such option. Here is my Makefile.inc

SCOTCHDIR  = /usr
ISCOTCH    = -I$(SCOTCHDIR)/include/scotch 
LSCOTCH    = -L$(SCOTCHDIR)/lib -lptesmumps -lptscotch -lptscotcherr

LPORDDIR = $(topdir)/PORD/lib/
IPORD    = -I$(topdir)/PORD/include/
LPORD    = -L$(LPORDDIR) -lpord

# (parmetis is supplied by aur mpi vers):
LMETISDIR = /usr/lib
IMETIS = -I/usr/include
LMETIS    = -L$(LMETISDIR) -lparmetis -lmetis

ORDERINGSF = -Dptscotch -Dparmetis -Dpord 
ORDERINGSC  = $(ORDERINGSF)

LORDERINGS = $(LMETIS) $(LPORD) $(LSCOTCH)
IORDERINGSF = $(ISCOTCH)
IORDERINGSC = $(IMETIS) $(IPORD) $(ISCOTCH)

PLAT    =
RM = /bin/rm -f
CC = mpicc
FC = mpif90
FL = mpif90
AR = ar vr
RANLIB  = echo
# I don't have any blacs libs there, although scalapack installer reports that blacs is contained in libscalapack
SCALAP  = /opt/scalapack/lib/libscalapack.a 
# INCPAR = -I/usr/local/mpich/include
INCPAR = -I/usr/include
LIBPAR = $(SCALAP) # -L/usr/lib/openmpi/ -lmpi 
# See point 17 in the FAQ to have more details on the compilation of mpich with gfortran
INCSEQ = -I$(topdir)/libseq
LIBSEQ  =  -L$(topdir)/libseq -lmpiseq
#LIBBLAS = -L/usr/lib/xmm/ -lf77blas -latlas
# SE: (using non optimized lapack /blas)
LIBBLAS = -L/opt/scalapack/lib -lreflapack -lrefblas
LIBOTHERS = -lpthread
#Preprocessor defs for calling Fortran from C (-DAdd_ or -DAdd__ or -DUPPER)
CDEFS   = -DAdd_

#Begin Optimized options
OPTF    = -O -Dintel_ -DALLOW_NON_INIT 
OPTL    = -O 
OPTC    = -O
#End Optimized options
INC = $(INCPAR)
LIB = $(LIBPAR)
LIBSEQNEEDED =

I have:

$ mpif90 --showme
/usr/bin/gfortran -I/usr/include -pthread -I/usr/lib/openmpi -L/usr/lib/openmpi -lmpi_f90 -lmpi_f77 -lmpi -ldl -lhwloc
$ mpicc --showme
gcc -pthread -L/usr/lib/openmpi -lmpi -ldl -lhwloc

Furthermore, I compiled parmetis, metis and (pt)scotch_esmumps (as indicated by the paths above). If you could take a look at that and the config I supplied in my last post I’d be very grateful.

I also tried to install aster with the aid of setup.py (fresh install), but that always when I wanted to link against the openmpi version of mumps. I may open a new thread on this.

merci!

PS: If I recompile mumps now, I do not need to rebuild aster, correct? If no, do I need to do `make clean` or will a `make` suffice? I should probably try to systematically deselect parmetis, ptscotch…

Last edited by Anna (2012-05-07 19:41:18)

«the document you have requested to translate is already in English»

Moderators: arango, robertson