cmor | Climate Model Output Rewriter | Code Quality library

You will need to explicitly rename dimension names to match, e.g., hs = hs.rename({'lat': 'latitude', 'longitude': 'longitude'}). You might also need to reindex with nearest-neighbor indexing, if the coordinate labels don't match exactly, e.g., hs.reindex_like(ws10, method='nearest', tolerance=0.01).

Or, less safely, you could strip out the coordinate labels from the second argument, and just pass in a unlabeled array instead, e.g., ws10.where(hs.data > 0). But I don't recommend this option, because nothing guarantees the consistency of the metadata.