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_sync_spin_to_orbit
21 :
22 : ! NOTE: remember to set true:
23 : ! use_other_sync_spin_to_orbit = .true.
24 :
25 : implicit none
26 :
27 : contains
28 :
29 0 : subroutine null_other_sync_spin_to_orbit(id, nz, osep, qratio, rl, dt_next, Ftid, sync_type, sync_mode, ierr)
30 : use const_def, only: dp, strlen
31 : use binary_def, only : binary_info, binary_ptr
32 : use star_def, only : star_info, star_ptr
33 : integer, intent(in) :: id
34 : integer, intent(in) :: nz
35 : real(dp), intent(in) :: osep ! orbital separation (cm)
36 : real(dp), intent(in) :: qratio ! mass_other_star/mass_this_star
37 : real(dp), intent(in) :: rl ! roche lobe radius (cm)
38 : real(dp), intent(in) :: dt_next ! next timestep
39 : real(dp), intent(in) :: Ftid ! efficiency of tidal synchronization. (time scale / Ftid ).
40 :
41 : character (len=strlen), intent(in) :: sync_type ! synchronization timescale
42 : character (len=strlen), intent(in) :: sync_mode ! where to put/take angular momentum
43 : integer, intent(out) :: ierr
44 :
45 : type (star_info), pointer :: s
46 : integer :: k
47 :
48 0 : call star_ptr(id, s, ierr)
49 0 : if (ierr /= 0) then
50 0 : write(*,*) 'failed in binary_ptr'
51 : return
52 : end if
53 :
54 0 : do k=1,nz
55 : s% extra_jdot(k) = s% extra_jdot(k) - 0d0
56 : end do
57 :
58 0 : write(*,*) "Warning: no implementation for other_sync_spin_to_orbit"
59 0 : write(*,*) "Not tidal torques are being included"
60 0 : end subroutine null_other_sync_spin_to_orbit
61 :
62 : end module mod_other_sync_spin_to_orbit
63 :
|