GET_MONITOR_VALUE.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_monitor_value (nb_nod, xq, val, xref, yref, zref, re
	Computes the mean value.

Function/Subroutine Documentation

get_monitor_value()

subroutine get_monitor_value	(	integer*4	nb_nod,
		real*8, dimension(nb_nod)	xq,
		real*8, dimension(nb_nod,nb_nod,nb_nod)	val,
		real*8	xref,
		real*8	yref,
		real*8	zref,
			re
	)

Computes the mean value.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_nod	number of 1D GLL nodes
[in]	xq	coor GLL node
[in]	val	value in each GLL node
[in]	xref	x-coordinate of the monitor
[in]	yref	y-coordinate of the monitor
[in]	zref	z-coordinate of the monitor
[out]	res	mean value (interpolation)

Definition at line 32 of file GET_MONITOR_VALUE.f90.

            
      implicit none
      
      integer*4 :: nb_nod
      integer*4 :: i,j,k,h
      
      real*8 :: xref,yref,zref,res
      real*8 :: fx,fy,fz
 
      real*8, dimension(nb_nod) :: xq
 
      real*8, dimension(nb_nod,nb_nod,nb_nod) :: val
 
      res = 0.0d0
      
      do k = 1,nb_nod
         fz = 1.0d0
         do h = 1,nb_nod
            if (h.ne.k) fz = fz * (zref - xq(h)) / (xq(k) - xq(h))
         enddo
         
         do j = 1,nb_nod
            fy = 1.0d0
            do h = 1,nb_nod
               if (h.ne.j) fy = fy * (yref - xq(h)) / (xq(j) - xq(h))
            enddo
            
            do i = 1,nb_nod
               fx = 1.0d0
               do h = 1,nb_nod
                  if (h.ne.i) fx = fx * (xref - xq(h)) / (xq(i) - xq(h))
               enddo
               
               res = res + val(i,j,k) * fx * fy * fz
            enddo
         enddo
      enddo
      

Referenced by compute_sdof_input(), get_max_values(), and write_output().

Here is the caller graph for this function: