READ_FILEXYZ.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_filexyz (filec, num_nodes, num_tria, x, y, z, node1, node2, node3, max_length)
	Reads mesh of triangular elements (topography or alluvial basin)

Function/Subroutine Documentation

read_filexyz()

subroutine read_filexyz	(	character*70	filec,
		integer*4	num_nodes,
		integer*4	num_tria,
		real*8, dimension(num_nodes)	x,
		real*8, dimension(num_nodes)	y,
		real*8, dimension(num_nodes)	z,
		integer*4, dimension(num_tria)	node1,
		integer*4, dimension(num_tria)	node2,
		integer*4, dimension(num_tria)	node3,
		real*8	max_length
	)

Reads mesh of triangular elements (topography or alluvial basin)

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	filec	file name for the mesh
[in]	num_nodes	number of nodes
[in]	num_tria	number of triangular fault
[out]	x	x-coordinate for the node
[out]	y	y-coordinate for the node
[out]	z	z-coordinate for the node
[out]	node1	1-vertex of the triangular fault
[out]	node2	2-vertex of the triangular fault
[out]	node3	3-vertex of the triangular fault
[out]	max_length	max diameter of the triangles

Definition at line 34 of file READ_FILEXYZ.f90.

 
      implicit none
 
      integer*4 :: num_nodes,num_tria
      integer*4 :: i
      integer*4 :: trash
      real*8, dimension(num_nodes) :: x,y,z
      integer*4, dimension(num_tria) :: node1,node2,node3
      character*70 :: filec
      integer*4 :: ileft,iright
      integer*4 :: status
      character*100000 :: input_line
      real*8 :: length_edge1,length_edge2,length_edge3
      real*8 :: max_edge_length,max_length
 
      max_length = 0.0d0
 
      open(20,file=filec)
      read(20,'(A)',iostat = status) input_line
      ileft = 1
      iright = len(input_line)
      read(input_line(ileft:iright),*) num_nodes,num_tria 
 
      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
 
      do i = 1,num_tria
        read(20,'(A)',iostat = status) input_line
        if (status.ne.0) exit
                ileft = 1
                iright = len(input_line)
                read(input_line(ileft:iright),*)trash,node1(i),node2(i),node3(i) 
      enddo 
                    
      close(20)
 
     ! compute max tria edge length
 
      do i = 1,num_tria
 
        length_edge1 = ( x(node1(i)) - x(node2(i)) )**2 + ( y(node1(i)) - y(node2(i)) )**2
        length_edge2 = ( x(node2(i)) - x(node3(i)) )**2 + ( y(node2(i)) - y(node3(i)) )**2
        length_edge3 = ( x(node3(i)) - x(node1(i)) )**2 + ( y(node3(i)) - y(node1(i)) )**2
 
        max_edge_length = max(length_edge1,length_edge2,length_edge3)
 
        if (max_length.lt.max_edge_length) then
                 max_length = max_edge_length
        endif
                
      enddo
       
        max_length = sqrt(max_length)
 
      return
 

Referenced by make_nothonoring().

Here is the caller graph for this function: