diff --git a/gwinc/__init__.py b/gwinc/__init__.py
index 7099de5c418c8b95ba68e9b4a4f530f6daab66fd..fcabced5271dfd79c041c0b2664fb7453618b34c 100644
--- a/gwinc/__init__.py
+++ b/gwinc/__init__.py
@@ -124,8 +124,7 @@ def gwinc(freq, ifo, source=None, plot=False, PRfixed=True):
     finesse = ifo.gwinc.finesse
     prfactor = ifo.gwinc.prfactor
     if ifo.Laser.Power * prfactor != pbs:
-        pass
-        #warning(['Thermal lensing limits input power to ' num2str(pbs/prfactor, 3) ' W']);
+        logging.warning("Thermal lensing limits input power to {} W".format(pbs/prfactor))
 
     # report astrophysical scores if so desired
     score = None