Line data    Source code

       1              : ! ***********************************************************************
       2              : !
       3              : !   Copyright (C) 2012  The MESA Team
       4              : !
       5              : !   This program is free software: you can redistribute it and/or modify
       6              : !   it under the terms of the GNU Lesser General Public License
       7              : !   as published by the Free Software Foundation,
       8              : !   either version 3 of the License, or (at your option) any later version.
       9              : !
      10              : !   This program is distributed in the hope that it will be useful,
      11              : !   but WITHOUT ANY WARRANTY; without even the implied warranty of
      12              : !   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
      13              : !   See the GNU Lesser General Public License for more details.
      14              : !
      15              : !   You should have received a copy of the GNU Lesser General Public License
      16              : !   along with this program. If not, see <https://www.gnu.org/licenses/>.
      17              : !
      18              : ! ***********************************************************************
      19              : 
      20              : module other_neu
      21              : 
      22              :    ! consult star/other/README for general usage instructions
      23              :    ! control name: use_other_neu = .true.
      24              :    ! procedure pointer: s% other_neu => my_routine
      25              : 
      26              :    implicit none
      27              : 
      28              : contains
      29              : 
      30            0 :    subroutine null_other_neu( &
      31              :       id, k, T, log10_T, Rho, log10_Rho, abar, zbar, log10_Tlim, flags, &
      32              :       loss, sources, ierr)
      33              :       use star_def
      34              :       use neu_lib, only: neu_get
      35              :       use neu_def
      36              :       integer, intent(in) :: id  ! id for star
      37              :       integer, intent(in) :: k  ! cell number or 0 if not for a particular cell
      38              :       real(dp), intent(in) :: T  ! temperature
      39              :       real(dp), intent(in) :: log10_T  ! log10 of temperature
      40              :       real(dp), intent(in) :: Rho  ! density
      41              :       real(dp), intent(in) :: log10_Rho  ! log10 of density
      42              :       real(dp), intent(in) :: abar  ! mean atomic weight
      43              :       real(dp), intent(in) :: zbar  ! mean charge
      44              :       real(dp), intent(in) :: log10_Tlim
      45              :       logical, intent(inout) :: flags(num_neu_types)  ! true if should include the type of loss
      46              :       real(dp), intent(inout) :: loss(num_neu_rvs)  ! total from all sources
      47              :       real(dp), intent(inout) :: sources(num_neu_types, num_neu_rvs)
      48              :       integer, intent(out) :: ierr
      49              :       call neu_get( &
      50              :          T, log10_T, Rho, log10_Rho, abar, zbar, log10_Tlim, flags, &
      51            0 :          loss, sources, ierr)
      52            0 :    end subroutine null_other_neu
      53              : 
      54              : end module other_neu
      55              :