daq-bbb-code/daq.c

/*
*/

#include<stdio.h>
#include<fcntl.h>
#include<unistd.h>
#include<string.h>
#include<stdlib.h>
#include<time.h>
#include"uart.c"

#define data_file_path  "/var/lib/cloud9/daq-serial-bbb/data.dat"

#define GC0017     4    //co2
#define CM31911    1    //co, t
#define OX0052     2    //o2

#define Rep_Dev_ID      "Y\r\n"
#define Pulling_Mode    "K 2\r\n"
#define OX_P_Mode       "M 1\r\n"
#define FILnUNFIL       "M 6\r\n"
#define unf_gas_con     "z\r\n"
#define fil_gas_con     "Z\r\n"
#define Temperature     "T\r\n"
#define get_readigns    "Q\r\n"
#define percent_oxigen  "%\r\n"
#define ppm_oxigen      "O\r\n"
#define M_zZTHBD        "M 14406\r\n"
#define Readings_OX     "A\r\n"


int sensConf(unsigned char uartNumber, int baudRate, unsigned char mode[], unsigned char response[], int tries);
int measure(unsigned char uartNumber, unsigned char command[], int multiplier, int excess);
char *gmeasures(char src[], char fval, int nchar);

int str2int(unsigned char *ptr);

int DAQ(int t_hrs, int tm_s);


int main(int argc, char *argv[]){

    //char bioxCsen[], OXsen[], monoxCsen[];

    int T,t;

    if(argc!=3){
        printf("Invalid number of arguments,\n usage: daq <time(h)> <sampling period(s)>\nexiting!\n");
        return -2;
    }

    //configuración de sensores
    printf("Configuring CO2 sensor\n");
    sensConf(GC0017, B9600, Pulling_Mode, " K 00002\r\n", 5);
    sensConf(GC0017, B9600, FILnUNFIL, " M 00006\r\n", 5);
    printf("Configuring CO sensor\n");
    sensConf(CM31911, B9600, Pulling_Mode, "K 00002\r\n", 5);
    sensConf(CM31911, B9600, M_zZTHBD, "M 14406\r\n", 5);
    printf("Configuring O2 sensor\n");
    sensConf(OX0052, B9600, OX_P_Mode, "M 01\r\n", 5);

    //adquisición
    T = atoi(argv[1]);//str to int
    t = atoi(argv[2]);

    printf("Starting data acquisition with duration of %dh and every %ds\n", T, t);

    DAQ(T, t);
    sleep(1);

    //finalización
    //deinit uart
    uartClose(GC0017);
    uartClose(CM31911);
    uartClose(OX0052);

    printf("Exiting of the program...\n");
    return 0;
}


int sensConf(unsigned char uartNumber, int baudRate, unsigned char mode[], unsigned char response[], int tries){

    int count;

    uartConf(uartNumber, baudRate);

    while(tries){
        uartTransmit(uartNumber, mode);
        tcdrain(uartFile[uartNumber]);  //wait all data has been sent
        printf("Command sended.\n");
        //usleep(100000);                  //give the sensor a chance to respond
        count = uartReceive(uartNumber);
        if (count == 0) printf("There was no data available to read!\n");
        else if (strcmp(receive[uartNumber], response) == 0) {
                printf("Sensor configurated.\n");
                return 0;
        } else {
            printf("The following was read in [%d]: %s\n",count,receive[uartNumber]);
            //char *c = receive[uartNumber];
            //while (*c != '\0'){
            //    printf("%d = '%c'\n",*c, *c);
            //    c++;
            //}
        }
        tries --;
    }
    printf("Sensor configuration failed.\n");
    return -1;
}

int measure(unsigned char uartNumber, unsigned char command[], int multiplier, int excess){
    int measure;
    uartTransmit(uartNumber, command);
    uartReceive(uartNumber);
    measure = (str2int(&receive[uartNumber]) - excess)*multiplier ;
    return measure;
}

int read_measure(unsigned char uartNumber, int multiplier, int excess){
    int measure;
    uartReceive(uartNumber);
    measure = (str2int(&receive[uartNumber]) - excess)*multiplier ;
    return measure;
}

int str2int(unsigned char *ptr){

  int number = 0;

  while (*ptr != '\0') {
    if ((*ptr >= '0') && (*ptr <= '9')) {
      int dig = (*ptr) - '0';
      number = number*10 + dig;
    }
    ptr++;
  }
  return number;
}

