add test for new polygons code (off until v0.4.1 release)

Author: gaelforget

test/runtests.jl

@@ -1,6 +1,7 @@
 using Test, Documenter, Suppressor, MeshArrays, CairoMakie
 import DataDeps, JLD2, Shapefile, GeoJSON, Proj, GeometryOps
-
+#import MITgcm
+    
 MeshArrays.Dataset("GRID_LL360")
 MeshArrays.Dataset("GRID_LLC90")
 MeshArrays.Dataset("GRID_LLC270")
@@ -349,12 +350,22 @@ end
     plot_examples(:baseproj,proj,lon0,pol=pol_shp)
 end
 
-@testset "polygon operattions" begin
+if false
+#@testset "polygon operations" begin
     name,rule=MeshArrays.within_pol(pol_json; ID=11)
     rule_vec = (x,y) -> rule.(x,y)
 
     np=10000; lo=-180 .+360*rand(np); la=-90 .+180*rand(np);
     sum(rule_vec(lo,la))
+
+    path_MITgcm=MITgcm.getdata("mitgcmsmallverif")
+    path_grid=joinpath(path_MITgcm,"MITgcm","verification","tutorial_held_suarez_cs","input")
+    pols,pols3D=MeshArrays.Polygons.polygons_demo(path_grid)
+    Depth=GridLoadVar("Depth",GridSpec(ID=:CS32))
+
+    MeshArrays.plot_examples(:polygons_plot,pols,color=Depth)
+    MeshArrays.plot_examples(:polygons_plot_dev1,pols,pols3D,sphere_view=true)
+    MeshArrays.plot_examples(:polygons_plot_dev1,pols,pols3D,sphere_view=false)
 end
 
 @testset "doctests" begin