initialize g1,g2 <> 0 to avoid division by zero in gotm_fabm::light() for calc_temp=F

ba9b62f0 · Knut · fe4eb016 · ba9b62f0
Commit ba9b62f0 authored May 13, 2018 by Knut
Hide whitespace changes
Inline Side-by-side

Showing with 4 additions and 0 deletions

src/3d/variables_3d.F90 src/3d/variables_3d.F90 +4 -0

No files found.
--- a/src/3d/variables_3d.F90
+++ b/src/3d/variables_3d.F90
@@ -213,6 +213,10 @@
   idpdy=_ZERO_
 #endif

+!  must be nonzero for gotm_fabm in case of calc_temp=F
+   g1 = -9999*_ONE_
+   g2 = -9999*_ONE_
+
 #ifdef DEBUG
   write(debug,*) 'Leaving init_variables_3d()'
   write(debug,*)