interp_FOcollision_p Subroutine

  1. korc_interp.f90
  2. korc_interp
  3. interp_FOcollision_p

public subroutine interp_FOcollision_p(Y_R, Y_PHI, Y_Z, ne, Te, Zeff, flag_cache)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in), DIMENSION(8):: Y_R
real(kind=rp), intent(in), DIMENSION(8):: Y_PHI
real(kind=rp), intent(in), DIMENSION(8):: Y_Z
real(kind=rp), intent(out), DIMENSION(8):: ne
real(kind=rp), intent(out), DIMENSION(8):: Te
real(kind=rp), intent(out), DIMENSION(8):: Zeff
integer(kind=is), intent(inout), DIMENSION(8):: flag_cache

Contents

Source Code

interp_FOcollision_p

Source Code

subroutine interp_FOcollision_p(Y_R,Y_PHI,Y_Z,ne,Te,Zeff,flag_cache)

  REAL(rp),DIMENSION(8),INTENT(IN)   :: Y_R,Y_PHI,Y_Z
  REAL(rp),DIMENSION(8),INTENT(OUT)   :: ne,Te,Zeff
  INTEGER(is),DIMENSION(8),INTENT(INOUT)   :: flag_cache

  call check_if_in_profiles_domain_p(Y_R,Y_PHI,Y_Z,flag_cache)
!  write(6,'("YR: ",E17.10)') Y_R(1)
!  write(6,'("YPHI: ",E17.10)') Y_PHI(1)
!  write(6,'("YZ: ",E17.10)') Y_Z(1)
  
!  write(6,'("Te_interp_R",E17.10)') profiles_2d%Te%x1
!  write(6,'("Te_interp_Z",E17.10)') profiles_2d%Te%x2
  
  call EZspline_interp(profiles_2d%ne,profiles_2d%Te, &
       profiles_2d%Zeff,8,Y_R,Y_Z,ne,Te,Zeff,ezerr)
  ! this will call PSPLINE routine EZspline_interp2_bmag_cloud_r8 as there
  ! is the same number of entries
  call EZspline_error(ezerr)
 
  
end subroutine interp_FOcollision_p