READ_FILE_MPGM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_file_mpgm (filec, nmonitors_pgm, n_monitor_pgm, el_monitor_pgm, xr_monitor_pgm, yr_monitor_pgm, zr_monitor_pgm)
	Reads files such as MLST.input, MLST.position or MPGM.input.

Function/Subroutine Documentation

read_file_mpgm()

subroutine read_file_mpgm	(	character*70	filec,
		integer*4	nmonitors_pgm,
		integer*4, dimension(nmonitors_pgm)	n_monitor_pgm,
		integer*4, dimension(nmonitors_pgm)	el_monitor_pgm,
		real*8, dimension(nmonitors_pgm)	xr_monitor_pgm,
		real*8, dimension(nmonitors_pgm)	yr_monitor_pgm,
		real*8, dimension(nmonitors_pgm)	zr_monitor_pgm
	)

Reads files such as MLST.input, MLST.position or MPGM.input.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	filec	file name to read
[in]	nmonitors_pgm	number of monitors
[out]	n_monitor_pgm	id monitor number (id spectral node)
[out]	el_monitor_pgm	id monitor element (id grid element)
[out]	xr_monitor_pgm	x-coordinate for monitor
[out]	yr_monitor_pgm	y-coordinate for monitor
[out]	zr_monitor_pgm	z-coordinate for monitor

Definition at line 31 of file READ_FILE_MPGM.f90.

 
      implicit none
 
      character*70 :: filec
      character*100000 :: input_line
 
      integer*4 :: nmonitors_pgm
      integer*4 :: i
      integer*4 :: trash
      integer*4 :: ileft,iright
      integer*4 :: status
 
      integer*4,dimension(nmonitors_pgm) :: n_monitor_pgm
      integer*4,dimension(nmonitors_pgm) :: el_monitor_pgm
 
      real*8,dimension(nmonitors_pgm) :: xr_monitor_pgm
      real*8,dimension(nmonitors_pgm) :: yr_monitor_pgm
      real*8,dimension(nmonitors_pgm) :: zr_monitor_pgm
 
 
              open(20,file=filec)
                read(20,'(A)',iostat = status) input_line
                ileft = 1
                iright = len(input_line)
                read(input_line(ileft:iright),*)nmonitors_pgm
 
                do i = 1,nmonitors_pgm
                        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)
 
                enddo
                    
        close(20)
 
      return
 

Referenced by find_monitor_position().

Here is the caller graph for this function: