Use default_num_points when available

`ClimaDiagnostics` introduced a function to compute the default number
of points for the NetCDFWriter. This commit uses this new function when
available.

Author: Sbozzolo

src/callbacks/get_callbacks.jl

@@ -33,26 +33,7 @@ function get_diagnostics(parsed_args, atmos_model, Y, p, sim_info, output_dir)
         num_netcdf_points =
             tuple(parsed_args["netcdf_interpolation_num_points"]...)
     else
-        # Estimate the number of points we need to cover the entire domain
-        # ncolumns is the number of local columns
-        tot_num_columns =
-            ClimaComms.nprocs(context) * Fields.ncolumns(axes(Y.c))
-        if parsed_args["config"] == "plane"
-            num1, num2 = tot_num_columns, 0
-        elseif parsed_args["config"] == "sphere"
-            num2 = round(Int, sqrt(tot_num_columns / 2))
-            num1 = 2num2
-        elseif parsed_args["config"] == "box"
-            num2 = round(Int, sqrt(tot_num_columns))
-            num1 = num2
-        elseif parsed_args["config"] == "column"
-            # We need at least two points horizontally because our column is
-            # actually a box
-            num1, num2 = 2, 2
-        else
-            error("Uncaught case")
-        end
-        num_netcdf_points = (num1, num2, Spaces.nlevels(axes(Y.c)))
+        num_netcdf_points = default_netcdf_points(axes(Y.c), parsed_args)
     end
 
     z_sampling_method =

src/compat.jl

@@ -1,6 +1,7 @@
 import ClimaCore
 import ClimaUtilities
 import ClimaCore: Domains, Spaces, Topologies
+import ClimaDiagnostics
 import RRTMGP
 
 # To allow for backwards compatibility of ClimaCore:
@@ -96,3 +97,32 @@ else
     WallTimeInfo = ClimaUtilities.OnlineLogging.WallTimeInfo
     report_walltime = ClimaUtilities.OnlineLogging.report_walltime
 end
+
+if pkgversion(ClimaDiagnostics) < v"0.2.13"
+    function default_netcdf_points(space, parsed_args)
+        # Estimate the number of points we need to cover the entire domain
+        # ncolumns is the number of local columns
+        tot_num_columns =
+            ClimaComms.nprocs(ClimaComms.context(space)) * Fields.ncolumns(space)
+        if parsed_args["config"] == "plane"
+            num1, num2 = tot_num_columns, 0
+        elseif parsed_args["config"] == "sphere"
+            num2 = round(Int, sqrt(tot_num_columns / 2))
+            num1 = 2num2
+        elseif parsed_args["config"] == "box"
+            num2 = round(Int, sqrt(tot_num_columns))
+            num1 = num2
+        elseif parsed_args["config"] == "column"
+            # We need at least two points horizontally because our column is
+            # actually a box
+            num1, num2 = 2, 2
+        else
+            error("Uncaught case")
+        end
+        return (num1, num2, Spaces.nlevels(space))
+    end
+else
+    function default_netcdf_points(space, _)
+        return ClimaDiagnostics.Writers.default_num_points(space)
+    end
+end