diff --git a/concore.py b/concore.py
index 2da1250..7a99b8f 100644
--- a/concore.py
+++ b/concore.py
@@ -7,9 +7,17 @@
 import zmq
 import numpy as np
 logging.basicConfig(
-    level=logging.INFO,
-    format='%(levelname)s - %(message)s'
-) # Added for ZeroMQ
+    level=logging.DEBUG, 
+    format='%(asctime)s - %(name)s - %(levelname)s - %(message)s', 
+    force=True 
+)
+
+#these lines mute the noisy library
+logging.getLogger('matplotlib').setLevel(logging.WARNING)
+logging.getLogger('PIL').setLevel(logging.WARNING)
+logging.getLogger('urllib3').setLevel(logging.WARNING) 
+logging.getLogger('requests').setLevel(logging.WARNING) 
+
 
 # if windows, create script to kill this process 
 # because batch files don't provide easy way to know pid of last command
diff --git a/tools/bangbang.py b/tools/bangbang.py
index eed72b8..3bf9e5a 100644
--- a/tools/bangbang.py
+++ b/tools/bangbang.py
@@ -1,5 +1,6 @@
 import numpy as np
 import concore
+import logging
 
 
 def bangbang_controller(ym):
@@ -19,6 +20,9 @@ def bangbang_controller(ym):
 init_simtime_u = "[0.0, 0.0,0.0]"
 init_simtime_ym = "[0.0, 70.0,91]"
 u = np.array([concore.initval(init_simtime_u)]).T
+
+logging.info("Starting Bang-Bang Controller")
+
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
@@ -26,5 +30,5 @@ def bangbang_controller(ym):
     
     ustar = bangbang_controller(ym)
     
-    print(str(concore.simtime) + " u="+str(ustar) + "ym="+str(ym))
+    logging.debug(f"{concore.simtime} u={ustar} ym={ym}")
     concore.write(1,"u",list(ustar),delta=0)
diff --git a/tools/cardiac_pm.py b/tools/cardiac_pm.py
index d3edb30..0fd5924 100644
--- a/tools/cardiac_pm.py
+++ b/tools/cardiac_pm.py
@@ -2,6 +2,7 @@
 import pulsatile_model_functions as pmf
 import healthy_params as K
 import concore
+import logging
 
 #x0 = np.loadtxt('pulsatile_steady.txt')
 
@@ -20,6 +21,9 @@ def cardiac_pm(x0,u):
 init_simtime_u = "[0.0, 0.0,0.0]"
 init_simtime_ym = "[0.0, 70,0]"
 ym = np.array([concore.initval(init_simtime_ym)]).T
+
+logging.info("Starting Cardiac PM Model")
+
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         u = concore.read(1,"u",init_simtime_u)
@@ -29,7 +33,7 @@ def cardiac_pm(x0,u):
     
     #dummy = np.array([hr,mapp])
     dummy = np.array([mapp,hr])
-    print(str(concore.simtime) + " u="+str(u) + "ym=" + str(dummy))
+    logging.debug(f"{concore.simtime} u={u} ym={dummy}")
     #concore.write(1,"ym",list(np.array([hr,mapp])),delta=1)
     concore.write(1,"ym",list(np.array([mapp,hr])),delta=1)
 
diff --git a/tools/cwrap.py b/tools/cwrap.py
index b8cfb4b..3523d00 100644
--- a/tools/cwrap.py
+++ b/tools/cwrap.py
@@ -4,6 +4,7 @@
 import time
 from ast import literal_eval
 import os
+import logging
 
 #time.sleep(7)
 timeout_max = 20
@@ -58,10 +59,10 @@
     except:
         init_simtime_ym = "[0.0, 0.0, 0.0]"
 
