sigma_v_calculation_2H(a, g)6Li.f90

program sigma_v_calculation_2H_a_g_6Li
    ! This program calculates sigma_v depending on given sigma or s_factor     
        implicit none
       
        ! Type declarations
           double precision :: new_sigma_v
           
           double precision, dimension(60) :: T9, dLi

           integer :: i
           
    
           T9 = (/0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009,0.010, 0.011, 0.012,0.013, 0.014, & 
           0.015, 0.016, 0.018, 0.020, 0.025, 0.030, 0.040, 0.050, 0.060, 0.070, 0.080, 0.090, 0.100, 0.110, 0.120, 0.130, &
           0.140, 0.150, 0.160, 0.180, 0.200, 0.250, 0.300, 0.350, 0.400, 0.450, 0.500, 0.600, 0.700, 0.800, 0.900, 1.000, &
           1.250, 1.500, 1.750, 2.000, 2.500, 3.000, 3.500, 4.000, 5.000, 6.000, 7.000, 8.000, 9.000, 10.00 /)
    
           dLi = (/ 0.9153e-29, 0.2610e-22,0.3458e-9, 0.3190e-17,  0.7929e-6, 0.9163e-15,0.7672e-14,0.4990e-13, 0.2100e-12, &
           0.6547e-12, 0.1655e-11, 0.3612e-11,0.7142e-11,0.1325e-10,0.2363e-10,0.4103e-10,0.1157e-09,0.2965e-09,0.2014e-08, &
           0.8452e-08, 0.6594e-07, 0.2827e-06, 0.8598e-06, 0.2094e-05, 0.4372e-05,0.8156e-05, 0.1397e-04, 0.2240e-04, & 
           0.3406e-04, 0.4959e-04, & 
           0.6967e-04, 0.9495e-04, 0.1261e-03, 0.2090e-03, 0.3237e-03, 0.7846e-03, 0.1557e-02, 0.2715e-02, 0.4325e-02, &
           0.6453e-02,0.9169e-02,0.1674e-01, 0.2813e-01, 0.4502e-01, 0.6944e-01, 0.1033e+00, 0.2359e+00, 0.4350e+00, &
           0.6839e+00, 0.9623e+00,0.1549e+01,0.2132e+01, 0.2705e+01, 0.3280e+01, 0.4476e+01, 0.5754e+01, 0.7088e+01, &
           0.8438e+01, 0.9773e+01,0.1107e+02 /)

           do i= 1, size(T9)
       
            open(1, file = 'sigma_v_2H(a,g)6Li.dat', status = 'unknown')  
            write(1,*) T9(i), dLi(i)


            new_sigma_v  = 3.01e+01*(T9(i))**(-(2.0/3.0))*EXP(-7.423/((T9(i))**(1.0/3.0))) &
            *(1.+.056*((T9(i))**(1.0/3.0))-4.85*(T9(i))**(2.0/3.0)+8.85*(T9(i))-.585*((T9(i))**(4.0/3.0)) &
            -.584*((T9(i))**(5.0/3.0))) &
            + 8.55e+1*(T9(i))**(-(3.0/2.0))*EXP(-8.228/(T9(i)))
   
            open(2, file = 'new_sigma_v_2H(a,g)6Li.dat', status = 'unknown')  
            write(2,*) T9(i), new_sigma_v
   

         end do
        
    end program sigma_v_calculation_2H_a_g_6Li