OM2-01 fix for grid topology

Author: briochemc

Project.toml

@@ -1,7 +1,7 @@
 name = "OceanTransportMatrixBuilder"
 uuid = "c2b4a04e-6049-4fc4-aa6a-5508a29a1e1c"
 authors = ["Benoit Pasquier <[email protected]> and contributors"]
-version = "0.5.3"
+version = "0.5.4"
 
 [deps]
 Distances = "b4f34e82-e78d-54a5-968a-f98e89d6e8f7"

src/gridtopology.jl

@@ -31,7 +31,7 @@ function getgridtopology(lon_vertices, lat_vertices, lev)
     if all(NPlat .== 90)
         return BipolarGridTopology(nx,ny,nz)
     # Otherwise check if the north pole is split in two
-    elseif NPlon == rot180(NPlon) && NPlat == rot180(NPlat)
+    elseif isapprox(NPlon, rot180(NPlon)) && isapprox(NPlat, rot180(NPlat))
         return TripolarGridTopology(nx,ny,nz)
     else
         return UnknownGridTopology(nx,ny,nz)