add logic to remove incrments from restart if outside IAU window

NOAA-GFDL · Jan 28, 2022 · 78b8d79 · 78b8d79
1 parent a5d9065
commit 78b8d79
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/src/ocean_data_assim/MOM_oda_incupd.F90 b/src/ocean_data_assim/MOM_oda_incupd.F90
@@ -554,6 +554,16 @@ subroutine apply_oda_incupd(h, tv, u, v, dt, G, GV, US, CS)
   if (CS%ncount >= CS%nstep_incupd) then
     if (is_root_pe()) call MOM_error(NOTE,"ended updating fields with increments. ")
     return
+    if (CS%ncount == CS%nstep_incupd) then
+       call register_restart_field_as_obsolete("oda_incupd_ncount", "none", CS)
+       call register_restart_field_as_obsolete("T_inc", "none", CS)
+       call register_restart_field_as_obsolete("S_inc", "none", CS)
+       call register_restart_field_as_obsolete("h_obs", "none", CS)
+       if (CS%uv_inc) then
+          call register_restart_field_as_obsolete("u_inc", "none", CS)
+          call register_restart_field_as_obsolete("v_inc", "none", CS)
+       endif
+     endif
   endif !ncount>CS%nstep_incupd
 
   ! update counter