int DAQ(int t_hrs, int tm_s)
{
  FILE* dfp;                                    // create a file pointer
  unsigned char co2_uf[10]="", co2_f[10]="",
              co_uf[10]="", co_f[10]="",
              o2_ppm[10]="", o2_xcent[10]="",
              co2_temp[10]="", o2_temp[10]="",
              co2_press[10]="", o2_press[10]="",
              co2_relH[10]="", DATA[100]="";    //measurements varirables
  clock_t start_t, end_t;                       //processing time measurements variables
  time_t new_time, prev_time, t0;               //time control variables
	double tm_ms = tm_s*1e6 - 0.20e6;             //tiempo de muestreo en milisegundos (-) margen
	double dif = tm_ms;                           //diferencia de tiempo(time to sleep)

  //registro datos de inicio
  dfp = fopen(data_file_path, "w");             // open file for writing
  time_t curtime;                               //current time (date)
  time(&curtime);                               //saving date in curtime
  fprintf(dfp, "%sStarting DAQ\ntime(s)\tCO2 unfil(ppm)\tCO2 fil(ppm)\tCO unfil(ppm)\tCO fil(ppm)"\
              "\tO2(ppm)\tO2(%c)\tTemperaure C02(ºC*10)\tTemperaure 02(ºC)"\
              "\tPressure C02(.mBar)\tPresure 02(mBar)\tRelative Humidity(.)\n",
              ctime(&curtime), '%');            // send the value to the file
  fclose(dfp);                                  // close the file using the file pointer

  //ciclo
  t0 = time(NULL) + tm_s;                       //intitial time = current time + tiempo de muestreo
  new_time = t0;                                //saving initial time in time control variable
  for(int i = 0; i < (t_hrs*3600); i=i+tm_s){   //cycle from 0 to adquisition time (seconds), incrementing sampling period

    if (dif <= 0) printf("T=%d Communication time exceded.\n", (int)(new_time-t0));
    else usleep((int)dif);                      //checking if there's time to sleep
    //waiting new second (time.h system) to start measurements
    if (time(NULL) < new_time) while(new_time != time(NULL));
    else printf("T=%d time error.\n", (int)(new_time-t0));

    start_t = clock();                          //saving start time
      /*new_time = time(NULL);                    //refresh new_time (next iteration)
        no necesario por:
        while(new_time != time(NULL));
        se avanza cuando new_time = time(NULL)
        por lo tanto "new_time = time(NULL);" es redundante
      */
      //transmiting commands to sensors
      uartTransmit(CM31911, get_readigns);
      uartTransmit(OX0052, Readings_OX);
      uartTransmit(GC0017, get_readigns
      //receiving replys from sensors
      uartReceive(OX0052);
      uartReceive(CM31911);
      uartReceive(GC0017);
      //saving gas measurements in measurements variables
      memcpy(co2_uf,   gmeasures(receive[GC0017],  'z', 5), 5);
      memcpy(co2_f,    gmeasures(receive[GC0017],  'Z', 5), 5);
      memcpy(co_uf,    gmeasures(receive[CM31911], 'z', 5), 5);
      memcpy(co_f,     gmeasures(receive[CM31911], 'Z', 5), 5);
      memcpy(o2_ppm,   gmeasures(receive[OX0052],  'O', 6), 6);
      memcpy(o2_xcent, gmeasures(receive[OX0052],  '%', 6), 6);
      memcpy(co2_temp, gmeasures(receive[CM31911], 'T', 5), 5);
      memcpy(o2_temp,  gmeasures(receive[OX0052],  'T', 5), 5);
      memcpy(co2_press,gmeasures(receive[CM31911], 'B', 5), 5);
      memcpy(o2_press, gmeasures(receive[OX0052],  'P', 4), 4);
      memcpy(co2_relH, gmeasures(receive[CM31911], 'H', 5), 5);
      //saving formated measurements in string DATA
      sprintf(DATA, "%d\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s",
              (int)(new_time-t0), co2_uf, co2_f, co_uf, co_f, o2_ppm, o2_xcent,
              co2_temp, o2_temp, co2_press, o2_press, co2_relH);
      //showing measurements on display
      printf("time(s)\tCO2 unfil(ppm)\tCO2 fil(ppm)\tCO unfil(ppm)\tCO fil(ppm)"\
          "\tO2(ppm)\tO2(%c)\tTemperaure C02(ºC*10)\tTemperaure 02(ºC)"\
          "\tPressure C02(.mBar)\tPresure 02(mBar)\tRelative Humidity(.)\n", '%');
      printf("%s\n", DATA);

      dfp = fopen(data_file_path, "a");         // open file for writing
          fprintf(dfp, "%s\n", DATA);           // saving measurements string to data file
      fclose(dfp);                              // close the file using the file pointer
      new_time += tm_s;                         //adding sampling period to new_time
    end_t = clock();                            //saving end time
    //calculate time to sleep
    dif =  tm_ms - ((double)(end_t - start_t)*1e6 / CLOCKS_PER_SEC);
  }

  return 0;
}

char *gmeasures(char src[], char fval, int nchar){
	char * ptr = &src[0];
	static char s[10]="";

	while(*ptr != fval) ptr++;
	ptr += 2;
	memcpy(s, ptr, nchar);
	return s;
}