From 833db6a612fd8ec1d78e92e595fdb308ccba7e05 Mon Sep 17 00:00:00 2001 From: Charles Stock Date: Thu, 4 Apr 2024 11:31:24 -0400 Subject: [PATCH] added felim to diazo iron limitation to address issue #27 --- generic_tracers/generic_COBALT.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/generic_tracers/generic_COBALT.F90 b/generic_tracers/generic_COBALT.F90 index c3eb4c7..da36db4 100644 --- a/generic_tracers/generic_COBALT.F90 +++ b/generic_tracers/generic_COBALT.F90 @@ -8276,7 +8276,7 @@ subroutine generic_COBALT_update_from_source(tracer_list,Temp,Salt,rho_dzt,dzt,h do k = 1, nk ; do j = jsc, jec ; do i = isc, iec !{ n=DIAZO phyto(n)%liebig_lim(i,j,k) = phyto(n)%o2lim(i,j,k)* & - min(phyto(n)%po4lim(i,j,k), phyto(n)%def_fe(i,j,k)) + min(phyto(n)%po4lim(i,j,k), max(phyto(n)%def_fe(i,j,k),phyto(n)%felim(i,j,k))) do n= 2, NUM_PHYTO !{ phyto(n)%liebig_lim(i,j,k) = min(phyto(n)%no3lim(i,j,k)+phyto(n)%nh4lim(i,j,k),& phyto(n)%po4lim(i,j,k), max(phyto(n)%def_fe(i,j,k),phyto(n)%felim(i,j,k)))