-print(apikey)
-print(yuyu)
-print(name1+'='+init_simtime_u)
-print(name2+'='+init_simtime_ym)
+logging.info(f"API Key: {apikey}")
+logging.info(f"Yuyu: {yuyu}")
+logging.info(f"{name1}={init_simtime_u}")
+logging.info(f"{name2}={init_simtime_ym}")
 
 while not os.path.exists(concore.inpath+'1/'+name1):
     time.sleep(concore.delay)
@@ -78,47 +79,45 @@
 
 
 while(concore.simtime<concore.maxtime):
-    print("CW outer loop")
+    logging.debug("CW outer loop")
     while concore.unchanged():
         u = concore.read(1,name1,init_simtime_u)
     f = {'file1': open(concore.inpath+'1/'+name1, 'rb')}
-    print("CW: before post u="+str(u))
-    print('http://www.controlcore.org/pm/'+yuyu+apikey+'&fetch='+name2)
+    logging.debug(f"CW: before post u={u}")
+    logging.debug(f'http://www.controlcore.org/pm/{yuyu}{apikey}&fetch={name2}')
     r = requests.post('http://www.controlcore.org/pm/'+yuyu+apikey+'&fetch='+name2, files=f,timeout=timeout_max)
     if r.status_code!=200:
-        print("bad POST request "+str(r.status_code))
+        logging.error(f"bad POST request {r.status_code}")
         quit()
     if len(r.text)!=0:
         try:
             t=literal_eval(r.text)[0]
         except:
-            print("bad eval "+r.text)
+            logging.error(f"bad eval {r.text}")
     timeout_count = 0
     t1 = time.perf_counter()
-    print("CW: after post status="+str(r.status_code)+" r.content="+str(r.content)+" t="+str(t))
+    logging.debug(f"CW: after post status={r.status_code} r.content={r.content} t={t}")
     #while r.text==oldym or len(r.content)==0:
     while oldt==t or len(r.content)==0:
         time.sleep(concore.delay)
-        print("CW waiting status="+str(r.status_code)+" content="+ r.content.decode('utf-8')+" t="+str(t))
+        logging.debug(f"CW waiting status={r.status_code} content={r.content.decode('utf-8')} t={t}")
         f = {'file1': open(concore.inpath+'1/'+name1, 'rb')}
         try:
             r = requests.post('http://www.controlcore.org/pm/'+yuyu+apikey+'&fetch='+name2, files=f,timeout=timeout_max)
         except:
-            print("CW: bad request")
+            logging.error("CW: bad request")
         timeout_count += 1
         if r.status_code!=200 or time.perf_counter()-t1 > 1.1*timeout_max: #timeout_count>100:
-            print("timeout or bad POST request "+str(r.status_code))
+            logging.error(f"timeout or bad POST request {r.status_code}")
             quit()
         if len(r.text)!=0:
             try:
                 t=literal_eval(r.text)[0]
             except:
-                print("bad eval "+r.text)
+                logging.error(f"bad eval {r.text}")
     oldt = t
     oldym = r.text
-    print("CW: oldym="+oldym+" t="+str(concore.simtime))
+    logging.debug(f"CW: oldym={oldym} t={concore.simtime}")
     concore.write(1,name2,oldym)
 #concore.write(1,"ym",init_simtime_ym)
-print("retry="+str(concore.retrycount))
-
-
+logging.info(f"retry={concore.retrycount}")
\ No newline at end of file
diff --git a/tools/learn.py b/tools/learn.py
index a29e1e6..53bb50e 100644
--- a/tools/learn.py
+++ b/tools/learn.py
@@ -2,6 +2,7 @@
 import numpy as np
 import matplotlib.pyplot as plt
 import time
+import logging
 GENERATE_PLOT = 1
 
 concore.delay = 0.002
@@ -21,7 +22,7 @@
         ut[int(concore.simtime)] = np.array(u).T
         ymt[int(concore.simtime)] = np.array(ym).T
     oldsimtime = concore.simtime
-print("retry="+str(concore.retrycount))
+logging.info(f"retry={concore.retrycount}")
 
 #################
 # plot inputs and outputs
