discretization.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/discretization.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/discretization.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module discretization
 
  use prec
  use results
 
  implicit none
  private
  
  integer,     public, protected :: &
    discretization_nip, &
    discretization_nelem
    
  integer,     public, protected, dimension(:),   allocatable :: &
    discretization_homogenizationat, &
    discretization_microstructureat   
 
  real(preal), public, protected, dimension(:,:), allocatable :: & 
    discretization_ipcoords0, &
    discretization_ipcoords, &
    discretization_nodecoords0, &
    discretization_nodecoords
    
  integer :: &
    discretization_sharednodesbegin
 
  public :: &
    discretization_init, &
    discretization_results, &
    discretization_setipcoords, &
    discretization_setnodecoords
 
contains
  
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine discretization_init(homogenizationAt,microstructureAt,&
                               IPcoords0,NodeCoords0,&
                               sharedNodesBegin)
 
  integer,     dimension(:),   intent(in) :: &
    homogenizationat, &
    microstructureat
  real(preal), dimension(:,:), intent(in) :: &
    ipcoords0, &
    nodecoords0
  integer, optional,           intent(in) :: &
    sharednodesbegin
 
  write(6,'(/,a)') ' <<<+-  discretization init  -+>>>'; flush(6)
 
  discretization_nelem = size(microstructureat,1)
  discretization_nip   = size(ipcoords0,2)/discretization_nelem
 
  discretization_homogenizationat = homogenizationat
  discretization_microstructureat = microstructureat  
 
  discretization_ipcoords0   = ipcoords0
  discretization_ipcoords    = ipcoords0
 
  discretization_nodecoords0 = nodecoords0
  discretization_nodecoords  = nodecoords0
  
  if(present(sharednodesbegin)) then
    discretization_sharednodesbegin = sharednodesbegin
  else
    discretization_sharednodesbegin = size(discretization_nodecoords0,2)
  endif
  
end subroutine discretization_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine discretization_results
 
  real(preal), dimension(:,:), allocatable :: u
  
  call results_closegroup(results_addgroup('current/geometry'))
  
  u = discretization_nodecoords(1:3,:discretization_sharednodesbegin) &
    - discretization_nodecoords0(1:3,:discretization_sharednodesbegin)
  call results_writedataset('current/geometry',u,'u_n','displacements of the nodes','m')
  
  u = discretization_ipcoords &
    - discretization_ipcoords0
  call results_writedataset('current/geometry',u,'u_p','displacements of the materialpoints','m')
 
end subroutine discretization_results
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine discretization_setipcoords(IPcoords)
 
  real(preal), dimension(:,:), intent(in) :: ipcoords
  
  discretization_ipcoords = ipcoords
 
end subroutine discretization_setipcoords
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine discretization_setnodecoords(NodeCoords)
 
  real(preal), dimension(:,:), intent(in) :: nodecoords
  
  discretization_nodecoords = nodecoords
 
end subroutine discretization_setnodecoords
 
 
end module discretization