diff --git a/src/physics/cam/check_energy.F90 b/src/physics/cam/check_energy.F90 index 1072211fea..0ba5ce11d1 100644 --- a/src/physics/cam/check_energy.F90 +++ b/src/physics/cam/check_energy.F90 @@ -1241,7 +1241,9 @@ subroutine enthalpy_adjustment(ncol, lchnk, state, cam_in, pbuf, ztodt, itim_old ! ! state%te_cur(:ncol,dyn_te_idx) holds te_endphys(:ncol) ! - state%te_cur(:ncol,dyn_te_idx) = state%te_cur(:ncol,dyn_te_idx)+ztodt*enthalpy_flux_tot(:ncol) + ! ocnfrac: for now apply fixer + ! + state%te_cur(:ncol,dyn_te_idx) = state%te_cur(:ncol,dyn_te_idx)+ztodt*enthalpy_flux_tot(:ncol)*cam_in%ocnfrac(:ncol) call pbuf_set_field(pbuf, teout_idx, state%te_cur(:,dyn_te_idx), (/1,itim_old/),(/pcols,1/)) ! ! the amount of total energy we need energy fixer to fix (associated with enthalpy flux)