diff --git a/src/MBPT/Bethe_Salpeter_Tmatrix.f90 b/src/MBPT/Bethe_Salpeter_Tmatrix.f90
index 869b7e4..94541b3 100644
--- a/src/MBPT/Bethe_Salpeter_Tmatrix.f90
+++ b/src/MBPT/Bethe_Salpeter_Tmatrix.f90
@@ -139,6 +139,7 @@ subroutine Bethe_Salpeter_Tmatrix(TDA_T,TDA,dBSE,dTDA,evDyn,singlet,triplet,eta,
 
       if(evDyn) then
 
+        print*, ' Iterative dynamical correction for BSE@GT NYI'
 !       call Bethe_Salpeter_dynamic_perturbation_iterative(dTDA,eta,nBas,nC,nO,nV,nR,nS,eGW,dipole_int,OmRPA,rho_RPA, &
 !                                                          OmBSE(:,ispin),XpY_BSE(:,:,ispin),XmY_BSE(:,:,ispin))
       else
@@ -178,6 +179,8 @@ subroutine Bethe_Salpeter_Tmatrix(TDA_T,TDA,dBSE,dTDA,evDyn,singlet,triplet,eta,
 
       if(evDyn) then
 
+    
+        print*, ' Iterative dynamical correction for BSE@GT NYI'
 !       call Bethe_Salpeter_dynamic_perturbation_iterative(dTDA,eta,nBas,nC,nO,nV,nR,nS,eGW,dipole_int,OmRPA,rho_RPA, &
 !                                                          OmBSE(:,ispin),XpY_BSE(:,:,ispin),XmY_BSE(:,:,ispin))
       else
diff --git a/src/MBPT/Bethe_Salpeter_Tmatrix_dynamic_perturbation.f90 b/src/MBPT/Bethe_Salpeter_Tmatrix_dynamic_perturbation.f90
index e0cab84..05b9fc2 100644
--- a/src/MBPT/Bethe_Salpeter_Tmatrix_dynamic_perturbation.f90
+++ b/src/MBPT/Bethe_Salpeter_Tmatrix_dynamic_perturbation.f90
@@ -99,6 +99,7 @@ subroutine Bethe_Salpeter_Tmatrix_dynamic_perturbation(dTDA,eta,nBas,nC,nO,nV,nR
 
     else
 
+      print*, ' Beyond-TDA dynamical correction for BSE@GT NYI'
       ! Resonant and anti-resonant part of the BSE correction
 
 !     call dynamic_Tmatrix_TAB(eta,nBas,nC,nO,nV,nR,nS,nOO,nVV,1d0,eGT,Omega1,Omega2,rho1,rho2,OmBSE(ia), &