utilities/RocfracPrep/meshdata.f90

Go to the documentation of this file.

!*********************************************************************
!* Illinois Open Source License                                      *
!*                                                                   *
!* University of Illinois/NCSA                                       * 
!* Open Source License                                               *
!*                                                                   *
!* Copyright@2008, University of Illinois.  All rights reserved.     *
!*                                                                   *
!*  Developed by:                                                    *
!*                                                                   *
!*     Center for Simulation of Advanced Rockets                     *
!*                                                                   *
!*     University of Illinois                                        *
!*                                                                   *
!*     www.csar.uiuc.edu                                             *
!*                                                                   *
!* Permission is hereby granted, free of charge, to any person       *
!* obtaining a copy of this software and associated documentation    *
!* files (the "Software"), to deal with the Software without         *
!* restriction, including without limitation the rights to use,      *
!* copy, modify, merge, publish, distribute, sublicense, and/or      *
!* sell copies of the Software, and to permit persons to whom the    *
!* Software is furnished to do so, subject to the following          *
!* conditions:                                                       *
!*                                                                   *
!*                                                                   *
!* @ Redistributions of source code must retain the above copyright  * 
!*   notice, this list of conditions and the following disclaimers.  *
!*                                                                   * 
!* @ Redistributions in binary form must reproduce the above         *
!*   copyright notice, this list of conditions and the following     *
!*   disclaimers in the documentation and/or other materials         *
!*   provided with the distribution.                                 *
!*                                                                   *
!* @ Neither the names of the Center for Simulation of Advanced      *
!*   Rockets, the University of Illinois, nor the names of its       *
!*   contributors may be used to endorse or promote products derived * 
!*   from this Software without specific prior written permission.   *
!*                                                                   *
!* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
!* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
!* OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
!* NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
!* COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       * 
!* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
!* ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
!* USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
!*********************************************************************
!* Please acknowledge The University of Illinois Center for          *
!* Simulation of Advanced Rockets in works and publications          *
!* resulting from this software or its derivatives.                  *
!*********************************************************************
MODULE meshdata


  USE generic_list, ONLY : link_ptr_type,link_type,list_type
  USE generic_list, ONLY : li_init_list,li_add_to_head,li_get_head, &
       li_remove_head,li_get_next,li_associated,li_get_len

  IMPLICIT NONE
  
  INTEGER :: io_input       ! unit id of control deck
  parameter(io_input = 10)

!  - DECK FILE PRAMETERS
  CHARACTER*20 :: prefx   ! I/O file prefx
  INTEGER :: prefx_lngth  ! length of prefx
!  - Fastest Dilatitial wave speed
  REAL*8 :: cd_fastest
!  - Derive Type for the boundary conditions
  TYPE bcvalues
     INTEGER :: b1, b2, b3
     REAL*8 :: bc1, bc2, bc3
  END TYPE bcvalues
!  - Stores the boundary conditions data values
  TYPE(bcvalues), DIMENSION(32) :: bc_conditions

! PRIMARY MESH DATA (i.e. mesh created with ansys, patran, or mesh3d):

!   - Number of nodes 
  INTEGER :: numnp_prmry
!   - Number of volumetric elements  
  INTEGER :: numelv_prmry
!   - meshing software flag
  INTEGER :: iansys         ! 0- no, 1-yes
  INTEGER :: ipatran        ! 0- no, 1-yes
  INTEGER :: itetmesh       ! 0- no, 1-yes
  INTEGER :: ipatcohin      ! 0- no, 1-yes
  INTEGER :: itetcohin      ! 0- no, 1-yes
  
  INTEGER :: numvertx
  
  REAL*8 :: convertunit
      

!   - Coordinate array
  REAL*8, ALLOCATABLE, DIMENSION(:,:) :: coor
!   - Coordinate array
  REAL*8, ALLOCATABLE, DIMENSION(:) :: press_nodal
!   - Connectivity array for volumetric elements
  INTEGER, ALLOCATABLE, DIMENSION(:,:) :: lmelv_prmry


!   - Nodal corner boundary flag
  INTEGER, ALLOCATABLE, DIMENSION(:) :: ibcaxi
!   - Partioned element's processor (from METIS)
  INTEGER, ALLOCATABLE, DIMENSION(:) :: epart

!   - Number of surface triangles with applied tractions
!  INTEGER, ALLOCATABLE, DIMENSION(:,:) :: numel_2d
!   - Number of surface nodes with applied mesh motion velocity
  INTEGER, ALLOCATABLE, DIMENSION(:) :: numnp_2d
  INTEGER, ALLOCATABLE, DIMENSION(:) :: numelperproc
  INTEGER, ALLOCATABLE, DIMENSION(:) :: numndperproc
