error_status.py


import time
import math
from sensors import *
#Import the firmware code for the Adafruit MAX31855 Amplifier for the Thermocouple
import Adafruit_MAX31855.MAX31855 as MAX31855
max_sensor = MAX31855.MAX31855(21, 20, 16)

#import smbus for use with the DHT20 sensor and setting up the bus with the address
import  smbus
import numpy as np
bus = smbus.SMBus(1)
#SHT20 address, 0x40(64)

#import MCP3008 (ADC) Library 
import Adafruit_MCP3008
mcp = Adafruit_MCP3008.MCP3008(clk=22, cs=25, miso=23, mosi=24)

#importing the DHT11 sensor library
import Adafruit_DHT
 
 
# a function that returns a string with the system status and stores that in the Store/status.csv file
def statusUpdate():
    try:
        while True:
            statusString = 0
            sensorString = 0
            
            #check for tank leaks
            tankLeak = checkTankLeak()
            if tankLeak["status"] == 1:
                statusString += 1
            
            #check for compressor overheat
            overHeat = compressorOverHeat()
            if overHeat["status"] == 1:
                statusString += 10
            
            #check for compressor current overload
            overLoad = compressorCurrentOverload()
            if overLoad["status"] == 1:
                statusString += 100
            
            #check for dryer failure
            dryerFailureState = dryerFailure()
            if dryerFailureState["status"] == 1:
                statusString += 1000
                
            #check for radiator failure
            radiatorFailureState = radiatorFailure()
            if radiatorFailureState["status"] == 1:
                statusString += 10000
            
            #check if AM2302 sensor is online
            if am2302Offline() == 1:
                sensorString += 1
            
            #check if SHT20 sensor is online
            if sht20Offline() == 1:
                sensorString += 10
                
            #check if Thermocouple is online
            if thermocoupleOffline() == 1:
                sensorString += 100
                
            #check if current sensor is online
            if currentSensorOffline() == 1:
                sensorString += 1000
                
            #check if pressure transducer is online
            if pressureTransducerOffline() == 1:
                sensorString += 10000

            
            with open("Store/status.csv", "w") as csvFile:
                fieldnames = ["status", "sensor"]
                writer = csv.DictWriter(csvFile, dialect='excel', fieldnames=fieldnames)
                writer.writerow({"status": statusString, "sensor": sensorString})
            
            
    except KeyboardInterrupt:
        print("we stopped")
        return
        
def checkTankLeak():
    array = readerCSV("pressure.csv",10)
    deltaArray = []
    for i in range(len(array)):
        try:
            deltaArray.append(abs(array[i+1]-array[i]))
        except LookupError:
            break     
    avgRate = round(sum(deltaArray)/len(deltaArray),3)
    if avgRate > 0.5 and checkCompressorState() == 0:
        return {"status":1, "value":avgRate}
    elif avgRate < 0.5 and checkCompressorState() != 0:
        return {"status":1, "value":avgRate}
    else:
        return {"status":0, "value":avgRate}

def compressorOverHeat():
    array = readerCSV("thermocouple.csv", 10)
    avgTemp = round(sum(array)/len(array),3)
    
    if avgTemp > 40:
        return {"status":1, "value":avgTemp}
    else:
        return {"status":0, "value":avgTemp}

def compressorCurrentOverload():
    array = readerCSV("current.csv", 50)
    avg = round(sum(array)/len(array),3)
    
    if avg > 6.8:
        return {"status":1, "value":avg}
    else:
        return {"status":0, "value":avg}

# formula for dewpoint 
def dewPoint(T, RH):
    d = 234.5
    if T < 0:
        b = 17.966
        c = 247.15
    else:
        b = 17.368
        c = 238.88
    val = (b-(T/d))*(T/(c+T))
    if val<= 0:
        ardBuck = RH/100
    else:
        ardBuck = math.log((RH/100)*(math.exp(val)))
    dp = (c*ardBuck)/(b-ardBuck)
    return dp
    
def dryerFailure():
    #Read last 1000 values from CSV file
    tankRH = readerCSV("sht20Humidity.csv", 1000)
    tankTemp = readerCSV("sht20Temperature.csv", 1000)
    avgTankRH = sum(tankRH)/len(tankRH)
    
    if avgTankRH > 40:
       data = 1
    else:
       data = 0
    
    dpArray = []
    for i in range(len(tankTemp)):
        dpArray.append(dewPoint(tankTemp[i], tankRH[i]))
    dp = round(sum(dpArray)/len(dpArray),3)
    return {"status":data, "value":dp}

def radiatorFailure():
    ambientTemp = readerCSV("am2302Temperature.csv",25)
    tankTemp = readerCSV("sht20Temperature.csv", 25)
    avgAmbientTemp = sum(ambientTemp)/len(ambientTemp)
    avgTankTemp = round(sum(tankTemp)/len(tankTemp),3)
    
    if avgTankTemp > avgAmbientTemp:
        data = 1
    else:
       data = 0
    return {"status":data, "value":avgTankTemp}

def am2302Offline():
    DHTsensor = Adafruit_DHT.DHT22
    am2302Humidity,am2302Temperature = Adafruit_DHT.read_retry(DHTsensor,18)
    try: 
        if am2302Humidity is None or am2302Temperature is None:
            data = 1
        else:
            data = 0
    except:
        data = 1
    finally:
        return data

def sht20Offline():
    try:
        if GetShtTemp_c() is None or GetShtHumid() is None:
            data = 1
        else:
            data = 0
    except:
        data = 1
    finally:
        return data

def thermocoupleOffline():
    try:
        if np.isnan(max_sensor.readTempC()) or max_sensor.readTempC() is None or not max_sensor.readTempC():
            data = 1
        else:
            data = 0
    except:
        data = 1
    finally:   
        return data

def currentSensorOffline():
    try:
        if mcp.read_adc(4) == 0:
            data = 1
        else:
            data = 0
    except:
        data = 1
    finally:   
        return data

def pressureTransducerOffline():
    try:
        if (mcp.read_adc(3)-73)*150.0/595.2 <= 0:
            data = 1
        else:
            data = 0
    except:
        data = 1
    finally:
        return data


if __name__ == "__main__":
    statusUpdate()