Obs at t0. Top charge for all t. Read not updated.

src/io.jl

@@ -171,7 +171,13 @@ function save_data()
         end
     end
 
-    BDIO_write!(fb,Eoft)
+
+    BDIO_write!(fb,Eoft0)
+    BDIO_write!(fb,Qt0)
+    for fl in 1:params["Frontflow"]["nsteps"]+1
+        BDIO_write!(fb,Eoft[:,fl])
+        BDIO_write!(fb,Qt[:,fl])
+    end
 
     for noi in 1:params["Backflow"]["N_noise"]
 

src/meas.jl

@@ -16,9 +16,19 @@ function load_fields()
     global pp_density = Array{Float64}(undef,lp.bsz,lp.rsz);
     global ap_density = Array{Complex{Float64}}(undef,lp.bsz,lp.rsz);
 
+    # Obs ending in 0 are only relevant if tzero > 0
+
     global pp_corr_t0 = fill(zero(Float64),(lp.iL[4],params["Frontflow"]["N_noise"]));
     global ap_corr_t0 = fill(zero(ComplexF64),(lp.iL[4],params["Frontflow"]["N_noise"]));
 
+    global pphat_t0 = Array{Float64}(undef,lp.iL[4],params["Frontflow"]["N_noise"]);
+    global pptilde_t0 = Array{Complex{Float64}}(undef,lp.iL[4],params["Frontflow"]["N_noise"]);
+
+    global Qt0 = Array{Complex{Float64}}(undef,lp.iL[4]);
+    global Eoft0 = Array{Complex{Float64}}(undef,lp.iL[4]);
+
+    # Frontflow
+
     global pp_corr_t = fill(zero(Float64),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1));
     global ap_corr_t = fill(zero(ComplexF64),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1));
 
@@ -28,9 +38,12 @@ function load_fields()
     global pphat_t = Array{Float64}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1);
     global pptilde_t = Array{Complex{Float64}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1);
 
-    global Eoft = Array{Complex{Float64}}(undef,2 + params["Frontflow"]["nsteps"],lp.iL[4]);
+    global Qt = Array{Complex{Float64}}(undef,params["Frontflow"]["nsteps"] + 1,lp.iL[4]);
+    global Eoft = Array{Complex{Float64}}(undef,params["Frontflow"]["nsteps"] + 1,lp.iL[4]);
     global Eofpla = Array{Complex{Float64}}(undef,lp.iL[4],lp.npls);
 
+    # Backflow
+
     global Quark_cond = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times));
     global Quark_cond_cfl = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times));
 
@@ -55,7 +68,9 @@ function Frontflow_pt() # Will be Frontflow
     flush(log_file)
 
     Eoft_plaq(Eofpla,U, gp, lp, ymws)
-    Eoft[1,:] .= sum(Eofpla,dims = 2)
+    Eoft0 .= sum(Eofpla,dims = 2)
+    Qtop(pp_density,U, gp, lp, ymws)
+    Qt0 .= pp_density
 
     for noi in 1:params["Frontflow"]["N_noise"]
 
@@ -84,7 +99,9 @@ function Frontflow_pt() # Will be Frontflow
 
             if noi ==1
                 Eoft_plaq(Eofpla,U, gp, lp, ymws)
-                Eoft[tstep+1,:] .= sum(Eofpla,dims = 2)
+                Eoft[tstep,:] .= sum(Eofpla,dims = 2)
+                Qtop(pp_density, U, gp, lp, ymws)
+                Qt[tstep,:] .= pp_density
             end
 
             pp_corr_t[:,noi,tstep] .= zero(Float64)
@@ -100,6 +117,8 @@ function Frontflow_pt() # Will be Frontflow
         if noi ==1
             Eoft_plaq(Eofpla,U, gp, lp, ymws)
             Eoft[end,:] .= sum(Eofpla,dims = 2)
+            Qtop(pp_density, U, gp, lp, ymws)
+            Qt[end,:] .= pp_density
         end
 
         pp_density .= Array(norm2.(psi))

src/utils.jl

@@ -3,7 +3,6 @@ using ADerrors
 using BDIO
 
 
-## NEDS UPDATE
 function read_ff(name::String)
 
     file = BDIO_open(name,"r","FerFlow correlators")