Functions/Subroutines
subroutine	make_nh_enhanced ()
	...Not-Honoring Enhanced (NHE) Implementation

Function/Subroutine Documentation

◆ make_nh_enhanced()

subroutine make_nh_enhanced

...Not-Honoring Enhanced (NHE) Implementation

Author: Srihari Sangaraju

Date: August, 2020

Version

1.0 @ Detailed Explanation

User Can specify that a single block/multiple blocks in mesh has to use material properties as per provided Tomography text file
For the each GLL node in NHE-specified blocks, this subroutine assigns Mechanical properties of the corresponding nearest point it can find in "Tomography text file"
For the nodes in other blocks It assigns Mechanical properties specified in *.mate file Note: In this version, This subroutine will be run, only if atleast 1 block is specified to use NHE in *.mate file

Parameters

[in]	loc_n_num.	Global node number of 'i'th local node is loc_n_num(i)
[in]	nn_loc.	No. of nodes in Local/Current Partition
[in]	nmat_nhe	No. of Blocks specified with NHE case
[in]	nhe_mat	Tag/Labels of Blocks where NHE has to be implemented
[in]	xs_loc	x-coordinate of spectral nodes
[in]	ys_loc	y-coordinate of spectral nodes
[in]	zs_loc	z-coordinate of spectral nodes
[out]	lambda_nhe	Lame coefficient lambda
[out]	mu_nhe	Lame coefficient mu
[out]	rho_nhe	material density
[out]	Qs_nhe	quality factor for S-waves (Standard Linear Solid damping)
[out]	Qp_nhe	quality factor for P-waves (Standard Linear Solid damping)
[out]	gamma_dmp_nhe	damping coefficient gamma (Kosloff&Kosloff)

Definition at line 52 of file MAKE_NH_Enhanced.f90.

 
        use speed_par
 
        implicit none
 
      include 'SPEED.MPI'
 
      integer*4 :: count
      integer*4, dimension(:), allocatable :: node_nhe_flag, NN_src_ind_loc
 
 
        if (mpi_id.eq.0) write(*,'(A)')
      if (mpi_id.eq.0) write(*,'(A)')'---------------Setup Not-Honoring Enhanced ---------------' 
 
 
      allocate(node_nhe_flag(nnod_loc))
      call make_nh_enhanced_initialise(nnod_loc, nmat_nhe, val_nhe, &
                             nmat, tag_mat, nelem_loc, con_nnz_loc, con_spx_loc, &
                             xx_spx_loc, yy_spx_loc, zz_spx_loc, &
                             count, &
                             node_nhe_flag, mpi_id, mpi_comm, mpi_file)
 
 
      allocate(nn_src_ind_loc(count))
      call make_nh_enhanced_nnsearch(nnod_loc, count, mpi_id, mpi_np, mpi_comm, &
                                    mpi_file, nn_src_ind_loc)
 
 
      allocate(lambda_nhe(nnod_loc), mu_nhe(nnod_loc), rho_nhe(nnod_loc))
      allocate(qs_nhe_el(nelem_loc), qp_nhe_el(nelem_loc)) ! gamma_nhe_el(nelem_loc))
 
      call make_nh_enhanced_assign_prop(nnod_loc, nmat, prop_mat, sdeg_mat, &
                                           nelem_loc, con_nnz_loc, con_spx_loc, &
                             node_nhe_flag, count, nn_src_ind_loc, qs, qp, &
                             lambda_nhe, mu_nhe, rho_nhe, qs_nhe_el, qp_nhe_el, &
                             mpi_id, mpi_comm)
 
 
      deallocate(node_nhe_flag, nn_src_ind_loc)
 
      if (mpi_id.eq.0) write(*,'(A)')'Completed.' 
 

References speed_par::con_nnz_loc, speed_par::con_spx_loc, speed_par::lambda_nhe, make_nh_enhanced_assign_prop(), make_nh_enhanced_initialise(), make_nh_enhanced_nnsearch(), speed_par::mpi_comm, speed_par::mpi_file, speed_par::mpi_id, speed_par::mpi_np, speed_par::mu_nhe, speed_par::nelem_loc, speed_par::nmat, speed_par::nmat_nhe, speed_par::nnod_loc, speed_par::prop_mat, speed_par::qp, speed_par::qp_nhe_el, speed_par::qs, speed_par::qs_nhe_el, speed_par::rho_nhe, speed_par::sdeg_mat, speed_par::tag_mat, speed_par::val_nhe, speed_par::xx_spx_loc, speed_par::yy_spx_loc, and speed_par::zz_spx_loc.

Here is the call graph for this function:

Functions/Subroutines

Function/Subroutine Documentation

◆ make_nh_enhanced()