add_interp_SCE_p_FS Subroutine

public subroutine add_interp_SCE_p_FS(params, F, PSIp, E_PHI)

Arguments

Type IntentOptional AttributesName
type(KORC_PARAMS), intent(in) :: params
type(FIELDS), intent(in) :: F
real(kind=rp), intent(in), DIMENSION(params%pchunk):: PSIp
real(kind=rp), intent(inout), DIMENSION(params%pchunk):: E_PHI

Calls

proc~~add_interp_sce_p_fs~~CallsGraph proc~add_interp_sce_p_fs add_interp_SCE_p_FS ezspline_interp ezspline_interp proc~add_interp_sce_p_fs->ezspline_interp ezspline_error ezspline_error proc~add_interp_sce_p_fs->ezspline_error

Called by

proc~~add_interp_sce_p_fs~~CalledByGraph proc~add_interp_sce_p_fs add_interp_SCE_p_FS proc~advance_gcinterp_psi_vars_fs advance_GCinterp_psi_vars_FS proc~advance_gcinterp_psi_vars_fs->proc~add_interp_sce_p_fs

Contents

Source Code


Source Code

subroutine add_interp_SCE_p_FS(params,F,PSIp,E_PHI)
  TYPE(KORC_PARAMS), INTENT(IN)                              :: params
  TYPE(FIELDS), INTENT(IN)                               :: F
  REAL(rp), DIMENSION(params%pchunk), INTENT(IN)      :: PSIp
  REAL(rp), DIMENSION(params%pchunk), INTENT(INOUT)      :: E_PHI

  REAL(rp),DIMENSION(params%pchunk) :: E_SC_PHI
  INTEGER :: cc,pchunk

  pchunk=params%pchunk

  call EZspline_interp(efield_SC1d%PHI,pchunk, PSIp, E_SC_PHI, ezerr)
  call EZspline_error(ezerr)

  !$OMP SIMD
  do cc=1_idef,pchunk
     E_PHI(cc)=E_PHI(cc)+E_SC_PHI(cc)
  end do
  !$OMP END SIMD  
  
end subroutine add_interp_SCE_p_FS