John Wise avatar John Wise committed b86d8ba

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

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);
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.