Lots of changes :)

Prototype SQUID practicum software PyRPL.md → archive/Prototype SQUID practicum software PyRPL.md

@@ -13,7 +13,9 @@ jupyter:
     name: python3
 ---
 
+<!-- #region hide_input=false -->
 # Load libraries + setup functions
+<!-- #endregion -->
 
 ```python hide_input=true
 import warnings
@@ -38,6 +40,9 @@ from scipy.signal import welch
 
 from datetime import datetime
 
+#Arduino Serial connection
+import serial
+
 # Red Pitaya Stuff
 
 HOSTNAME = "rp-f06897.local"
@@ -87,11 +92,13 @@ def get_IV_t(decimation = 1024):
     return I*1e6,V*1e6,t
 ```
 
+<!-- #region hide_input=false -->
 # Your information
 
 Please update the practicum round
+<!-- #endregion -->
 
-```python
+```python hide_input=false
 academic_year = "2020_2021"
 practicum_round = 4
 ```
@@ -104,7 +111,7 @@ if not os.path.exists(folder):
     os.makedirs(folder)
 ```
 
-<!-- #region hide_input=true -->
+<!-- #region hide_input=false -->
 # IV vs Time Live Display
 <!-- #endregion -->
 
@@ -335,14 +342,34 @@ def update_iv():
     I,V,_ = get_IV_t()
     source1.data = dict(x=I, y=V)
 
-def get_temperature():
-    # For Jasper to update
-    return 300
+def get_temperature(serial_port, baudrate, timeout):
+    # serial_port='COM7', baudrate=9600, timeout=1
+    
+    #variables (could be adjusted to inputs?)
+    alpha = 4.1139e-3; R0 = 1000; R1=1081; Vin=3.272        #should be according to literature: alpha = 3.85e-3 & R0 = 1000
+    
+    #open Serial port to Arduino + flush what's still on it
+    ser = serial.Serial(serial_port,baudrate,timeout=0.5)
+    ser.flush()
+    
+    #sent command to Arduino so we can read it out:
+    time.sleep(0.5)                       #needs to sleep, otherwise we get a: 'can not convert float... at the Vout-line'
+    ser.reset_input_buffer()
+    ser.write(b'g')
+    
+    #Serial read + Temperature calculations:
+    Vout = float(str(ser.readline().decode("utf-8")))*(Vin/1023.0)          #calculates the voltage from the voltage divider (over the PT1000)
+    Rpt = (Vout*R1)/(Vin-Vout)                                              #calculates the resistance of the PT1000
+    Temp_C = ((Rpt/R0)-1)/alpha                                             #with the above calculated resistance and 'known' variable we calculate the temperature in Celsius
+    Temp_K = Temp_C + 273.15                                                #Temperature in Kelvin
+    
+    return Temp_K
+
 
 def take_measurement():
     global timer_t0 
     timer_t0 = time.time()
-    T.append(get_temperature()) 
+    T.append(get_temperature('COM7',9600,1)) 
     I,V,_ = get_IV_t()
     R.append(np.polyfit(I,V,1)[0])
     Iavg.append(np.average(I))

instruction_manual/figs/BaseDoc.blg

+This is BibTeX, Version 0.99dThe top-level auxiliary file: BaseDoc.aux
+I found no \citation commands---while reading file BaseDoc.aux
+I found no \bibdata command---while reading file BaseDoc.aux
+I found no \bibstyle command---while reading file BaseDoc.aux
+(There were 3 error messages)