Commits

John Wise committed b86d8ba

Missing the scale factor when calculating the timestep, using safety
velocities, for the regular and RT timestep.

  • Participants
  • Parent commits b13093a
  • Branches euler_improve
  • Tags Wise10_GalaxyBirth

Comments (0)

Files changed (2)

src/enzo/Grid_ComputePhotonTimestep.C

 
   float dtSafetyVelocity = huge_number;
   if (TimestepSafetyVelocity > 0)
-    dtSafetyVelocity = CellWidth[0][0] / 
+    dtSafetyVelocity = a*CellWidth[0][0] / 
       (TimestepSafetyVelocity*1e5 / VelocityUnits);    // parameter in km/s
 
   dt = min(dt, dtSafetyVelocity);
   float dtPhotonSafety = tiny_number;
 
   if (RadiativeTransferTimestepVelocityLimit > 0)
-    dtPhotonSafety = CellWidth[0][0] / 
+    dtPhotonSafety = a*CellWidth[0][0] / 
       (RadiativeTransferTimestepVelocityLimit*1e5 / VelocityUnits);
   dt = max(dt, dtPhotonSafety);
 

src/enzo/Grid_ComputeTimeStep.C

 
   float dtSafetyVelocity = huge_number;
   if (TimestepSafetyVelocity > 0)
-    dtSafetyVelocity = CellWidth[0][0] / 
+    dtSafetyVelocity = a*CellWidth[0][0] / 
       (TimestepSafetyVelocity*1e5 / VelocityUnits);    // parameter in km/s
 
   dt = min(dt, dtSafetyVelocity);