--- MITgcm_contrib/bling/pkg/bling_dvm.F	2016/05/19 16:30:00	1.3
+++ MITgcm_contrib/bling/pkg/bling_dvm.F	2016/05/19 20:00:53	1.4
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/bling/pkg/bling_dvm.F,v 1.3 2016/05/19 16:30:00 mmazloff Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/bling/pkg/bling_dvm.F,v 1.4 2016/05/19 20:00:53 mmazloff Exp $
 C $Name:  $
 
 #include "BLING_OPTIONS.h"
@@ -340,15 +340,16 @@
 
 #else
 
-        z_dvm_regr = 398. _d 0 - 0.56 _d 0*min(300. _d 0,
-     & max(-10. _d 0,(o2_upper - o2_lower))) - 
-     & 115. _d 0*min(0.85 _d 0,max(-1.80 _d 0,log10(chl(i,j,1,bi,bj))))
-     &  + 0.36 _d 0*min(500. _d 0,max(epsln,mld(i,j))) - 
-     &  2.4 _d 0*min(20. _d 0,max(-3. _d 0, (temp_upper - temp_lower)))
-
-        ! Limit the depth of migration in polar winter.
-        ! Use irr_mem since this is averaged over multiple days, dampening the diurnal cycle.
-        ! Tapers Z_DVM to the minimum when surface irradince is below a given threshold (here 10 W/m2).
+        z_dvm_regr = 398. _d 0 
+     &   - 0.56 _d 0*min(300. _d 0,max(-10. _d 0,(o2_upper - o2_lower)))
+     &   - 115. _d 0*min(0.85 _d 0,max(-1.80 _d 0,
+     &                            log10(max(chl(i,j,1,bi,bj),chl_min))))
+     &   + 0.36 _d 0*min(500. _d 0,max(epsln,mld(i,j)))
+     &   - 2.40 _d 0*min(20. _d 0,max(-3. _d 0,(temp_upper-temp_lower)))
+
+c        ! Limit the depth of migration in polar winter.
+c        ! Use irr_mem since this is averaged over multiple days, dampening the diurnal cycle.
+c        ! Tapers Z_DVM to the minimum when surface irradince is below a given threshold (here 10 W/m2).
 
         if ( irr_mem(i,j,1,bi,bj) .lt. 10. ) then                       
           z_dvm_regr = 150. _d 0 + (z_dvm_regr - 150. _d 0) *