src/service/itm_grid_examples.f90

Go to the documentation of this file.

module itm_grid_examples

  use itm_types
  use euITM_schemas ! IGNORE
  use euITM_routines ! IGNORE
  use itm_grid
  use itm_grid_structured
  use itm_grid_helper
  
  implicit none

  integer, parameter :: TEST1_SHOT  = 9001 
  integer, parameter :: TEST1_RUN   = 1
  integer, parameter :: TEST1_TIME  = 0.0

contains

  ! Example 1: 2d structured grid

  !> Example routine that assembles a 2d grid consisting of two 1d spaces
  function example_2d() result( grid )
    type(type_complexgrid) :: grid
    
    ! internal
    integer, parameter :: NPOINTR = 6
    integer, parameter :: NPOINTZ = 5    

    integer :: i

    call gridSetupStructuredSep( grid, &
         & ndim = 2, &
         & c1 = COORDTYPE_R, &
         & x1 = (/ ( 1.0_R8 * i, i=0,NPOINTR-1) /), &
         & c2 = COORDTYPE_Z, &
         & x2 = (/ ( 0.5_R8 * i, i=0,NPOINTZ-1) /),&
         & id = '2d_structured' )

    call gridSanityCheck( grid, doStop = .false. ) 
    
  end function example_2d


!!$  subroutine itm_grid_write_example_2d
!!$    
!!$    ! variables
!!$    type (type_edge),pointer :: edgecpo(:) => null()
!!$    integer :: idx
!!$
!!$    ! internal
!!$    integer, parameter :: NSPC = 2
!!$    integer, parameter :: NDIM = 1
!!$    integer, parameter :: NPOINTR = 6
!!$    integer, parameter :: NPOINTZ = 5   
!!$   
!!$    integer :: ir, iz
!!$    real(R8) :: nodeData(NPOINTR, NPOINTZ), cellData(NPOINTR - 1, NPOINTZ - 1)
!!$
!!$    ! set up cpo
!!$    allocate( edgecpo(1) )    
!!$    call setupEdgeCPO( edgecpo(1) )
!!$
!!$    ! set up grid
!!$    edgecpo(1)%grid = example_2d()
!!$
!!$    ! Write some fake data to CPO.
!!$    allocate(edgecpo(1)%fluid%ne%value(4))
!!$
!!$    do ir = 1, NPOINTR
!!$        nodeData(ir,:) = ir
!!$    end do
!!$    call gridStructWriteData2d( edgecpo(1)%grid, edgecpo(1)%fluid%ne%value(1), GRID_STRUCT_NODES, nodeData )
!!$
!!$    do ir = 1, NPOINTR - 1 
!!$        cellData(ir, :) = ir
!!$    enddo
!!$    call gridStructWriteData2d( edgecpo(1)%grid, edgecpo(1)%fluid%ne%value(2), GRID_STRUCT_FACES, cellData )
!!$
!!$    do iz = 1, NPOINTZ - 1 
!!$        cellData(:, iz) = iz
!!$    enddo
!!$    call gridStructWriteData2d( edgecpo(1)%grid, edgecpo(1)%fluid%ne%value(3), GRID_STRUCT_FACES, cellData )
!!$
!!$    do ir = 1, NPOINTR - 1 
!!$       do iz = 1, NPOINTZ - 1
!!$           cellData(ir, iz) = ir * iz          
!!$       enddo
!!$    enddo
!!$    call gridStructWriteData2d( edgecpo(1)%grid, edgecpo(1)%fluid%ne%value(4), GRID_STRUCT_FACES, cellData )
!!$
!!$
!!$    ! Write CPO to ual
!!$    write (*,*) "Example 1: writing to shot ", TEST1_SHOT, ", run ", TEST1_RUN
!!$    call euitm_create( 'euitm', TEST1_SHOT, TEST1_RUN, 0, 0, idx)
!!$    call euitm_put(idx,"edge", edgecpo)
!!$    call euitm_close(idx)
!!$    call euitm_deallocate(edgecpo)   
!!$
!!$  end subroutine itm_grid_write_example_2d


  subroutine itm_grid_read_example_2d( edgecpo )
    type(type_edge), intent(out) :: edgecpo
      
    ! variables
    !type(type_complexgrid) :: grid
    type(type_edge),pointer :: ledgecpo(:) => null()
    integer :: idx
    
    ! read cpo
    ! treename, shot, run
    call euitm_open( 'euitm', TEST1_SHOT, TEST1_RUN, idx )
    call euitm_get( idx,"edge", ledgecpo ) 

    edgecpo = ledgecpo(1)

  end subroutine itm_grid_read_example_2d


end module itm_grid_examples