READ_FILEMESH_LOC.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_filemesh_loc (file_grid, nn_loc, xx_loc, yy_loc, zz
	Reads the local mesh.

Function/Subroutine Documentation

read_filemesh_loc()

subroutine read_filemesh_loc	(	character*70	file_grid,
		integer*4	nn_loc,
		real*8, dimension(nn_loc), intent(inout)	xx_loc,
		real*8, dimension(nn_loc), intent(inout)	yy_loc,
		real*8	zz
	)

Reads the local mesh.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	file_grid	file name
[in]	nn_loc	number of local nodes
[in]	loc_n_num	local node numbering
[out]	xx_loc	x-coordinate of local nodes
[out]	yy_loc	y-coordinate of local nodes
[out]	zz_loc	z-coordinate of local nodes

Definition at line 30 of file READ_FILEMESH_LOC.f90.

                            loc_n_num)
 
      implicit none
      
      character*70 :: file_grid
      character*110 :: input_line
 
      integer*4 :: nn_loc
      integer*4 :: inode
      integer*4 :: nnode,nelem,i,ic,status
 
      integer*4, dimension(nn_loc) :: loc_n_num
 
      real*8 :: xx,yy,zz      
      real*8, dimension(nn_loc), intent(inout) :: xx_loc(nn_loc), yy_loc(nn_loc), zz_loc(nn_loc)
      
      inode = 0
      status = 0 
      
      open(23,file=file_grid)
      
      do 
         read(23,'(A)') input_line
         if (input_line(1:1) .ne. '#') exit
      enddo
      
      read(input_line,*) nnode, nelem
      
      do i = 1,nnode
        read(23,*) inode, xx, yy, zz
        if (inode.ne.i) then
          status = 1
        endif
        
        call get_indloc_from_indglo(loc_n_num, nn_loc, inode, ic)
        
        
        if (ic .ne. 0) then
          xx_loc(ic) = xx
          yy_loc(ic) = yy
          zz_loc(ic) = zz
        endif
        
      enddo
      
      
      close(23)
      
      return
      

References get_indloc_from_indglo().

Referenced by make_spx_grid_with_loc_numeration().

Here is the call graph for this function:

Here is the caller graph for this function: