Add SOPAC, SOIND, and SOATL functions

Author: briochemc

Project.toml

@@ -1,7 +1,7 @@
 name = "OceanBasins"
 uuid = "d1bb7020-b2be-4340-9d18-d24ca645bddb"
 authors = ["Benoit Pasquier <briochemc@gmail.com>"]
-version = "0.1.3"
+version = "0.1.4"
 
 [deps]
 DataDeps = "124859b0-ceae-595e-8997-d05f6a7a8dfe"

src/OceanBasins.jl

@@ -57,9 +57,13 @@ isantarctic(P, oceans) = P.lat ≤ -40
 isatlantic2(P, oceans) = isatlantic(P, oceans) && !isantarctic(P, oceans)
 ispacific2(P, oceans) = ispacific(P, oceans) && !isantarctic(P, oceans)
 isindian2(P, oceans) = isindian(P, oceans) && !isantarctic(P, oceans)
+convert_65S_to_90S(P::Point2D{T}) where T = Point2D(P.lat < -65 ? T(-65) : P.lat, P.lon)
+isSOatlantic(P, oceans) = isatlantic(convert_65S_to_90S(P), oceans) && isantarctic(P, oceans)
+isSOpacific(P, oceans) = ispacific(convert_65S_to_90S(P), oceans) && isantarctic(P, oceans)
+isSOindian(P, oceans) = isindian(convert_65S_to_90S(P), oceans) && isantarctic(P, oceans)
 isindonesian(P, oceans) = foldl(|, P ∈ oceans[i] for i in 11:37) # Indo throughflow (TODO fix name?)
 for ocn in (:pacific, :atlantic, :indian, :arctic, :mediterranean, :antarctic,
-            :pacific2, :atlantic2, :indian2, :indonesian)
+            :pacific2, :atlantic2, :indian2, :indonesian, :SOatlantic, :SOpacific, :SOindian)
     f = Symbol(:is, ocn)
     @eval begin
         """

test/runtests.jl

@@ -8,6 +8,12 @@ const OCEANS = oceanpolygons()
     @test isarctic(89, 0, OCEANS)
     @test isatlantic(0, -30, OCEANS)
     @test isindian(0, 90, OCEANS)
+    @test ispacific2(0, -160, OCEANS)
+    @test isatlantic2(0, -30, OCEANS)
+    @test isindian2(0, 90, OCEANS)
+    @test isSOpacific(-80, -160, OCEANS)
+    @test isSOatlantic(-80, -30, OCEANS)
+    @test isSOindian(-80, 90, OCEANS)
     @test isantarctic(-40, 0, OCEANS)
     @test ismediterranean(31, 15, OCEANS)
     @test isindonesian(0, 120, OCEANS)