READ_FILEPG.f90

Go to the documentation of this file.

!    Copyright (C) 2012 The SPEED FOUNDATION
!    Author: Ilario Mazzieri
!
!    This file is part of SPEED.
!
!    SPEED is free software; you can redistribute it and/or modify it
!    under the terms of the GNU Affero General Public License as
!    published by the Free Software Foundation, either version 3 of the
!    License, or (at your option) any later version.
!
!    SPEED is distributed in the hope that it will be useful, but
!    WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
!    Affero General Public License for more details.
!
!    You should have received a copy of the GNU Affero General Public License
!    along with SPEED.  If not, see <http://www.gnu.org/licenses/>.
 
 
      subroutine read_filepg(filec,num_nodes, x,y,z)
 
      implicit none
 
      character*70 :: filec
      character*100000 :: input_line
 
      integer*4 :: i,trash, num_nodes
      integer*4 :: ileft,iright, status
      
      real*8 :: length_edge1,length_edge2,length_edge3
 
      real*8, dimension(num_nodes) :: x,y,z
 
 
      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),*)trash,x(i),y(i),z(i) 
      enddo
                    
      close(20)
 
      return
 
      end subroutine read_filepg