#include <cstdio>
// Robert Grumbine
// 4 Febuary 2009

#ifndef SSMICLASS_INCLUDE
  #define SSMICLASS_INCLUDE
  
  #ifndef POINTSH
    #include "points.h"
  #endif
  #ifndef PARAMSH
    #include "params.h"
  #endif
  #ifndef SSMI_INCLUDE
    #include "ssmi.h"
  #endif
  #ifndef ICESSMI
    #include "icessmi.h"
  #endif

// ssmiclass begun 5/5/2004  Robert Grumbine
//   Purpose is to transition to having a class for ssmi processing,
//   rather than requiring C-only codes 
// 4 Feb 2009 RG note: want to do something stable about the flag values
//    used in processing -- LAND, ...
// 4 February 2009: add member functions polarization, gradient, differences, qc, show


// Rely on ssmi.h for the elementary definitions of functions,
//   types, and flags.  Leave just the classlike behaviors here.
// Rely on icessmi.h for others
// 7 May 2013


class ssmipt {
  public:
    ssmi obs;
    ssmipt();
    operator double();
    ssmipt(ssmipt &);
    ssmipt & operator=(ssmipt &);
// Added 4 February 2009, after versions in 28 Jul 2000 experiment
    point3<float> polarization();
    point3<float> gradient();
    point3<float> differences();
    int qc();
    void show();
// Future: To add, after methods in ssmiclass.operations.h 
//    ssmipt & operator=(int );
};


ssmipt::ssmipt() {
  obs.t19v = 0;
  obs.t19h = 0;
  obs.t22v = 0;
  obs.t37v = 0;
  obs.t37h = 0;
  obs.t85v = 0;
  obs.t85h = 0;
  obs.conc_bar      = 0;
  obs.bar_conc      = 0;
  obs.count         = 0;
  obs.hires_conc    = 0;
  obs.weather_count = 0;
}
ssmipt::operator double() {
  return (double) obs.conc_bar;
}
ssmipt::ssmipt(ssmipt &x) {
  obs.t19v = x.obs.t19v;
  obs.t19h = x.obs.t19h;
  obs.t22v = x.obs.t22v;
  obs.t37v = x.obs.t37v;
  obs.t37h = x.obs.t37h;
  obs.t85v = x.obs.t85v;
  obs.t85h = x.obs.t85h;
  obs.conc_bar      = x.obs.conc_bar;
  obs.bar_conc      = x.obs.bar_conc;
  obs.count         = x.obs.count;
  obs.hires_conc    = x.obs.hires_conc;
  obs.weather_count = x.obs.weather_count;
}
ssmipt & ssmipt::operator=(ssmipt &x) {
  obs.t19v = x.obs.t19v;
  obs.t19h = x.obs.t19h;
  obs.t22v = x.obs.t22v;
  obs.t37v = x.obs.t37v;
  obs.t37h = x.obs.t37h;
  obs.t85v = x.obs.t85v;
  obs.t85h = x.obs.t85h;
  obs.conc_bar      = x.obs.conc_bar;
  obs.bar_conc      = x.obs.bar_conc;
  obs.count         = x.obs.count;
  obs.hires_conc    = x.obs.hires_conc;
  obs.weather_count = x.obs.weather_count;
  return *this;
}
// Added 4 February 2009, after versions in 28 Jul 2000 experiment
void ssmipt::show() {
  printf("ssmipt %5d %5d %5d %5d %5d %5d %5d   %3d %3d\n",
    obs.t19v, obs.t19h,obs.t22v, obs.t37v, obs.t37h, obs.t85v, obs.t85h,
    obs.conc_bar, obs.bar_conc);
}
int ssmipt::qc() {
  int ret=0;
// If not a data point, don't worry
  if (obs.conc_bar == LAND ||
      obs.conc_bar == WEATHER ||
      obs.conc_bar == COAST ||
      obs.conc_bar == NO_DATA ||
      obs.bar_conc == LAND ||
      obs.bar_conc == WEATHER ||
      obs.bar_conc == COAST ||
      obs.bar_conc == NO_DATA ) {
    return ret;
  }
  if (obs.t19v < 50*100 || obs.t19v > 300*100) {
    ret += 1;
  }
  if (obs.t19h < 50*100 || obs.t19h > 300*100) {
    ret += 1;
  }
  if (obs.t22v < 50*100 || obs.t22v > 300*100) {
    ret += 1;
  }
  if (obs.t37v < 50*100 || obs.t37v > 300*100) {
    ret += 1;
  }
  if (obs.t37h < 50*100 || obs.t37h > 300*100) {
    ret += 1;
  }
  if (obs.t85v < 50*100 || obs.t85v > 300*100) {
    ret += 1;
  }
  if (obs.t85h < 50*100 || obs.t85h > 300*100) {
    ret += 1;
  }
  
  return ret;
}

point3<float> ssmipt::gradient() {
  point3<float> x;
  x.i = (float) obs.t22v / (float) obs.t19v;
  x.j = (float) obs.t37v / (float) obs.t19v;
  x.k = (float) obs.t85v / (float) obs.t19v;
  return x;
}
point3<float> ssmipt::polarization() {
  point3<float> x;
  x.i = (float) obs.t19v / (float) obs.t19h;
  x.j = (float) obs.t37v / (float) obs.t37h;
  x.k = (float) obs.t85v / (float) obs.t85h;
  return x;
}
point3<float> ssmipt::differences() {
  point3<float> x;
  x.i = (float) obs.t19v - (float) obs.t19h;
  x.j = (float) obs.t37v - (float) obs.t37h;
  x.k = (float) obs.t85v - (float) obs.t85h;
  return x;
}


#endif