READ_FILE_MPGM_SYS.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_file_mpgm_sys(filec,&
                        SDOFnum,n_monitor_pgm,el_monitor_pgm,&
                        xr_monitor_pgm,yr_monitor_pgm,zr_monitor_pgm,axis_rot)
 
      implicit none
 
      character*70 :: filec
      character*100000 :: input_line
 
      integer*4 :: SDOFnum
      integer*4 :: i
      integer*4 :: trash
      integer*4 :: ileft,iright
      integer*4 :: status
 
      integer*4,dimension(SDOFnum) :: n_monitor_pgm
      integer*4,dimension(SDOFnum) :: el_monitor_pgm
 
      real*8,dimension(SDOFnum) :: xr_monitor_pgm
      real*8,dimension(SDOFnum) :: yr_monitor_pgm
      real*8,dimension(SDOFnum) :: zr_monitor_pgm
      real*8,dimension(SDOFnum) :: axis_rot
 
 
        open(20,file=filec)
        read(20,'(A)',iostat = status) input_line
        ileft = 1
        iright = len(input_line)
        read(input_line(ileft:iright),*)sdofnum
 
        do i = 1,sdofnum
                read(20,'(A)',iostat = status) input_line
                if (status.ne.0) exit
 
                ileft = 1
                iright = len(input_line)
                read(input_line(ileft:iright),*)trash,n_monitor_pgm(i),el_monitor_pgm(i),&
                            xr_monitor_pgm(i),yr_monitor_pgm(i),zr_monitor_pgm(i),axis_rot(i)
 
        enddo
                    
        close(20)
 
      return
 
      end subroutine read_file_mpgm_sys