Index: /issm/trunk-jpl/src/m/classes/sealevelmodel.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/sealevelmodel.py	(revision 27529)
+++ /issm/trunk-jpl/src/m/classes/sealevelmodel.py	(revision 27530)
@@ -241,5 +241,5 @@
             self.eltransitions.append(meshintersect3d(xe, ye, ze, xei, yei, zei, 'force', force))
 
-            self.earth.solidearth.transfercount[self.transitions[i]] = self.earth.solidearth/transfercount[self.transitions[i]] + 1
+            self.earth.solidearth.transfercount[self.transitions[i]] = self.earth.solidearth.transfercount[self.transitions[i]] + 1
 
         for i in range(len(self.icecaps)):
Index: /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py	(revision 27529)
+++ /issm/trunk-jpl/src/m/coordsystems/epsg2proj.py	(revision 27530)
@@ -20,23 +20,10 @@
 
     #First, get GDAL version
-    subproc_args = "gdalsrsinfo --version | awk '{print $2}' | cut -d '.' -f1"
-    subproc = subprocess.Popen(subproc_args, shell=True, stdout=subprocess.PIPE, stderr=subprocess.PIPE, universal_newlines=True)
-    outs, errs = subproc.communicate()
-    outs = outs.decode()
-    errs = errs.decode()
-    print(errs != '', errs == '', errs)
-    if errs != '':
-        raise RuntimeError("epsg2proj: call to gdalsrsinfo failed: {}".format(errs))
-
-    version_major = int(outs)
-
-    subproc_args = "gdalsrsinfo epsg:{} | command grep PROJ.4 | tr -d '\n' | sed 's/PROJ.4 : //'".format(epsg)
+    #subproc_args = "gdalsrsinfo --version | awk '{print $2}' | cut -d '.' -f1"
+    subproc_args = "projinfo -o PROJ -q epsg:{}".format(epsg)
     subproc = subprocess.Popen(subproc_args, shell=True, stdout=subprocess.PIPE, stderr=subprocess.PIPE, universal_newlines=True)
     outs, errs = subproc.communicate()
     if errs != '':
-        raise RuntimeError("epsg2proj: call to gdalsrsinfo failed: {}".format(errs))
-
-    if version_major == 1:
-        outs = outs[1:-1]
+        raise RuntimeError("epsg2proj: call to projinfo failed: {}".format(errs))
 
     return outs