Fix a bug with the swredist code in ktherm=2. (#332)

* Bug fix for dt_rhoi_hlyr * Fix bug with swredist in ktherm 2
CICE-Consortium · Aug 19, 2020 · ddec41e · ddec41e
1 parent 4c42a82
commit ddec41e
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 0 deletions.
diff --git a/columnphysics/icepack_therm_mushy.F90 b/columnphysics/icepack_therm_mushy.F90
@@ -219,6 +219,8 @@ subroutine temperature_changes_salinity(dt,                 &
 
     if (sw_redist) then
 
+    dt_rhoi_hlyr = dt / (rhoi*hilyr)
+
     do k = 1, nilyr
 
        Iswabs_tmp = c0 ! all Iswabs is moved into fswsfc

diff --git a/configuration/scripts/options/set_nml.swredist b/configuration/scripts/options/set_nml.swredist
@@ -0,0 +1,7 @@
+    tr_pond_lvl  = .true.
+    tr_pond_topo = .false.
+    tr_pond_cesm = .false.
+    ktherm            = 2
+    sw_redist         = .true.
+    sw_frac           = 0.9d0
+    sw_dtemp          = 0.02d0
diff --git a/configuration/scripts/tests/base_suite.ts b/configuration/scripts/tests/base_suite.ts
@@ -5,6 +5,7 @@ smoke          col     1x1        debug,bgcISPOL
 smoke          col     1x1        debug,bgcNICE
 smoke          col     1x1        debug,bgcsklNICE
 smoke          col     1x1        debug,run1year,thermo1
+smoke          col     1x1        debug,run1year,swredist
 smoke          col     1x1        debug,run1year,swccsm3
 smoke          col     1x1        debug,run1year,alt01
 smoke          col     1x1        debug,run1year,alt02