test_template.f

Go to the documentation of this file.

C> @file tests/test_template.f
C> Unit tests of subroutines from formats/template.f and convert/template.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_template

      USE gridman
      USE gridman_lib
      IMPLICIT NONE

      INTEGER(GRIDMAN_SP) :: m,l
      INTEGER :: ierr,res
      INTEGER(GRIDMAN_SP),ALLOCATABLE :: n(:)
      REAL(GRIDMAN_DP),ALLOCATABLE :: x(:),y(:)  
      INTEGER(GRIDMAN_SP),PARAMETER :: m0=3
      INTEGER(GRIDMAN_SP),PARAMETER :: l0=11
      INTEGER(GRIDMAN_SP) :: n0(m0)
      REAL(GRIDMAN_DP) :: x0(l0),y0(l0)
      TYPE(gridman_grid) :: grid1,grid2

      gridman_check=.true.
      gridman_dbg=.false. 

      n0(1)=4
      x0(1)=1.
      y0(1)=1.
      x0(2)=5.
      y0(2)=1.
      x0(3)=5.
      y0(3)=5.
      x0(4)=1.
      y0(4)=1.

      n0(2)=2
      x0(5)=3.
      y0(5)=5.
      x0(6)=1.
      y0(6)=3.

      n0(3)=5
      x0(7)=0.
      y0(7)=0.
      x0(8)=0.
      y0(8)=6.
      x0(9)=6.
      y0(9)=6.
      x0(10)=6.
      y0(10)=0.              
      x0(11)=0.
      y0(11)=0

      CALL gridman_template2grid(grid1,m0,n0,l0,x0,y0,ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_grid2template(grid1,m,n,l,x,y,ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"
  
cc      IL=0
cc      WRITE(*,*)
cc      DO IP=1,M
cc       DO IV=1,N(IP)
cc        IL=IL+1
cc        WRITE(*,*) X(IL),Y(IL)
cc       END DO
cc       WRITE(*,*)
cc      END DO  

      CALL gridman_template_write('tmp.txt',m0,n0,l0,x0,y0,ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_template_read(grid2,'tmp.txt',ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_grid_compare(grid1,grid2,res,ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"
      IF(res.NE.0) THEN
       WRITE(*,*) "ERROR after GRIDMAN_GRID_COMPARE"
       WRITE(*,*) "Grids are different, RES ",res
       stop "TEST_TEMPLATE TERMINATED"
      END IF 

      CALL gridman_template_read(grid2,'./input/54001.wall',ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

cc      IL=0
cc      WRITE(*,*)
cc      DO IP=1,M
cc       DO IV=1,N(IP)
cc        IL=IL+1
cc        WRITE(*,*) X(IL),Y(IL)
cc       END DO
cc       WRITE(*,*)
cc      END DO  

      CALL gridman_template_write(grid2,'tmp2.txt',ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_grid_write(grid1,'tmp1.grd',ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_grid_write(grid2,'tmp2.grd',ierr)
      IF(ierr.NE.0) stop "TEST_TEMPLATE TERMINATED"

      CALL gridman_grid_deallocate(grid1,ierr)
      CALL gridman_grid_deallocate(grid2,ierr)

      WRITE(*,*) "TEST_TEMPLATE COMPLETED"

      END PROGRAM test_template