Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Browse files Browse the repository at this point in the history 
Scaled the quantities in symba_frag_pos to keep things of O(~1)
  • Loading branch information
@daminton
daminton committed May 15, 2021
1 parent 4041b4d commit 1416a5f
Show file tree
Hide file tree
Showing 7 changed files with 98 additions and 104 deletions.
18 changes: 9 additions & 9 deletions src/modules/module_interfaces.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -723,8 +723,8 @@ function symba_casedisruption (symba_plA, family, nmergeadd, mergeadd_list, x, v
integer(I4B), dimension(:), intent(in) :: family
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
real(DP), dimension(:,:), intent(in) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(in) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(inout) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, mass_res
type(user_input_parameters),intent(inout) :: param
real(DP), intent(in) :: Qloss
integer(I4B) :: status
Expand All @@ -741,8 +741,8 @@ function symba_casehitandrun (symba_plA, family, nmergeadd, mergeadd_list, name,
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
integer(I4B), dimension(:), intent(in) :: name
real(DP), dimension(:,:), intent(in) :: x, v, Lspin, Ip
real(DP), dimension(:), intent(in) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(inout) :: x, v, Lspin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, mass_res
type(user_input_parameters),intent(inout) :: param
real(DP), intent(in) :: Qloss
integer(I4B) :: status
Expand All @@ -757,8 +757,8 @@ function symba_casemerge (symba_plA, family, nmergeadd, mergeadd_list, x, v, mas
integer(I4B), dimension(:), intent(in) :: family
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
real(DP), dimension(:,:), intent(in) :: x, v, lspin, Ip
real(DP), dimension(:), intent(in) :: mass, radius
real(DP), dimension(:,:), intent(inout) :: x, v, lspin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius
type(user_input_parameters),intent(inout) :: param
integer(I4B) :: status
end function symba_casemerge
Expand Down Expand Up @@ -930,9 +930,9 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi
type(symba_pl), intent(inout) :: symba_plA
integer(I4B), dimension(:), intent(in) :: family
real(DP), intent(in) :: Qloss
real(DP), dimension(:,:), intent(in) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(in) :: mass, radius, m_frag, rad_frag
real(DP), dimension(:,:), intent(in) :: Ip_frag
real(DP), dimension(:,:), intent(inout) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, m_frag, rad_frag
real(DP), dimension(:,:), intent(inout) :: Ip_frag
real(DP), dimension(:,:), intent(out) :: xb_frag, vb_frag, rot_frag
logical, intent(out) :: lmerge
end subroutine symba_frag_pos
Expand Down
4 changes: 2 additions & 2 deletions src/symba/symba_casedisruption.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,8 +16,8 @@ function symba_casedisruption (symba_plA, family, nmergeadd, mergeadd_list, x, v
integer(I4B), dimension(:), intent(in) :: family
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
real(DP), dimension(:), intent(in) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(in) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(inout) :: x, v, L_spin, Ip
type(user_input_parameters),intent(inout) :: param
real(DP), intent(inout) :: Qloss
! Result
Expand Down
4 changes: 2 additions & 2 deletions src/symba/symba_casehitandrun.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,8 +16,8 @@ function symba_casehitandrun(symba_plA, family, nmergeadd, mergeadd_list, id, x,
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
integer(I4B), dimension(:), intent(in) :: id
real(DP), dimension(:,:), intent(in) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(in) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(inout) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, mass_res
type(user_input_parameters),intent(inout) :: param
real(DP), intent(inout) :: Qloss
! Result
Expand Down
4 changes: 2 additions & 2 deletions src/symba/symba_casesupercatastrophic.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,8 +15,8 @@ function symba_casesupercatastrophic (symba_plA, family, nmergeadd, mergeadd_lis
integer(I4B), dimension(:), intent(in) :: family
integer(I4B), intent(inout) :: nmergeadd
type(symba_merger), intent(inout) :: mergeadd_list
real(DP), dimension(:,:), intent(in) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(in) :: mass, radius, mass_res
real(DP), dimension(:,:), intent(inout) :: x, v, L_spin, Ip
real(DP), dimension(:), intent(inout) :: mass, radius, mass_res
type(user_input_parameters),intent(inout) :: param
real(DP), intent(inout) :: Qloss
! Result
Expand Down
Loading

0 comments on commit 1416a5f

Please sign in to comment.