TL: debug (againg ...)

Author: tlunet

pySDC/playgrounds/dedalus/problems/rbc.py

@@ -256,6 +256,19 @@ def runSimulation(cls, runDir, tEnd, baseDt, tBeg=0, logEvery=100,
 
         p = cls(**pParams)
 
+        dt = baseDt/p.resFactor
+        nSteps = round(float(tEnd-tBeg)/dt, ndigits=3)
+        if float(tEnd-tBeg) != round(nSteps*dt, ndigits=3):
+            raise ValueError(f"{tEnd=} is not divisible by timestep {dt=} ({nSteps=})")
+        nSteps = int(nSteps)
+        p.infos.update(tEnd=tEnd, dt=dt, nSteps=nSteps)
+
+        if os.path.isfile(f"{runDir}/01_finalized.txt"):
+            cls.log(" -- simulation already finalized, skipping !")
+            return p
+        os.makedirs(runDir, exist_ok=True)
+        p.infos.update(dirName=runDir)
+
         if writeDecomposition:
             decompFile = f"{runDir}/decomp.txt"
             if MPI_RANK == 0:
@@ -281,19 +294,6 @@ def runSimulation(cls, runDir, tEnd, baseDt, tBeg=0, logEvery=100,
                     )
             COMM_WORLD.Barrier()
 
-        dt = baseDt/p.resFactor
-        nSteps = round(float(tEnd-tBeg)/dt, ndigits=3)
-        if float(tEnd-tBeg) != round(nSteps*dt, ndigits=3):
-            raise ValueError(f"{tEnd=} is not divisible by timestep {dt=} ({nSteps=})")
-        nSteps = int(nSteps)
-        p.infos.update(tEnd=tEnd, dt=dt, nSteps=nSteps)
-
-        if os.path.isfile(f"{runDir}/01_finalized.txt"):
-            cls.log(" -- simulation already finalized, skipping !")
-            return p
-        os.makedirs(runDir, exist_ok=True)
-        p.infos.update(dirName=runDir)
-
         # Solver
         cls.log(" -- building dedalus solver ...")
         solver = p.problem.build_solver(TimeStepper)