Merge pull request #49 from MichaelBarryBMT/IN-IP-lims-upper-bound

Upper limit to internal nutrient limit functions
AquaticEcoDynamics · Mar 7, 2022 · c8dc4bd · c8dc4bd
2 parents 13fdd43 + c1324dd
commit c8dc4bd
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/aed_bio_utils.F90 b/src/aed_bio_utils.F90
@@ -336,6 +336,8 @@ FUNCTION phyto_fN(phytos, group, IN, din, don) RESULT(fN)
  ENDIF
 
  IF ( fN < zero_ ) fN = zero_
+ IF ( fN > 1.000 ) fN = 1.000
+
 END FUNCTION phyto_fN
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
@@ -366,7 +368,8 @@ FUNCTION phyto_fP(phytos, group, IP, frp) RESULT(fP)
  (phytos(group)%X_pmax-phytos(group)%X_pmin)
  ENDIF
 
- IF( fP<zero_ ) fP=zero_
+ IF( fP < zero_ ) fP = zero_ 
+ IF( fP > 1.000 ) fP = 1.000 
 
 END FUNCTION phyto_fP
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++