moved getting om_cal only if ISWCAL not equal to 0

when ISWCAL=0 the dependency is not defined (because is not needed)
pmlmodelling · Nov 9, 2022 · 596d710 · 596d710
1 parent 447efb9
commit 596d710
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/benthic_calcite.F90 b/src/benthic_calcite.F90
@@ -82,13 +82,14 @@ subroutine do_bottom(self,_ARGUMENTS_DO_BOTTOM_)
       _HORIZONTAL_LOOP_BEGIN_
 
          _GET_HORIZONTAL_(self%id_c, bL2c)
-         _GET_(self%id_om_cal, om_cal)
 
          if (self%iswcal==0) then  ! NB select case would be cleaner but makes vectorization impossible for ifort 14
             fdiss = 0._rk
          elseif (self%iswcal==1) then
+            _GET_(self%id_om_cal, om_cal)
             fdiss = (max(1._rk-om_cal,0._rk))**self%ndiss
          else
+            _GET_(self%id_om_cal, om_cal)
             fdiss = max(0._rk,(1._rk-om_cal)/(1._rk-om_cal+self%KcalomX))
          end if