test_f90.f90

Go to the documentation of this file.

program minimal
  use itm_types
  use amns_types
  use amns_module

  implicit none

  type (amns_handle_type) :: amns                            ! AMNS global handle
  type (amns_handle_rx_type)  :: amns_rx                     ! AMNS table handle
  type (amns_reaction_type) :: xx_rx
  type (amns_reactants_type) :: species
  real (kind=R8) :: te=100.0_r8, ne=1e20_r8, rate

  call itm_amns_setup(amns)                                  ! set up the AMNS system
  allocate(species%components(4))                            ! set up reactants
  species%components = (/ amns_reactant_type(6, 1, 12, 0), amns_reactant_type(1, 0, 2, 0), amns_reactant_type(6, 0, 12, 1), amns_reactant_type(1, 1, 2, 1) /)
  xx_rx%string='CX'                                          ! set up reaction
  call itm_amns_setup_table(amns, xx_rx, species, amns_rx)   ! set up table
  call itm_amns_rx(amns_rx, rate, te, ne)                    ! get results
  write(*,*) 'Rate = ', rate
  call itm_amns_finish_table(amns_rx)                        ! finish with table
  call itm_amns_finish(amns)                                 ! finish with amns

end program minimal