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

Commit

Permalink
Added in placeholder select statement to switch between fragmenation …
Browse files Browse the repository at this point in the history
…regimes
  • Loading branch information
daminton committed Aug 5, 2021
1 parent 663ad8f commit b274576
Showing 1 changed file with 13 additions and 1 deletion.
14 changes: 13 additions & 1 deletion src/symba/symba_collision.f90
Original file line number Diff line number Diff line change
Expand Up @@ -466,7 +466,19 @@ module subroutine symba_collision_resolve_fragmentations(self, system, param)
mass_res(:) = (mass_res(:) / param%MU2KG) * param%GU
Qloss = Qloss * (param%GU / param%MU2KG) * (param%TU2S / param%DU2M)**2

!status = symba_fragmentation_casemerge(system, param, family, x, v, mass, radius, L_spin, Ip)
select case (regime)
case (COLLRESOLVE_REGIME_DISRUPTION)
!status = symba_fragmentation_casedisruption(system, param, family, x, v, mass, radius, L_spin, Ip, mass_res, Qloss)
case (COLLRESOLVE_REGIME_SUPERCATASTROPHIC)
!status = symba_fragmentation_casesupercatastrophic(system, param, family, x, v, mass, radius, L_spin, Ip, mass_res, Qloss)
case (COLLRESOLVE_REGIME_HIT_AND_RUN)
!status = symba_fragmentation_casehitandrun(system, param, family, x, v, mass, radius, L_spin, Ip, mass_res, Qloss)
case (COLLRESOLVE_REGIME_MERGE, COLLRESOLVE_REGIME_GRAZE_AND_MERGE)
status = symba_fragmentation_casemerge(system, param, family, x, v, mass, radius, L_spin, Ip)
case default
write(*,*) "Error in symba_collision, unrecognized collision regime"
call util_exit(FAILURE)
end select
end do
end select
end associate
Expand Down

0 comments on commit b274576

Please sign in to comment.