! sub
  INTEGER, DIMENSION(:,:), ALLOCATABLE :: lmtri
  INTEGER, DIMENSION(:), ALLOCATABLE :: elm_2d
  INTEGER, DIMENSION(:), ALLOCATABLE :: elm_2d_flag
  INTEGER, ALLOCATABLE, DIMENSION(:,:) :: neigh_2d
  INTEGER, DIMENSION(:), ALLOCATABLE :: epart_2d
! sub

!   - Node tracking
  INTEGER :: numnodeio
  INTEGER, ALLOCATABLE, DIMENSION(:) :: nodeio

!   - Node 

  INTEGER :: ioformat     ! 0 = unformatted (default), 1 = formatted 
        ! 0 = no output, 1 = output
  INTEGER :: iopmvis  ! PMVIS software to view the partition 0=no, 1 = yes
  INTEGER :: ipress  ! flag for if we have pressure loading 0=no, 1 = yes

  INTEGER :: intfaceflag

!
! -- derived surface type
!
!$$$        TYPE request
!$$$            INTEGER, DIMENSION(1:6) :: Conn ! list of nodes
!$$$            INTEGER :: NumEl
!$$$            INTEGER :: NumNp
!$$$            TYPE(request), POINTER :: next
!$$$        END TYPE request
!$$$      
!$$$        TYPE(request), POINTER :: head, tail
!$$$        TYPE(request), POINTER :: item, first

  INTEGER, ALLOCATABLE, DIMENSION(:) :: indsburnflg
  
  INTEGER :: nummat
  
  INTEGER, DIMENSION(:), ALLOCATABLE :: matid
  INTEGER, DIMENSION(:), ALLOCATABLE :: eltypeid

      ! limit of 10 nodes to monitor history

  INTEGER, DIMENSION(1:10) :: ndhistory
  INTEGER :: numndhistory

  INTEGER :: numbc_prmry
  INTEGER :: numbc_prmry_mm
  INTEGER :: numbc_prmry_ht
!----- dimension:  <1>  old node number   <2> processor
  INTEGER, ALLOCATABLE, DIMENSION(:,:) :: ik1

! User-defined list element
! The Link_Type field MUST be the FIRST in the user-defined list element 
! Note pointer to data so as to easily create sublists

  TYPE user_type_bc
     TYPE(link_type) :: link
     TYPE(user_data_type_bc), POINTER :: data  
  END TYPE user_type_bc
  
  TYPE user_data_type_bc
     INTEGER :: bc_nodeglb
     INTEGER :: bc_flagglb
  END TYPE user_data_type_bc
  
! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_bc
     TYPE(user_type_bc), POINTER :: p
  END TYPE user_ptr_type_bc
  
  TYPE(user_ptr_type_bc)  :: user_bc

  TYPE(list_type) :: bc_structural
  TYPE(list_type) :: bc_meshmotion
  TYPE(list_type) :: bc_thermal
  
  TYPE user_type_mapnd
     TYPE(link_type) :: link
     TYPE(user_data_type_mapnd), POINTER :: data  
  END TYPE user_type_mapnd
  
  TYPE user_data_type_mapnd
     INTEGER :: locnd
     INTEGER :: proc
  END TYPE user_data_type_mapnd

! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_mapnd
     TYPE(user_type_mapnd), POINTER :: p
  END TYPE user_ptr_type_mapnd
  
  TYPE(user_ptr_type_mapnd)  :: user_mapnd
  
  TYPE(list_type), ALLOCATABLE, DIMENSION(:) :: mapnd_glb2locproc
  
  INTEGER, ALLOCATABLE, DIMENSION(:) :: numbc_structural, numbc_meshmotion, numbc_thermal
  
  
  INTEGER, PARAMETER :: maxnumberofprocstosharenode = 8
  
  INTEGER, ALLOCATABLE, DIMENSION(:,:) :: procndlist
  INTEGER, ALLOCATABLE, DIMENSION(:) ::  numprocpernd
  
  INTEGER, ALLOCATABLE, DIMENSION(:) :: nodeflag

! User-defined list element
! The Link_Type field MUST be the FIRST in the user-defined list element 
! Note pointer to data so as to easily create sublists

  TYPE user_type_surfmesh_tri3
     TYPE(link_type) :: link
     TYPE(user_data_type_surfmesh_tri3), POINTER :: data  
  END TYPE user_type_surfmesh_tri3

  TYPE user_data_type_surfmesh_tri3
     INTEGER, DIMENSION(1:4) :: elemdata
  END TYPE user_data_type_surfmesh_tri3
  
! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_surfmesh_tri3
     TYPE(user_type_surfmesh_tri3), POINTER :: p
  END TYPE user_ptr_type_surfmesh_tri3
  
  TYPE(user_ptr_type_surfmesh_tri3)  :: user_surfmesh_tri3
  
  TYPE(list_type) :: surfmesh_tri3_s
  TYPE(list_type) :: surfmesh_tri3_sf
  
! User-defined list element
! The Link_Type field MUST be the FIRST in the user-defined list element 
! Note pointer to data so as to easily create sublists

  TYPE user_type_surfmesh_tri6
     TYPE(link_type) :: link
     TYPE(user_data_type_surfmesh_tri6), POINTER :: data  
  END TYPE user_type_surfmesh_tri6
  
  TYPE user_data_type_surfmesh_tri6
     INTEGER, DIMENSION(1:7) :: elemdata
  END TYPE user_data_type_surfmesh_tri6

! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_surfmesh_tri6
     TYPE(user_type_surfmesh_tri6), POINTER :: p
  END TYPE user_ptr_type_surfmesh_tri6
  
  TYPE(user_ptr_type_surfmesh_tri6)  :: user_surfmesh_tri6
  
  TYPE(list_type) :: surfmesh_tri6_s
  TYPE(list_type) :: surfmesh_tri6_sf
  
! User-defined list element
! The Link_Type field MUST be the FIRST in the user-defined list element 
! Note pointer to data so as to easily create sublists

  TYPE user_type_surfmesh_hex8
     TYPE(link_type) :: link
     TYPE(user_data_type_surfmesh_hex8), POINTER :: data  
  END TYPE user_type_surfmesh_hex8
  
  TYPE user_data_type_surfmesh_hex8
     INTEGER, DIMENSION(1:5) :: elemdata
  END TYPE user_data_type_surfmesh_hex8

! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_surfmesh_hex8
     TYPE(user_type_surfmesh_hex8), POINTER :: p
  END TYPE user_ptr_type_surfmesh_hex8
  
  TYPE(user_ptr_type_surfmesh_hex8)  :: user_surfmesh_hex8

  TYPE(list_type) :: surfmesh_hex8_s
  TYPE(list_type) :: surfmesh_hex8_sf

! User-defined list element
! The Link_Type field MUST be the FIRST in the user-defined list element 
! Note pointer to data so as to easily create sublists

  TYPE user_type_procnodelist
     TYPE(link_type) :: link
     TYPE(user_data_type_procnodelist), POINTER :: data  
  END TYPE user_type_procnodelist

  TYPE user_data_type_procnodelist
     INTEGER :: locndnum
  END TYPE user_data_type_procnodelist

! Auxilliary data type required for the transfer function
  TYPE user_ptr_type_procnodelist
     TYPE(user_type_procnodelist), POINTER :: p
  END TYPE user_ptr_type_procnodelist
  
  TYPE(user_ptr_type_procnodelist)  :: user_procnodelist
  
  TYPE(list_type), POINTER, DIMENSION(:) :: procnodelist

! type of surface mesh
      ! = 3 for all tri
      ! = 4 for all quad
      ! = 5 for mixed

  INTEGER :: meshtype2d
  
  INTEGER, DIMENSION(:), POINTER :: elflag
  
  INTEGER, DIMENSION(:,:),POINTER :: numelhex2d
  INTEGER, DIMENSION(:,:),POINTER :: numeltet2d

  INTEGER, DIMENSION(:,:), POINTER :: elconntable
  INTEGER, DIMENSION(:), POINTER :: mattype

  
  INTEGER, DIMENSION(:), POINTER :: numelvolmat
  INTEGER, DIMENSION(:), POINTER :: numelpartbndrymat


  INTEGER, DIMENSION(:,:),POINTER :: nodeflag_str
  INTEGER, DIMENSION(:,:),POINTER :: nodeflag_mm
  INTEGER, DIMENSION(:,:),POINTER :: nodeflag_th

  INTEGER, DIMENSION(:), POINTER :: numneighprocs_list

  INTEGER, DIMENSION(:), POINTER :: nodestocommunicate,id_sendto_list
  INTEGER :: nodestocommunicate_cnt

  INTEGER, DIMENSION(:,:), ALLOCATABLE :: bc_flag
  
  INTEGER, DIMENSION(:), ALLOCATABLE :: numbc_flag
  INTEGER, POINTER, DIMENSION(:) :: pconn_comm

  INTEGER :: maxnumbc_str,maxnumbc_mm,maxnumbc_th, numserbc

  INTEGER, POINTER, DIMENSION(:,:) :: bc_values_mm, bc_values_str
  INTEGER, POINTER, DIMENSION(:) :: bc_values_th
  REAL*8, POINTER, DIMENSION(:) :: bcvalue

  LOGICAL :: interactmesh
  LOGICAL :: overlaymesh


END MODULE meshdata