diff --git a/.gitignore b/.gitignore
index 1cfddae..b0229c9 100644
--- a/.gitignore
+++ b/.gitignore
@@ -9,12 +9,14 @@ EventBuilder*
 *.seq
 *.png
 *.gif
+*.msh
 Mapper
 AnasenMS
 
 data/
 data_proton/
 Sudarshan/
+wires2d/
 Analyzer_C_ACLiC_dict0713aaa966_dictContent.h
 .gitignore
 Analyzer_C_ACLiC_dict5411fecd5c_dictUmbrella.h
diff --git a/anasen_fem/paraview_plotter.py b/anasen_fem/paraview_plotter.py
index 1a3150f..8c83940 100755
--- a/anasen_fem/paraview_plotter.py
+++ b/anasen_fem/paraview_plotter.py
@@ -1,4 +1,4 @@
-#!~/ParaView-6.1.0-RC1-MPI-Linux-Python3.12-x86_64/bin/pvbatch
+#!/home/vsitaraman/ParaView-6.1.0-RC1-MPI-Linux-Python3.12-x86_64/bin/pvbatch
 import numpy as np
 import sys
 from paraview.simple import *
@@ -9,7 +9,7 @@ contour_filter = Contour(Input=reader,ContourBy = 'potential')
 contour_filter.Isosurfaces = [i for i in np.arange(0,660,650/24.)]
 
 renderView = GetActiveViewOrCreate('RenderView')
-renderView.ViewSize = [800,800]
+renderView.ViewSize = [2000,2000]
 renderView.OrientationAxesVisibility = 0 # Hide axis
 renderView.UseColorPaletteForBackground=0
 renderView.Background = [0.1, 0.1, 0.1] # Set background to dark gray (RGB 0-1)
@@ -19,13 +19,14 @@ renderView.MultiSamples = 8  # 0 disables it, 4-8 is usually sufficient
 ResetCamera()
 
 contour_display = Show(contour_filter, renderView)
-
+contour_display.LineWidth = 3.0          # Increase this for thicker lines
+contour_display.RenderLinesAsTubes = 1    # Makes lines look smoother at high res
 #colorbar
 contour_display_potentialLUT = GetColorTransferFunction('potential', contour_display, separate=True)
 contour_display_potentialLUT.ApplyPreset('Cool to Warm', True)
 contour_display.SetScalarBarVisibility(renderView, True)
 
-#axesGrid = renderView.AxesGrid
+#axesGrid = renderView.AxesGridrfcxgdtv
 #axesGrid.Visibility = 1
 #axesGrid.XTitle = "x (mm)"
 #axesGrid.YTitle = "y (mm)"
diff --git a/anasen_fem/run.py b/anasen_fem/run.py
index 2b4309c..8034e1c 100755
--- a/anasen_fem/run.py
+++ b/anasen_fem/run.py
@@ -2,7 +2,7 @@ import os
 
 #val=-178.3
 val=17.83
-count=11
+count=11    
 while val<178.3+0.1:
     print(val)
     os.system("python3 wires_gmsh2d_bc.py "+str(val))
diff --git a/anasen_fem/scalars.dat.names b/anasen_fem/scalars.dat.names
index b9970f2..947f957 100755
--- a/anasen_fem/scalars.dat.names
+++ b/anasen_fem/scalars.dat.names
@@ -1,8 +1,8 @@
 Metadata for SaveScalars file: ./scalars.dat
 Elmer version: 26.1
-Elmer compilation date: 2026-03-05
+Elmer compilation date: 2026-03-15
 Solver input file: wires2d.sif
-File started at: 2026/03/11 23:33:58
+File started at: 2026/03/23 17:44:06
  
 Variables in columns of matrix: 
    1: res: potential difference
diff --git a/anasen_fem/wires2d.sif b/anasen_fem/wires2d.sif
index 6bb38f6..185c6ac 100755
--- a/anasen_fem/wires2d.sif
+++ b/anasen_fem/wires2d.sif
@@ -87,17 +87,17 @@ End
 
 Boundary Condition 1
   Target Boundaries = 10
-  Potential = 650
-  Calculate Electric Force = True
+  Potential = 0
 End
 
 
 Boundary Condition 2
   Target Boundaries = 20
-  Potential = 0
+  Potential = 650
+  Calculate Electric Force = True
 End
 
