READ_FILESYS.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_filesys (filec, num_nodes, label, x, y, z, axis_rot)
	Reads SYS.input file.

Function/Subroutine Documentation

read_filesys()

subroutine read_filesys	(	character*70	filec,
		integer*4	num_nodes,
		integer*4, dimension(num_nodes)	label,
		real*8, dimension(num_nodes)	x,
		real*8, dimension(num_nodes)	y,
		real*8, dimension(num_nodes)	z,
		real*8, dimension(num_nodes)	axis_rot
	)

Reads SYS.input file.

Author

Aline Herlin, Srihari Sangaraju

Date

February, 2021

Version

1.0

Parameters

[in]	filec	file name
[in]	num_nodes	number of oscillators
[out]	label	system ID
[out]	x,y,z	coord. of the structure

Definition at line 28 of file READ_FILESYS.f90.

 
  implicit none
 
  character*70 :: filec
  character*100000 :: input_line
  integer*4 :: i, num_nodes, ileft, iright, status
  real*8, dimension(num_nodes) :: x, y, z, axis_rot
  integer*4, dimension(num_nodes) :: label
  real*8 :: dummy, pi
 
  pi = 4.0*atan(1.d0);
 
  open(20,file=filec)
  read(20,'(A)',iostat = status) input_line
  ileft = 1
  iright = len(input_line)
  read(input_line(ileft:iright),*) num_nodes
 
  do i = 1,num_nodes
    read(20,'(A)',iostat = status) input_line
    if (status.ne.0) exit
    ileft = 1
    iright = len(input_line)
    read(input_line(ileft:iright),*) label(i), x(i), y(i), z(i), dummy
    axis_rot(i) = dummy*pi/180.d0
  enddo
 
  close(20)
 
  return
 

Referenced by read_system_position().

Here is the caller graph for this function: