fixed another bug (DIR in fb.solve)

scan_fixed_boundaries.py

@@ -42,8 +42,8 @@ parameters = {
 
 # parameters to scan
 lamda_list = np.arange(1.0, 61.0, 10.0)
-turnover_rho_list = np.arange(0, 1.0, 0.01)
-elasticity_list = np.arange(0, 1.0, 0.01)
+turnover_rho_list = np.arange(0, 1.0, 0.1)
+elasticity_list = np.arange(0, 1.0, 0.1)
 
 plist = list(itertools.product(lamda_list, turnover_rho_list, elasticity_list))
 
@@ -66,12 +66,13 @@ for i in range(len(plist)):
 
     try:
         fb = FixedBoundaries(params)
-        fb.solve(extend=False)
         DIR = args.outputdir
+
+        fb.solve(extend=False, DIR=DIR)
         with open(os.path.join(DIR, params['timestamp']+'_fixed_boundaries.json'), "w") as fp:
             json.dump(params, fp, indent=4)
             
-        visualize(params, DIR, interactive=False)
+        visualize(params, DIR=DIR, interactive=False)
     except:
         print('*** FAILED ****')
         os.system("rm -f %s*" % timestamp)