diff --git a/outputdipoleamplitude100vertical b/outputdipoleamplitude100vertical
new file mode 100644
index 0000000000000000000000000000000000000000..fe03adf40d8791d1cb773553c7f1f169177a53fd
--- /dev/null
+++ b/outputdipoleamplitude100vertical
@@ -0,0 +1,4160 @@
+#Current working directory is /n/work01/pyjorunk/JIMWLK/wilsonline_reader
+#Running on compute-1-35
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-0
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00043926
+0.03000000, 0.00217743
+0.04000000, 0.00217743
+0.05000000, 0.00217743
+0.06000000, 0.00221980
+0.07000000, 0.00221980
+0.08000000, 0.00320051
+0.09000000, 0.00320051
+
+0.10000000, 0.00320051
+0.20000000, 0.01184412
+0.30000000, 0.03873750
+0.40000000, 0.05839374
+0.50000000, 0.08598646
+0.60000000, 0.11529871
+0.70000000, 0.10355687
+0.80000000, 0.08995341
+0.90000000, 0.08255312
+
+1.00000000, 0.09409725
+2.00000000, 0.26095787
+3.00000000, 0.35175660
+4.00000000, 0.21235294
+5.00000000, 0.08404967
+6.00000000, 0.02914264
+7.00000000, 0.00908423
+8.00000000, 0.00324735
+9.00000000, 0.00116483
+
+10.00000000, 0.00041937
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-1
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00033807
+0.03000000, 0.00344073
+0.04000000, 0.00344073
+0.05000000, 0.00344073
+0.06000000, 0.00671055
+0.07000000, 0.00671055
+0.08000000, 0.00993143
+0.09000000, 0.00993143
+
+0.10000000, 0.00993143
+0.20000000, 0.03714395
+0.30000000, 0.06266321
+0.40000000, 0.07256698
+0.50000000, 0.09750225
+0.60000000, 0.12510604
+0.70000000, 0.16297871
+0.80000000, 0.18715515
+0.90000000, 0.20362877
+
+1.00000000, 0.20376835
+2.00000000, 0.30115586
+3.00000000, 0.16978037
+4.00000000, 0.10506653
+5.00000000, 0.04681729
+6.00000000, 0.02117845
+7.00000000, 0.00722118
+8.00000000, 0.00273217
+9.00000000, 0.00101788
+
+10.00000000, 0.00037528
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-10
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00029225
+0.03000000, 0.00120184
+0.04000000, 0.00120184
+0.05000000, 0.00120184
+0.06000000, 0.00213095
+0.07000000, 0.00213095
+0.08000000, 0.00532670
+0.09000000, 0.00532670
+
+0.10000000, 0.00532670
+0.20000000, 0.02538620
+0.30000000, 0.03875759
+0.40000000, 0.05062893
+0.50000000, 0.04617985
+0.60000000, 0.05683408
+0.70000000, 0.07126105
+0.80000000, 0.08912326
+0.90000000, 0.10785096
+
+1.00000000, 0.15729950
+2.00000000, 0.53166511
+3.00000000, 0.29586967
+4.00000000, 0.13018850
+5.00000000, 0.04194406
+6.00000000, 0.01322771
+7.00000000, 0.00480096
+8.00000000, 0.00180026
+9.00000000, 0.00066339
+
+10.00000000, 0.00024247
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-11
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00046189
+0.03000000, 0.00362302
+0.04000000, 0.00362302
+0.05000000, 0.00362302
+0.06000000, 0.00873407
+0.07000000, 0.00873407
+0.08000000, 0.01145051
+0.09000000, 0.01145051
+
+0.10000000, 0.01145051
+0.20000000, 0.03672023
+0.30000000, 0.04288684
+0.40000000, 0.06614074
+0.50000000, 0.10447356
+0.60000000, 0.13029985
+0.70000000, 0.14387247
+0.80000000, 0.17369493
+0.90000000, 0.21925646
+
+1.00000000, 0.26122133
+2.00000000, 0.37932503
+3.00000000, 0.14416460
+4.00000000, 0.06090030
+5.00000000, 0.02208654
+6.00000000, 0.00800221
+7.00000000, 0.00290695
+8.00000000, 0.00111083
+9.00000000, 0.00041674
+
+10.00000000, 0.00015433
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-12
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00030043
+0.03000000, 0.00181720
+0.04000000, 0.00181720
+0.05000000, 0.00181720
+0.06000000, 0.00403361
+0.07000000, 0.00403361
+0.08000000, 0.00773041
+0.09000000, 0.00773041
+
+0.10000000, 0.00773041
+0.20000000, 0.02986811
+0.30000000, 0.06335811
+0.40000000, 0.08597103
+0.50000000, 0.11182885
+0.60000000, 0.15907236
+0.70000000, 0.21515775
+0.80000000, 0.22189262
+0.90000000, 0.23329635
+
+1.00000000, 0.24664267
+2.00000000, 0.25783247
+3.00000000, 0.12912279
+4.00000000, 0.08644545
+5.00000000, 0.02343666
+6.00000000, 0.00656294
+7.00000000, 0.00243518
+8.00000000, 0.00098082
+9.00000000, 0.00037645
+
+10.00000000, 0.00014083
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-13
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00029828
+0.03000000, 0.00166811
+0.04000000, 0.00166811
+0.05000000, 0.00166811
+0.06000000, 0.00553725
+0.07000000, 0.00553725
+0.08000000, 0.00830215
+0.09000000, 0.00830215
+
+0.10000000, 0.00830215
+0.20000000, 0.02112589
+0.30000000, 0.04576964
+0.40000000, 0.09466658
+0.50000000, 0.15091200
+0.60000000, 0.14798388
+0.70000000, 0.15874835
+0.80000000, 0.16104144
+0.90000000, 0.19304413
+
+1.00000000, 0.20240668
+2.00000000, 0.38077359
+3.00000000, 0.28316232
+4.00000000, 0.11308796
+5.00000000, 0.02861761
+6.00000000, 0.01332912
+7.00000000, 0.00486978
+8.00000000, 0.00181446
+9.00000000, 0.00066849
+
+10.00000000, 0.00024485
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-14
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00059506
+0.03000000, 0.00319132
+0.04000000, 0.00319132
+0.05000000, 0.00319132
+0.06000000, 0.00845223
+0.07000000, 0.00845223
+0.08000000, 0.01396711
+0.09000000, 0.01396711
+
+0.10000000, 0.01396711
+0.20000000, 0.05766103
+0.30000000, 0.11008018
+0.40000000, 0.13761962
+0.50000000, 0.15303367
+0.60000000, 0.18614253
+0.70000000, 0.22605476
+0.80000000, 0.27519521
+0.90000000, 0.32294200
+
+1.00000000, 0.33830731
+2.00000000, 0.58469222
+3.00000000, 0.30836079
+4.00000000, 0.14107800
+5.00000000, 0.06113022
+6.00000000, 0.02166142
+7.00000000, 0.00820999
+8.00000000, 0.00319499
+9.00000000, 0.00120199
+
+10.00000000, 0.00044459
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-15
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00018760
+0.03000000, 0.00163835
+0.04000000, 0.00163835
+0.05000000, 0.00163835
+0.06000000, 0.00448933
+0.07000000, 0.00448933
+0.08000000, 0.00841534
+0.09000000, 0.00841534
+
+0.10000000, 0.00841534
+0.20000000, 0.02017047
+0.30000000, 0.02234207
+0.40000000, 0.02978363
+0.50000000, 0.04277830
+0.60000000, 0.09799130
+0.70000000, 0.13317608
+0.80000000, 0.18231543
+0.90000000, 0.20442526
+
+1.00000000, 0.23481455
+2.00000000, 0.59559893
+3.00000000, 0.28870708
+4.00000000, 0.11068279
+5.00000000, 0.05061721
+6.00000000, 0.02201769
+7.00000000, 0.00849678
+8.00000000, 0.00329769
+9.00000000, 0.00123971
+
+10.00000000, 0.00045890
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-16
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00031845
+0.03000000, 0.00143662
+0.04000000, 0.00143662
+0.05000000, 0.00143662
+0.06000000, 0.00363707
+0.07000000, 0.00363707
+0.08000000, 0.00629869
+0.09000000, 0.00629869
+
+0.10000000, 0.00629869
+0.20000000, 0.02531700
+0.30000000, 0.06029913
+0.40000000, 0.08678076
+0.50000000, 0.14156660
+0.60000000, 0.16714496
+0.70000000, 0.20614557
+0.80000000, 0.24839519
+0.90000000, 0.31515689
+
+1.00000000, 0.32095358
+2.00000000, 0.21091823
+3.00000000, 0.14860206
+4.00000000, 0.08769605
+5.00000000, 0.04368018
+6.00000000, 0.01672331
+7.00000000, 0.00563490
+8.00000000, 0.00209499
+9.00000000, 0.00077046
+
+10.00000000, 0.00028164
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-17
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00031086
+0.03000000, 0.00301814
+0.04000000, 0.00301814
+0.05000000, 0.00301814
+0.06000000, 0.00726105
+0.07000000, 0.00726105
+0.08000000, 0.00847822
+0.09000000, 0.00847822
+
+0.10000000, 0.00847822
+0.20000000, 0.02981365
+0.30000000, 0.05532029
+0.40000000, 0.07051891
+0.50000000, 0.08910289
+0.60000000, 0.14739921
+0.70000000, 0.19738940
+0.80000000, 0.22150046
+0.90000000, 0.25015089
+
+1.00000000, 0.25525874
+2.00000000, 0.27536812
+3.00000000, 0.12778513
+4.00000000, 0.04429886
+5.00000000, 0.01785186
+6.00000000, 0.00908279
+7.00000000, 0.00350085
+8.00000000, 0.00136070
+9.00000000, 0.00051415
+
+10.00000000, 0.00019128
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-18
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00065590
+0.03000000, 0.00399660
+0.04000000, 0.00399660
+0.05000000, 0.00399660
+0.06000000, 0.00663099
+0.07000000, 0.00663099
+0.08000000, 0.00917842
+0.09000000, 0.00917842
+
+0.10000000, 0.00917842
+0.20000000, 0.02763802
+0.30000000, 0.05883117
+0.40000000, 0.08491567
+0.50000000, 0.09561065
+0.60000000, 0.11799672
+0.70000000, 0.14166682
+0.80000000, 0.16496758
+0.90000000, 0.21487570
+
+1.00000000, 0.28320634
+2.00000000, 0.42491454
+3.00000000, 0.11724093
+4.00000000, 0.06969390
+5.00000000, 0.03453736
+6.00000000, 0.01262016
+7.00000000, 0.00412751
+8.00000000, 0.00143595
+9.00000000, 0.00050053
+
+10.00000000, 0.00017621
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-19
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00034271
+0.03000000, 0.00172605
+0.04000000, 0.00172605
+0.05000000, 0.00172605
+0.06000000, 0.00311224
+0.07000000, 0.00311224
+0.08000000, 0.00387484
+0.09000000, 0.00387484
+
+0.10000000, 0.00387484
+0.20000000, 0.01527554
+0.30000000, 0.05414006
+0.40000000, 0.13424887
+0.50000000, 0.19119214
+0.60000000, 0.22367071
+0.70000000, 0.24719421
+0.80000000, 0.27189150
+0.90000000, 0.31941271
+
+1.00000000, 0.37691771
+2.00000000, 0.33603234
+3.00000000, 0.26204175
+4.00000000, 0.18972739
+5.00000000, 0.05714852
+6.00000000, 0.01697103
+7.00000000, 0.00570146
+8.00000000, 0.00213360
+9.00000000, 0.00078712
+
+10.00000000, 0.00028813
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-2
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00032317
+0.03000000, 0.00124830
+0.04000000, 0.00124830
+0.05000000, 0.00124830
+0.06000000, 0.00296341
+0.07000000, 0.00296341
+0.08000000, 0.00475114
+0.09000000, 0.00475114
+
+0.10000000, 0.00475114
+0.20000000, 0.01808923
+0.30000000, 0.05328454
+0.40000000, 0.07969171
+0.50000000, 0.11976971
+0.60000000, 0.17373354
+0.70000000, 0.22347779
+0.80000000, 0.25914454
+0.90000000, 0.28748137
+
+1.00000000, 0.35330857
+2.00000000, 0.43721085
+3.00000000, 0.29098210
+4.00000000, 0.12141782
+5.00000000, 0.03434624
+6.00000000, 0.01079785
+7.00000000, 0.00395842
+8.00000000, 0.00150308
+9.00000000, 0.00055612
+
+10.00000000, 0.00020358
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-20
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00049797
+0.03000000, 0.00273279
+0.04000000, 0.00273279
+0.05000000, 0.00273279
+0.06000000, 0.00854437
+0.07000000, 0.00854437
+0.08000000, 0.01470828
+0.09000000, 0.01470828
+
+0.10000000, 0.01470828
+0.20000000, 0.04854286
+0.30000000, 0.10159423
+0.40000000, 0.15521276
+0.50000000, 0.20944762
+0.60000000, 0.24118208
+0.70000000, 0.33918452
+0.80000000, 0.35586151
+0.90000000, 0.35877702
+
+1.00000000, 0.34526569
+2.00000000, 0.27970155
+3.00000000, 0.24510382
+4.00000000, 0.10796002
+5.00000000, 0.03504857
+6.00000000, 0.01242780
+7.00000000, 0.00401126
+8.00000000, 0.00146686
+9.00000000, 0.00053487
+
+10.00000000, 0.00019462
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-21
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00079516
+0.03000000, 0.00323207
+0.04000000, 0.00323207
+0.05000000, 0.00323207
+0.06000000, 0.00686287
+0.07000000, 0.00686287
+0.08000000, 0.01036327
+0.09000000, 0.01036327
+
+0.10000000, 0.01036327
+0.20000000, 0.01830218
+0.30000000, 0.03597655
+0.40000000, 0.05353723
+0.50000000, 0.06916465
+0.60000000, 0.08622121
+0.70000000, 0.12137315
+0.80000000, 0.13265699
+0.90000000, 0.10931610
+
+1.00000000, 0.12394322
+2.00000000, 0.14177924
+3.00000000, 0.15581494
+4.00000000, 0.07789386
+5.00000000, 0.03376997
+6.00000000, 0.01013385
+7.00000000, 0.00343713
+8.00000000, 0.00126362
+9.00000000, 0.00045787
+
+10.00000000, 0.00016544
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-22
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00061988
+0.03000000, 0.00317439
+0.04000000, 0.00317439
+0.05000000, 0.00317439
+0.06000000, 0.00858469
+0.07000000, 0.00858469
+0.08000000, 0.01686996
+0.09000000, 0.01686996
+
+0.10000000, 0.01686996
+0.20000000, 0.05642422
+0.30000000, 0.09851344
+0.40000000, 0.12986093
+0.50000000, 0.14517065
+0.60000000, 0.19400850
+0.70000000, 0.24408092
+0.80000000, 0.27641299
+0.90000000, 0.34005516
+
+1.00000000, 0.37027430
+2.00000000, 0.32216069
+3.00000000, 0.23014979
+4.00000000, 0.09805837
+5.00000000, 0.01651247
+6.00000000, 0.00465509
+7.00000000, 0.00187506
+8.00000000, 0.00074267
+9.00000000, 0.00028320
+
+10.00000000, 0.00010593
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-23
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00043731
+0.03000000, 0.00275503
+0.04000000, 0.00275503
+0.05000000, 0.00275503
+0.06000000, 0.00540513
+0.07000000, 0.00540513
+0.08000000, 0.00815699
+0.09000000, 0.00815699
+
+0.10000000, 0.00815699
+0.20000000, 0.02895444
+0.30000000, 0.05489269
+0.40000000, 0.08123691
+0.50000000, 0.08191445
+0.60000000, 0.06319660
+0.70000000, 0.04821727
+0.80000000, 0.05017167
+0.90000000, 0.07086450
+
+1.00000000, 0.08561407
+2.00000000, 0.16874793
+3.00000000, 0.07665472
+4.00000000, 0.03267769
+5.00000000, 0.01809785
+6.00000000, 0.01023892
+7.00000000, 0.00346082
+8.00000000, 0.00124805
+9.00000000, 0.00045139
+
+10.00000000, 0.00016367
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-24
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00063710
+0.03000000, 0.00491023
+0.04000000, 0.00491023
+0.05000000, 0.00491023
+0.06000000, 0.01078955
+0.07000000, 0.01078955
+0.08000000, 0.01944608
+0.09000000, 0.01944608
+
+0.10000000, 0.01944608
+0.20000000, 0.06779088
+0.30000000, 0.11420135
+0.40000000, 0.19321195
+0.50000000, 0.20976422
+0.60000000, 0.19877727
+0.70000000, 0.18760391
+0.80000000, 0.19041373
+0.90000000, 0.19045918
+
+1.00000000, 0.21505017
+2.00000000, 0.32294855
+3.00000000, 0.17608625
+4.00000000, 0.06509687
+5.00000000, 0.02159254
+6.00000000, 0.00712226
+7.00000000, 0.00294759
+8.00000000, 0.00121425
+9.00000000, 0.00047460
+
+10.00000000, 0.00017983
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-25
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00024600
+0.03000000, 0.00173959
+0.04000000, 0.00173959
+0.05000000, 0.00173959
+0.06000000, 0.00335290
+0.07000000, 0.00335290
+0.08000000, 0.00661821
+0.09000000, 0.00661821
+
+0.10000000, 0.00661821
+0.20000000, 0.03073681
+0.30000000, 0.03281048
+0.40000000, 0.03700376
+0.50000000, 0.04097871
+0.60000000, 0.07346036
+0.70000000, 0.10109708
+0.80000000, 0.11351340
+0.90000000, 0.14084683
+
+1.00000000, 0.15630899
+2.00000000, 0.26443007
+3.00000000, 0.14715278
+4.00000000, 0.03657385
+5.00000000, 0.01535761
+6.00000000, 0.00750220
+7.00000000, 0.00244870
+8.00000000, 0.00089696
+9.00000000, 0.00032938
+
+10.00000000, 0.00012072
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-26
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00016816
+0.03000000, 0.00134635
+0.04000000, 0.00134635
+0.05000000, 0.00134635
+0.06000000, 0.00296808
+0.07000000, 0.00296808
+0.08000000, 0.00530493
+0.09000000, 0.00530493
+
+0.10000000, 0.00530493
+0.20000000, 0.03128993
+0.30000000, 0.07115999
+0.40000000, 0.09049646
+0.50000000, 0.12402681
+0.60000000, 0.17007501
+0.70000000, 0.21178804
+0.80000000, 0.22021509
+0.90000000, 0.20430453
+
+1.00000000, 0.18695665
+2.00000000, 0.14824179
+3.00000000, 0.13954467
+4.00000000, 0.08342412
+5.00000000, 0.04184282
+6.00000000, 0.01330178
+7.00000000, 0.00496134
+8.00000000, 0.00186177
+9.00000000, 0.00068536
+
+10.00000000, 0.00025035
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-27
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00107917
+0.03000000, 0.00716467
+0.04000000, 0.00716467
+0.05000000, 0.00716467
+0.06000000, 0.01643345
+0.07000000, 0.01643345
+0.08000000, 0.03074089
+0.09000000, 0.03074089
+
+0.10000000, 0.03074089
+0.20000000, 0.09423046
+0.30000000, 0.14197301
+0.40000000, 0.20158376
+0.50000000, 0.27433612
+0.60000000, 0.36017405
+0.70000000, 0.37873337
+0.80000000, 0.39583650
+0.90000000, 0.40586678
+
+1.00000000, 0.42858685
+2.00000000, 0.29577997
+3.00000000, 0.12179075
+4.00000000, 0.04027246
+5.00000000, 0.01695624
+6.00000000, 0.00662055
+7.00000000, 0.00236651
+8.00000000, 0.00088001
+9.00000000, 0.00032144
+
+10.00000000, 0.00011686
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-28
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00039123
+0.03000000, 0.00319487
+0.04000000, 0.00319487
+0.05000000, 0.00319487
+0.06000000, 0.00711650
+0.07000000, 0.00711650
+0.08000000, 0.01142627
+0.09000000, 0.01142627
+
+0.10000000, 0.01142627
+0.20000000, 0.03758632
+0.30000000, 0.06998579
+0.40000000, 0.11323080
+0.50000000, 0.13091232
+0.60000000, 0.12959496
+0.70000000, 0.12570681
+0.80000000, 0.12780756
+0.90000000, 0.17419107
+
+1.00000000, 0.19679461
+2.00000000, 0.44627598
+3.00000000, 0.20228525
+4.00000000, 0.08611092
+5.00000000, 0.04491541
+6.00000000, 0.01796008
+7.00000000, 0.00640636
+8.00000000, 0.00247831
+9.00000000, 0.00093758
+
+10.00000000, 0.00034931
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-29
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00018794
+0.03000000, 0.00193884
+0.04000000, 0.00193884
+0.05000000, 0.00193884
+0.06000000, 0.00464121
+0.07000000, 0.00464121
+0.08000000, 0.00658761
+0.09000000, 0.00658761
+
+0.10000000, 0.00658761
+0.20000000, 0.01704651
+0.30000000, 0.03122461
+0.40000000, 0.03171355
+0.50000000, 0.05054940
+0.60000000, 0.08027455
+0.70000000, 0.10329362
+0.80000000, 0.10263105
+0.90000000, 0.08428555
+
+1.00000000, 0.08551459
+2.00000000, 0.22438423
+3.00000000, 0.15231016
+4.00000000, 0.09030431
+5.00000000, 0.06143154
+6.00000000, 0.02509450
+7.00000000, 0.00809597
+8.00000000, 0.00288496
+9.00000000, 0.00102590
+
+10.00000000, 0.00036639
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-3
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00103570
+0.03000000, 0.00655955
+0.04000000, 0.00655955
+0.05000000, 0.00655955
+0.06000000, 0.01057478
+0.07000000, 0.01057478
+0.08000000, 0.01775879
+0.09000000, 0.01775879
+
+0.10000000, 0.01775879
+0.20000000, 0.04338533
+0.30000000, 0.08086703
+0.40000000, 0.12731430
+0.50000000, 0.17043241
+0.60000000, 0.23940732
+0.70000000, 0.27593774
+0.80000000, 0.32675214
+0.90000000, 0.39168485
+
+1.00000000, 0.37703666
+2.00000000, 0.25946285
+3.00000000, 0.13256494
+4.00000000, 0.04754561
+5.00000000, 0.02122370
+6.00000000, 0.00922670
+7.00000000, 0.00295678
+8.00000000, 0.00107646
+9.00000000, 0.00039306
+
+10.00000000, 0.00014349
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-30
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00020610
+0.03000000, 0.00134360
+0.04000000, 0.00134360
+0.05000000, 0.00134360
+0.06000000, 0.00341127
+0.07000000, 0.00341127
+0.08000000, 0.00569046
+0.09000000, 0.00569046
+
+0.10000000, 0.00569046
+0.20000000, 0.02516216
+0.30000000, 0.06080448
+0.40000000, 0.10207029
+0.50000000, 0.15499532
+0.60000000, 0.18236498
+0.70000000, 0.20535603
+0.80000000, 0.22393385
+0.90000000, 0.26164040
+
+1.00000000, 0.25860409
+2.00000000, 0.20391323
+3.00000000, 0.15302697
+4.00000000, 0.04763085
+5.00000000, 0.01509246
+6.00000000, 0.00442966
+7.00000000, 0.00188464
+8.00000000, 0.00076766
+9.00000000, 0.00029828
+
+10.00000000, 0.00011270
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-31
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00099298
+0.03000000, 0.00549007
+0.04000000, 0.00549007
+0.05000000, 0.00549007
+0.06000000, 0.01147596
+0.07000000, 0.01147596
+0.08000000, 0.01690882
+0.09000000, 0.01690882
+
+0.10000000, 0.01690882
+0.20000000, 0.04207629
+0.30000000, 0.07702428
+0.40000000, 0.10048250
+0.50000000, 0.11450438
+0.60000000, 0.15659456
+0.70000000, 0.19704647
+0.80000000, 0.20730042
+0.90000000, 0.21049831
+
+1.00000000, 0.22286746
+2.00000000, 0.11573390
+3.00000000, 0.02825458
+4.00000000, 0.03168374
+5.00000000, 0.00899328
+6.00000000, 0.00377181
+7.00000000, 0.00158611
+8.00000000, 0.00065966
+9.00000000, 0.00026178
+
+10.00000000, 0.00010056
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-32
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00042315
+0.03000000, 0.00116252
+0.04000000, 0.00116252
+0.05000000, 0.00116252
+0.06000000, 0.00272319
+0.07000000, 0.00272319
+0.08000000, 0.00684105
+0.09000000, 0.00684105
+
+0.10000000, 0.00684105
+0.20000000, 0.03329983
+0.30000000, 0.06502656
+0.40000000, 0.09285002
+0.50000000, 0.14682511
+0.60000000, 0.17946243
+0.70000000, 0.23511005
+0.80000000, 0.27231916
+0.90000000, 0.28129382
+
+1.00000000, 0.30967247
+2.00000000, 0.16723577
+3.00000000, 0.15140111
+4.00000000, 0.07993870
+5.00000000, 0.03537220
+6.00000000, 0.01233737
+7.00000000, 0.00404247
+8.00000000, 0.00143720
+9.00000000, 0.00051193
+
+10.00000000, 0.00018330
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-33
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00020206
+0.03000000, 0.00070291
+0.04000000, 0.00070291
+0.05000000, 0.00070291
+0.06000000, 0.00165314
+0.07000000, 0.00165314
+0.08000000, 0.00271527
+0.09000000, 0.00271527
+
+0.10000000, 0.00271527
+0.20000000, 0.00862383
+0.30000000, 0.03650707
+0.40000000, 0.08035929
+0.50000000, 0.11787654
+0.60000000, 0.16671237
+0.70000000, 0.21992536
+0.80000000, 0.28475309
+0.90000000, 0.34014986
+
+1.00000000, 0.39903599
+2.00000000, 0.37503813
+3.00000000, 0.24343353
+4.00000000, 0.15370444
+5.00000000, 0.06245299
+6.00000000, 0.02646702
+7.00000000, 0.00954564
+8.00000000, 0.00354636
+9.00000000, 0.00129231
+
+10.00000000, 0.00046838
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-34
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00031365
+0.03000000, 0.00191386
+0.04000000, 0.00191386
+0.05000000, 0.00191386
+0.06000000, 0.00432382
+0.07000000, 0.00432382
+0.08000000, 0.00595384
+0.09000000, 0.00595384
+
+0.10000000, 0.00595384
+0.20000000, 0.02214034
+0.30000000, 0.03892977
+0.40000000, 0.04464123
+0.50000000, 0.06132989
+0.60000000, 0.08931146
+0.70000000, 0.12795750
+0.80000000, 0.15577216
+0.90000000, 0.19238795
+
+1.00000000, 0.19878634
+2.00000000, 0.26621419
+3.00000000, 0.28359166
+4.00000000, 0.09888727
+5.00000000, 0.03683285
+6.00000000, 0.01181178
+7.00000000, 0.00401492
+8.00000000, 0.00146692
+9.00000000, 0.00053234
+
+10.00000000, 0.00019308
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-35
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00016601
+0.03000000, 0.00124531
+0.04000000, 0.00124531
+0.05000000, 0.00124531
+0.06000000, 0.00312269
+0.07000000, 0.00312269
+0.08000000, 0.00531322
+0.09000000, 0.00531322
+
+0.10000000, 0.00531322
+0.20000000, 0.01302189
+0.30000000, 0.02056641
+0.40000000, 0.03019819
+0.50000000, 0.05580047
+0.60000000, 0.07536566
+0.70000000, 0.09111148
+0.80000000, 0.11333969
+0.90000000, 0.11485270
+
+1.00000000, 0.11787794
+2.00000000, 0.17647364
+3.00000000, 0.08772164
+4.00000000, 0.08715773
+5.00000000, 0.03366474
+6.00000000, 0.01497921
+7.00000000, 0.00450085
+8.00000000, 0.00162818
+9.00000000, 0.00058982
+
+10.00000000, 0.00021365
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-36
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00046802
+0.03000000, 0.00423751
+0.04000000, 0.00423751
+0.05000000, 0.00423751
+0.06000000, 0.01138904
+0.07000000, 0.01138904
+0.08000000, 0.01719317
+0.09000000, 0.01719317
+
+0.10000000, 0.01719317
+0.20000000, 0.05534026
+0.30000000, 0.09805150
+0.40000000, 0.15827612
+0.50000000, 0.14573494
+0.60000000, 0.13509320
+0.70000000, 0.13821180
+0.80000000, 0.15562917
+0.90000000, 0.15563010
+
+1.00000000, 0.17144090
+2.00000000, 0.31575592
+3.00000000, 0.11246610
+4.00000000, 0.06205141
+5.00000000, 0.03235302
+6.00000000, 0.01249619
+7.00000000, 0.00420119
+8.00000000, 0.00151029
+9.00000000, 0.00054193
+
+10.00000000, 0.00019498
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-37
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00025195
+0.03000000, 0.00155682
+0.04000000, 0.00155682
+0.05000000, 0.00155682
+0.06000000, 0.00307666
+0.07000000, 0.00307666
+0.08000000, 0.00451793
+0.09000000, 0.00451793
+
+0.10000000, 0.00451793
+0.20000000, 0.02206913
+0.30000000, 0.03925907
+0.40000000, 0.05291215
+0.50000000, 0.05507548
+0.60000000, 0.04710550
+0.70000000, 0.06988884
+0.80000000, 0.09237954
+0.90000000, 0.11717600
+
+1.00000000, 0.15597036
+2.00000000, 0.24502098
+3.00000000, 0.10746780
+4.00000000, 0.05657258
+5.00000000, 0.03910675
+6.00000000, 0.01299273
+7.00000000, 0.00388168
+8.00000000, 0.00138034
+9.00000000, 0.00049412
+
+10.00000000, 0.00017763
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-38
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00058426
+0.03000000, 0.00303040
+0.04000000, 0.00303040
+0.05000000, 0.00303040
+0.06000000, 0.00761468
+0.07000000, 0.00761468
+0.08000000, 0.01255985
+0.09000000, 0.01255985
+
+0.10000000, 0.01255985
+0.20000000, 0.03448351
+0.30000000, 0.06597632
+0.40000000, 0.07710761
+0.50000000, 0.09530219
+0.60000000, 0.10417510
+0.70000000, 0.10998728
+0.80000000, 0.14556079
+0.90000000, 0.16434258
+
+1.00000000, 0.16622568
+2.00000000, 0.15320601
+3.00000000, 0.05012356
+4.00000000, 0.05708221
+5.00000000, 0.03325750
+6.00000000, 0.00956137
+7.00000000, 0.00309133
+8.00000000, 0.00111023
+9.00000000, 0.00039811
+
+10.00000000, 0.00014311
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-39
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00060165
+0.03000000, 0.00172839
+0.04000000, 0.00172839
+0.05000000, 0.00172839
+0.06000000, 0.00210936
+0.07000000, 0.00210936
+0.08000000, 0.00306789
+0.09000000, 0.00306789
+
+0.10000000, 0.00306789
+0.20000000, 0.01228612
+0.30000000, 0.03228729
+0.40000000, 0.08694724
+0.50000000, 0.10175770
+0.60000000, 0.12494063
+0.70000000, 0.14073961
+0.80000000, 0.17315416
+0.90000000, 0.17493309
+
+1.00000000, 0.18609998
+2.00000000, 0.25061664
+3.00000000, 0.12878595
+4.00000000, 0.05664194
+5.00000000, 0.03018505
+6.00000000, 0.01238711
+7.00000000, 0.00454558
+8.00000000, 0.00174762
+9.00000000, 0.00065438
+
+10.00000000, 0.00024158
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-4
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00021551
+0.03000000, 0.00224554
+0.04000000, 0.00224554
+0.05000000, 0.00224554
+0.06000000, 0.00661799
+0.07000000, 0.00661799
+0.08000000, 0.01160674
+0.09000000, 0.01160674
+
+0.10000000, 0.01160674
+0.20000000, 0.03035110
+0.30000000, 0.05531434
+0.40000000, 0.06615644
+0.50000000, 0.06562372
+0.60000000, 0.09413180
+0.70000000, 0.11547767
+0.80000000, 0.15118615
+0.90000000, 0.18360606
+
+1.00000000, 0.19295633
+2.00000000, 0.17859269
+3.00000000, 0.11393600
+4.00000000, 0.11572037
+5.00000000, 0.06150515
+6.00000000, 0.02844778
+7.00000000, 0.00869475
+8.00000000, 0.00302657
+9.00000000, 0.00106323
+
+10.00000000, 0.00037688
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-40
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00007681
+0.03000000, 0.00104577
+0.04000000, 0.00104577
+0.05000000, 0.00104577
+0.06000000, 0.00190948
+0.07000000, 0.00190948
+0.08000000, 0.00440781
+0.09000000, 0.00440781
+
+0.10000000, 0.00440781
+0.20000000, 0.01830199
+0.30000000, 0.00930771
+0.40000000, 0.01460408
+0.50000000, 0.02502207
+0.60000000, 0.05064099
+0.70000000, 0.04209834
+0.80000000, 0.05902855
+0.90000000, 0.08278442
+
+1.00000000, 0.10502961
+2.00000000, 0.30932523
+3.00000000, 0.27207080
+4.00000000, 0.09503117
+5.00000000, 0.02322673
+6.00000000, 0.00818649
+7.00000000, 0.00290842
+8.00000000, 0.00115485
+9.00000000, 0.00044503
+
+10.00000000, 0.00016786
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-41
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00023856
+0.03000000, 0.00138278
+0.04000000, 0.00138278
+0.05000000, 0.00138278
+0.06000000, 0.00330911
+0.07000000, 0.00330911
+0.08000000, 0.00594384
+0.09000000, 0.00594384
+
+0.10000000, 0.00594384
+0.20000000, 0.02405435
+0.30000000, 0.06167601
+0.40000000, 0.11437778
+0.50000000, 0.14925404
+0.60000000, 0.20870395
+0.70000000, 0.26989612
+0.80000000, 0.33518829
+0.90000000, 0.40242950
+
+1.00000000, 0.47315553
+2.00000000, 0.69451508
+3.00000000, 0.31392019
+4.00000000, 0.06823992
+5.00000000, 0.02688265
+6.00000000, 0.01092475
+7.00000000, 0.00358305
+8.00000000, 0.00133232
+9.00000000, 0.00049093
+
+10.00000000, 0.00017982
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-42
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00033800
+0.03000000, 0.00126940
+0.04000000, 0.00126940
+0.05000000, 0.00126940
+0.06000000, 0.00233733
+0.07000000, 0.00233733
+0.08000000, 0.00388290
+0.09000000, 0.00388290
+
+0.10000000, 0.00388290
+0.20000000, 0.01618749
+0.30000000, 0.03356920
+0.40000000, 0.06358528
+0.50000000, 0.07988384
+0.60000000, 0.10919107
+0.70000000, 0.14224151
+0.80000000, 0.15531738
+0.90000000, 0.17516841
+
+1.00000000, 0.18767215
+2.00000000, 0.11215741
+3.00000000, 0.05297956
+4.00000000, 0.10610819
+5.00000000, 0.03847118
+6.00000000, 0.01090498
+7.00000000, 0.00414386
+8.00000000, 0.00159398
+9.00000000, 0.00059731
+
+10.00000000, 0.00022094
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-43
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00058069
+0.03000000, 0.00572290
+0.04000000, 0.00572290
+0.05000000, 0.00572290
+0.06000000, 0.01150127
+0.07000000, 0.01150127
+0.08000000, 0.01728139
+0.09000000, 0.01728139
+
+0.10000000, 0.01728139
+0.20000000, 0.05059752
+0.30000000, 0.12544900
+0.40000000, 0.19684207
+0.50000000, 0.25611018
+0.60000000, 0.26871311
+0.70000000, 0.31724362
+0.80000000, 0.38514943
+0.90000000, 0.40513448
+
+1.00000000, 0.41372486
+2.00000000, 0.27894912
+3.00000000, 0.11009050
+4.00000000, 0.08026537
+5.00000000, 0.03798018
+6.00000000, 0.01296127
+7.00000000, 0.00429001
+8.00000000, 0.00153389
+9.00000000, 0.00055107
+
+10.00000000, 0.00019899
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-44
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00039676
+0.03000000, 0.00213163
+0.04000000, 0.00213163
+0.05000000, 0.00213163
+0.06000000, 0.00528338
+0.07000000, 0.00528338
+0.08000000, 0.01070260
+0.09000000, 0.01070260
+
+0.10000000, 0.01070260
+0.20000000, 0.02402567
+0.30000000, 0.02503696
+0.40000000, 0.03397484
+0.50000000, 0.04116517
+0.60000000, 0.05109567
+0.70000000, 0.08087052
+0.80000000, 0.12756143
+0.90000000, 0.14614018
+
+1.00000000, 0.18487367
+2.00000000, 0.24024465
+3.00000000, 0.07904517
+4.00000000, 0.06310015
+5.00000000, 0.02526644
+6.00000000, 0.00859481
+7.00000000, 0.00341846
+8.00000000, 0.00136202
+9.00000000, 0.00051873
+
+10.00000000, 0.00019343
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-45
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00040893
+0.03000000, 0.00221951
+0.04000000, 0.00221951
+0.05000000, 0.00221951
+0.06000000, 0.00497420
+0.07000000, 0.00497420
+0.08000000, 0.00826578
+0.09000000, 0.00826578
+
+0.10000000, 0.00826578
+0.20000000, 0.01261189
+0.30000000, 0.01683114
+0.40000000, 0.01694151
+0.50000000, 0.02851971
+0.60000000, 0.04587074
+0.70000000, 0.06638988
+0.80000000, 0.10589113
+0.90000000, 0.14115750
+
+1.00000000, 0.17936953
+2.00000000, 0.37345128
+3.00000000, 0.14302137
+4.00000000, 0.04911909
+5.00000000, 0.03191004
+6.00000000, 0.01283469
+7.00000000, 0.00450360
+8.00000000, 0.00171255
+9.00000000, 0.00063697
+
+10.00000000, 0.00023415
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-46
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00023975
+0.03000000, 0.00297088
+0.04000000, 0.00297088
+0.05000000, 0.00297088
+0.06000000, 0.00651951
+0.07000000, 0.00651951
+0.08000000, 0.01104461
+0.09000000, 0.01104461
+
+0.10000000, 0.01104461
+0.20000000, 0.02393119
+0.30000000, 0.04396861
+0.40000000, 0.04771098
+0.50000000, 0.04603519
+0.60000000, 0.05771251
+0.70000000, 0.08907210
+0.80000000, 0.13463829
+0.90000000, 0.16920974
+
+1.00000000, 0.17879289
+2.00000000, 0.22418029
+3.00000000, 0.22217176
+4.00000000, 0.15186315
+5.00000000, 0.07385494
+6.00000000, 0.02573448
+7.00000000, 0.00900027
+8.00000000, 0.00332338
+9.00000000, 0.00120875
+
+10.00000000, 0.00043806
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-47
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00019155
+0.03000000, 0.00225743
+0.04000000, 0.00225743
+0.05000000, 0.00225743
+0.06000000, 0.00691450
+0.07000000, 0.00691450
+0.08000000, 0.01166792
+0.09000000, 0.01166792
+
+0.10000000, 0.01166792
+0.20000000, 0.04180395
+0.30000000, 0.07265757
+0.40000000, 0.09369890
+0.50000000, 0.16679035
+0.60000000, 0.17536885
+0.70000000, 0.17577436
+0.80000000, 0.20627450
+0.90000000, 0.19926381
+
+1.00000000, 0.18268000
+2.00000000, 0.11680791
+3.00000000, 0.16839339
+4.00000000, 0.07573513
+5.00000000, 0.04546650
+6.00000000, 0.02248870
+7.00000000, 0.00787638
+8.00000000, 0.00301106
+9.00000000, 0.00112524
+
+10.00000000, 0.00041513
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-48
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00007954
+0.03000000, 0.00046739
+0.04000000, 0.00046739
+0.05000000, 0.00046739
+0.06000000, 0.00219974
+0.07000000, 0.00219974
+0.08000000, 0.00358893
+0.09000000, 0.00358893
+
+0.10000000, 0.00358893
+0.20000000, 0.00830257
+0.30000000, 0.02597832
+0.40000000, 0.04738509
+0.50000000, 0.07273911
+0.60000000, 0.09520412
+0.70000000, 0.12135333
+0.80000000, 0.15441627
+0.90000000, 0.20632291
+
+1.00000000, 0.24635273
+2.00000000, 0.22366808
+3.00000000, 0.05642008
+4.00000000, 0.03837508
+5.00000000, 0.01789558
+6.00000000, 0.00675554
+7.00000000, 0.00260250
+8.00000000, 0.00098831
+9.00000000, 0.00036725
+
+10.00000000, 0.00013528
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-49
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00048310
+0.03000000, 0.00140501
+0.04000000, 0.00140501
+0.05000000, 0.00140501
+0.06000000, 0.00302532
+0.07000000, 0.00302532
+0.08000000, 0.00526125
+0.09000000, 0.00526125
+
+0.10000000, 0.00526125
+0.20000000, 0.02862632
+0.30000000, 0.06256233
+0.40000000, 0.08923162
+0.50000000, 0.10773323
+0.60000000, 0.11603681
+0.70000000, 0.14115306
+0.80000000, 0.17620384
+0.90000000, 0.21323434
+
+1.00000000, 0.25688673
+2.00000000, 0.35769890
+3.00000000, 0.14767975
+4.00000000, 0.10219618
+5.00000000, 0.05742952
+6.00000000, 0.02653704
+7.00000000, 0.00991575
+8.00000000, 0.00375206
+9.00000000, 0.00138512
+
+10.00000000, 0.00050666
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-5
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00035606
+0.03000000, 0.00224551
+0.04000000, 0.00224551
+0.05000000, 0.00224551
+0.06000000, 0.00405798
+0.07000000, 0.00405798
+0.08000000, 0.00828472
+0.09000000, 0.00828472
+
+0.10000000, 0.00828472
+0.20000000, 0.02794044
+0.30000000, 0.07896036
+0.40000000, 0.17195891
+0.50000000, 0.26351770
+0.60000000, 0.31970350
+0.70000000, 0.36252141
+0.80000000, 0.42596216
+0.90000000, 0.51206231
+
+1.00000000, 0.56764746
+2.00000000, 0.44861719
+3.00000000, 0.16829697
+4.00000000, 0.08112137
+5.00000000, 0.03731441
+6.00000000, 0.00977531
+7.00000000, 0.00335245
+8.00000000, 0.00122115
+9.00000000, 0.00044144
+
+10.00000000, 0.00015961
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-50
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00027612
+0.03000000, 0.00306060
+0.04000000, 0.00306060
+0.05000000, 0.00306060
+0.06000000, 0.00765728
+0.07000000, 0.00765728
+0.08000000, 0.01272641
+0.09000000, 0.01272641
+
+0.10000000, 0.01272641
+0.20000000, 0.04950321
+0.30000000, 0.09138903
+0.40000000, 0.10136138
+0.50000000, 0.10126915
+0.60000000, 0.13202277
+0.70000000, 0.18325527
+0.80000000, 0.23857281
+0.90000000, 0.32665115
+
+1.00000000, 0.39604977
+2.00000000, 0.42367198
+3.00000000, 0.26300087
+4.00000000, 0.09466751
+5.00000000, 0.03554137
+6.00000000, 0.01518868
+7.00000000, 0.00504168
+8.00000000, 0.00188622
+9.00000000, 0.00069946
+
+10.00000000, 0.00025746
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-51
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00003547
+0.03000000, 0.00124689
+0.04000000, 0.00124689
+0.05000000, 0.00124689
+0.06000000, 0.00267384
+0.07000000, 0.00267384
+0.08000000, 0.00438088
+0.09000000, 0.00438088
+
+0.10000000, 0.00438088
+0.20000000, 0.01169854
+0.30000000, 0.02576251
+0.40000000, 0.04459882
+0.50000000, 0.06107642
+0.60000000, 0.08937920
+0.70000000, 0.10475070
+0.80000000, 0.10393342
+0.90000000, 0.12497223
+
+1.00000000, 0.16205751
+2.00000000, 0.45922166
+3.00000000, 0.36921682
+4.00000000, 0.22575608
+5.00000000, 0.09703175
+6.00000000, 0.03182379
+7.00000000, 0.01099894
+8.00000000, 0.00411390
+9.00000000, 0.00151289
+
+10.00000000, 0.00055243
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-52
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00051033
+0.03000000, 0.00308554
+0.04000000, 0.00308554
+0.05000000, 0.00308554
+0.06000000, 0.00646052
+0.07000000, 0.00646052
+0.08000000, 0.01269048
+0.09000000, 0.01269048
+
+0.10000000, 0.01269048
+0.20000000, 0.05696741
+0.30000000, 0.10861213
+0.40000000, 0.15602429
+0.50000000, 0.19437902
+0.60000000, 0.23744222
+0.70000000, 0.31686771
+0.80000000, 0.39476500
+0.90000000, 0.43541495
+
+1.00000000, 0.46886233
+2.00000000, 0.40954359
+3.00000000, 0.23634797
+4.00000000, 0.08696404
+5.00000000, 0.03938226
+6.00000000, 0.01743951
+7.00000000, 0.00531294
+8.00000000, 0.00188340
+9.00000000, 0.00067112
+
+10.00000000, 0.00024032
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-53
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00098511
+0.03000000, 0.00806230
+0.04000000, 0.00806230
+0.05000000, 0.00806230
+0.06000000, 0.01583426
+0.07000000, 0.01583426
+0.08000000, 0.02557612
+0.09000000, 0.02557612
+
+0.10000000, 0.02557612
+0.20000000, 0.09179562
+0.30000000, 0.15406945
+0.40000000, 0.28103813
+0.50000000, 0.36544053
+0.60000000, 0.38509232
+0.70000000, 0.38268018
+0.80000000, 0.35555981
+0.90000000, 0.35817425
+
+1.00000000, 0.36733585
+2.00000000, 0.23481523
+3.00000000, 0.11365606
+4.00000000, 0.04094189
+5.00000000, 0.02681471
+6.00000000, 0.00877579
+7.00000000, 0.00267532
+8.00000000, 0.00093301
+9.00000000, 0.00032741
+
+10.00000000, 0.00011606
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-54
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00103855
+0.03000000, 0.00540361
+0.04000000, 0.00540361
+0.05000000, 0.00540361
+0.06000000, 0.01105754
+0.07000000, 0.01105754
+0.08000000, 0.01862196
+0.09000000, 0.01862196
+
+0.10000000, 0.01862196
+0.20000000, 0.05829738
+0.30000000, 0.13213130
+0.40000000, 0.21201254
+0.50000000, 0.26066592
+0.60000000, 0.29640792
+0.70000000, 0.32890243
+0.80000000, 0.32862962
+0.90000000, 0.34523026
+
+1.00000000, 0.36258289
+2.00000000, 0.29895121
+3.00000000, 0.14898750
+4.00000000, 0.10688406
+5.00000000, 0.06813457
+6.00000000, 0.01971138
+7.00000000, 0.00696399
+8.00000000, 0.00257578
+9.00000000, 0.00093629
+
+10.00000000, 0.00033883
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-55
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00021477
+0.03000000, 0.00118745
+0.04000000, 0.00118745
+0.05000000, 0.00118745
+0.06000000, 0.00311650
+0.07000000, 0.00311650
+0.08000000, 0.00384830
+0.09000000, 0.00384830
+
+0.10000000, 0.00384830
+0.20000000, 0.02103616
+0.30000000, 0.04520517
+0.40000000, 0.06981078
+0.50000000, 0.10763590
+0.60000000, 0.13360331
+0.70000000, 0.14627328
+0.80000000, 0.19315261
+0.90000000, 0.28472746
+
+1.00000000, 0.33788664
+2.00000000, 0.35749485
+3.00000000, 0.45387971
+4.00000000, 0.22592810
+5.00000000, 0.06299886
+6.00000000, 0.02418928
+7.00000000, 0.00863462
+8.00000000, 0.00320456
+9.00000000, 0.00117152
+
+10.00000000, 0.00042622
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-56
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00093879
+0.03000000, 0.00480523
+0.04000000, 0.00480523
+0.05000000, 0.00480523
+0.06000000, 0.01185623
+0.07000000, 0.01185623
+0.08000000, 0.01871339
+0.09000000, 0.01871339
+
+0.10000000, 0.01871339
+0.20000000, 0.04496247
+0.30000000, 0.08822107
+0.40000000, 0.09789580
+0.50000000, 0.13659957
+0.60000000, 0.18486653
+0.70000000, 0.21472577
+0.80000000, 0.23532734
+0.90000000, 0.27069391
+
+1.00000000, 0.27790579
+2.00000000, 0.30571453
+3.00000000, 0.15555520
+4.00000000, 0.05144647
+5.00000000, 0.02198854
+6.00000000, 0.00723104
+7.00000000, 0.00262127
+8.00000000, 0.00099229
+9.00000000, 0.00036874
+
+10.00000000, 0.00013581
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-57
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00012774
+0.03000000, 0.00269532
+0.04000000, 0.00269532
+0.05000000, 0.00269532
+0.06000000, 0.00899968
+0.07000000, 0.00899968
+0.08000000, 0.01853102
+0.09000000, 0.01853102
+
+0.10000000, 0.01853102
+0.20000000, 0.05859345
+0.30000000, 0.10051346
+0.40000000, 0.15538554
+0.50000000, 0.20858157
+0.60000000, 0.25908285
+0.70000000, 0.30063010
+0.80000000, 0.35038881
+0.90000000, 0.40808196
+
+1.00000000, 0.46544003
+2.00000000, 0.57945522
+3.00000000, 0.43119986
+4.00000000, 0.19944554
+5.00000000, 0.06897897
+6.00000000, 0.02502719
+7.00000000, 0.00894608
+8.00000000, 0.00339316
+9.00000000, 0.00125553
+
+10.00000000, 0.00045944
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-58
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00056136
+0.03000000, 0.00366859
+0.04000000, 0.00366859
+0.05000000, 0.00366859
+0.06000000, 0.00707414
+0.07000000, 0.00707414
+0.08000000, 0.01436265
+0.09000000, 0.01436265
+
+0.10000000, 0.01436265
+0.20000000, 0.03116785
+0.30000000, 0.06882938
+0.40000000, 0.10010854
+0.50000000, 0.12877462
+0.60000000, 0.13586528
+0.70000000, 0.12192188
+0.80000000, 0.13007977
+0.90000000, 0.12353922
+
+1.00000000, 0.13136095
+2.00000000, 0.31354214
+3.00000000, 0.29355078
+4.00000000, 0.18687690
+5.00000000, 0.08947194
+6.00000000, 0.03428247
+7.00000000, 0.01174939
+8.00000000, 0.00434122
+9.00000000, 0.00158304
+
+10.00000000, 0.00057490
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-59
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00025979
+0.03000000, 0.00166523
+0.04000000, 0.00166523
+0.05000000, 0.00166523
+0.06000000, 0.00389786
+0.07000000, 0.00389786
+0.08000000, 0.00579291
+0.09000000, 0.00579291
+
+0.10000000, 0.00579291
+0.20000000, 0.01590012
+0.30000000, 0.03913995
+0.40000000, 0.06966625
+0.50000000, 0.10848242
+0.60000000, 0.14313932
+0.70000000, 0.15772038
+0.80000000, 0.19305409
+0.90000000, 0.25063975
+
+1.00000000, 0.32932212
+2.00000000, 0.39278339
+3.00000000, 0.10515472
+4.00000000, 0.07677066
+5.00000000, 0.02345515
+6.00000000, 0.00764628
+7.00000000, 0.00210012
+8.00000000, 0.00072825
+9.00000000, 0.00025447
+
+10.00000000, 0.00008962
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-6
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00029756
+0.03000000, 0.00286892
+0.04000000, 0.00286892
+0.05000000, 0.00286892
+0.06000000, 0.00980902
+0.07000000, 0.00980902
+0.08000000, 0.01849286
+0.09000000, 0.01849286
+
+0.10000000, 0.01849286
+0.20000000, 0.06946752
+0.30000000, 0.13572957
+0.40000000, 0.22580272
+0.50000000, 0.28968489
+0.60000000, 0.34883173
+0.70000000, 0.38803280
+0.80000000, 0.39282118
+0.90000000, 0.37764953
+
+1.00000000, 0.35215724
+2.00000000, 0.31974328
+3.00000000, 0.22162262
+4.00000000, 0.11546990
+5.00000000, 0.03910396
+6.00000000, 0.00965783
+7.00000000, 0.00367604
+8.00000000, 0.00141245
+9.00000000, 0.00052639
+
+10.00000000, 0.00019374
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-60
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00049911
+0.03000000, 0.00333756
+0.04000000, 0.00333756
+0.05000000, 0.00333756
+0.06000000, 0.00578474
+0.07000000, 0.00578474
+0.08000000, 0.00722605
+0.09000000, 0.00722605
+
+0.10000000, 0.00722605
+0.20000000, 0.02458061
+0.30000000, 0.04572145
+0.40000000, 0.07798331
+0.50000000, 0.11172813
+0.60000000, 0.13545957
+0.70000000, 0.14845684
+0.80000000, 0.17994198
+0.90000000, 0.21212447
+
+1.00000000, 0.21352425
+2.00000000, 0.20476726
+3.00000000, 0.16563578
+4.00000000, 0.12744227
+5.00000000, 0.05094379
+6.00000000, 0.01819931
+7.00000000, 0.00643555
+8.00000000, 0.00241018
+9.00000000, 0.00088607
+
+10.00000000, 0.00032326
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-61
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00028092
+0.03000000, 0.00223059
+0.04000000, 0.00223059
+0.05000000, 0.00223059
+0.06000000, 0.00631606
+0.07000000, 0.00631606
+0.08000000, 0.01161108
+0.09000000, 0.01161108
+
+0.10000000, 0.01161108
+0.20000000, 0.03092378
+0.30000000, 0.04722237
+0.40000000, 0.05316468
+0.50000000, 0.05223254
+0.60000000, 0.05190446
+0.70000000, 0.06191705
+0.80000000, 0.07842446
+0.90000000, 0.08626182
+
+1.00000000, 0.09986281
+2.00000000, 0.26548420
+3.00000000, 0.24033655
+4.00000000, 0.12380974
+5.00000000, 0.06290249
+6.00000000, 0.02026846
+7.00000000, 0.00718829
+8.00000000, 0.00265170
+9.00000000, 0.00096366
+
+10.00000000, 0.00034886
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-62
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00036125
+0.03000000, 0.00133702
+0.04000000, 0.00133702
+0.05000000, 0.00133702
+0.06000000, 0.00290278
+0.07000000, 0.00290278
+0.08000000, 0.00328575
+0.09000000, 0.00328575
+
+0.10000000, 0.00328575
+0.20000000, 0.01280963
+0.30000000, 0.02792423
+0.40000000, 0.06010806
+0.50000000, 0.08011378
+0.60000000, 0.08256097
+0.70000000, 0.10002715
+0.80000000, 0.11364393
+0.90000000, 0.13488343
+
+1.00000000, 0.14257682
+2.00000000, 0.11598887
+3.00000000, 0.12725415
+4.00000000, 0.10983488
+5.00000000, 0.06545753
+6.00000000, 0.02321335
+7.00000000, 0.00866787
+8.00000000, 0.00333026
+9.00000000, 0.00124419
+
+10.00000000, 0.00045841
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-63
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00017752
+0.03000000, 0.00163663
+0.04000000, 0.00163663
+0.05000000, 0.00163663
+0.06000000, 0.00463343
+0.07000000, 0.00463343
+0.08000000, 0.00934003
+0.09000000, 0.00934003
+
+0.10000000, 0.00934003
+0.20000000, 0.02010128
+0.30000000, 0.04725794
+0.40000000, 0.07876488
+0.50000000, 0.14036701
+0.60000000, 0.18790975
+0.70000000, 0.24074274
+0.80000000, 0.26908579
+0.90000000, 0.30431087
+
+1.00000000, 0.32200797
+2.00000000, 0.45991147
+3.00000000, 0.22063416
+4.00000000, 0.13927227
+5.00000000, 0.05384444
+6.00000000, 0.02009056
+7.00000000, 0.00756976
+8.00000000, 0.00288153
+9.00000000, 0.00106839
+
+10.00000000, 0.00039183
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-64
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00025824
+0.03000000, 0.00269791
+0.04000000, 0.00269791
+0.05000000, 0.00269791
+0.06000000, 0.00624686
+0.07000000, 0.00624686
+0.08000000, 0.01175296
+0.09000000, 0.01175296
+
+0.10000000, 0.01175296
+0.20000000, 0.04352788
+0.30000000, 0.07626763
+0.40000000, 0.10596032
+0.50000000, 0.14262971
+0.60000000, 0.19063981
+0.70000000, 0.20515913
+0.80000000, 0.22919022
+0.90000000, 0.22923774
+
+1.00000000, 0.23615862
+2.00000000, 0.32052491
+3.00000000, 0.17704332
+4.00000000, 0.09128647
+5.00000000, 0.02914677
+6.00000000, 0.00830269
+7.00000000, 0.00276652
+8.00000000, 0.00101677
+9.00000000, 0.00036965
+
+10.00000000, 0.00013392
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
+# Reading Wilson line: binary format, size is 1024, Nc 3, length is [fm] 5.12, a is [fm]0.005
+# Dipole amplitudes for b=0
+# 1 rotations 
+# Filename /n/work01/pyjorunk/JIMWLK/ipglasma/outputs//V-65
+# -------------------------
+# r, D_p(r,0)
+
+0.00010000, 0.00000000
+0.00020000, 0.00000000
+0.00030000, 0.00000000
+0.00040000, 0.00000000
+0.00050000, 0.00000000
+0.00060000, 0.00000000
+0.00070000, 0.00000000
+0.00080000, 0.00000000
+0.00090000, 0.00000000
+
+0.00100000, 0.00000000
+0.00200000, 0.00000000
+0.00300000, 0.00000000
+0.00400000, 0.00000000
+0.00500000, 0.00000000
+0.00600000, 0.00000000
+0.00700000, 0.00000000
+0.00800000, 0.00000000
+0.00900000, 0.00000000
+
+0.01000000, 0.00000000
+0.02000000, 0.00022247
+0.03000000, 0.00076631
+0.04000000, 0.00076631
+0.05000000, 0.00076631
+0.06000000, 0.00062860
+0.07000000, 0.00062860
+0.08000000, 0.00162917
+0.09000000, 0.00162917
+
+0.10000000, 0.00162917
+0.20000000, 0.00655856
+0.30000000, 0.01362686
+0.40000000, 0.01669954
+0.50000000, 0.02102393
+0.60000000, 0.04134764
+0.70000000, 0.07485351
+0.80000000, 0.09763002
+0.90000000, 0.15298641
+
+1.00000000, 0.21205807
+2.00000000, 0.52749953
+3.00000000, 0.26724634
+4.00000000, 0.18283752
+5.00000000, 0.07228655
+6.00000000, 0.02288751
+7.00000000, 0.00774401
+8.00000000, 0.00289826
+9.00000000, 0.00106912
+
+10.00000000, 0.00039139
+20.00000000, 0.00000000
+30.00000000, 0.00000000
+40.00000000, 0.00000000
+50.00000000, 0.00000000
+60.00000000, 0.00000000
+70.00000000, 0.00000000
+80.00000000, 0.00000000
+90.00000000, 0.00000000
diff --git a/src/ipglasma.cpp b/src/ipglasma.cpp
index 1b4a818913a4ebd869d600474e2455d82b8e400a..a683f5dec1b0de7281004ca4da6d004c6c5bcd67 100644
--- a/src/ipglasma.cpp
+++ b/src/ipglasma.cpp
@@ -470,3 +470,207 @@ int IPGlasma::WilsonLineCoordinate(int  xind, int yind)
 {
     return xcoords.size()*xind + yind; 
 }
