Added max cnfg run adn output tzero

input/FerFlow.in

@@ -8,13 +8,14 @@ size       = [8,8,8,8]
 blocks     = [4,4,4,4]
 
 [HMC]
-nth   = 100
+nth      = 10
 beta     = 6.0
 c0       = 1.0
 cG       = 1.0
 dmeas    = 5
 eps      = 0.01
 ns       = 30
+mclength = -1 # Number of measurements. Negative to ignore
 
 [Fermion]
 kappa   = 0.1348
@@ -30,10 +31,10 @@ N_noise = 6
 t_zero = 2.0
 epsilon = 0.01
 nsteps = 100
-tsource = 8
+tsource = 4
 
 [Backflow]
 N_noise = 6
-Flow_times = [2.47,3.86]
+Flow_times = [2.0,3.0]
 Nsaves = 25
-tsource = 8
+tsource = 4

main.jl

@@ -38,6 +38,8 @@ while RUN_ON
 
     @timeit "Saving" save_data()
 
+    check_run_status()
+
 end
 
 print_timer(log_file)

src/io.jl

@@ -99,6 +99,7 @@ function write_log()
     println(log_file,"tolerance          = ", params["Fermion"]["tolerance"])
     println(log_file,"maxiter            = ", params["Fermion"]["maxiter"])
     println(log_file,"N noise frontflow  = ", params["Frontflow"]["N_noise"])
+    println(log_file,"t_zero             = ", params["Frontflow"]["t_zero"])
     println(log_file,"eps frontflw       = ", params["Frontflow"]["epsilon"])
     println(log_file,"t_source frontflow = ", params["Frontflow"]["tsource"])
     println(log_file,"Flow times         = ", params["Backflow"]["Flow_times"])

src/mc.jl

@@ -25,6 +25,20 @@ function therm_or_readlist()
     return nothing
 end
 
+function check_run_status()
+
+    if parsed_args["c"] == nothing
+        global RUN_ON = isfile("RUN_ON_"*params["Run"]["name"])
+    else
+        global RUN_ON = (isfile("RUN_ON_"*params["Run"]["name"]) && CONFIG_LIST[MCid+1] != "./cnfg/")
+    end
+
+    if MCid == params["HMC"]["mclength"]
+        global RUN_ON = false
+    end
+
+    return nothing
+end
 
 function Gauge_update()
     if parsed_args["c"] == nothing