diff --git a/Source/Fortran/FermiOperatorModule.F90 b/Source/Fortran/FermiOperatorModule.F90
index 61f52cac..7acfcfbb 100644
--- a/Source/Fortran/FermiOperatorModule.F90
+++ b/Source/Fortran/FermiOperatorModule.F90
@@ -123,7 +123,6 @@ SUBROUTINE ComputeDenseFOE(H, ISQ, nel, K, inv_temp_in, &
        ELSE
           sval = tlist%DATA(II)%point_value - chemical_potential
           occ = 1.0 / (1 + EXP(inv_temp * sval))
-          WRITE(*,*) II, occ
           energy_value = energy_value + occ * tlist%DATA(II)%point_value
           tlist%DATA(II)%point_value = occ
        END IF