+
+
+// attempt to overload amplitude with one wilson line and two indices
+double IPGlasma::Amplitude(double xpom, WilsonLine quark, int q2[2])
+{
+    //TODO: make work with int?
+    //ApplyPeriodicBoundaryConditions(q1);
+    //ApplyPeriodicBoundaryConditions(q2);
+
+    // TODO: make work with indices (xMax?
+    // Out of grid? Return 0 (probably very large dipole)
+
+    //if (q1[0] < xcoords[0] or q1[0] > xcoords[xcoords.size()-1]
+    //    or q1[1] < ycoords[0] or q1[1] > ycoords[ycoords.size()-1]
+    //    etc...
+    // ?TODO?: (double?) check if dipole is less than lattice spacing   
+
+    // First find corresponding grid indeces
+    WilsonLine antiquark = GetWilsonLine(WilsonLineCoordinate(q2[0], q2[1]));
+
+    WilsonLine prod;
+    try {
+        prod = quark.MultiplyByHermitianConjugate(antiquark);
+    } catch (...) {
+        cerr << "Matrix multiplication failed!" << endl;
+        cout << quark << endl;
+        cout << "Antiquark: " << q2[0] << ", " << q2[1] << endl;
+        cout << antiquark << endl;
+        exit(1);
+    }
+    std::complex<double > amp =  1.0 - 1.0/NC * prod.Trace();
+
+    double result = amp.real();
+    if (result < 0) return 0;
+    return result;
+    if (result > 1)
+        return 1;
+    if (result < 0)
+        return 0;
+
+    if (isnan(result))
+    {
+	cerr << "Wilson line trance NaN, antiquark coords " << q2[0] << ", " << q2[1] << "and quark wline" << quark << endl;
+        exit(1);
+    }
+
+    return result;
+}
+
+// attempt to overload amplitude when given coordinates are integers i.e. indices
+double IPGlasma::Amplitude(double xpom, int q1[2], int q2[2])
+{
+    //TODO: make work with int?
+    //ApplyPeriodicBoundaryConditions(q1);
+    //ApplyPeriodicBoundaryConditions(q2);
+
+    // TODO: make work with indices (xMax?
+    // Out of grid? Return 0 (probably very large dipole)
+
+    //if (q1[0] < xcoords[0] or q1[0] > xcoords[xcoords.size()-1]
+    //    or q1[1] < ycoords[0] or q1[1] > ycoords[ycoords.size()-1]
+    //    etc...
+    // ?TODO?: (double?) check if dipole is less than lattice spacing   
+
+    // First find corresponding grid indeces
+    WilsonLine quark = GetWilsonLine(WilsonLineCoordinate(q1[0], q1[1]));
+    WilsonLine antiquark = GetWilsonLine(WilsonLineCoordinate(q2[0], q2[1]));
+
+    //antiquark = antiquark.HermitianConjugate();
+
+    WilsonLine prod;
+    try {
+        prod = quark.MultiplyByHermitianConjugate(antiquark);
+        //prod =  quark*antiquark;
+    } catch (...) {
+        cerr << "Matrix multiplication failed!" << endl;
+        cout << "Quark: " << q1[0] << ", " << q1[1] << endl;
+        cout << quark << endl;
+        cout << "Antiquark: " << q2[0] << ", " << q2[1] << endl;
+        cout << antiquark << endl;
+        exit(1);
+    }
+    std::complex<double > amp =  1.0 - 1.0/NC * prod.Trace();
+
+    double result = amp.real();
+    if (result < 0) return 0;
+    return result;
+    if (result > 1)
+        return 1;
+    if (result < 0)
+        return 0;
+
+    if (isnan(result))
+    {
+        cerr << "Wilson line trance NaN, quark coords " << q1[0] << ", " << q1[1] << " and " << q2[0] << ", " << q2[1] << endl;
+        exit(1);
+    }
+
+    return result;
+}
+
+//TODO: check
+double IPGlasma::DistanceToOrigin(int q[2]){
+    // Get x and y coordinates corresponding to these indices
+    double x = X(q[0]);
+    double y = Y(q[1]);
+    // calculate length of the coordinate vector
+    return sqrt(x*x+y*y);
+
+}
+
+
+//TODO: check
+double IPGlasma::SumAmplitudes(int l)
+{
+    double sum = 0.0;
+    int L = xcoords.size();
+    cout << "# Summing with lattice size L = " << L << endl;
+    WilsonLine quark;
+    for(int i =0; i < L; i += 1){
+        for(int j =0; j < L; j += 1){
+            
+            int q1[2] = {i,j};
+            //TODO: pointer instead of copy?
+            quark = GetWilsonLine(WilsonLineCoordinate(i, j));
+            for(int k =0; k < L; k +=1){
+                int q2[2] = {k,l};
+                double r1 = DistanceToOrigin(q1);
+                double r2 = DistanceToOrigin(q2);
+                if (r1 > 1. && r2 > 1.){
+                    //cout << "#!!Throwing away large dipole separation!" << endl;
+                    //cout << "#Quark coordinate indices " << q1[0] << "," << q1[1] << " and " << q2[0] << "," << q2[1] << endl;
+                    //cout << "#Dipole size is " << r1-r2 << endl;
+                    //cout << "#Dipole amplitude is " << Amplitude(0, q1, q2) << endl;
+                    continue;
+                }
+
+                // Haw w2 = WilsonLine(k,l)
+                // Laske w1.MultiplbyHermitianConjugate(w2);
+                // sum += amplitude(0, q1, q2, ipglasma);
+                sum += Amplitude(0, quark, q2);
+            }
+        }
+    }
+    cout << "# Intermediate sum" << endl;
+    cout << sum << endl;
+    cout << "#---------" << endl;
+    
+    return sum;
+}
+
+//TODO: check
+double IPGlasma::SumAmplitudes()
+{
+    double sum = 0.0;
+    int L = xcoords.size();
+    cout << "# Summing with lattice size L = " << L << endl;
+    WilsonLine quark;
+    int thrownAway = 0;
+    int iWasHere = 0;
+    #pragma omp parallel for reduction(+:sum,thrownAway,iWasHere) collapse(4)
+    for(int i =0; i < L; i += 1){
+        for(int j =0; j < L; j += 1){
+            
+            //#pragma omp parallel for reduction(+:sum) collapse(2)
+            for(int k =0; k < L; k +=1){
+//                cout << "# Thread " << omp_get_thread_num() << " is processing (" << i << ", " << j << ", " << k <<")" << endl;
+                for(int l =0; l < L; l +=1){
+                    cout << "# Hello from inside the l loop!" << endl;
+            //TODO: pointer instead of copy?
+                    //TODO: faster to get these two layers above?
+                    iWasHere += 1;
+                    int q1[2] = {i,j};
+                    //quark = GetWilsonLine(WilsonLineCoordinate(i, j));
+                    int q2[2] = {k,l};
+		    double r1 = DistanceToOrigin(q1);
+                    double r2 = DistanceToOrigin(q2);
+                    if (r1 > 1. && r2 > 1.){
+                       // cout << "#!!Throwing away large dipole separation!" << endl;
+                        //cout << "#Quark coordinate indices " << q1[0] << "," << q1[1] << " and " << q2[0] << "," << q2[1] << endl;
+                        //cout << "#Quark coordinates " << X(q1[0]) << "," << Y(q1[1]) << " and " << X(q2[0]) << "," << Y(q2[1]) << endl;
+                        //cout << "#Dipole size is " << r1-r2 << endl;
+                        //cout << "#Dipole amplitude is " << Amplitude(0, q1, q2) << endl;
+                        thrownAway += 1;
+                        continue;
+                }
+
+                    // Haw w2 = WilsonLine(k,l)
+                    // Laske w1.MultiplbyHermitianConjugate(w2);
+
+                   sum += Amplitude(0, q1, q2);
+                   //sum += Amplitude(0, quark, q2);
+                   iWasHere += 1;  
+                }
+            }
+        }
+    }
+    cout << "# total sum = " << sum << endl;
+    cout << "# i counted " << iWasHere << " values" << endl;
+    cout << "# threw away " << thrownAway << " values (" << 100.*thrownAway*pow(L,-4) << "%)" << endl;
+    return sum;
+}
+
+
diff --git a/src/ipglasma.hpp b/src/ipglasma.hpp
index 2f8c52288dccec5930d2d271d3465d7fcf97de07..1a979b92fcc70346ef06b5ae3daae602aba72aad 100644
--- a/src/ipglasma.hpp
+++ b/src/ipglasma.hpp
@@ -54,6 +54,14 @@ public:
 
     double X(int ix) { return xcoords[ix]; }
     double Y(int iy) { return ycoords[iy]; }
