lebedev.f90 Source File


Source Code

!-----------------------------------------------------------------------------------------!
! Copyright (c) 2018 Peter Grünberg Institut, Forschungszentrum Jülich, Germany           !
! This file is part of Jülich KKR code and available as free software under the conditions!
! of the MIT license as expressed in the LICENSE.md file in more detail.                  !
!-----------------------------------------------------------------------------------------!

!------------------------------------------------------------------------------------
!> Summary: This subroutine includes the data from the old mesh-file lebedev
!> Author: M. Kosuth 
!> This subroutine includes the data from the old mesh-file lebedev
!------------------------------------------------------------------------------------
module mod_lebedev
  use :: mod_datatypes, only: dp
  private :: dp

contains
  !-------------------------------------------------------------------------------
  !> Summary: This subroutine includes the data from the old mesh-file lebedev
  !> Author: M. Kosuth
  !> Category: xc-potential, KKRhost 
  !> Deprecated: False 
  !> This subroutine includes the data from the old mesh-file lebedev
  !------------------------------------------------------------------------------- 
  subroutine lebedev(ij, rij1, rij2, rij3, wij)

    implicit none
    integer :: ijd
    parameter (ijd=434)
    integer :: i, ij

    real (kind=dp), dimension(ijd) :: r1, r2, r3, w
    real (kind=dp) :: rij1, rij2, rij3, wij

    data (r1(i), i=1, 30)/1.000000000000000000e00_dp, -1.000000000000000000e00_dp, 0.000000000000000000e00_dp, 0.000000000000000000e00_dp, 0.000000000000000000e00_dp, &
      0.000000000000000000e00_dp, 0.577350269189625842e00_dp, 0.577350269189625842e00_dp, 0.577350269189625842e00_dp, 0.577350269189625842e00_dp, -0.577350269189625842e00_dp, &
      -0.577350269189625842e00_dp, -0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.707106781186547462e00_dp, 0.000000000000000000e00_dp, 0.707106781186547462e00_dp, &
      0.707106781186547462e00_dp, 0.000000000000000000e00_dp, -0.707106781186547462e00_dp, -0.707106781186547462e00_dp, 0.000000000000000000e00_dp, 0.707106781186547462e00_dp, &
      -0.707106781186547462e00_dp, 0.000000000000000000e00_dp, -0.707106781186547462e00_dp, 0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, &
      0.212646824707551779e00_dp/
    data (r1(i), i=31, 60)/0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, &
      0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, &
      -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, &
      -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, &
      0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, &
      0.407712664897769694e00_dp/
    data (r1(i), i=61, 90)/0.645666470742425558e00_dp, 0.645666470742425558e00_dp, -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, &
      -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, &
      -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, 0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, &
      0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, &
      0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, &
      -0.719016501040843359e00_dp/
    data (r1(i), i=91, 120)/ -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, &
      -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, &
      0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, &
      0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, &
      -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, &
      -0.831584400419232317e00_dp/
    data (r1(i), i=121, 150)/ -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, &
      0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, &
      0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, &
      -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, &
      -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, 0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, &
      0.967987158791472790e00_dp/
    data (r1(i), i=151, 180)/0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, &
      0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, &
      -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, &
      -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, &
      0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, &
      0.994255912631277883e00_dp/
    data (r1(i), i=181, 210)/0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, &
      -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, &
      -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, 0.000000000000000000e00_dp, 0.977642811118264898e00_dp, 0.210272522857306965e00_dp, &
      0.000000000000000000e00_dp, 0.210272522857306965e00_dp, 0.977642811118264898e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, 0.210272522857306965e00_dp, &
      0.000000000000000000e00_dp, 0.210272522857306965e00_dp, -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.977642811118264898e00_dp, -0.210272522857306965e00_dp, &
      0.000000000000000000e00_dp/
    data (r1(i), i=211, 240)/ -0.210272522857306965e00_dp, 0.977642811118264898e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, -0.210272522857306965e00_dp, &
      0.000000000000000000e00_dp, -0.210272522857306965e00_dp, -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, 0.471598691151315919e00_dp, &
      0.000000000000000000e00_dp, 0.471598691151315919e00_dp, 0.881813287779428800e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, 0.471598691151315919e00_dp, &
      0.000000000000000000e00_dp, 0.471598691151315919e00_dp, -0.881813287779428800e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, -0.471598691151315919e00_dp, &
      0.000000000000000000e00_dp, -0.471598691151315919e00_dp, 0.881813287779428800e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, -0.471598691151315919e00_dp, &
      0.000000000000000000e00_dp/
    data (r1(i), i=241, 270)/ -0.471598691151315919e00_dp, -0.881813287779428800e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, &
      0.334436314534345469e00_dp, 0.334436314534345469e00_dp, 0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, &
      0.334436314534345469e00_dp, 0.334436314534345469e00_dp, -0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, &
      -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, 0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, &
      -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, &
      0.334436314534345469e00_dp/
    data (r1(i), i=271, 300)/0.334436314534345469e00_dp, 0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, &
      0.334436314534345469e00_dp, 0.334436314534345469e00_dp, -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, &
      -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, 0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, &
      -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, -0.937180985855372239e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, 0.868946032287241210e00_dp, &
      0.450233038258262552e00_dp, 0.450233038258262552e00_dp, 0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, -0.868946032287241210e00_dp, &
      0.450233038258262552e00_dp/
    data (r1(i), i=301, 330)/0.450233038258262552e00_dp, -0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, 0.868946032287241210e00_dp, &
      -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, 0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, -0.868946032287241210e00_dp, &
      -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, 0.868946032287241210e00_dp, &
      0.450233038258262552e00_dp, 0.450233038258262552e00_dp, 0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, &
      0.450233038258262552e00_dp, 0.450233038258262552e00_dp, -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, 0.868946032287241210e00_dp, &
      -0.450233038258262552e00_dp/
    data (r1(i), i=331, 360)/ -0.450233038258262552e00_dp, 0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, &
      -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, -0.868946032287241210e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, 0.771746262691590079e00_dp, &
      0.555015236107680665e00_dp, 0.555015236107680665e00_dp, 0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, -0.771746262691590079e00_dp, &
      0.555015236107680665e00_dp, 0.555015236107680665e00_dp, -0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, 0.771746262691590079e00_dp, &
      -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, 0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, -0.771746262691590079e00_dp, &
      -0.555015236107680665e00_dp/
    data (r1(i), i=361, 390)/ -0.555015236107680665e00_dp, -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, 0.771746262691590079e00_dp, &
      0.555015236107680665e00_dp, 0.555015236107680665e00_dp, 0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, &
      0.555015236107680665e00_dp, 0.555015236107680665e00_dp, -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, 0.771746262691590079e00_dp, &
      -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, 0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, &
      -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, -0.771746262691590079e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, 0.799927854385728554e00_dp, &
      0.590515704892527138e00_dp/
    data (r1(i), i=391, 420)/0.590515704892527138e00_dp, 0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, -0.799927854385728554e00_dp, &
      0.590515704892527138e00_dp, 0.590515704892527138e00_dp, -0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, 0.799927854385728554e00_dp, &
      -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, 0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, -0.799927854385728554e00_dp, &
      -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, -0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, 0.799927854385728554e00_dp, &
      0.590515704892527138e00_dp, 0.590515704892527138e00_dp, 0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, -0.799927854385728554e00_dp, &
      0.590515704892527138e00_dp/
    data (r1(i), i=421, ijd)/0.590515704892527138e00_dp, -0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, 0.799927854385728554e00_dp, &
      -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, 0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, -0.799927854385728554e00_dp, &
      -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, -0.799927854385728554e00_dp/

    data (r2(i), i=1, 30)/0.000000000000000000e00_dp, 0.000000000000000000e00_dp, 1.000000000000000000e00_dp, -1.000000000000000000e00_dp, 0.000000000000000000e00_dp, &
      0.000000000000000000e00_dp, 0.577350269189625842e00_dp, 0.577350269189625842e00_dp, -0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.577350269189625842e00_dp, &
      0.577350269189625842e00_dp, -0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.707106781186547462e00_dp, 0.707106781186547462e00_dp, 0.000000000000000000e00_dp, &
      -0.707106781186547462e00_dp, 0.707106781186547462e00_dp, 0.000000000000000000e00_dp, 0.707106781186547462e00_dp, -0.707106781186547462e00_dp, 0.000000000000000000e00_dp, &
      -0.707106781186547462e00_dp, -0.707106781186547462e00_dp, 0.000000000000000000e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, 0.690934630750911105e00_dp, &
      0.690934630750911105e00_dp/
    data (r2(i), i=31, 60)/0.212646824707551779e00_dp, 0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, &
      -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, 0.690934630750911105e00_dp, &
      0.690934630750911105e00_dp, 0.212646824707551779e00_dp, 0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, &
      -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, -0.690934630750911105e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, &
      0.645666470742425558e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, &
      -0.645666470742425558e00_dp/
    data (r2(i), i=61, 90)/ -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, &
      0.645666470742425558e00_dp, 0.407712664897769694e00_dp, 0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, &
      -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, -0.645666470742425558e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, 0.491434263778474656e00_dp, &
      0.491434263778474656e00_dp, 0.719016501040843359e00_dp, 0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, &
      -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, 0.491434263778474656e00_dp, &
      0.491434263778474656e00_dp/
    data (r2(i), i=91, 120)/0.719016501040843359e00_dp, 0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, &
      -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, -0.491434263778474656e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, &
      0.392725976336800231e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, &
      -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, &
      0.392725976336800231e00_dp, 0.831584400419232317e00_dp, 0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, &
      -0.392725976336800231e00_dp/
    data (r2(i), i=121, 150)/ -0.831584400419232317e00_dp, -0.392725976336800231e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, &
      0.286128901030763827e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, &
      -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, &
      0.286128901030763827e00_dp, 0.914472801120872480e00_dp, 0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, &
      -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, -0.286128901030763827e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, 0.177483605460915767e00_dp, &
      0.177483605460915767e00_dp/
    data (r2(i), i=151, 180)/0.967987158791472790e00_dp, 0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, &
      -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, 0.177483605460915767e00_dp, &
      0.177483605460915767e00_dp, 0.967987158791472790e00_dp, 0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, &
      -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, -0.177483605460915767e00_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, &
      0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, &
      -0.756808436717801847e-1_dp/
    data (r2(i), i=181, 210)/ -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, &
      0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, &
      -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, -0.756808436717801847e-1_dp, 0.210272522857306965e00_dp, 0.000000000000000000e00_dp, 0.977642811118264898e00_dp, &
      0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.210272522857306965e00_dp, 0.210272522857306965e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, &
      -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.210272522857306965e00_dp, -0.210272522857306965e00_dp, 0.000000000000000000e00_dp, 0.977642811118264898e00_dp, &
      0.977642811118264898e00_dp/
    data (r2(i), i=211, 240)/0.000000000000000000e00_dp, -0.210272522857306965e00_dp, -0.210272522857306965e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, &
      -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, -0.210272522857306965e00_dp, 0.471598691151315919e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, &
      0.881813287779428800e00_dp, 0.000000000000000000e00_dp, 0.471598691151315919e00_dp, 0.471598691151315919e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, &
      -0.881813287779428800e00_dp, 0.000000000000000000e00_dp, 0.471598691151315919e00_dp, -0.471598691151315919e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, &
      0.881813287779428800e00_dp, 0.000000000000000000e00_dp, -0.471598691151315919e00_dp, -0.471598691151315919e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, &
      -0.881813287779428800e00_dp/
    data (r2(i), i=241, 270)/0.000000000000000000e00_dp, -0.471598691151315919e00_dp, 0.334436314534345469e00_dp, 0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, &
      0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, -0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, &
      0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, 0.334436314534345469e00_dp, -0.334436314534345469e00_dp, 0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, &
      0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, -0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, &
      0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, -0.334436314534345469e00_dp, 0.334436314534345469e00_dp, 0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, &
      -0.992176963642923682e-1_dp/
    data (r2(i), i=271, 300)/0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, &
      -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, 0.334436314534345469e00_dp, -0.334436314534345469e00_dp, 0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, &
      -0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, &
      -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, -0.334436314534345469e00_dp, 0.450233038258262552e00_dp, 0.868946032287241210e00_dp, 0.205482369640304363e00_dp, &
      0.205482369640304363e00_dp, 0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, -0.868946032287241210e00_dp, 0.205482369640304363e00_dp, &
      0.205482369640304363e00_dp/
    data (r2(i), i=301, 330)/ -0.868946032287241210e00_dp, 0.450233038258262552e00_dp, -0.450233038258262552e00_dp, 0.868946032287241210e00_dp, 0.205482369640304363e00_dp, &
      0.205482369640304363e00_dp, 0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, -0.868946032287241210e00_dp, 0.205482369640304363e00_dp, &
      0.205482369640304363e00_dp, -0.868946032287241210e00_dp, -0.450233038258262552e00_dp, 0.450233038258262552e00_dp, 0.868946032287241210e00_dp, -0.205482369640304363e00_dp, &
      -0.205482369640304363e00_dp, 0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, &
      -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, 0.450233038258262552e00_dp, -0.450233038258262552e00_dp, 0.868946032287241210e00_dp, -0.205482369640304363e00_dp, &
      -0.205482369640304363e00_dp/
    data (r2(i), i=331, 360)/0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, &
      -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, -0.450233038258262552e00_dp, 0.555015236107680665e00_dp, 0.771746262691590079e00_dp, 0.310428403516654128e00_dp, &
      0.310428403516654128e00_dp, 0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, -0.771746262691590079e00_dp, 0.310428403516654128e00_dp, &
      0.310428403516654128e00_dp, -0.771746262691590079e00_dp, 0.555015236107680665e00_dp, -0.555015236107680665e00_dp, 0.771746262691590079e00_dp, 0.310428403516654128e00_dp, &
      0.310428403516654128e00_dp, 0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, -0.771746262691590079e00_dp, 0.310428403516654128e00_dp, &
      0.310428403516654128e00_dp/
    data (r2(i), i=361, 390)/ -0.771746262691590079e00_dp, -0.555015236107680665e00_dp, 0.555015236107680665e00_dp, 0.771746262691590079e00_dp, -0.310428403516654128e00_dp, &
      -0.310428403516654128e00_dp, 0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, &
      -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, 0.555015236107680665e00_dp, -0.555015236107680665e00_dp, 0.771746262691590079e00_dp, -0.310428403516654128e00_dp, &
      -0.310428403516654128e00_dp, 0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, &
      -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, -0.555015236107680665e00_dp, 0.590515704892527138e00_dp, 0.799927854385728554e00_dp, 0.106801826075804823e00_dp, &
      0.106801826075804823e00_dp/
    data (r2(i), i=391, 420)/0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, -0.799927854385728554e00_dp, 0.106801826075804823e00_dp, &
      0.106801826075804823e00_dp, -0.799927854385728554e00_dp, 0.590515704892527138e00_dp, -0.590515704892527138e00_dp, 0.799927854385728554e00_dp, 0.106801826075804823e00_dp, &
      0.106801826075804823e00_dp, 0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, -0.799927854385728554e00_dp, 0.106801826075804823e00_dp, &
      0.106801826075804823e00_dp, -0.799927854385728554e00_dp, -0.590515704892527138e00_dp, 0.590515704892527138e00_dp, 0.799927854385728554e00_dp, -0.106801826075804823e00_dp, &
      -0.106801826075804823e00_dp, 0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, -0.799927854385728554e00_dp, -0.106801826075804823e00_dp, &
      -0.106801826075804823e00_dp/
    data (r2(i), i=421, ijd)/ -0.799927854385728554e00_dp, 0.590515704892527138e00_dp, -0.590515704892527138e00_dp, 0.799927854385728554e00_dp, -0.106801826075804823e00_dp, &
      -0.106801826075804823e00_dp, 0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, -0.799927854385728554e00_dp, -0.106801826075804823e00_dp, &
      -0.106801826075804823e00_dp, -0.799927854385728554e00_dp, -0.590515704892527138e00_dp/

    data (r3(i), i=1, 30)/0.000000000000000000e00_dp, 0.000000000000000000e00_dp, 0.000000000000000000e00_dp, 0.000000000000000000e00_dp, 1.000000000000000000e00_dp, &
      -1.000000000000000000e00_dp, 0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.577350269189625842e00_dp, &
      -0.577350269189625842e00_dp, 0.577350269189625842e00_dp, -0.577350269189625842e00_dp, 0.000000000000000000e00_dp, 0.707106781186547462e00_dp, 0.707106781186547462e00_dp, &
      0.000000000000000000e00_dp, -0.707106781186547462e00_dp, 0.707106781186547462e00_dp, 0.000000000000000000e00_dp, 0.707106781186547462e00_dp, -0.707106781186547462e00_dp, &
      0.000000000000000000e00_dp, -0.707106781186547462e00_dp, -0.707106781186547462e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, &
      -0.690934630750911105e00_dp/
    data (r3(i), i=31, 60)/ -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, &
      -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, &
      -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, 0.690934630750911105e00_dp, 0.690934630750911105e00_dp, 0.212646824707551779e00_dp, &
      -0.690934630750911105e00_dp, -0.690934630750911105e00_dp, -0.212646824707551779e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, &
      -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, &
      -0.645666470742425558e00_dp/
    data (r3(i), i=61, 90)/ -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, &
      -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, 0.645666470742425558e00_dp, 0.645666470742425558e00_dp, 0.407712664897769694e00_dp, &
      -0.645666470742425558e00_dp, -0.645666470742425558e00_dp, -0.407712664897769694e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, &
      -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, &
      -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, &
      -0.491434263778474656e00_dp/
    data (r3(i), i=91, 120)/ -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, 0.491434263778474656e00_dp, 0.491434263778474656e00_dp, 0.719016501040843359e00_dp, &
      -0.491434263778474656e00_dp, -0.491434263778474656e00_dp, -0.719016501040843359e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, &
      -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, &
      -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, &
      -0.392725976336800231e00_dp, -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, 0.392725976336800231e00_dp, 0.392725976336800231e00_dp, 0.831584400419232317e00_dp, &
      -0.392725976336800231e00_dp/
    data (r3(i), i=121, 150)/ -0.392725976336800231e00_dp, -0.831584400419232317e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, &
      -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, &
      -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, &
      -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, 0.286128901030763827e00_dp, 0.286128901030763827e00_dp, 0.914472801120872480e00_dp, &
      -0.286128901030763827e00_dp, -0.286128901030763827e00_dp, -0.914472801120872480e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, &
      -0.177483605460915767e00_dp/
    data (r3(i), i=151, 180)/ -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, &
      -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, &
      -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, 0.177483605460915767e00_dp, 0.177483605460915767e00_dp, 0.967987158791472790e00_dp, &
      -0.177483605460915767e00_dp, -0.177483605460915767e00_dp, -0.967987158791472790e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, &
      -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, &
      -0.756808436717801847e-1_dp/
    data (r3(i), i=181, 210)/ -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, &
      -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, 0.756808436717801847e-1_dp, 0.756808436717801847e-1_dp, 0.994255912631277883e00_dp, &
      -0.756808436717801847e-1_dp, -0.756808436717801847e-1_dp, -0.994255912631277883e00_dp, 0.977642811118264898e00_dp, 0.210272522857306965e00_dp, 0.000000000000000000e00_dp, &
      0.210272522857306965e00_dp, 0.977642811118264898e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, 0.210272522857306965e00_dp, 0.000000000000000000e00_dp, &
      0.210272522857306965e00_dp, -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.977642811118264898e00_dp, -0.210272522857306965e00_dp, 0.000000000000000000e00_dp, &
      -0.210272522857306965e00_dp/
    data (r3(i), i=211, 240)/0.977642811118264898e00_dp, 0.000000000000000000e00_dp, -0.977642811118264898e00_dp, -0.210272522857306965e00_dp, 0.000000000000000000e00_dp, &
      -0.210272522857306965e00_dp, -0.977642811118264898e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, 0.471598691151315919e00_dp, 0.000000000000000000e00_dp, &
      0.471598691151315919e00_dp, 0.881813287779428800e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, 0.471598691151315919e00_dp, 0.000000000000000000e00_dp, &
      0.471598691151315919e00_dp, -0.881813287779428800e00_dp, 0.000000000000000000e00_dp, 0.881813287779428800e00_dp, -0.471598691151315919e00_dp, 0.000000000000000000e00_dp, &
      -0.471598691151315919e00_dp, 0.881813287779428800e00_dp, 0.000000000000000000e00_dp, -0.881813287779428800e00_dp, -0.471598691151315919e00_dp, 0.000000000000000000e00_dp, &
      -0.471598691151315919e00_dp/
    data (r3(i), i=241, 270)/ -0.881813287779428800e00_dp, 0.000000000000000000e00_dp, 0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, &
      0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, &
      -0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, &
      0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, &
      -0.937180985855372239e00_dp, 0.992176963642923682e-1_dp, 0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, &
      0.937180985855372239e00_dp/
    data (r3(i), i=271, 300)/ -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, 0.334436314534345469e00_dp, 0.334436314534345469e00_dp, &
      -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, 0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, &
      0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, -0.937180985855372239e00_dp, -0.334436314534345469e00_dp, -0.334436314534345469e00_dp, &
      -0.937180985855372239e00_dp, -0.992176963642923682e-1_dp, -0.992176963642923682e-1_dp, 0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, &
      0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, -0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, &
      -0.868946032287241210e00_dp/
    data (r3(i), i=301, 330)/0.205482369640304363e00_dp, 0.205482369640304363e00_dp, 0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, &
      0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, -0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, &
      -0.868946032287241210e00_dp, 0.205482369640304363e00_dp, 0.205482369640304363e00_dp, 0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, &
      0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, 0.450233038258262552e00_dp, 0.450233038258262552e00_dp, &
      -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, 0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, &
      0.868946032287241210e00_dp/
    data (r3(i), i=331, 360)/ -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, -0.868946032287241210e00_dp, -0.450233038258262552e00_dp, -0.450233038258262552e00_dp, &
      -0.868946032287241210e00_dp, -0.205482369640304363e00_dp, -0.205482369640304363e00_dp, 0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, &
      0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, -0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, &
      -0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, 0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, &
      0.771746262691590079e00_dp, 0.310428403516654128e00_dp, 0.310428403516654128e00_dp, -0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, &
      -0.771746262691590079e00_dp/
    data (r3(i), i=361, 390)/0.310428403516654128e00_dp, 0.310428403516654128e00_dp, 0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, &
      0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, 0.555015236107680665e00_dp, 0.555015236107680665e00_dp, &
      -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, 0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, &
      0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, -0.771746262691590079e00_dp, -0.555015236107680665e00_dp, -0.555015236107680665e00_dp, &
      -0.771746262691590079e00_dp, -0.310428403516654128e00_dp, -0.310428403516654128e00_dp, 0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, &
      0.799927854385728554e00_dp/
    data (r3(i), i=391, 420)/0.106801826075804823e00_dp, 0.106801826075804823e00_dp, -0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, &
      -0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, 0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, &
      0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, -0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, &
      -0.799927854385728554e00_dp, 0.106801826075804823e00_dp, 0.106801826075804823e00_dp, 0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, &
      0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, -0.799927854385728554e00_dp, 0.590515704892527138e00_dp, 0.590515704892527138e00_dp, &
      -0.799927854385728554e00_dp/
    data (r3(i), i=421, ijd)/ -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, 0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, &
      0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp, -0.799927854385728554e00_dp, -0.590515704892527138e00_dp, -0.590515704892527138e00_dp, &
      -0.799927854385728554e00_dp, -0.106801826075804823e00_dp, -0.106801826075804823e00_dp/

    data (w(i), i=1, 30)/0.526589796822443581e-3_dp, 0.526589796822443581e-3_dp, 0.526589796822443581e-3_dp, 0.526589796822443581e-3_dp, 0.526589796822443581e-3_dp, &
      0.526589796822443581e-3_dp, 0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, &
      0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, 0.251231741892730736e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, &
      0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, &
      0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.254821997200260731e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, &
      0.253040380118635519e-2_dp/
    data (w(i), i=31, 60)/0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, &
      0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, &
      0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, &
      0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.253040380118635519e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, &
      0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, &
      0.251326717459756429e-2_dp/
    data (w(i), i=61, 90)/0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, &
      0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, &
      0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.251326717459756429e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, &
      0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, &
      0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, &
      0.250172516840293591e-2_dp/
    data (w(i), i=91, 120)/0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, &
      0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.250172516840293591e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, &
      0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, &
      0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, &
      0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, &
      0.244537343731297992e-2_dp/
    data (w(i), i=121, 150)/0.244537343731297992e-2_dp, 0.244537343731297992e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, &
      0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, &
      0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, &
      0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, &
      0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.230269478222741575e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, &
      0.201427902091852809e-2_dp/
    data (w(i), i=151, 180)/0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, &
      0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, &
      0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, &
      0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.201427902091852809e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, &
      0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, &
      0.146249562159461374e-2_dp/
    data (w(i), i=181, 210)/0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, &
      0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, &
      0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.146249562159461374e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, &
      0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, &
      0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, &
      0.191095128217953227e-2_dp/
    data (w(i), i=211, 240)/0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, &
      0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.191095128217953227e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, &
      0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, &
      0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, &
      0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, &
      0.241744237563898097e-2_dp/
    data (w(i), i=241, 270)/0.241744237563898097e-2_dp, 0.241744237563898097e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp/
    data (w(i), i=271, 300)/0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, &
      0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.223660776043784881e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp/
    data (w(i), i=301, 330)/0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp/
    data (w(i), i=331, 360)/0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, &
      0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.241693004432477548e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp/
    data (w(i), i=361, 390)/0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, &
      0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.249664405455308613e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp/
    data (w(i), i=391, 420)/0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp/
    data (w(i), i=421, ijd)/0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, &
      0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp, 0.251223685456349521e-2_dp/

    if (ij>ijd) then
      write (*, *) 'IN < LEBEDEV >: IJ = ', ij, ' > IJD = ', ijd
      stop
    end if

    rij1 = r1(ij)
    rij2 = r2(ij)
    rij3 = r3(ij)
    wij = w(ij)

    return
  end subroutine lebedev

end module mod_lebedev