READ_FACES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_faces (filename, fac, nodes, n, mpi_id, mpi_com)
	Reads FACES.input.

Function/Subroutine Documentation

read_faces()

subroutine read_faces	(	character*70	filename,
		integer*4, dimension(3,n), intent(inout)	fac,
		real*8, dimension(25,n), intent(inout)	nodes,
		integer*4	n,
		integer*4	mpi_id,
		integer*4	mpi_com
	)

Reads FACES.input.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	filename	file name
[in]	n	number of DG faces
[in]	mpi_id	mpi processor identity
[in]	mpi_com	mpi communicator
[out]	fac	material, element, face for a DG face
[out]	nodes	constants for the bilienar map from (-1,1)^3 to the current element

Definition at line 30 of file READ_FACES.f90.

 
     implicit none
 
     character*70 :: filename
 
     integer*4 :: j, n, mpi_id, mpi_com, error, unitname
 
     integer*4, dimension(3,n), intent(inout) :: fac
 
     real*8, dimension(25,n), intent(inout) :: nodes
    
     unitname = 400 + mpi_id
     open(unitname, file=filename)
     do j = 1,n
         read(unitname,*) &
                         fac(1,j), fac(2,j), fac(3,j), &
                         nodes(1,j), nodes(2,j), nodes(3,j), & 
                         nodes(4,j), nodes(5,j), nodes(6,j), & 
                         nodes(7,j), nodes(8,j), nodes(9,j), & 
                         nodes(10,j), nodes(11,j), nodes(12,j), & 
                         nodes(13,j), nodes(14,j), nodes(15,j), & 
                         nodes(16,j), nodes(17,j), nodes(18,j), & 
                         nodes(19,j), nodes(20,j), nodes(21,j), & 
                         nodes(22,j), nodes(23,j), nodes(24,j), & 
                         nodes(25,j) 
 
      enddo
      
      close(unitname)
 
 
      return