check_if_in_LCFS Subroutine

  1. korc_interp.f90
  2. korc_interp
  3. check_if_in_LCFS

private subroutine check_if_in_LCFS(F, Y, inLCFS)

External fields and interpolants can have different spatial domains where they are defined. Therefore, it is necessary to check if a given particle has left these spatial domains to stop following it, otherwise this will cause an error in the simulation.

Arguments

Type IntentOptional AttributesName
type(FIELDS), intent(in) :: F
real(kind=rp), intent(in), DIMENSION(:,:), ALLOCATABLE:: Y

Particles' position in cylindrical coordinates, Y(1,:) = , Y(2,:) = , and Y(3,:) = .
integer(kind=is), intent(inout), DIMENSION(:), ALLOCATABLE:: inLCFS

Flag that determines whether particles are followed in the simulation (flag=1), or not (flag=0).

Called by

proc~~check_if_in_lcfs~~CalledByGraph proc~check_if_in_lcfs check_if_in_LCFS proc~interp_fields interp_fields proc~interp_fields->proc~check_if_in_lcfs proc~get_fields get_fields proc~get_fields->proc~interp_fields proc~mh_psi MH_psi proc~mh_psi->proc~get_fields proc~sample_hollmann_distribution_3d_psi sample_Hollmann_distribution_3D_psi proc~sample_hollmann_distribution_3d_psi->proc~get_fields proc~unitvectors unitVectors proc~unitvectors->proc~get_fields proc~get_hollmann_distribution_3d_psi get_Hollmann_distribution_3D_psi proc~get_hollmann_distribution_3d_psi->proc~sample_hollmann_distribution_3d_psi proc~gyro_distribution gyro_distribution proc~gyro_distribution->proc~unitvectors

Contents

Source Code

check_if_in_LCFS

Source Code

  subroutine check_if_in_LCFS(F,Y,inLCFS)
    !! @note Subrotuine that checks if particles in the simulation are within
    !! the spatial domain where interpolants and fields are known. @endnote
    !! External fields and interpolants can have different spatial domains where
    !! they are defined. Therefore, it is necessary to
    !! check if a given particle has left these spatial domains to stop
    !! following it, otherwise this will cause an error in the simulation.
    TYPE(FIELDS), INTENT(IN)                                   :: F
    REAL(rp), DIMENSION(:,:), ALLOCATABLE, INTENT(IN)      :: Y
    !! Particles' position in cylindrical coordinates,
    !! Y(1,:) = \(R\), Y(2,:) = \(\phi\), and Y(3,:) = \(Z\).
    INTEGER(is), DIMENSION(:), ALLOCATABLE, INTENT(INOUT)  :: inLCFS
    !! Flag that determines whether particles are followed in the
    !! simulation (flag=1), or not (flag=0).
    INTEGER                                                :: IR
    !! Variable used to localize the grid cell in the \((R,\phi,Z)\)
    !! or \((R,Z)\) grid containing the fields data that corresponds
    !! to the radial position of the particles.
    INTEGER                                                :: IPHI
    !! Variable used to localize the grid cell in the \((R,\phi,Z)\)
    !! or \((R,Z)\) grid containing the fields data that corresponds
    !! to the azimuthal position of the particles.
    INTEGER                                                :: IZ
    !! Variable used to localize the grid cell in the \((R,\phi,Z)\)
    !! or \((R,Z)\) grid containing the fields data that corresponds
    !! to the vertical position of the particles.
    INTEGER(ip)                                            :: pp
    !! Particle iterator.
    INTEGER(ip)                                            :: ss
    !! Species iterator.
    REAL(rp) :: Rwall

    if (size(Y,1).eq.1) then
       ss = size(Y,1)
    else
       if (Y(2,1).eq.0) then
          ss=1_idef
       else
          ss = size(Y,1)
       end if
    endif

!    write(output_unit_write,'("R: ",E15.10)') Y(1,1)
!    write(output_unit_write,'("PHI: ",E15.10)') Y(2,1)
!    write(output_unit_write,'("Z: ",E15.10)') Y(1,3)

    !    write(output_unit_write,*) 'Flag',flag(1)

    !write(6,*) 'ss',ss
    !write(6,*) 'inLCFS',inLCFS
    

       
    !$OMP PARALLEL DO FIRSTPRIVATE(ss) PRIVATE(pp,IR,IZ) &
    !$OMP& SHARED(Y,inLCFS,fields_domain,bfield_2d)
    do pp=1_idef,ss

       IR = INT(FLOOR((Y(pp,1)  - fields_domain%Ro + 0.5_rp* &
            fields_domain%DR)/fields_domain%DR) + 1.0_rp,idef)
       IZ = INT(FLOOR((Y(pp,3)  + ABS(fields_domain%Zo) + 0.5_rp* &
            fields_domain%DZ)/fields_domain%DZ) + 1.0_rp,idef)

#if DBG_CHECK    
       if ((IR.lt.0).or.(IZ.lt.0).or.(IR.GT. &
            bfield_2d%NR).OR.(IZ.GT.bfield_2d%NZ)) then
          write(6,'("YR:",E17.10)') Y(1,1)
          write(6,'("YZ:",E17.10)') Y(1,3)
          write(6,'("IR: ",I16)') IR
          write(6,'("IZ: ",I16)') IZ
          !call KORC_ABORT(23)
       end if
#endif

       !write(output_unit_write,'("IR: ",I16)') IR
       !write(output_unit_write,'("IZ: ",I16)') IZ

       if (((IR.lt.0).or.(IZ.lt.0).or. &
            (IR.GT.bfield_2d%NR).OR.(IZ.GT.bfield_2d%NZ)).or. &
            (fields_domain%LCFS2D(IR,IZ).NE.1_is)) then
          !write(output_unit_write,*) 'here'

          inLCFS(pp) = 0_is

       end if
    end do
    !$OMP END PARALLEL DO

  
  end subroutine check_if_in_LCFS