hecmw_adapt_active.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
 
subroutine hecmw_adapt_active (hecMESH)
  !C
  !C***
  !C*** hecmw_adapt_active
  !C***
  !C
  !C    FIND active NODEs, EDGEs, ELEMs
  !C
  use  hecmw_util
  type (hecmwST_local_mesh) :: hecMESH
  !C
  !C-- arrays
  allocate (hecmesh%adapt_edge_home(hecmesh%n_adapt_edge))
  allocate (hecmesh%adapt_act_edge (hecmesh%n_adapt_edge))
 
  hecmesh%adapt_edge_home= hecmesh%my_rank
  hecmesh%adapt_act_edge = 0
 
  !C
  !C-- EDGE home
  icou= 0
  do ie= 1, hecmesh%n_adapt_edge
    in1= hecmesh%adapt_edge_node(2*ie-1)
    in2= hecmesh%adapt_edge_node(2*ie  )
    ip1= hecmesh%node_ID(2*in1)
    ip2= hecmesh%node_ID(2*in2)
    hecmesh%adapt_edge_home(ie)= min(ip1,ip2)
 
    if ((ip1.eq.hecmesh%my_rank.and.ip2.eq.hecmesh%my_rank).or.     &
        &      (ip1.eq.hecmesh%my_rank.and.ip1.lt.ip2            ).or.     &
        &      (ip2.eq.hecmesh%my_rank.and.ip2.lt.ip1    )) then
      icou = icou + 1
      hecmesh%adapt_act_edge(icou)= ie
    endif
  enddo
  hecmesh%n_adapt_act_edge= icou
 
  return
end