Line data    Source code

       1              : ! ***********************************************************************
       2              : !
       3              : !   Copyright (C) 2012  Bill Paxton, Pablo Marchant & 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 mod_other_adjust_mdots
      21              : 
      22              :       ! NOTE: remember to set true:
      23              :       ! use_other_adjust_mdots = .true.
      24              : 
      25              :       implicit none
      26              : 
      27              :       contains
      28              : 
      29            0 :       subroutine null_other_adjust_mdots(binary_id, ierr)
      30              :          use binary_def, only : binary_info, binary_ptr
      31              :          use const_def, only: dp
      32              :          integer, intent(in) :: binary_id
      33              :          integer, intent(out) :: ierr
      34              :          type (binary_info), pointer :: b
      35              :          ierr = 0
      36            0 :          call binary_ptr(binary_id, b, ierr)
      37            0 :          if (ierr /= 0) then
      38            0 :             write(*,*) 'failed in binary_ptr'
      39              :             return
      40              :          end if
      41            0 :          b% fixed_xfer_fraction = 0d0
      42              :          ! note that you should also eval mdot_edd here and the wind mass transfer eff
      43              :          ! and angular momentum accretion. You can call the default functions using
      44              :          ! the ones provided through binary_lib
      45            0 :          b% mdot_edd = 0d0
      46            0 :          b% mdot_edd_eta = 0d0
      47            0 :          b% mdot_wind_transfer(:) = 0d0
      48            0 :          b% wind_xfer_fraction(:) = 0d0
      49            0 :          b% s_donor% mstar_dot = 0d0
      50            0 :          if (b% point_mass_i == 0) then
      51            0 :             b% component_mdot(b% a_i) = 0d0
      52              :          else
      53            0 :             b% s_accretor% mstar_dot = 0d0
      54              :          end if
      55            0 :          b% accretion_luminosity = 0d0
      56            0 :          b% mdot_system_transfer(:) = 0d0
      57            0 :          b% mdot_system_cct = 0d0
      58              : 
      59            0 :       end subroutine null_other_adjust_mdots
      60              : 
      61              :       end module mod_other_adjust_mdots
      62              :