src/examples/itm_grid_example8_separatrix.f90

Go to the documentation of this file.

program itm_grid_example8_separatrix


  use itm_types         ! ITM standard Fortran data types, e.g. real(R8)
  use itm_constants
  use euITM_schemas     ! ITM data structure/CPO definitions
  use euITM_routines    ! ITM UAL routines
  use itm_grid_common   ! ITM grid service library: constant definitions like COORDTYPE_R, COORDTYPE_Z
  use itm_grid_structured ! ITM grid service libary, routines for handling structured grids

  implicit none

  ! parameters
  integer, parameter :: NPOINTR = 6
  integer, parameter :: NPOINTZ = 5   

  integer, parameter :: SHOTNUM = 17151
  integer, parameter :: RUNNUM = 20407

  ! variables
  type (type_edge) :: edgecpo
  integer :: idx
  integer :: ir, iz, i, j, k
  real(R8) :: cellData(NPOINTR - 1, NPOINTZ - 1)
  real(R8) :: nodeData(NPOINTR, NPOINTZ)

  integer :: sepId,ompId,n_sep_objs,n_omp_objs,match
  type(GridObject) :: obj,face,node
  type(GridObject),dimension(:),allocatable :: nodes
  integer,dimension(:),allocatable :: sepNodes, ompNodes
  
  real(R8) :: ne, Te, Ti

  ! === 1. Set up CPO ===

  write (*,*) "Example 8: reading from shot ", SHOTNUM, ", run ", RUNNUM
  call euitm_open_env( 'euitm', SHOTNUM, RUNNUM, idx, 'coster', 'aug', '4.09a')
  call euitm_get_slice( idx, "edge", edgecpo, 2.5d0, UAL_PREVIOUS_SAMPLE ) 

  ! === 

  ! figure out subgrid id for separatrix
  sepId = gridFindSubgridByName(edgecpo%grid, "Separatrix")
  ompId = gridFindSubgridByName(edgecpo%grid, "Outer midplane")

  !print *,sepId,ompId

  n_sep_objs = gridSubGridSize(edgecpo%grid%subgrids(sepId))
  n_omp_objs = gridSubGridSize(edgecpo%grid%subgrids(ompId))

  allocate(sepNodes(2*n_sep_objs),ompNodes(n_omp_objs))
  sepNodes = -9999
  ompNodes = -9999

  k = 0
  do i = 1,n_sep_objs
     face = subGridGetObject(edgecpo%grid%subgrids(sepId),i)
     !call gridWriteGridObject(face)
     call getComposingObjects(edgecpo%grid,face,nodes)
     !call writeObjectList(nodes)
     do j = 1,size(nodes)
        !call gridWriteGridObject(nodes(j))
        !if (k .eq. 1) then
        !   sepNodes(k) = nodes(j)%ind
        !   k = k + 1
        !else
           if (all(sepNodes .ne. nodes(j)%ind(1),1)) then
              k = k + 1
              !print *,k,nodes(j)%ind(1)
              sepNodes(k) = nodes(j)%ind(1)
           end if
        !end if
     end do
     deallocate(nodes)
  end do

  !print *, sepNodes(:k)

  do i = 1,n_omp_objs
     node = subGridGetObject(edgecpo%grid%subgrids(ompId),i)
     !call gridWriteGridObject(node)
     if (any(sepNodes(:k) .eq. node%ind(1))) then
        match = i !node%ind(1)
        !print *, "Node ", match, " at position ", i, " is in outer midplane and separatrix"
    end if
  end do  

  do i = 1,size(edgecpo%fluid%ne%value)
     !print *,ompId, edgecpo%fluid%ne%value(i)%subgrid
     if (edgecpo%fluid%ne%value(i)%subgrid .eq. ompId) then
         ne = edgecpo%fluid%ne%value(i)%scalar(match)
         te = edgecpo%fluid%te%value(i)%scalar(match)
         ti = edgecpo%fluid%ti(1)%value(i)%scalar(match)
         print *, ne, te, ti
    else 
       !print *, "Cannot find the outer midplane values in the edge cpo."
    end if
  end do

end program itm_grid_example8_separatrix