default version of enthalpy fixer

PeterHjortLauritzen · Jul 29, 2024 · 8ac1716 · 8ac1716
1 parent 7e73102
commit 8ac1716
Showing 1 changed file with 3 additions and 1 deletion.
diff --git 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)