Added crater_soften routine for inside crater softening

MintonGroup · Mar 7, 2017 · be20dc9 · be20dc9
1 parent b0c96e2
commit be20dc9
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90
@@ -172,8 +172,11 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
          crater%maxinc = 0
          ! Do seismic shaking
          if (user%doseismic) call seismic_shake(user,surf,crater,domain)
+
+         ! Generate interior anomalous diffusion
+         if (user%dosoftening) call crater_soften(user,surf,crater,domain)
 
-         ! Generate dynamic diffusion
+         ! Generate distal anomalous diffusion
          if (user%dosoftening) call crater_soften_accumulate(user,surf,crater,domain,kdiff)
 
          ! find the average height and slope at crater location