protect merge of meltsliq when tr_snow is off (#376)

CICE-Consortium · Nov 12, 2021 · 152bd70 · 152bd70
1 parent 27b872c
commit 152bd70
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/columnphysics/icepack_flux.F90 b/columnphysics/icepack_flux.F90
@@ -12,7 +12,7 @@ module icepack_flux
       use icepack_parameters, only: c1, emissivity
       use icepack_warnings, only: warnstr, icepack_warnings_add
       use icepack_warnings, only: icepack_warnings_setabort, icepack_warnings_aborted
-      use icepack_tracers, only: tr_iso
+      use icepack_tracers, only: tr_iso, tr_snow
 
       implicit none
       private
@@ -217,7 +217,9 @@ subroutine merge_fluxes (aicen,                &
       meltt     = meltt     + melttn    * aicen
       meltb     = meltb     + meltbn    * aicen
       melts     = melts     + meltsn    * aicen
-      meltsliq  = meltsliq  + meltsliqn * aicen
+      if (tr_snow) then
+         meltsliq  = meltsliq  + meltsliqn * aicen
+      endif
       dsnow     = dsnow     + dsnown    * aicen
       congel    = congel    + congeln   * aicen
       snoice    = snoice    + snoicen   * aicen