test_grid3d.f

Go to the documentation of this file.

C> @file tests/test_grid3d.f
C> Unit tests of subroutines reletaed to the 3D grid type
C GRIDMAN, grid managment library. Author: Vladislav Kotov, v.kotov@fz-juelich.de

!    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
!    along with GRIDMAN.  If not, see <http://www.gnu.org/licenses/>.

      PROGRAM test_grid3d
      USE gridman
      USE gridman_lib
      IMPLICIT NONE

      TYPE(gridman_grid) :: grid
      INTEGER :: ierr,st
      REAL(GRIDMAN_DP),ALLOCATABLE :: data(:,:)
      INTEGER(GRIDMAN_SP) :: ie,iedge,isrf

      gridman_check=.true.
 
      CALL gridman_addsurf3d_read(grid,'input/input.eir.3D',ierr) 
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'

      CALL gridman_grid_write(grid,'grid3d.grd',ierr)
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'
     
      WRITE(*,*) 'eirene3D.vtk'
      CALL gridman_vtk_grid3d_write('eirene3D.vtk','VTK 3D test',
     s                              grid,grid%NEDGES,0,0,0,0,ierr)
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'

C     cell scalar added 
      WRITE(*,*) 'eirene3D-2.vtk'
      ALLOCATE(DATA(grid%NEDGES,1),stat=st)
      IF(st.NE.0) THEN
       WRITE(*,*) "Cannot allocate memory"
       stop 'TEST_GRID3D TERMINATED'
      END IF       
      DATA(:,1)=0
      DO ie=1,grid%EDGEINDEX(1)%NELEMENTS
       iedge=grid%EDGEINDEX(1)%INDEXES(0,ie)
       isrf=grid%EDGEINDEX(1)%INDEXES(1,ie)
       DATA(iedge,1)=isrf
      END DO

      CALL gridman_vtk_grid3d_write('eirene3D_2.vtk','VTK 3D test',
     s                              grid,grid%NEDGES,
     s                              1,0,0,0,ierr, 
     s                              cell_scalar=DATA,csname='ADDSURF')
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'

C     index removed
      WRITE(*,*) 'eirene3D-3.vtk'
      CALL gridman_index_deallocate(grid%EDGEINDEX(1),ierr)
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'
      DEALLOCATE(grid%EDGEINDEX)
      grid%NEDGEINDEX=0

      CALL gridman_vtk_grid3d_write('eirene3D_3.vtk','VTK 3D test',
     s                              grid,grid%NEDGES,
     s                              1,0,0,0,ierr, 
     s                              cell_scalar=DATA,csname='ADDSURF')
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'

C     w/o cell scalar
      WRITE(*,*) 'eirene3D-4.vtk'
      CALL gridman_vtk_grid3d_write('eirene3D_4.vtk','VTK 3D test',
     s                              grid,grid%NEDGES,
     s                              0,0,0,0,ierr)
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'


      DEALLOCATE(DATA,stat=st)
      IF(st.NE.0) THEN
       WRITE(*,*) "Cannot deallocate memory"
       stop 'TEST_GRID3D ERMINATED'
      END IF

      CALL gridman_grid_deallocate(grid,ierr)
      IF(ierr.NE.0) stop 'TEST_GRID3D TERMINATED'

      WRITE(*,*) "TEST_GRID3D COMPLETED"

      END PROGRAM test_grid3d