interpol1D_coeff Subroutine

  1. fidasim.f90
  2. libfida
  3. interpol1D_coeff

public subroutine interpol1D_coeff(xmin, dx, nx, xout, c, err)

Linear interpolation coefficients and index for a 1D grid y(x)

Arguments

Type IntentOptional AttributesName
real(kind=Float64), intent(in) :: xmin

Minimum abscissa value
real(kind=Float64), intent(in) :: dx

Absissa spacing
integer, intent(in) :: nx

Number of abscissa
real(kind=Float64), intent(in) :: xout

Abscissa value to interpolate
type(InterpolCoeffs1D), intent(out) :: c

Interpolation Coefficients
integer, intent(out), optional :: err

Error code

Called by

proc~~interpol1d_coeff~~CalledByGraph proc~interpol1d_coeff interpol1D_coeff interface~interpol_coeff interpol_coeff interface~interpol_coeff->proc~interpol1d_coeff proc~interpol1d_coeff_arr interpol1D_coeff_arr interface~interpol_coeff->proc~interpol1d_coeff_arr proc~interpol1d_coeff_arr->proc~interpol1d_coeff proc~in_plasma in_plasma proc~in_plasma->interface~interpol_coeff proc~get_rate_matrix get_rate_matrix proc~get_rate_matrix->interface~interpol_coeff proc~interpol2d_arr interpol2D_arr proc~interpol2d_arr->interface~interpol_coeff proc~interpol3d_arr interpol3D_arr proc~interpol3d_arr->interface~interpol_coeff proc~interpol2d_2d_arr interpol2D_2D_arr proc~interpol2d_2d_arr->interface~interpol_coeff proc~get_neutron_rate get_neutron_rate proc~get_neutron_rate->interface~interpol_coeff proc~bt_cx_rates bt_cx_rates proc~bt_cx_rates->interface~interpol_coeff proc~bb_cx_rates bb_cx_rates proc~bb_cx_rates->interface~interpol_coeff proc~interpol1d_arr interpol1D_arr proc~interpol1d_arr->interface~interpol_coeff proc~interpol3d_2d_arr interpol3D_2D_arr proc~interpol3d_2d_arr->interface~interpol_coeff proc~get_fields get_fields proc~get_fields->proc~in_plasma proc~make_beam_grid make_beam_grid proc~make_beam_grid->proc~in_plasma proc~pnpa_mc pnpa_mc proc~pnpa_mc->proc~bt_cx_rates proc~pnpa_mc->proc~get_fields proc~get_plasma get_plasma proc~pnpa_mc->proc~get_plasma proc~store_npa store_npa proc~pnpa_mc->proc~store_npa proc~attenuate attenuate proc~pnpa_mc->proc~attenuate proc~track_cylindrical track_cylindrical proc~track_cylindrical->proc~in_plasma proc~track_cylindrical->proc~get_fields proc~mc_nbi mc_nbi proc~mc_nbi->proc~in_plasma proc~fida_weights_los fida_weights_los proc~fida_weights_los->proc~in_plasma proc~fida_weights_los->proc~bt_cx_rates proc~fida_weights_los->proc~bb_cx_rates proc~fida_weights_los->proc~get_fields proc~fida_weights_los->proc~get_plasma proc~track track proc~fida_weights_los->proc~track proc~colrad colrad proc~fida_weights_los->proc~colrad proc~read_mc read_mc proc~read_mc->proc~in_plasma proc~pfida_mc pfida_mc proc~pfida_mc->proc~bt_cx_rates proc~pfida_mc->proc~get_fields proc~pfida_mc->proc~track_cylindrical proc~pfida_mc->proc~get_plasma proc~get_plasma->proc~in_plasma proc~dcx_spec dcx_spec proc~dcx_spec->proc~in_plasma proc~dcx_spec->proc~get_plasma proc~store_photons store_photons proc~dcx_spec->proc~store_photons proc~neutron_mc neutron_mc proc~neutron_mc->proc~get_neutron_rate proc~neutron_mc->proc~get_fields proc~neutron_mc->proc~get_plasma proc~track->proc~in_plasma proc~track->proc~get_fields proc~halo_spec halo_spec proc~halo_spec->proc~in_plasma proc~halo_spec->proc~get_plasma proc~halo_spec->proc~store_photons proc~nbi_spec nbi_spec proc~nbi_spec->proc~in_plasma proc~nbi_spec->proc~store_photons proc~get_plasma_extrema get_plasma_extrema proc~get_plasma_extrema->proc~in_plasma interface~interpol interpol interface~interpol->proc~interpol2d_arr interface~interpol->proc~interpol3d_arr interface~interpol->proc~interpol2d_2d_arr interface~interpol->proc~interpol1d_arr interface~interpol->proc~interpol3d_2d_arr proc~neutron_f neutron_f proc~neutron_f->proc~get_neutron_rate proc~neutron_f->proc~get_fields proc~neutron_f->proc~get_plasma proc~pfida_f pfida_f proc~pfida_f->proc~bt_cx_rates proc~pfida_f->proc~track_cylindrical proc~pfida_f->proc~get_plasma proc~mc_fastion_pass_grid mc_fastion_pass_grid proc~pfida_f->proc~mc_fastion_pass_grid proc~get_beam_cx_rate get_beam_cx_rate proc~get_beam_cx_rate->proc~bt_cx_rates proc~get_beam_cx_rate->proc~bb_cx_rates proc~get_beam_cx_rate->proc~get_plasma proc~colrad->proc~get_rate_matrix proc~pnpa_f pnpa_f proc~pnpa_f->proc~bt_cx_rates proc~pnpa_f->proc~get_plasma proc~pnpa_f->proc~mc_fastion_pass_grid proc~pnpa_f->proc~store_npa proc~pnpa_f->proc~attenuate proc~bremsstrahlung bremsstrahlung proc~bremsstrahlung->proc~get_plasma proc~read_distribution read_distribution proc~read_distribution->proc~read_mc proc~mc_fastion_pass_grid->proc~get_fields proc~get_distribution get_distribution proc~mc_fastion_pass_grid->proc~get_distribution proc~get_distribution->interface~interpol proc~npa_mc npa_mc proc~npa_mc->proc~get_fields proc~npa_mc->proc~get_beam_cx_rate proc~npa_mc->proc~store_npa proc~npa_mc->proc~attenuate proc~ndmc ndmc proc~ndmc->proc~get_fields proc~ndmc->proc~mc_nbi proc~ndmc->proc~track program~fidasim fidasim program~fidasim->proc~make_beam_grid program~fidasim->proc~pnpa_mc program~fidasim->proc~fida_weights_los program~fidasim->proc~pfida_mc program~fidasim->proc~dcx_spec program~fidasim->proc~neutron_mc program~fidasim->proc~halo_spec program~fidasim->proc~nbi_spec program~fidasim->proc~neutron_f program~fidasim->proc~pfida_f program~fidasim->proc~pnpa_f program~fidasim->proc~bremsstrahlung program~fidasim->proc~read_distribution program~fidasim->proc~npa_mc program~fidasim->proc~ndmc proc~read_equilibrium read_equilibrium program~fidasim->proc~read_equilibrium proc~npa_weights npa_weights program~fidasim->proc~npa_weights proc~fida_mc fida_mc program~fidasim->proc~fida_mc proc~halo halo program~fidasim->proc~halo proc~dcx dcx program~fidasim->proc~dcx proc~fida_f fida_f program~fidasim->proc~fida_f proc~npa_f npa_f program~fidasim->proc~npa_f proc~fida_weights_mc fida_weights_mc program~fidasim->proc~fida_weights_mc proc~cold_spec cold_spec program~fidasim->proc~cold_spec proc~make_diagnostic_grids make_diagnostic_grids program~fidasim->proc~make_diagnostic_grids proc~read_equilibrium->proc~colrad proc~npa_weights->proc~get_fields proc~npa_weights->proc~get_beam_cx_rate proc~npa_weights->proc~attenuate proc~fida_mc->proc~get_fields proc~fida_mc->proc~track proc~mc_halo mc_halo proc~mc_halo->proc~get_plasma proc~halo->proc~get_plasma proc~halo->proc~track proc~halo->proc~mc_halo proc~dcx->proc~get_plasma proc~dcx->proc~track proc~dcx->proc~mc_halo proc~store_npa->proc~get_fields proc~fida_f->proc~get_plasma proc~fida_f->proc~track proc~mc_fastion mc_fastion proc~fida_f->proc~mc_fastion proc~store_photons->proc~get_fields proc~npa_f->proc~get_plasma proc~npa_f->proc~get_beam_cx_rate proc~npa_f->proc~store_npa proc~npa_f->proc~attenuate proc~npa_f->proc~mc_fastion proc~attenuate->proc~get_plasma proc~attenuate->proc~colrad proc~store_fw_photons store_fw_photons proc~store_fw_photons->proc~get_fields proc~fida_weights_mc->proc~get_fields proc~fida_weights_mc->proc~get_plasma proc~fida_weights_mc->proc~track proc~get_ep_denf get_ep_denf proc~fida_weights_mc->proc~get_ep_denf proc~cold_spec->proc~get_plasma proc~cold_spec->proc~store_photons proc~make_diagnostic_grids->proc~get_fields proc~make_diagnostic_grids->proc~track_cylindrical proc~make_diagnostic_grids->proc~track proc~make_passive_grid make_passive_grid proc~make_diagnostic_grids->proc~make_passive_grid proc~get_ep_denf->interface~interpol proc~make_passive_grid->proc~get_plasma_extrema proc~mc_fastion->proc~get_fields proc~mc_fastion->proc~get_distribution proc~store_bes_photons store_bes_photons proc~store_bes_photons->proc~store_photons proc~store_fida_photons store_fida_photons proc~store_fida_photons->proc~store_photons

Contents

Source Code

interpol1D_coeff

Source Code

subroutine interpol1D_coeff(xmin,dx,nx,xout,c,err)
    !+ Linear interpolation coefficients and index for a 1D grid y(x)
    real(Float64), intent(in)           :: xmin
        !+ Minimum abscissa value
    real(Float64), intent(in)           :: dx
        !+ Absissa spacing
    integer, intent(in)                 :: nx
        !+ Number of abscissa
    real(Float64), intent(in)           :: xout
        !+ Abscissa value to interpolate
    type(InterpolCoeffs1D), intent(out) :: c
        !+ Interpolation Coefficients
    integer, intent(out), optional      :: err
        !+ Error code

    real(Float64) :: x1, xp, b1, b2
    integer :: i, err_status

    err_status = 1
    xp = max(xout,xmin)
    i = floor((xp - xmin)/dx)+1

    if ((i.gt.0).and.(i.le.(nx-1))) then
        x1 = xmin + (i-1)*dx

        b2 = (xp - x1)/dx
        b1 = (1.0 - b2)

        c%i = i
        c%b1 = b1
        c%b2 = b2
        err_status = 0
    endif

    if(present(err)) err = err_status

end subroutine interpol1D_coeff