Fixing bugs in Fraggle

examples/Fragmentation/Fragmentation_Movie.py

@@ -206,7 +206,7 @@ def data_stream(self, frame=0):
         # Set fragmentation parameters
         minimum_fragment_gmass = 0.2 * body_Gmass[style][1] # Make the minimum fragment mass a fraction of the smallest body
         gmtiny = 0.99 * body_Gmass[style][1] # Make GMTINY just smaller than the smallest original body. This will prevent runaway collisional cascades
-        sim.set_parameter(collision_model="bounce", encounter_save="both", gmtiny=gmtiny, minimum_fragment_gmass=minimum_fragment_gmass, verbose=False)
+        sim.set_parameter(collision_model="fraggle", encounter_save="both", gmtiny=gmtiny, minimum_fragment_gmass=minimum_fragment_gmass, verbose=False)
         sim.run(dt=1e-3, tstop=1.0e-3, istep_out=1, dump_cadence=0)
 
         print("Generating animation")

src/fraggle/fraggle_generate.f90

@@ -252,7 +252,6 @@ module subroutine fraggle_generate_fragments(collider, nbody_system, param, lfai
             return
          end if
 
-
          ! This is a factor that will "distort" the shape of the frgment cloud in the direction of the impact velocity 
          f_spin= .mag. (runit(:) .cross. vunit(:)) 
 

src/fraggle/fraggle_setup.f90

@@ -22,8 +22,8 @@ module subroutine fraggle_setup_fragments_system(self, nfrag)
       integer(I4B),          intent(in)    :: nfrag !! Number of fragments to create
 
       if (allocated(self%fragments)) deallocate(self%fragments)
-      allocate(fraggle_fragments(nfrag) :: self%fragments)
-
+      allocate(fraggle_fragments(nbody=nfrag) :: self%fragments)
+      self%fragments%nbody = nfrag
       return
    end subroutine fraggle_setup_fragments_system