diff --git a/tools/pid2.py b/tools/pid2.py
index 6f88285..1c720f8 100644
--- a/tools/pid2.py
+++ b/tools/pid2.py
@@ -1,5 +1,7 @@
 import numpy as np
 import concore
+import logging
+
 setpoint = 67.5
 setpointF = 75.0
 KpF = 0.1
@@ -51,7 +53,9 @@ def  pid_controller(ym):
 init_simtime_u = "[0.0, 0.0,0.0]"
 init_simtime_ym = "[0.0, 70.0,91]"
 u = np.array([concore.initval(init_simtime_u)]).T
-print("Shannon's PID controller: setpoint is "+str(setpoint))
+
+logging.info(f"Shannon's PID controller: setpoint is {setpoint}")
+
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
@@ -59,6 +63,6 @@ def  pid_controller(ym):
     
     ustar =  pid_controller(ym)
     
-    print(str(concore.simtime) + " u="+str(ustar) + "ym="+str(ym))
+    logging.debug(f"{concore.simtime} u={ustar} ym={ym}")
     concore.write(1,"u",list(ustar),delta=0)
 
diff --git a/tools/pidmayuresh.py b/tools/pidmayuresh.py
index 900d32f..519c750 100644
--- a/tools/pidmayuresh.py
+++ b/tools/pidmayuresh.py
@@ -1,6 +1,7 @@
 import numpy as np
 import math
 import concore
+import logging
 dT = 0.1
 
 sp = concore.tryparam('sp', 67.5)
@@ -19,7 +20,7 @@ def  pid_controller(state, ym, sp, Kp, Ki, Kd, sigout, cin, low, up):
     elif cin == 'map':
         Error = sp - ym[0]
     else:
-        print('invalid control input '+cin)
+        logging.error(f'invalid control input {cin}')
         quit()
     P = Error
     I = I + Error*dT 
@@ -37,15 +38,15 @@ def  pid_controller(state, ym, sp, Kp, Ki, Kd, sigout, cin, low, up):
 init_simtime_ym = "[0.0, 70.0,91]"
 ym = np.array(concore.initval(init_simtime_ym))
 state = [0.0, 0.0]
-print("Mayuresh's PID controller: sp is "+str(sp))
-print(concore.params)
+logging.info(f"Mayuresh's PID controller: sp is {sp}")
+logging.info(concore.params)
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
     ym = np.array(ym)
     (state,amp) =  pid_controller(state,ym,sp,Kp,Ki,Kd,sigout,cin,0,3)
     u = np.array([amp,freq])    
-    print(str(concore.simtime) + " u="+str(u) + "ym="+str(ym))
+    logging.debug(f"{concore.simtime} u={u} ym={ym}")
     concore.write(1,"u",list(u),delta=0)
 
 
diff --git a/tools/pidsig.py b/tools/pidsig.py
index b89a917..42b3dda 100644
--- a/tools/pidsig.py
+++ b/tools/pidsig.py
@@ -1,6 +1,7 @@
 import numpy as np
 import math
 import concore
+import logging
 dT = 0.1
 global Prev_Error, I, freq
 Prev_Error = 0
@@ -42,7 +43,7 @@ def  pid_controller(ym):
     elif cin == 'map':
         Error = sp - ym[0]
     else:
-        print('invalid control input '+cin)
+        logging.error(f'invalid control input {cin}')
         quit()
     P = Error
     I = I + Error*dT 
@@ -60,8 +61,8 @@ def  pid_controller(ym):
 init_simtime_u = "[0.0, 0.0,0.0]"
 init_simtime_ym = "[0.0, 70.0,91]"
 u = np.array([concore.initval(init_simtime_u)]).T
-print("Shannon's PID controller: sp is "+str(sp))
-print(concore.params)
+logging.info(f"Shannon's PID controller: sp is {sp}")
+logging.info(concore.params)
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
@@ -72,7 +73,7 @@ def  pid_controller(ym):
     else:
         ustar =  pid_controller(ym)
     
