validation dieSwell

validation/dieSwell/dieSwell.py

@@ -55,6 +55,8 @@ class DieSwell(unittest.TestCase):
             return tension
 
         # Define output directory
+        dir_path = os.path.dirname(os.path.realpath(__file__))
+        os.chdir(dir_path)
         outputdir = "output"
         if not os.path.isdir(outputdir) :
             os.makedirs(outputdir)
@@ -125,4 +127,5 @@ class DieSwell(unittest.TestCase):
                                                     # The shape of low-speed capillary jets of Newtonian liquids. 
                                                     # Journal of Fluid Mechanics, 25(1), 185-198."
         error = np.abs(swelling - swelling_experiment)
+        # print(f"swelling = {swelling}\t error = {error}")
         self.assertLess(error, 2*10-2, "error is too large in dieSwell")
\ No newline at end of file