Avoid chl=0 in lookup_ohlmann_opacity()

previous commit only avoiding log10(0) in lookup_ohlmann_swpen() but we have another log10(chl) in the opacity function
NCAR · Oct 1, 2024 · 8fb8aaf · 8fb8aaf
1 parent 5726d94
commit 8fb8aaf
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/parameterizations/vertical/MOM_opacity.F90 b/src/parameterizations/vertical/MOM_opacity.F90
@@ -1377,7 +1377,11 @@ function lookup_ohlmann_opacity(chl,optics) result(B)
   integer :: n
 
   ! Make sure we are in the table
-  log10chl = max(optics%log10chl_min,min(log10(chl),optics%log10chl_max))
+  if (chl > optics%chl_min) then
+    log10chl = min(log10(chl),optics%log10chl_max)
+  else
+    log10chl = optics%log10chl_min
+  endif
   ! Do a nearest neighbor lookup
   n = nint( (log10chl - optics%log10chl_min)/optics%dlog10chl ) + 1