test_addsurf.f

Go to the documentation of this file.

C> @file tests/test_addsurf.f
C> Unit tests of subroutines from formats/addsurf.f and convert/addsurf.f
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_addsurf
      USE gridman
      USE gridman_lib
      IMPLICIT NONE

      INTEGER(GRIDMAN_SP) :: nsurf,ip1,ip2,i
      INTEGER :: res,ierr
      REAL(GRIDMAN_DP),ALLOCATABLE :: x(:,:),y(:,:),s(:),vn(:,:)
      INTEGER(GRIDMAN_SP),ALLOCATABLE :: isrf(:)   
      TYPE(gridman_grid) :: grid,grid2
      REAL(GRIDMAN_DP) :: vol(1),xcn(2,1),x1,y1,x2,y2
      LOGICAL :: ltmp(1)

      gridman_check=.true.

      CALL gridman_addsurf2d_read(grid,'./input/iterm.carre.105',ierr) 
      IF(ierr.NE.100) THEN
        WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ_ARRAY"
        WRITE(*,*) "Expected value IERR=100, IERR ",ierr
        stop "TEST_ADDSURF TERMINATED"
      END IF

      CALL gridman_addsurf2d_read(grid,'./input/input.eir.corrupt',ierr) 
      IF(ierr.NE.300) THEN
        WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ_ARRAY"
        WRITE(*,*) "Expected value IERR=300, IERR ",ierr
        stop "TEST_ADDSURF TERMINATED"
      END IF

      CALL gridman_addsurf2d_read('./input/input.eir',
     c                          nsurf,x,y,isrf,ierr) 
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"    
      CALL gridman_addsurf2grid_2d(grid,nsurf,x,y,isrf,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_addsurf2d_read(grid,'./input/input.eir',ierr) 
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"    

      IF(grid%NEDGES.NE.180_gridman_sp) THEN
       WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ"
       WRITE(*,*) "Expected value NEDGES=180, NEDGES ",grid%NEDGES
       stop "TEST_ADDSURF TERMINATED"      
      END IF

      ip1=grid%POINTS(1,180)
      ip2=grid%POINTS(2,180)
      x1=grid%X(1,ip1)
      y1=grid%X(2,ip1)
      x2=grid%X(1,ip2)
      y2=grid%X(2,ip2)
      IF(abs(x1-491.389000000000).GT.gridman_tol*abs(x1).OR.
     f   abs(y1+394.193000000000).GT.gridman_tol*abs(y1).OR.
     f   abs(x2-481.061000000000).GT.gridman_tol*abs(x2).OR.
     f   abs(y2+390.124000000000).GT.gridman_tol*abs(y2)) THEN
        WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ: ",
     w             "wrong coordinates"
        WRITE(*,*) " X1, Y1, X2, Y2 ",x1,y1,x2,y2
        WRITE(*,*) "Expected ", 491.389,-394.193,481.061,-390.124
        stop "TEST_ADDSURF TERMINATED"      
      END IF

      IF(grid%EDGEINDEX(1)%NELEMENTS.NE.grid%NEDGES) THEN
        WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ: ",
     w              "wrong index"
        WRITE(*,*) " NELEMENETS, NEDGES ",
     w               grid%EDGEINDEX(1)%NELEMENTS,grid%NEDGES
        stop "TEST_ADDSURF TERMINATED"      
      END IF
      DO i=1,grid%EDGEINDEX(1)%NELEMENTS
       ip1=grid%EDGEINDEX(1)%INDEXES(0,i)
       ip2=grid%EDGEINDEX(1)%INDEXES(1,i)
       IF(ip1.NE.i.OR.ip2.NE.i) THEN
        WRITE(*,*) "ERROR after GRIDMAN_ADDSURF2D_READ: ",
     w              "wrong index"
        WRITE(*,*) " IE, IEDGE, IADD ",i,ip1,ip2
        stop "TEST_ADDSURF TERMINATED"  
       END IF
      END DO

cc      DO I=1,GRID%NEDGES
cc       WRITE(*,*) " ISRF ",I
cc       IP1=GRID%POINTS(1,1,I)
cc       IP2=GRID%POINTS(1,2,I)
cc       WRITE(*,*) " X1, Y1 ",GRID%X1(IP1),GRID%X2(IP1)
cc       WRITE(*,*) " X2, Y2 ",GRID%X1(IP2),GRID%X2(IP2)
cc     END DO

      ALLOCATE(s(grid%NEDGES)) 
      CALL gridman_grid2d_cylareas(grid,s,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid2d_cylvolumes(grid,vol,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid2d_center(grid,xcn,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      ALLOCATE(vn(2,grid%NEDGES)) 
      CALL gridman_grid2d_norm(grid,vn,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid2d_isconvex(grid,ltmp,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid2d_crossect(grid,vol,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid_write(grid,'addsrf.grd',ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid_copy(grid2,grid,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid_write(grid,'addsrf.grd',ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid_read(grid,'addsrf.grd',ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      CALL gridman_grid_compare(grid,grid2,res,ierr)
      IF(ierr.NE.0) stop "TEST_ADDSURF TERMINATED"

      IF(res.NE.0) THEN
       WRITE(*,*) "ERROR after GRIDMAN_GRID_READ"
       WRITE(*,*) "Expected value RES=0, RES ",res
       stop "TEST_ADDSURF TERMINATED"
      END IF

      WRITE(*,*) "TEST_ADDSURF COMPLETED"

      END  PROGRAM test_addsurf