diff --git a/src/tc_scf/fock_three_hermit.irp.f b/src/tc_scf/fock_three_hermit.irp.f
index 6c132189..00d47fae 100644
--- a/src/tc_scf/fock_three_hermit.irp.f
+++ b/src/tc_scf/fock_three_hermit.irp.f
@@ -95,7 +95,12 @@ BEGIN_PROVIDER [double precision, diag_three_elem_hf]
 
   if(.not. three_body_h_tc) then
 
-    diag_three_elem_hf = 0.d0
+   if(noL_standard) then
+      PROVIDE noL_0e
+      diag_three_elem_hf = noL_0e
+    else
+      diag_three_elem_hf = 0.d0
+    endif
 
   else
 
diff --git a/src/tc_scf/rh_tcscf_diis.irp.f b/src/tc_scf/rh_tcscf_diis.irp.f
index 0504373c..66fc83bd 100644
--- a/src/tc_scf/rh_tcscf_diis.irp.f
+++ b/src/tc_scf/rh_tcscf_diis.irp.f
@@ -71,10 +71,7 @@ subroutine rh_tcscf_diis()
   etc_tot = TC_HF_energy
   etc_1e  = TC_HF_one_e_energy
   etc_2e  = TC_HF_two_e_energy
-  etc_3e  = 0.d0
-  if(three_body_h_tc) then
-    etc_3e = diag_three_elem_hf
-  endif
+  etc_3e  = diag_three_elem_hf
   !tc_grad = grad_non_hermit
   er_DIIS = maxval(abs(FQS_SQF_mo))
   e_delta = dabs(etc_tot - e_save)
@@ -202,10 +199,7 @@ subroutine rh_tcscf_diis()
     etc_tot = TC_HF_energy
     etc_1e  = TC_HF_one_e_energy
     etc_2e  = TC_HF_two_e_energy
-    etc_3e  = 0.d0
-    if(three_body_h_tc) then
-      etc_3e = diag_three_elem_hf
-    endif
+    etc_3e  = diag_three_elem_hf
     !tc_grad  = grad_non_hermit
     er_DIIS  = maxval(abs(FQS_SQF_mo))
     e_delta  = dabs(etc_tot - e_save)
@@ -245,7 +239,7 @@ subroutine rh_tcscf_diis()
     write(json_unit, json_real_fmt) ' delta Energy   ', e_delta
     write(json_unit, json_real_fmt) ' DIIS error     ', er_DIIS
     write(json_unit, json_real_fmt) ' level_shift    ', level_shift_tcscf
-    write(json_unit, json_real_fmt) ' DIIS           ', dim_DIIS
+    write(json_unit, json_int_fmt)  ' DIIS           ', dim_DIIS
     write(json_unit, json_real_fmt) ' Wall time (min)', (t1-t0)/60.d0
     call unlock_io