Use numpy rad2deg

Author: cbegeman

polaris/tasks/ocean/customizable_viz/viz_horiz_field.py

@@ -36,8 +36,8 @@ def run(self):  # noqa:C901
         max_latitude = section.getfloat('max_latitude')
         min_longitude = section.getfloat('min_longitude')
         max_longitude = section.getfloat('max_longitude')
-        lat_cell = ds_mesh['latCell'] * 180.0 / np.pi
-        lon_cell = ds_mesh['lonCell'] * 180.0 / np.pi
+        lat_cell = np.rad2deg(ds_mesh['latCell'])
+        lon_cell = np.rad2deg(ds_mesh['lonCell'])
         if min_longitude < 0.0 and lon_cell.min().values > 0.0:
             max_longitude_copy = max_longitude
             max_longitude = 360.0 - min_longitude

polaris/viz/spherical.py

@@ -148,10 +148,10 @@ def plot_global_mpas_field(
     )
 
     if enforce_aspect_ratio:
-        min_latitude = mesh_ds.latCell.min().values * 180.0 / np.pi
-        max_latitude = mesh_ds.latCell.max().values * 180.0 / np.pi
-        min_longitude = mesh_ds.lonCell.min().values * 180.0 / np.pi
-        max_longitude = mesh_ds.lonCell.max().values * 180.0 / np.pi
+        min_latitude = np.rad2deg(mesh_ds.latCell.min().values)
+        max_latitude = np.rad2deg(mesh_ds.latCell.max().values)
+        min_longitude = np.rad2deg(mesh_ds.lonCell.min().values)
+        max_longitude = np.rad2deg(mesh_ds.lonCell.max().values)
         geod = Geodesic()
         x_distance = geod.inverse(
             [min_longitude, min_latitude], [max_longitude, min_latitude]