static_echo.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
module m_static_echo
contains
  !C
  !C***
  !C***
  !C
  subroutine fstr_echo (hecMESH)
    use m_fstr
    use m_hecmw2fstr_mesh_conv
    implicit none
    type (hecmwST_local_mesh):: hecMESH
 
    !** Local variables
    integer(kind=kint) :: i,j,is,iE,nn,ic_type
    integer(kind=kint) :: icel,isect,ig1,iS0,iE0,ik
    integer(kind=kint) :: in,nid, itype
    real(kind=kreal) :: x,y,z
    integer(kind=kint):: nids(20)
    !C====
    !C +-------------------------------+
    !C | NODE                          |
    !C +-------------------------------+
    !C===
    !C** nodal coordinate
    write(ilog,*) '### Number of nodes',hecmesh%n_node
    write(ilog,*) 'ID X Y Z'
    do i=1,hecmesh%n_node
      nid = hecmesh%global_node_ID(i)
      x = hecmesh%node(3*i-2)
      y = hecmesh%node(3*i-1)
      z = hecmesh%node(3*i)
      write(ilog,'(i8,e15.5,e15.5,e15.5)') nid,x,y,z
    enddo
    !C
    !C +-------------------------------+
    !C | ELEMENT                       |
    !C +-------------------------------+
    !C===
    call fstr2hecmw_mesh_conv( hecmesh )
    write(ilog,*) '### Elements', hecmesh%n_elem
    do itype= 1, hecmesh%n_elem_type
      is= hecmesh%elem_type_index(itype-1) + 1
      ie= hecmesh%elem_type_index(itype  )
      ic_type= hecmesh%elem_type_item(itype)
      !C** Set number of nodes
      nn = hecmw_get_max_node(ic_type)
      !C** element loop
      do icel= is, ie
        !C** node ID
        is= hecmesh%elem_node_index(icel-1)
        do j=1,nn
          if( hecmesh%n_refine > 0 ) then
            nids(j)= hecmesh%elem_node_item (is+j)
          else
            nids(j)= hecmesh%global_node_ID( hecmesh%elem_node_item (is+j))
          endif
        enddo
        !C** section  ID
        isect= hecmesh%section_ID(icel)
        write(ilog,*) '### Element ID=',ic_type,isect,hecmesh%global_elem_id(icel)
        write(ilog,*) (nids(j),j=1,nn)
      enddo
    enddo
    call hecmw2fstr_mesh_conv( hecmesh )
    !C +-------------------------------+
    !C | NODE GROUP                    |
    !C +-------------------------------+
    write(ilog,*) '### Ngroup'
    do ig1= 1, hecmesh%node_group%n_grp
      write(ilog,*)
      write(ilog,'(a80)') hecmesh%node_group%grp_name(ig1)
      is0= hecmesh%node_group%grp_index(ig1-1) + 1
      ie0= hecmesh%node_group%grp_index(ig1  )
      do ik= is0, ie0
        in   = hecmesh%node_group%grp_item(ik)
        write(ilog,*) hecmesh%global_node_ID(in)
      enddo
    enddo
    !C +-------------------------------+
    !C | ELEMENT GROUP                 |
    !C +-------------------------------+
    write(ilog,*) '### Egroup'
    do ig1= 1, hecmesh%elem_group%n_grp
      write(ilog,*)
      write(ilog,'(a80)') hecmesh%elem_group%grp_name(ig1)
      is0= hecmesh%elem_group%grp_index(ig1-1) + 1
      ie0= hecmesh%elem_group%grp_index(ig1  )
      do ik= is0, ie0
        in   = hecmesh%elem_group%grp_item(ik)
        write(ilog,*) hecmesh%global_elem_ID(in)
      enddo
    enddo
    write(ilog,*) '### Reftemp',ref_temp
    !C====
  end subroutine fstr_echo
end module m_static_echo