indlist.f

Go to the documentation of this file.

C> @file indlist.f
C> Methods of data-type GRIDMAN_INDLIS, see gridman.f
C>
C> List of elements with variable number of indices for each element

!    Copyright (c) 2017  Forschungszentrum Juelich GmbH
!                        Vladislav Kotov 
!
!    This file is part of GRIDMAN.
!
!    GRIDMAN is free software: you can redistribute it and/or modify
!    it under the terms of the GNU General Public License as published by
!    the Free Software Foundation, either version 3 of the License, or
!    (at your option) any later version.
!
!    GRIDMAN is distributed in the hope that it will be useful,
!    but WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!    GNU General Public License for more details.
!
!    You should have received a copy of the GNU General Public License

C
C GRIDMAN, grid managment library. Author: Vladislav Kotov, v.kotov@fz-juelich.de
C
C *****************************************************************************************
C> Allocate list of elements
C>
C> WARNING: if INDLIST elready exists it will be re-allocated
C *****************************************************************************************
       SUBROUTINE gridman_indlist_allocate(INDLIST,N,L,IERR)
       USE gridman,ONLY:gridman_indlist,gridman_sp,
     u                  gridman_dbg,gridman_unit
       IMPLICIT NONE
C> Resulting list of elements object 
       TYPE(gridman_indlist) :: INDLIST
C> Number of elements
       INTEGER(GRIDMAN_SP),INTENT(IN) :: N
C> Total number of indices
       INTEGER(GRIDMAN_SP),INTENT(IN) :: L
C> Error code
       INTEGER,INTENT(OUT) :: IERR
             
       INTEGER :: ST

       IF(gridman_dbg) 
     w   WRITE(gridman_unit,*) "Starting GRIDMAN_INDLIST_ALLOCATE"

       ierr=0

       IF(n.LT.1.OR.l.LT.1.OR.n.GT.l) THEN
        ierr=100
        WRITE(gridman_unit,*) "ERROR in GRIDMAN_INDLIST_ALLOCATE: ",
     w                        "incorrect array size" 
        WRITE(gridman_unit,*) " N, L ",n,l
        RETURN
       END IF 
       
       CALL gridman_indlist_deallocate(indlist,ierr)
       IF(ierr.NE.0) THEN
        WRITE(gridman_unit,*) "GRIDMAN_INDLIST_ALLOCATE terminated"
        RETURN
       END IF

       indlist%N=n
       indlist%L=l

       ALLOCATE(indlist%IFIRST(n),stat=st)
       IF(st.NE.0) THEN
        ierr=200
        WRITE(gridman_unit,*) "ERROR in GRIDMAN_INDLIST_ALLOCATE: ",
     w                        "cannot allocate array IFIRST"
        WRITE(gridman_unit,*) " N ",n
        RETURN
       END IF

       ALLOCATE(indlist%ILAST(n),stat=st)
       IF(st.NE.0) THEN
        ierr=200
        WRITE(gridman_unit,*) "ERROR in GRIDMAN_INDLIST_ALLOCATE: ",
     w                        "cannot allocate array ILAST"
        WRITE(gridman_unit,*) " N ",n
        RETURN
       END IF

       ALLOCATE(indlist%IND(l),stat=st)
       IF(st.NE.0) THEN
        ierr=200
        WRITE(gridman_unit,*) "ERROR in GRIDMAN_INDLIST_ALLOCATE: ",
     w                        "cannot allocate array IND"
        WRITE(gridman_unit,*) " L ",l
        RETURN
       END IF


       IF(gridman_dbg) 
     w   WRITE(gridman_unit,*) "GRIDMAN_INDLIST_ALLOCATE finished"


       END SUBROUTINE gridman_indlist_allocate

C *****************************************************************************************
C> Deallocate list of indices
C *****************************************************************************************
      SUBROUTINE gridman_indlist_deallocate(INDLIST,IERR)
      USE gridman,ONLY:gridman_indlist,gridman_dbg,gridman_unit
      IMPLICIT NONE
C> List of indices to be deallocated
      TYPE(gridman_indlist) :: INDLIST
C> Error code
      INTEGER,INTENT(OUT) :: IERR
            
      INTEGER :: ST

      IF(gridman_dbg) 
     w   WRITE(gridman_unit,*) "Starting GRIDMAN_INDLIST_DEALLOCATE"

      ierr=0 

      IF(ALLOCATED(indlist%IND)) THEN
       DEALLOCATE(indlist%IND,stat=st)
       IF(st.NE.0) THEN
        ierr=-200
         WRITE(gridman_unit,*) 
     w            "WARNING from GRIDMAN_INDLIST_DEALLOCATE: ",
     w            "can not deallocate IND"       
       END IF
      END IF

      IF(ALLOCATED(indlist%IFIRST)) THEN
       DEALLOCATE(indlist%IFIRST,stat=st)
       IF(st.NE.0) THEN
        ierr=-200
         WRITE(gridman_unit,*) 
     w            "WARNING from GRIDMAN_INDLIST_DEALLOCATE: ",
     w            "can not deallocate IEDGESS"       
       END IF
      END IF

      IF(ALLOCATED(indlist%ILAST)) THEN
       DEALLOCATE(indlist%ILAST,stat=st)
       IF(st.NE.0) THEN
        ierr=-200
         WRITE(gridman_unit,*) 
     w            "WARNING from GRIDMAN_INDLIST_DEALLOCATE: ",
     w            "can not deallocate ILAST"       
       END IF
      END IF

      indlist%N=0
      indlist%L=0

      IF(gridman_dbg) 
     w   WRITE(gridman_unit,*) "GRIDMAN_CELLS_DEALLOCATE finished"

      END SUBROUTINE gridman_indlist_deallocate