diff --git a/bgcval2/bgcvaltools/bv2tools.py b/bgcval2/bgcvaltools/bv2tools.py
index 1fd050ef..a5054f58 100644
--- a/bgcval2/bgcvaltools/bv2tools.py
+++ b/bgcval2/bgcvaltools/bv2tools.py
@@ -953,10 +953,10 @@ def robinPlotQuad(lons,
 
             if scatter:
                 if doLogs[i] and spl in [221, 222]:
-                    rbmi = np.int(np.log10(rbmi))
+                    rbmi = np.int32(np.log10(rbmi))
                     rbma = np.log10(rbma)
-                    if rbma > np.int(rbma): rbma += 1
-                    rbma = np.int(rbma)
+                    if rbma > np.int32(rbma): rbma += 1
+                    rbma = np.int32(rbma)
 
                 if doLogs[i]:
                     ims.append(bms[i].scatter(
@@ -1199,10 +1199,10 @@ def HovPlotQuad(
         axs.append(fig.add_subplot(spl))
         if scatter:
             if doLogs[i] and spl in [221, 222]:
-                rbmi = np.int(np.log10(rbmi))
+                rbmi = np.int32(np.log10(rbmi))
                 rbma = np.log10(rbma)
-                if rbma > np.int(rbma): rbma += 1
-                rbma = np.int(rbma)
+                if rbma > np.int32(rbma): rbma += 1
+                rbma = np.int32(rbma)
 
             if doLogs[i]:
                 ims.append(
@@ -1463,10 +1463,10 @@ def ArcticTransectPlotQuad(
         axs.append(fig.add_subplot(spl))
         if scatter:
             if doLogs[i] and spl in [221, 222]:
-                rbmi = np.int(np.log10(rbmi))
+                rbmi = np.int32(np.log10(rbmi))
                 rbma = np.log10(rbma)
-                if rbma > np.int(rbma): rbma += 1
-                rbma = np.int(rbma)
+                if rbma > np.int32(rbma): rbma += 1
+                rbma = np.int32(rbma)
 
             if doLogs[i]:
                 ims.append(
@@ -1989,7 +1989,7 @@ def histsPlot(datax,
 
     if logx:
         maxd = np.ma.power(10.,
-                           np.int(np.ma.max(np.ma.abs(np.ma.log10(d))) + 1))
+                           np.int32(np.ma.max(np.ma.abs(np.ma.log10(d))) + 1))
         print(maxd, 1 / maxd)
         n, bins, patchesx = pyplot.hist(d,
                                         histtype='stepfilled',
diff --git a/bgcval2/functions/applyLandMask.py b/bgcval2/functions/applyLandMask.py
index 311a0db7..a44458a4 100644
--- a/bgcval2/functions/applyLandMask.py
+++ b/bgcval2/functions/applyLandMask.py
@@ -40,7 +40,7 @@
 def loadDataMask(areafile, maskname):
     global tmask
     nc = dataset(areafile, 'r')        
-    tmask[(areafile, maskname)] = np.array(nc.variables[maskname][:].squeeze(), dtype=np.bool)
+    tmask[(areafile, maskname)] = np.array(nc.variables[maskname][:].squeeze(), dtype=bool)
     nc.close()
     return tmask[(areafile, maskname)]
 
diff --git a/bgcval2/timeseries/timeseriesTools.py b/bgcval2/timeseries/timeseriesTools.py
index ef99eff6..38ff1ca2 100644
--- a/bgcval2/timeseries/timeseriesTools.py
+++ b/bgcval2/timeseries/timeseriesTools.py
@@ -227,7 +227,7 @@ def getHorizontalSlice(nc, coords, details, layer, data=''):
     if type(layer) in [
             type(0),
             np.int64,
-            np.int,
+            int,
     ]:
         k = layer
         try: