test_carre_DN.f

Go to the documentation of this file.

C> @file tests/test_carre_DN.f
C> Unit tests of subroutines from formats/carre.f and convert/carre.f 
C> for Double Null topolgy
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_carre_dn
      USE gridman
      USE gridman_lib
      IMPLICIT NONE

      TYPE(gridman_grid) :: sonnet_grid,fort30_grid
      REAL(GRIDMAN_DP),ALLOCATABLE :: vol1(:),vol2(:)
      REAL(GRIDMAN_DP) :: v1,v2
      INTEGER :: ierr  

      gridman_check=.true.

      WRITE(gridman_unit,*) "READ30"
      CALL gridman_carre_read30_grid(fort30_grid,
     c       './input/fort.30.kstar',ierr,.true.)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"

      WRITE(gridman_unit,*) "READSONNET"
      CALL gridman_carre_readsonnet_grid(sonnet_grid,
     c       './input/kstar.carre.082_2',ierr,.true.)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"

      WRITE(gridman_unit,*) "CYLVOLUMES"
      ALLOCATE(vol1(sonnet_grid%NCELLS),vol2(fort30_grid%NCELLS))
      CALL gridman_grid2d_cylvolumes(sonnet_grid,vol1,ierr)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"
      CALL gridman_grid2d_cylvolumes(fort30_grid,vol2,ierr)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"

      v1=sum(vol1)*sonnet_grid%UNIT2SI**3
      v2=sum(vol2)*fort30_grid%UNIT2SI**3
      IF(abs(v1-v2).GT.gridman_tol*(abs(v1)+abs(v2))) THEN
       WRITE(gridman_unit,*) "ERROR detected in TEST_CARRE_DDN: ",
     w       "mismatch of grid volumes between SONNET and FORT.30"              
       WRITE(gridman_unit,*) "V_SONNET, V_FORT30 ",v1,v2
       stop "TEST_CARRE_DDN TERMINATED"
      END IF
    
      WRITE(gridman_unit,*) "WRITE"
      CALL gridman_grid_write(sonnet_grid,'sonnetDN.grd',ierr)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"
      CALL gridman_grid_write(fort30_grid,'fort30DN.grd',ierr)
      IF(ierr.NE.0) stop "TEST_CARRE_DN TERMINATED"
     
      WRITE(*,*) "TEST_CARRE_DN COMPLETED"
     
      END PROGRAM test_carre_dn