/*******************************************************************************
*
*    Projeto 89 - Comparando sensor DHT11 X DHT22 na prática
*    Autor: Angelo Luis Ferreira
*    Data: 17/08/2020
*            http://squids.com.br/arduino
*
*******************************************************************************/
#include <dht.h> // Inclui a biblioteca no seu código
#include <LiquidCrystal_I2C.h>//Carrega a biblioteca lcd com I2C

//Define os pinos que serão utilizados para ligação ao display
LiquidCrystal_I2C lcd(0x27,2,1,0,4,5,6,7,3, POSITIVE);
 
dht DHT; // Cria um objeto da classe dht
uint32_t timer = 0;
 
void setup()
{
  Serial.begin(9600); // Inicializa serial com taxa de transmissão de 9600 bauds
  lcd.begin(20, 4); // Define o display com 16 colunas e 2 linhas
  lcd.clear(); // limpa a tela do display
}
 
void loop() {
  // Executa 1 vez a cada 2 segundos
  if(millis() - timer>= 2000)  {

    // le DHT 22
    // Exibe no display LCD o valor da temperatura
    DHT.read22(A1); // uso do sensor DHT22
    lcd.setCursor(0,0); // Define o cursor na posição de início
    lcd.print("DHT22");
    lcd.setCursor(0,1); // Define o cursor na posição de início
    lcd.print("T= ");
    lcd.print(DHT.temperature,1);
    lcd.write(B11011111); // Imprime o símbolo de grau
    lcd.print("C");

    // Exibe no display LCD o valor da humidade
    lcd.setCursor(0,2); // Define o cursor na posição de início
    lcd.print("U= ");
    lcd.print(DHT.humidity,1);
    lcd.print(" %");

    // Exibe no display o valor do ponto de orvalho (dew point)
    lcd.setCursor(0,3); // Define o cursor na posição de início
    lcd.print("DP ");
    lcd.print(dewPoint(DHT.temperature, DHT.humidity),1);
    lcd.write(B11011111); // Imprime o símbolo de grau
    lcd.print("C");

     // lê DHT11 
     // Exibe no display LCD o valor da temperatura
    DHT.read11(A0); // uso do sensor DHT11
    lcd.setCursor(10,0); // Define o cursor na posição de início
    lcd.print("DHT11");
    lcd.setCursor(10,1); // Define o cursor na posição de início
    lcd.print("T= ");
    lcd.print(DHT.temperature,1);
    lcd.write(B11011111); // Imprime o símbolo de grau
    lcd.print("C");

    // Exibe no display LCD o valor da humidade
    lcd.setCursor(10,2); // Define o cursor na posição de início
    lcd.print("U= ");
    lcd.print(DHT.humidity,1);
    lcd.print(" %");
    
    // Exibe no display o valor do ponto de orvalho (dew point)
    lcd.setCursor(10,3); // Define o cursor na posição de início
    lcd.print("DP ");
    lcd.print(dewPoint(DHT.temperature, DHT.humidity),1);
    lcd.write(B11011111); // Imprime o símbolo de grau
    lcd.print("C");

    timer = millis(); // Atualiza a referência
  }
}

// CALCULO DO PONTO DE ORVALHO (DEW POINT)
// Referência: https://playground.arduino.cc/Main/DHT11Lib/ 
// dewPoint function NOAA
// reference (1) : https://wahiduddin.net/calc/density_algorithms.htm
// reference (2) : https://www.colorado.edu/geography/weather_station/Geog_site/about.htm
//
double dewPoint(double celsius, double humidity) {
  // (1) Saturation Vapor Pressure = ESGG(T)
  double RATIO = 373.15 / (273.15 + celsius);
  double RHS = -7.90298 * (RATIO - 1);
  RHS += 5.02808 * log10(RATIO);
  RHS += -1.3816e-7 * (pow(10, (11.344 * (1 - 1/RATIO ))) - 1) ;
  RHS += 8.1328e-3 * (pow(10, (-3.49149 * (RATIO - 1))) - 1) ;
  RHS += log10(1013.246);

   // factor -3 is to adjust units - Vapor Pressure SVP * humidity
  double VP = pow(10, RHS - 3) * humidity;

  // (2) DEWPOINT = F(Vapor Pressure)
  double T = log(VP/0.61078);   // temp var
  return (241.88 * T) / (17.558 - T);
}