Merge pull request #74 from ACCESS-NRI/pearseb-patch-chl

Update generic_WOMBATlite.F90

Author: pearseb

generic_tracers/generic_WOMBATlite.F90

@@ -2638,7 +2638,7 @@ subroutine generic_WOMBATlite_update_from_source(tracer_list, Temp, Salt,  &
       phi = wombat%radmld(i,j,1) / (wombat%radmld(i,j,1) + wombat%chlkWm2)
 
       pchl_pisl = phy_pisl / ( wombat%phy_mumax(i,j,k) * 86400.0 * & 
-                  (1. - min(wombat%phy_lnit(i,j,k), wombat%phy_lfer(i,j,k))) + epsi )
+                  (1. - max(0.33, min(wombat%phy_lnit(i,j,k), wombat%phy_lfer(i,j,k)))) + epsi )
       wombat%pchl_lpar(i,j,k) = (1. - exp(-pchl_pisl * wombat%radmld(i,j,k)))
       pchl_mumin = wombat%phyminqc * wombat%phy_mu(i,j,k) * biophy * 12.0   ![mg/m3/s] 
       pchl_muopt = wombat%phyoptqc * wombat%phy_mu(i,j,k) * biophy * 12.0   ![mg/m3/s]