Remove non-finite fill values for pf.ice_ext() calculation

With the fill value change from 0 to nan in fesom2.7 the pf.ice_ext() did return incorrect values.

Author: FinnHeu

pyfesom2/diagnostics.py

@@ -76,18 +76,21 @@ def ice_ext(data, mesh, hemisphere="N", threshhold=0.15, attrs={}):
         data = data.where(data < threshhold, 1)
         data = data.where(data > threshhold)
 
-        ext = (data[:, hemis_mask] * mesh.lump2[hemis_mask]).sum(axis=1)
+        finite_mask = data.squeeze().values
+        
+        ext = (data[:, hemis_mask & finite_mask] * mesh.lump2[hemis_mask & finite_mask]).sum(axis=1)
         da = xr.DataArray(
             ext, dims=["time"], coords={"time": data.time}, name=varname, attrs=attrs
         )
         return da
 
     else:
         logger.debug(data)
+        finite_mask = data.squeeze()
         i, j = np.where(data < 0.15)
         data[:] = 1
         data[i, j] = 0
-        ext = (data[:, hemis_mask] * mesh.lump2[hemis_mask]).sum(axis=1)
+        ext = (data[:, hemis_mask & finite_mask] * mesh.lump2[hemis_mask & finite_mask]).sum(axis=1)
         return ext