Index: /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py	(revision 25345)
+++ /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py	(revision 25346)
@@ -1,7 +1,8 @@
+import shlex
 import subprocess
 
 
 def epsg2proj(epsg): #{{{
-    '''
+    """
     FUNCTION EPSG2PROJ - uses gdalsrsinfo to provide PROJ.4 compatible string 
     from EPSG code
@@ -18,5 +19,5 @@
         - In case of Python 2, except socket.timeout: https://docs.python.org/3/library/socket.html?highlight=socket%20timeout#socket.timeout
         - In case of Python 3, except TimeoutExpired: https://docs.python.org/3/library/subprocess.html#subprocess.SubprocessError
-    '''
+    """
 
     #First, get GDAL version
Index: /issm/trunk-jpl/src/m/io/loadmodel.py
===================================================================
--- /issm/trunk-jpl/src/m/io/loadmodel.py	(revision 25345)
+++ /issm/trunk-jpl/src/m/io/loadmodel.py	(revision 25346)
@@ -33,5 +33,5 @@
     #       try:
     #recover model on file and name it md
-    struc = loadvars(path, OL=onlylast)
+    struc = loadvars(path, onlylast=onlylast)
     name = [key for key in list(struc.keys())]
     if len(name) > 1:
Index: /issm/trunk-jpl/src/m/io/loadvars.py
===================================================================
--- /issm/trunk-jpl/src/m/io/loadvars.py	(revision 25345)
+++ /issm/trunk-jpl/src/m/io/loadvars.py	(revision 25346)
@@ -15,6 +15,6 @@
 
 
-def loadvars(*args, OL):
-    '''
+def loadvars(*args, **kwargs):
+    """
     LOADVARS - function to load variables from a file.
 
@@ -30,5 +30,5 @@
         nvdict = loadvars('shelve.dat', {'a':None, 'b':None})
         nvdict = loadvars('shelve.dat')
-    '''
+    """
 
     filename = ''
@@ -55,4 +55,10 @@
     else:
         raise TypeError("Unrecognized input arguments.")
+
+    onlylast = False
+
+    for key, value in kwargs.items():
+        if key == 'onlylast':
+            onlylast = value
 
     if whichdb(filename):   #We used python pickle for the save
@@ -110,5 +116,5 @@
                             nvdict['md'].__dict__[classtree[mod][0]].__dict__[classtree[mod][1]] = [OrderedDict() for i in range(max(1, len(curclass.groups)))]
                         else:
-                            if OL:   #we load only the last result to save on time and memory
+                            if onlylast:   #we load only the last result to save on time and memory
                                 nvdict['md'].__dict__[classtree[mod][0]].__dict__[classtree[mod][1]] = [getattr(__import__(listtype), listtype)()]
                                 Tree = nvdict['md'].__dict__[classtree[mod][0]].__dict__[classtree[mod][1]]
@@ -126,5 +132,5 @@
                         if AllHaveTime:
                             #Time dimension is in all the variables so we take that as stepnumber for the results
-                            if OL:   #we load only the last result to save on time and memory
+                            if onlylast:   #we load only the last result to save on time and memory
                                 nvdict['md'].__dict__[classtree[mod][0]].__dict__[classtree[mod][1]] = [getattr(__import__(listtype), listtype)()]
                                 Tree = nvdict['md'].__dict__[classtree[mod][0]].__dict__[classtree[mod][1]]
@@ -166,5 +172,5 @@
                         NewFormat = 'Time' in NCFile.dimensions
                         if type(Tree) == list and NewFormat:
-                            if OL:
+                            if onlylast:
                                 if vardim == 1:
                                     Tree[0].__dict__[str(var)] = varval[-1].data