Added catch for unusual case of (hicen_init(n+1) - hicen_init(n))>0 (#…

…337) * Added catch for unusual case of (hicen_init(n+1) - hicen_init(n))>0 * fix conditional Co-authored-by: Andrew Roberts <[email protected]>
CICE-Consortium · Mar 17, 2021 · 7223366 · 7223366
1 parent ee7cf8e
commit 7223366
Showing 1 changed file with 16 additions and 3 deletions.
diff --git a/columnphysics/icepack_therm_itd.F90 b/columnphysics/icepack_therm_itd.F90
@@ -305,11 +305,24 @@ subroutine linear_itd (ncat,        hin_max,     &
 
          if (hicen_init(n)   > puny .and. &
              hicen_init(n+1) > puny) then
-             ! interpolate between adjacent category growth rates
-             slope = (dhicen(n+1) - dhicen(n)) / &
+
+            if ((hicen_init(n+1) - hicen_init(n))>0) then
+
+              ! interpolate between adjacent category growth rates
+              slope = (dhicen(n+1) - dhicen(n)) / &
                  (hicen_init(n+1) - hicen_init(n))
-             hbnew(n) = hin_max(n) + dhicen(n) &
+              hbnew(n) = hin_max(n) + dhicen(n) &
                       + slope * (hin_max(n) - hicen_init(n))
+
+            else
+
+              write(warnstr,*) subname, &
+                 'ITD Thermodynamics: hicen_init(n+1) <= hicen_init(n)'
+              call icepack_warnings_setabort(.true.)
+              call icepack_warnings_add(warnstr)
+
+            endif
+
          elseif (hicen_init(n) > puny) then ! hicen_init(n+1)=0
              hbnew(n) = hin_max(n) + dhicen(n)
          elseif (hicen_init(n+1) > puny) then ! hicen_init(n)=0