Skip to content

Commit 440fea5

Browse files
Disable checkIniTemp/checkIniSalt for repeat-year restart ICs
The model produces theta=0 and salt=0 at 63 wet cells on the northern OBC boundary (j=423, mid-depth). These are legitimate model state, not uninitialized data, so the INI_THETA check must be disabled for pickup-derived initial conditions. Also remove the NaN/Inf fill from pickup_to_init — the big-endian fix eliminated the overflow, so the fill was no longer needed. Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
1 parent bd01a8d commit 440fea5

2 files changed

Lines changed: 3 additions & 5 deletions

File tree

simulations/glorysv12-curvilinear/input/data

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -73,6 +73,6 @@
7373
hydrogSaltFile='S.init.bin',
7474
uVelInitFile='U.init.bin',
7575
vVelInitFile='V.init.bin',
76-
# checkIniTemp=.FALSE.
77-
# checkIniSalt=.FALSE.
76+
checkIniTemp=.FALSE.,
77+
checkIniSalt=.FALSE.,
7878
&end

spectre_utils/pickup_to_init.py

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -123,9 +123,7 @@ def pickup_to_init(pickup_prefix: str, output_dir: str, nx: int, ny: int, nr: in
123123
init_path = out / init_name
124124
if init_path.is_symlink() or init_path.exists():
125125
init_path.unlink()
126-
data32 = data.astype(">f4")
127-
data32[~np.isfinite(data32)] = 0.0
128-
data32.tofile(init_path)
126+
data.astype(">f4").tofile(init_path)
129127
size_mb = init_path.stat().st_size / 1e6
130128
print(f" Wrote {init_path} ({size_mb:.1f} MB)")
131129

0 commit comments

Comments
 (0)