From 4f2f4b37e01a4dac7c127cc0d1ded577c5ec3ab2 Mon Sep 17 00:00:00 2001 From: Stefano Piani Date: Thu, 19 Mar 2026 01:45:14 +0100 Subject: [PATCH] Now nudging indices contain only 0 and 1 --- src/mitgcm_inputs/ob_indices/ob_indices.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/src/mitgcm_inputs/ob_indices/ob_indices.py b/src/mitgcm_inputs/ob_indices/ob_indices.py index b33692e..0922fb3 100644 --- a/src/mitgcm_inputs/ob_indices/ob_indices.py +++ b/src/mitgcm_inputs/ob_indices/ob_indices.py @@ -178,6 +178,9 @@ def generate_ob_indices( side, ) + # In the nudging indices, the only allowed indices are 0 and 1 + nudging_cells = np.where(sponge_cells == closed_cell, 0, 1) + open_bc_cells = np.copy(sponge_cells) for river in rivers_positions: if river["model"] != "stem_flux": @@ -205,7 +208,7 @@ def generate_ob_indices( ob_indices += "\n" ob_sponge += f"nudgOB_{direction_char}{side}=" - for i in rewrite_in_mitgcm_format(sponge_cells): + for i in rewrite_in_mitgcm_format(nudging_cells): ob_sponge += f"{i[0]}*{i[1]}," ob_sponge += "\n"