-    print(str(concore.simtime) + " u="+str(ustar) + "ym="+str(ym))
+    logging.debug(f"{concore.simtime} u={ustar} ym={ym}")
     concore.write(1,"u",list(ustar),delta=0)
 
 
diff --git a/tools/plotu.py b/tools/plotu.py
index 2507702..73edc6e 100644
--- a/tools/plotu.py
+++ b/tools/plotu.py
@@ -1,8 +1,9 @@
 import concore
+import logging
 import numpy as np
 import matplotlib.pyplot as plt
 import time
-print("plot u")
+logging.info("plot u")
 
 concore.delay = 0.005
 concore.default_maxtime(150)
@@ -15,9 +16,9 @@
     while concore.unchanged():
         u = concore.read(1,"u",init_simtime_u)
     concore.write(1,"u",u)
-    print("u="+str(u))
+    logging.debug(f"u={u}")
     ut.append(np.array(u).T)
-print("retry="+str(concore.retrycount))
+logging.info(f"retry={concore.retrycount}")
 
 #################
 
diff --git a/tools/plotym.py b/tools/plotym.py
index 9f0deb9..cd9bfae 100644
--- a/tools/plotym.py
+++ b/tools/plotym.py
@@ -1,8 +1,9 @@
 import concore
+import logging
 import numpy as np
 import matplotlib.pyplot as plt
 import time
-print("plot ym")
+logging.info("plot ym")
 
 concore.delay = 0.005
 concore.default_maxtime(150)
@@ -15,9 +16,9 @@
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
     concore.write(1,"ym",ym)
-    print(" ym="+str(ym))
+    logging.debug(f" ym={ym}")
     ymt.append(np.array(ym).T)
-print("retry="+str(concore.retrycount))
+logging.info(f"retry={concore.retrycount}")
 
 #################
 
diff --git a/tools/plotymlag.py b/tools/plotymlag.py
index 8df1586..e81fca6 100644
--- a/tools/plotymlag.py
+++ b/tools/plotymlag.py
@@ -1,11 +1,12 @@
 import concore
+import logging
 import numpy as np
 import matplotlib.pyplot as plt
 import time
 
 size = 10
 lag = concore.tryparam("lag", 0) 
-print("plot ym with lag="+str(lag))
+logging.info(f"plot ym with lag={lag}")
 
 concore.delay = 0.005
 concore.default_maxtime(150)
@@ -21,10 +22,10 @@
     while concore.unchanged():
         ym[cur] = concore.read(1,"ym",init_simtime_ym)
     concore.write(1,"ym",ym[cur])
-    print(" ym="+str(ym[cur]))
+    logging.debug(f" ym={ym[cur]}")
     ymt.append(np.array(ym[(cur-lag) % size]).T)
     cur = (cur+1) % size
-print("retry="+str(concore.retrycount))
+logging.info(f"retry={concore.retrycount}")
 
 #################
 
diff --git a/tools/pwrap.py b/tools/pwrap.py
index beb9011..73f1651 100644
--- a/tools/pwrap.py
+++ b/tools/pwrap.py
@@ -4,6 +4,7 @@
 import time
 from ast import literal_eval
 import os
+import logging
 
 #time.sleep(7)
 timeout_max=20
@@ -61,10 +62,10 @@
     except:
         init_simtime_ym = "[0.0, 0.0, 0.0]"
 
-print(apikey)
-print(yuyu)
-print(name1+'='+init_simtime_u)
-print(name2+'='+init_simtime_ym)
+logging.info(f"API Key: {apikey}")
+logging.info(f"Yuyu: {yuyu}")
+logging.info(f"{name1}={init_simtime_u}")
+logging.info(f"{name2}={init_simtime_ym}")
 
 while not os.path.exists(concore.inpath+'1/'+name2):
     time.sleep(concore.delay)
