Merge pull request #1005 from NCAR/roms-ssh-required

ROMS Rutgers assert sea surface height is in the state for model_interpolate

Author: hkershaw-brown

CHANGELOG.rst

@@ -22,6 +22,18 @@ individual files.
 
 The changes are now listed with the most recent at the top.
 
+**December 3 2025 :: BAMS article. Tag 11.19.1**
+
+Documentation updates:
+
+  - How to cite DART
+  - Bulletin of the American Meteorological Society 2025 article
+
+Bug-fix:
+
+  - Assert sea surface height is present in the state for ROMS Rutgers 
+    model_interpolate
+
 **November 13 2025 :: HF Radar converter. Tag 11.19.0**
 
 New features:

conf.py

@@ -21,7 +21,7 @@
 author = 'Data Assimilation Research Section'
 
 # The full version, including alpha/beta/rc tags
-release = '11.19.0'
+release = '11.19.1'
 root_doc = 'index'
 
 # -- General configuration ---------------------------------------------------

models/ROMS_rutgers/model_mod.f90

@@ -402,6 +402,7 @@ subroutine model_interpolate(state_handle, ens_size, location, qty, expected_obs
 ! because we need to compute the depth using it.
 ! Recall that SSH can be both -ve and +ve
 sshid = get_varid_from_kind(domid, QTY_SEA_SURFACE_HEIGHT)
+if (sshid < 0 ) call error_handler(E_ERR, 'model_interpolate', 'requires sea surface height in the state')
 do i = 1, Nc
    dartidx       = get_dart_vector_index(lon_c(i), lat_c(i), 1, domid, sshid)
    SSHcorn(i, :) = get_state(dartidx, state_handle)