-!Boundary Condition 2
-!  Target Boundaries = 30
-!  Potential = 0
-!End
+Boundary Condition 3
+  Target Boundaries = 30
+  Potential = 0
+End
diff --git a/anasen_fem/wires_gmsh2d_bc.py b/anasen_fem/wires_gmsh2d_bc.py
index 65e9446..7ec4144 100755
--- a/anasen_fem/wires_gmsh2d_bc.py
+++ b/anasen_fem/wires_gmsh2d_bc.py
@@ -3,7 +3,7 @@ import gmsh,sys
 
 gmsh.initialize()
 gmsh.model.add("adaptive_mesh")
-gmsh.option.setNumber('General.NumThreads', 4)
+gmsh.option.setNumber('General.NumThreads', 10)
 #gmsh.option.setNumber("Mesh.Adapt.MaxNumberOfElements", 200000)
 #gmsh.option.setNumber("Mesh.Adapt.MaxNumberOfNodes", 200000)
 #gmsh.option.setNumber("Mesh.Adapt.MaxIter",5)
@@ -23,6 +23,12 @@ z_loc = float(sys.argv[1])
 
 k=(2*np.pi/24.)
 
+#guard wires, plane 1 at -zmax/2
+kg=2*np.pi/24.
+offsetg = -4*kg + 2*kg -  np.pi/24 #-pi/4
+xarrg_1 = np.array([32*np.cos(kg*i+offsetg) for i in np.arange(0,24)])
+yarrg_1 = np.array([32*np.sin(kg*i+offsetg) for i in np.arange(0,24)])
+
 #anodes, plane 1 at -zmax/2
 k=-2*np.pi/24.
 offset = 6*k + 3*k #-pi/2
@@ -35,6 +41,11 @@ offsetc = -4*kc + 2*kc -  np.pi/24 #-pi/4
 xarrc_1 = np.array([42*np.cos(kc*i+offsetc) for i in np.arange(0,24)])
 yarrc_1 = np.array([42*np.sin(kc*i+offsetc) for i in np.arange(0,24)])
 
+#guard wires, plane 2 at +zmax/2
+offsetg = offsetg-3*kg
+xarrg_2 = np.array([32*np.cos(kg*i+offsetg) for i in np.arange(0,24)])
+yarrg_2 = np.array([32*np.sin(kg*i+offsetg) for i in np.arange(0,24)])
+
 #anodes, plane 2 at +zmax/2
 offset = offset-3*k
 xarra_2 = np.array([37*np.cos(k*i+offset) for i in np.arange(0,24)])
@@ -45,6 +56,9 @@ offsetc = offsetc-3*kc
 xarrc_2 = np.array([42*np.cos(kc*i+offsetc) for i in np.arange(0,24)])
 yarrc_2 = np.array([42*np.sin(kc*i+offsetc) for i in np.arange(0,24)])
 
+direction_guard_x = xarrg_2 - xarrg_1
+direction_guard_y = yarrg_2 - yarrg_1
+
 direction_anodes_x = xarra_2 - xarra_1
 direction_anodes_y = yarra_2 - yarra_1
 
@@ -52,23 +66,31 @@ direction_cathodes_x = xarrc_2 - xarrc_1
 direction_cathodes_y = yarrc_2 - yarrc_1
 
 t = (z_loc+178.3)/(2*178.3) #z=-178.3 is 0, z=+178.3 is 1
+xloc_g = xarrg_1 + t*direction_guard_x
+yloc_g = yarrg_1 + t*direction_guard_y
 xloc_a = xarra_1 + t*direction_anodes_x
 yloc_a = yarra_1 + t*direction_anodes_y
 xloc_c = xarrc_1 + t*direction_cathodes_x
 yloc_c = yarrc_1 + t*direction_cathodes_y
 
-wire_radius_a = 0.018 #mm
-wire_radius_c = 0.0762 #mm
+# wire_radius_a = 0.018 #mm
+# wire_radius_c = 0.0762 #mm
+# wire_radius_g = 0.0762 #mm
+wire_radius = 0.254 #mm
+guard_wires = []
 anode_wires = []
 cathode_wires = []
+gw_tags = [(3,i+24) for i in range(24)]
 aw_tags = [(3,i) for i in range(24)]
 cw_tags = [(3,i+24) for i in range(24)]
 
 #for i,[xa,ya,xc,yc] in enumerate(zip(xarra_1,yarra_1,xarrc_1,yarrc_1)):
-for i,[xa,ya,xc,yc] in enumerate(zip(xloc_a,yloc_a,xloc_c,yloc_c)):
-    print(i,xa,ya,-178.3,xc,yc,-178.3)
-    adisk = gmsh.model.occ.addDisk(xa,ya,0,wire_radius_a,wire_radius_a)
-    cdisk = gmsh.model.occ.addDisk(xc,yc,0,wire_radius_c,wire_radius_c)
+for i,[xg,yg,xa,ya,xc,yc] in enumerate(zip(xloc_g,yloc_g,xloc_a,yloc_a,xloc_c,yloc_c)):
+    print(i,xg,yg,-178.3,xa,ya,-178.3,xc,yc,-178.3)
+    gdisk = gmsh.model.occ.addDisk(xg,yg,0,wire_radius,wire_radius)
+    adisk = gmsh.model.occ.addDisk(xa,ya,0,wire_radius,wire_radius)
+    cdisk = gmsh.model.occ.addDisk(xc,yc,0,wire_radius,wire_radius)
+    guard_wires.append(gdisk)
     anode_wires.append(adisk)
     cathode_wires.append(cdisk)
 
@@ -80,6 +102,10 @@ gmsh.option.setNumber("Geometry.Tolerance", 1e-6)
 gmsh.option.setNumber("Geometry.OCCFixDegenerated", 1)
 gmsh.model.occ.synchronize()
 
+gwire_surfs = []
+for w in guard_wires:
+    gwire_surfs += [s[1] for s in gmsh.model.getBoundary([(2,w)], oriented=False) if s[0] == 1]
+
 awire_surfs = []
 for w in anode_wires:
     awire_surfs += [s[1] for s in gmsh.model.getBoundary([(2,w)], oriented=False) if s[0] == 1]
@@ -87,7 +113,7 @@ for w in anode_wires:
 cwire_surfs = []
 for w in cathode_wires:
     cwire_surfs += [s[1] for s in gmsh.model.getBoundary([(2,w)], oriented=False) if s[0] == 1]
-gmsh.model.mesh.embed(1,cwire_surfs+awire_surfs,2,anasen_barrel)
+gmsh.model.mesh.embed(1,cwire_surfs+awire_surfs+gwire_surfs,2,anasen_barrel)
 
 for s in gmsh.model.getBoundary([(2,w)],oriented=False):
     if s[0] == 1:
@@ -95,7 +121,7 @@ for s in gmsh.model.getBoundary([(2,w)],oriented=False):
 
 
 f1 = gmsh.model.mesh.field.add("Distance")
-gmsh.model.mesh.field.setNumbers(f1,"CurvesList",cwire_surfs+awire_surfs)
+gmsh.model.mesh.field.setNumbers(f1,"CurvesList",cwire_surfs+awire_surfs+gwire_surfs)
 
 f2 = gmsh.model.mesh.field.add("Threshold")
 gmsh.model.mesh.field.setNumber(f2,"InField",f1)
@@ -106,11 +132,14 @@ gmsh.model.mesh.field.setNumber(f2,"DistMax",20)
 
 gmsh.model.mesh.field.setAsBackgroundMesh(f2)
 
-gmsh.model.addPhysicalGroup(1, awire_surfs, tag=10)
-gmsh.model.setPhysicalName(1,10,"anode_wires")
+gmsh.model.addPhysicalGroup(1, gwire_surfs, tag=10)
+gmsh.model.setPhysicalName(1,10,"guard_wires")
 
-gmsh.model.addPhysicalGroup(1, cwire_surfs, tag=20)
-gmsh.model.setPhysicalName(1,20,"cathode_wires")
+gmsh.model.addPhysicalGroup(1, awire_surfs, tag=20)
+gmsh.model.setPhysicalName(1,20,"anode_wires")
+
+gmsh.model.addPhysicalGroup(1, cwire_surfs, tag=30)
+gmsh.model.setPhysicalName(1,30,"cathode_wires")
 
 #gmsh.model.addPhysicalGroup(1, [anasen_bdry], tag=30)
 #gmsh.model.setPhysicalName(1,30,"barrel_boundary")