READ_FILEOUT.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_fileout (filename, counter, procs, nb_vec, vec)
	Reads files *.out for the restart of the simulation.

Function/Subroutine Documentation

read_fileout()

subroutine read_fileout	(	character*70	filename,
		integer*4	counter,
		integer*4	procs,
		integer*4	nb_vec,
		real*8, dimension(nb_vec)	vec
	)

Reads files *.out for the restart of the simulation.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	filename	file name
[in]	counter	index for snapshot
[in]	procs	MPI process id
[in]	nb_vec	number of values
[out]	vec	values read

Definition at line 29 of file READ_FILEOUT.f90.

      
      implicit none
      
      integer*4 :: counter,procs,nb_vec
      real*8, dimension(nb_vec) :: vec
      character*70 :: filename
      
      character*70 :: out_file
      integer*4 :: i,lname
      
      lname = len_trim(filename)
      out_file = filename(1:lname) // '000000_000000.out'
      
      if (procs .lt. 10) then
         write(out_file(lname+6:lname+6),'(i1)') procs
      else if (procs .lt. 100) then
         write(out_file(lname+5:lname+6),'(i2)') procs
      else if (procs .lt. 1000) then
         write(out_file(lname+4:lname+6),'(i3)') procs     
      else if (procs .lt. 10000) then
         write(out_file(lname+3:lname+6),'(i4)') procs      
      else if (procs .lt. 100000) then
         write(out_file(lname+2:lname+6),'(i5)') procs
      else
         write(out_file(lname+1:lname+6),'(i6)') procs
      endif
      
      if (counter .lt. 10) then
         write(out_file(lname+13:lname+13),'(i1)') counter
      else if (counter .lt. 100) then
         write(out_file(lname+12:lname+13),'(i2)') counter
      else if (counter .lt. 1000) then
         write(out_file(lname+11:lname+13),'(i3)') counter    
      else if (counter .lt. 10000) then
         write(out_file(lname+10:lname+13),'(i4)') counter    
      else if (counter .lt. 100000) then
         write(out_file(lname+9:lname+13),'(i5)') counter    
      else
         write(out_file(lname+8:lname+13),'(i6)') counter
      endif
      
      open(20+procs, file=out_file)
      
      
      
      
      do i = 1,nb_vec
         read(20+procs,*) vec(i)
      enddo
      
      close(20+procs)
      
      return
 
      
      return