+    
+    //Pyry added:
+    double DistanceToOrigin(int q[2]);
+    double Amplitude(double xpom, int q1[2], int q2[2]);
+    double Amplitude(double xpom, WilsonLine quark, int q2[2]);
+    double SumAmplitudes();
+    double SumAmplitudes(int l);
+    
 
 private:
     
diff --git a/src/main.cpp b/src/main.cpp
index 9f9423f0cee7498e33fc35e988e920f6def38e52..74866c703d91e361975336bdaae08e6a230318f3 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -1,19 +1,187 @@
 #include <iostream>
 #include <string>
 #include "ipglasma.hpp"
+#include <cmath>
+#include <iostream>
+#include <fstream>
+#include <iomanip>
+#include <omp.h>
 
 using namespace std;
 
+int printDipoleAmplitudes(string fname, IPGlasma ipglasma)
+{
+    //ofstream output0;
+    //Initialize params
+    double q1[2] = {0,0};
+    double q2[2] = {0,0};
+    double r = 0;
+    double dipoleAmp = 0;
+    //TODO: set number of different rotations for dipole!
+    int N = 100;
+
+    cout << fixed << setprecision(8);
+    
+    cout << "# Dipole amplitudes for b=0"<< endl;
+    cout << "# " << N << " rotations " << endl;
+    cout << "# Filename " << fname << endl;
+    cout << "# -------------------------"<< endl;
+    cout << "# r, D_p(r,0)" << endl;
+
+    // loop to calculate dipole amplitude at radii r=j*10^i
+    for (int i=-4; i < 2; i += 1){
+        cout <<""<< endl;
+
+        // use index j to get radius
+        for (double j=1; j < 10; j += 1){
+            r = j*pow(10, i);
+            // go through different orientations of the dipole via rotations
+            for (double k=0; k < N; k += 1){
+            // inside cos/sin is k*1/N*2*3.141592
+                q1[0] = 0.5*r*cos(k*1./N*2*3.141592);
+                q1[1] = 0.5*r*sin(k*1./N*2*3.141592);
+                q2[0] = -q1[0];
+                q2[1] = -q1[1];
+    
+                dipoleAmp += ipglasma.Amplitude(0, q1, q2);
+            }
+
+            dipoleAmp = dipoleAmp/N;
+            cout << r << ", " << dipoleAmp << endl;
+            dipoleAmp = 0;
+        }
+    }
+
+    return 0; 
+}
+
+void integrateDipoleAmplitudes(IPGlasma* ipglasma, int task_id)
+{
+    
+    //TODO: does data type matter?
+    double sum = ipglasma->SumAmplitudes(task_id);
+
+    cout << "# Task " << task_id << " partial sum = " << sum << endl; 
+    //cout << "# Dipole amplitude integral = " << sum << endl;
+    
+    // Save partial results
+    ofstream outfile("partial_sum_" + to_string(task_id) + ".txt");
+    outfile << sum;
+    outfile.close();
+}
+
+//TODO: check
+double sumAmplitudes(IPGlasma *ipglasma)
+{
+    #pragma omp parallel
+    {
+        int thread_id = omp_get_thread_num();
+        int num_threads = omp_get_num_threads();
+    
+        // Print the number of threads (only from the master thread to avoid too many prints)
+        if (thread_id == 0) {
+            cout << "# Number of threads: " << num_threads << endl;
+        }
+
+        // Print the thread id for debugging (can be limited to a few threads for easier reading)
+        if (thread_id < 5) {  // Change 5 to a larger number if needed
+             cout << "# Thread " << thread_id << " is starting" << endl;
+        }
+
+    } 
+    cout << "# Hello from main.cpp sumAmplitudes!" << endl;
+    double sum = 0.0;
+    int L = 512; //ipglasma->xcoords.size();
+    cout << "# Summing with lattice size L = " << L << endl;
+    //WilsonLine quark;
+    //int thrownAway = 0;
+    //int iWasHere = 0;
+    #pragma omp parallel for reduction(+:sum)//,thrownAway,iWasHere) //collapse(4)
+
+    for(int i =0; i < L; i += 1){ 
+        for(int j =0; j < L; j += 1){ 
+    
+            //#pragma omp parallel for reduction(+:sum) collapse(2)
+            for(int k =0; k < L; k +=1){
+//                cout << "# Thread " << omp_get_thread_num() << " is processing (" << i << ", " << j << ", " << k <<")" << endl;
+                for(int l =0; l < L; l +=1){
+                    int q1[2] = {i,j};
+                    double r1 = ipglasma->DistanceToOrigin(q1);
+                    //#pragma omp critical
+                    //{
+                    //    cout << "# Hello from inside the l loop!" << endl;
+                    // }
+            //TODO: pointer instead of copy?
+                    //TODO: faster to get these two layers above?
+                    //iWasHere += 1;
+                    //quark = GetWilsonLine(WilsonLineCoordinate(i, j));
+                    int q2[2] = {k,l};
+                    //double r1 = ipglasma->DistanceToOrigin(q1);
+                    double r2 = ipglasma->DistanceToOrigin(q2);
+                    if (r1 > 1. && r2 > 1.){
+                       // cout << "#!!Throwing away large dipole separation!" << endl;
+                        //cout << "#Quark coordinate indices " << q1[0] << "," << q1[1] << " and " << q2[0] << "," << q2[1] << endl;
+                        //cout << "#Quark coordinates " << X(q1[0]) << "," << Y(q1[1]) << " and " << X(q2[0]) << "," << Y(q2[1]) << endl;
+                        //cout << "#Dipole size is " << r1-r2 << endl;
+                        //cout << "#Dipole amplitude is " << Amplitude(0, q1, q2) << endl;
+                        //thrownAway += 1;
+                        continue;
+                }
+
+                    // Haw w2 = WilsonLine(k,l)
+                    // Laske w1.MultiplbyHermitianConjugate(w2);
+
+                   sum += ipglasma->Amplitude(0, q1, q2);
+                   //sum += Amplitude(0, quark, q2);
+                   //iWasHere += 1;  
+                }
+            }
+        }
+    }   
+    cout << "# total sum = " << sum << endl;
+    //cout << "# i counted " << iWasHere << " values" << endl;
+    //cout << "# threw away " << thrownAway << " values (" << 100.*thrownAway*pow(L,-4) << "%)" << endl;
+    return sum;
+}
+
+    //TODO: does data type matter?
+double integrateDipoleAmplitudes(IPGlasma* ipglasma)
+{
+   //TODO: use ipglasma->SumAmplitudes(); 
+    return sumAmplitudes(ipglasma);
+
+}
+
+
 int main(int argc, char* argv[])
 {
-    // Arguments: filaname
+    if (argc < 3 || argc > 4){
+        cerr << "Usage: " << argv[0] << " filename 1(integrate)/0(print amplitudes) [task_id]" << endl;
+    }
+    
     string fname = argv[1];
     IPGlasma ipglasma(fname, 0, BINARY);
 
-    // Get Wilson line at x=1 GeV^-1, y=3 GeV^-1
-    WilsonLine w = ipglasma.GetWilsonLine(1,3);
-    cout << w.Element(1,2) << endl; // Print element [1,2], note that it is a complex number
+    int task_id = 0;
+
+    if (argc > 3){ 
+        task_id= atoi(argv[3]);
+        cout << "#Integrating array job #" << task_id << endl;
+        integrateDipoleAmplitudes(&ipglasma, task_id);   
+    }
 
+    else if (!argv[2]){ 
+        printDipoleAmplitudes(fname, ipglasma);
+    }
+    
+    else{ 
+        double sum = integrateDipoleAmplitudes(&ipglasma);
+        // Save results
+        //std::filesystem::create_directory("output_"+fname);
+        ofstream outfile("da_sum_ +" + fname + ".txt");
+        outfile << sum << endl;
+	outfile.close();
+    }
 
     return 0;
-}   
\ No newline at end of file
+}