!--------------------------------------------------------------------
! This program computes Rate Constants for transitions in H2O + H2O system.
! For more details and citing this work, please refer to:
! Bikramaditya Mandal and Dmitri Babikov, Astronomy and Astrphysics, 2023.
! This is the main driver program.
! User should change the value of variable "req_temp" to the desired one.
! For example, the current value is set to 300 K.
! The Rate Constants for all transitions are printed in a file named "Rates_for_T_300.00000_K.out".
! The unit of input Temterature is K and output TACS is cm^3/s.
!--------------------------------------------------------------------

!--------------------------------------------------------------------
! A driver program to execute the above mentioned modules
!--------------------------------------------------------------------

program obtain_sigma
  use module_engines
  implicit none
  integer i
  real*8 :: req_temp
  character (len = 1) :: input_data_type

  req_temp = 300.d0
  input_data_type = 'a'

  call interpolate(req_temp, input_data_type)
end program