UPDATE_DAMPING_TENSOR.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	update_damping_tensor (nnod_loc, n_sls, frequency_range,
	Update the damping tensor with the Heun method.

Function/Subroutine Documentation

update_damping_tensor()

subroutine update_damping_tensor	(	integer*4	nnod_loc,
		integer*4	n_sls,
		real*8, dimension(n_sls)	frequency_range
	)

Update the damping tensor with the Heun method.

Author

Ilario Mazzieri

Date

November, 2014

Version

1.0

Parameters

[in]	nnod_loc	number of local nodes
[in]	N_SLS	number of standard linear solids
[in]	frequency_range	frequency range for SLS model
[in]	strain	strain tensor
[in]	deltat	time step
[out]	strain_visc	viscoelastic tensor

Definition at line 31 of file UPDATE_DAMPING_TENSOR.f90.

                                  strain_visc, strain, deltat)
 
      
      implicit none
                                                
 
      integer*4 :: nnod_loc, N_SLS, i
                         
      
      real*8 :: deltat, frequency_range(n_sls), strain_visc(6*nnod_loc,n_sls), &
                strain(6*nnod_loc), temp(6*nnod_loc, n_sls)
 
 
      
      do i = 1, n_sls
         !EXPLICIT EULER
         !temp(:,i) = strain_visc(:,i) + deltat*frequency_range(i)*(strain - strain_visc(:,i))
         !EXPLICIT HEUN
         temp(:,i) = strain_visc(:,i) + &