@@ -86,44 +87,44 @@
 
 
 while(concore.simtime<concore.maxtime):
-    print("PW outer loop")
+    logging.debug("PW outer loop")
     while concore.unchanged():
         ym = concore.read(1,name2,init_simtime_ym)
     f = {'file1': open(concore.inpath+'1/'+name2, 'rb')}
-    print("PW: before post ym="+str(ym))
+    logging.debug(f"PW: before post ym={ym}")
     r = requests.post('http://www.controlcore.org/ctl/'+yuyu+apikey+'&fetch='+name1, files=f,timeout=timeout_max)
     if r.status_code!=200:
-        print("bad POST request "+str(r.status_code))
+        logging.error(f"bad POST request {r.status_code}")
         quit()
     if len(r.text)!=0:
         try:
             t=literal_eval(r.text)[0]
         except:
-            print("bad eval "+r.text)
+            logging.error(f"bad eval {r.text}")
     timeout_count = 0
     t1 = time.perf_counter()
-    print("PW: after post status="+str(r.status_code)+" r.content="+str(r.content)+"/"+r.text)
+    logging.debug(f"PW: after post status={r.status_code} r.content={r.content}/{r.text}")
     #while r.text==oldu or len(r.content)==0:
     while oldt==t or len(r.content)==0:
         time.sleep(concore.delay)
-        print("PW waiting status="+str(r.status_code)+" content="+ r.content.decode('utf-8')+" t="+str(t))
+        logging.debug(f"PW waiting status={r.status_code} content={r.content.decode('utf-8')} t={t}")
         f = {'file1': open(concore.inpath+'1/'+name2, 'rb')}
         try:
             r = requests.post('http://www.controlcore.org/ctl/'+yuyu+apikey+'&fetch='+name1, files=f,timeout=timeout_max)
         except:
-            print("PW: bad requests")
+            logging.error("PW: bad requests")
         timeout_count += 1
         if r.status_code!=200 or time.perf_counter()-t1 > 1.1*timeout_max: #timeout_count>200:
-            print("timeout or bad POST request "+str(r.status_code))
+            logging.error(f"timeout or bad POST request {r.status_code}")
             quit()
         if len(r.text)!=0:
             try:
                 t=literal_eval(r.text)[0]
             except:
-                print("bad eval "+r.text)
+                logging.error(f"bad eval {r.text}")
     oldt = t
     oldu = r.text
-    print("PW: oldu="+oldu+" t="+str(concore.simtime))
+    logging.debug(f"PW: oldu={oldu} t={concore.simtime}")
     concore.write(1,name1,oldu)
 #concore.write(1,"u",init_simtime_u)
-print("retry="+str(concore.retrycount))
+logging.info(f"retry={concore.retrycount}")
\ No newline at end of file
diff --git a/tools/shannon.py b/tools/shannon.py
index ba371af..9da6f51 100644
--- a/tools/shannon.py
+++ b/tools/shannon.py
@@ -1,5 +1,6 @@
 import numpy as np
 import concore
+import logging
 setpoint = 67.5
 Kp = 0.1
 Ki = 0.01
@@ -27,7 +28,7 @@ def  pid_controller(ym):
 init_simtime_u = "[0.0, 0.0,0.0]"
 init_simtime_ym = "[0.0, 70.0,91]"
 u = np.array([concore.initval(init_simtime_u)]).T
-print("Shannon's PID controller: setpoint is "+str(setpoint))
+logging.info(f"Shannon's PID controller: setpoint is {setpoint}")
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
@@ -35,6 +36,6 @@ def  pid_controller(ym):
     
     ustar =  pid_controller(ym)
     
-    print(str(concore.simtime) + " u="+str(ustar) + "ym="+str(ym))
+    logging.debug(f"{concore.simtime} u={ustar} ym={ym}")
     concore.write(1,"u",list(ustar),delta=0)