Bug in io. Meas last Eoft.

parent 69abc912
...@@ -145,7 +145,7 @@ function save_data() ...@@ -145,7 +145,7 @@ function save_data()
end end
BDIO_start_record!(fb, BDIO_BIN_GENERIC, 8) BDIO_start_record!(fb, BDIO_BIN_GENERIC, 8, true)
for noi in 1:params["Frontflow"]["N_noise"] for noi in 1:params["Frontflow"]["N_noise"]
BDIO_write!(fb,pp_corr_t0[:,noi]) BDIO_write!(fb,pp_corr_t0[:,noi])
......
...@@ -105,6 +105,11 @@ function Frontflow_pt() # Will be Frontflow ...@@ -105,6 +105,11 @@ function Frontflow_pt() # Will be Frontflow
flw(U, psi, int, 1, params["Frontflow"]["epsilon"], gp, dpar, lp, ymws, dws) flw(U, psi, int, 1, params["Frontflow"]["epsilon"], gp, dpar, lp, ymws, dws)
end end
if noi ==1
Eoft_plaq(Eofpla,U, gp, lp, ymws)
Eoft[end,:] .= sum(Eofpla,dims = 2)
end
pp_density .= Array(norm2.(psi)) pp_density .= Array(norm2.(psi))
ap_density .= Array(dot.(psi,dmul.(Gamma{4},psi))) ap_density .= Array(dot.(psi,dmul.(Gamma{4},psi)))
pp_corr_t[:,noi,end] .= zero(Float64) pp_corr_t[:,noi,end] .= zero(Float64)
......
...@@ -14,7 +14,7 @@ function read_ff(name::String) ...@@ -14,7 +14,7 @@ function read_ff(name::String)
ihdr = Vector{Int32}(undef,1) ihdr = Vector{Int32}(undef,1)
BDIO_read(file, ihdr) BDIO_read(file, ihdr)
if ihdr[1] == 1730280201 if ihdr[1] == 1730280201
print("Reading file ",name) println("Reading file ",name)
else else
error("Wrong IHDR in file ", name) error("Wrong IHDR in file ", name)
end end
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment