diff --git a/LoadCase/LoadsAndBCs.py b/LoadCase/LoadsAndBCs.py
index dc01fe5d944bbf07e540a6e8f88d2376ae35768b..6fad3122379732cc54fbb4a0b21741465597e1d3 100644
--- a/LoadCase/LoadsAndBCs.py
+++ b/LoadCase/LoadsAndBCs.py
@@ -32,10 +32,9 @@ def create_load(**kwargs):
     LoadVal = kwargs['CompressiveLoad']
     for i in range(len(ipMTnames)):
         if i % 2 == 0:  # Right pole
-            # TODO: Resolve issue with selecting the proper end of ipMT to apply loading
             v = a.instances[ipMTnames[i]].vertices
             loc = v[-1].pointOn
-            verts = v.findAt((loc, ), )
+            verts = v.findAt(loc, )
             region = a.Set(vertices=verts, name='ipMTtip-' + str(i))
             mdb.models[modelname].ConcentratedForce(name='Load-' + str(i),
                                                     createStepName=StepName,
@@ -47,7 +46,7 @@ def create_load(**kwargs):
         else:  # Left pole
             v = a.instances[ipMTnames[i]].vertices
             loc = v[-1].pointOn
-            verts = v.findAt((loc,), )
+            verts = v.findAt(loc, )
             region = a.Set(vertices=verts, name='ipMTtip-' + str(i))
             mdb.models[modelname].ConcentratedForce(name='Load-' + str(i),
                                                     createStepName=StepName,
diff --git a/LoadCase/Step.py b/LoadCase/Step.py
index a125dcd88dc37abd1af221dce2e2690d972a8346..a3cde5ff6f6ceb8b791107469aaf51af9d57e28d 100644
--- a/LoadCase/Step.py
+++ b/LoadCase/Step.py
@@ -1,5 +1,3 @@
-# -*- coding: mbcs -*-
-# Do not delete the following import lines
 from abaqus import *
 from abaqusConstants import *
 import __main__
@@ -26,10 +24,11 @@ def Step(**kwargs):
     # Define a standard buckling step with Subspace solver
     modelname = kwargs['modelname']
     StepName = kwargs['StepName']
+    NumberOfEigs = kwargs['NumberOfEigs']
     mdb.models[modelname].BuckleStep(name=StepName,
                                      previous='Initial',
-                                     numEigen=20,
-                                     vectors=40,
+                                     numEigen=NumberOfEigs,
+                                     vectors=2*NumberOfEigs,
                                      maxIterations=300)
     # Create a list of output variables
     mdb.models['test'].FieldOutputRequest(