Set oceanFracObserved to 0 where MALI's bed elevation is > 0

Author: cbegeman

compass/ocean/tests/global_ocean/mesh/remap_mali_topography/__init__.py

@@ -293,7 +293,9 @@ def _combine_topo(self):
 
         alpha = ds_out.maliFrac
         # NOTE: MALI's ocean fraction is already scaled by the MALI fraction
-        ds_out['oceanFracObserved'] = (
+        ds_out['oceanFracObserved'] = xr.where(
+            ds_mali.bed_elevation > 0.,
+            0.,
             ds_mali.oceanFrac +
             (1.0 - alpha) * ds_bedmachine.oceanFracObserved)