MAKE_NLE.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_nle (depth, nn, yon, nnt, zs_elev, cs_nnz, cs, ielem, ncase, tcase, zs_all, vs_tria)
	Find non-linear nodes.

Function/Subroutine Documentation

make_nle()

subroutine make_nle	(	real*8	depth,
		integer*4	nn,
		integer*4, dimension(nn,nn,nn)	yon,
		integer*4	nnt,
		real*8, dimension(nnt)	zs_elev,
		integer*4	cs_nnz,
		integer*4, dimension(0:cs_nnz)	cs,
		integer*4	ielem,
		integer*4	ncase,
		integer*4	tcase,
		real*8, dimension(nnt)	zs_all,
		real*8, dimension(nnt)	vs_tria
	)

Find non-linear nodes.

Author

Ilario Mazzieri

Date

September, 2015

Version

1.0

Parameters

[in]	Depth	depth for searching nodes
[in]	nn	polynomial degree
[in]	nnt	number of nodes
[in]	zs_elv	nodes elevation from surface
[in]	cs_nnz	length of cs
[in]	cs	connectivity vector
[in]	ielem	element index
[in]	ncase	id for case
[in]	tcase	label for case
[in]	zs_all	nodes elevation from alluvial
[in]	vs_tria	vs30 values
[out]	yon	1 if the node is non-linear elastic, 0 otherwise

Definition at line 39 of file MAKE_NLE.f90.

      
   
     
      implicit none
                                                              
      integer*4 :: nn
      integer*4 :: p,q,r
      integer*4 :: nnt
      integer*4 :: cs_nnz                                
      integer*4 :: ncase
      integer*4 :: tcase                                                                        
      integer*4 :: is,in,ielem
 
      integer*4, dimension(0:cs_nnz) :: cs                
 
      integer*4, dimension(nn,nn,nn) :: yon
 
      real*8 :: depth
          
      real*8, dimension(nnt) :: zs_elev
      real*8, dimension(nnt) :: zs_all
      real*8, dimension(nnt) :: vs_tria                 
      
!     STRESS CALCULATION
      
      do r = 1,nn
         do q = 1,nn
            do p = 1,nn
 
                yon(p,q,r) = 1
 
                is = nn*nn*(r -1) +nn*(q -1) +p
               in = cs(cs(ielem -1) +is)
               if (ncase.gt.0 .and. tcase .ne. 16) then
                        if (depth.lt.zs_elev(in) .or. zs_elev(in) .lt. 0.d0 ) yon(p,q,r) = 0
                        if ((tcase .gt. 0) .and. (zs_all(in) .lt. 0.0d0)) yon(p,q,r) = 0
               elseif (ncase.gt.0 .and. tcase .eq. 16) then
                        if (depth .lt. zs_elev(in) .or. zs_elev(in) .lt. 0.d0 ) yon(p,q,r) = 0
                        if (vs_tria(in) .ge. 450.0d0) yon(p,q,r) = 0
               endif        
 
            enddo
         enddo
      enddo
      
      return