Added option to not measure flow observables

main/ym.jl

@@ -168,10 +168,11 @@ for i in 1:nthm
     end
 end
 
-
-Ecl = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
-Epl = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
-Qt  = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
+if flw_dtr != 0
+    Ecl = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
+    Epl = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
+    Qt  = Array{T,2}(undef, flw_ns+1, div(nt,flw_dtr))
+end
 println(flog, "\n # Measurements")
 kfl = 0
 for i in 1:nt
@@ -189,31 +190,39 @@ for i in 1:nt
                 sfc[k,1], sfc[k,2])
     end
 
-    ymws.U1 .= U
-    if i%flw_dtr == 0
-        global kfl = kfl + 1
-        for j in 1:flw_ns
-            Epl[j,kfl] = Eoft_plaq(U, gp, lp, ymws)
-            Ecl[j,kfl] = Eoft_clover(U, gp, lp, ymws)
-            Qt[j,kfl]  = Qtop(U, gp, lp, ymws)
-
-            flw(U, flwint, 1, gp, lp, ymws)
-            if (j-1)%div(flw_ns,10) == 0
-                t = j*flwint.eps
-                @printf(flog, "     FLOW   t=%8.4f:  %20.12e  %20.12e      %10.4e \n",
-                        t, t^2*Epl[j,kfl], t^2*Ecl[j,kfl], Qt[j,kfl])
+    if flw_dtr != 0
+        ymws.U1 .= U
+        if i%flw_dtr == 0
+            global kfl = kfl + 1
+            for j in 1:flw_ns
+                Epl[j,kfl] = Eoft_plaq(U, gp, lp, ymws)
+                Ecl[j,kfl] = Eoft_clover(U, gp, lp, ymws)
+                Qt[j,kfl]  = Qtop(U, gp, lp, ymws)
+                
+                flw(U, flwint, 1, gp, lp, ymws)
+                if (j-1)%div(flw_ns,10) == 0
+                    t = j*flwint.eps
+                    @printf(flog, "     FLOW   t=%8.4f:  %20.12e  %20.12e      %10.4e \n",
+                            t, t^2*Epl[j,kfl], t^2*Ecl[j,kfl], Qt[j,kfl])
+                end
             end
+            Epl[end,kfl] = Eoft_plaq(U, gp, lp, ymws)
+            Ecl[end,kfl] = Eoft_clover(U, gp, lp, ymws)
+            Qt[end,kfl]  = Qtop(U, gp, lp, ymws)
         end
-        Epl[end,kfl] = Eoft_plaq(U, gp, lp, ymws)
-        Ecl[end,kfl] = Eoft_clover(U, gp, lp, ymws)
-        Qt[end,kfl]  = Qtop(U, gp, lp, ymws)
-    end
-    U .= ymws.U1
+        U .= ymws.U1
 
-    if i%10 == 0
-        write(flog, TranscodingStreams.TOKEN_END)
-        flush(flog)
+        if i%10 == 0
+            write(flog, TranscodingStreams.TOKEN_END)
+            flush(flog)
+        end
+    else
+        if i%100 == 0
+            write(flog, TranscodingStreams.TOKEN_END)
+            flush(flog)
+        end
     end
+        
 end
 
 print_timer(flog)
@@ -266,10 +275,12 @@ if BCS == BC_SF_AFWB
     write_record(fb, 8, sfc[:,1])
     write_record(fb, 10, sfc[:,2])
 end
-write_flow_record(fb, 12, Epl)
-write_flow_record(fb, 13, Ecl)
-write_flow_record(fb, 14, Qt)
 
+if flw_dtr != 0
+    write_flow_record(fb, 12, Epl)
+    write_flow_record(fb, 13, Ecl)
+    write_flow_record(fb, 14, Qt)
+end
 
 BDIO_close!(fb)