#ifndef GLM_POLITICAL_INPUTS_H
#define GLM_POLITICAL_INPUTS_H

#include <vector>
#include <math.h>
#include "config.h"
#include "grid_cell.h"

namespace glm {

   inline
     double fround2(double n, unsigned d){
     return floor(n*pow(10., d) + .5) / pow(10., d); 
   }

   class Zone;
   typedef std::vector<Zone> ZoneList;
   
   template<class T> class GriddedInput;

   

   class GridInfo {

   public:
      GridInfo (ZoneList& zones, Config& config);

      void readCellInfo (Config& config);
      void readCountryCodes (ZoneList& countries, Config& config);
      void readShiftCult (Config& config);
      void readContinentCodes (ZoneList& countries, Config& config);
      void readRegionCodes (ZoneList& countries, Config& config);
      void readOtherData (Config& config);
      void readCropTypesData (Config& config);
      void readFertilizerData (ZoneList& countries, Config& config);
      void readFloodedFractionData (Config& config);
      
      double lat[NY];
      double lon[NX];
      double cellArea[NY][NX];
      double cellArea2Deg[AGG_NY][AGG_NX];
      
      GriddedInput<double>* wtr_;
      GriddedInput<double>* ice_;
      double wtr2Deg_[AGG_NY][AGG_NX];
      double ice2Deg_[AGG_NY][AGG_NX];
      
      double forested2Deg_[AGG_NY][AGG_NX];

      int shiftCult[NY][NX];
      int countryCode[NY][NX];
      int countryIndex[NY][NX];
      int continentCode[NY][NX];
      double vrgnBiomassAboveGround[NY][NX];
      double vrgnNPPAboveGround[NY][NX];
      double future_pasture_fracs[NY][NX];
      int forested[NY][NX];
      double forestedFraction[NY][NX];
      double C3_annual_input_2000[NY][NX];
      double C4_annual_input_2000[NY][NX];
      double C3_perennial_input_2000[NY][NX];
      double C4_perennial_input_2000[NY][NX];
      double N_fixing_input_2000[NY][NX];
      double flooded_input_2000[NY][NX];
      double flooded_ratio[NY][NX];
      double crop2015[NY][NX];
      double crop2100[NY][NX];
      double SC2099[NY][NX];
      double rs_forest_loss_2012[NY][NX];

      GriddedInput<double>* restartC3Ann_;
      GriddedInput<double>* restartC4Ann_;
      GriddedInput<double>* restartC3Per_;
      GriddedInput<double>* restartC4Per_;
      GriddedInput<double>* restartC3NFx_;
      GriddedInput<double>* restartC3Ann2015_;
      GriddedInput<double>* restartC4Ann2015_;
      GriddedInput<double>* restartC3Per2015_;
      GriddedInput<double>* restartC4Per2015_;
      GriddedInput<double>* restartC3NFx2015_;
      GriddedInput<double>* restartPastr_;
      GriddedInput<double>* restartRange_;
      GriddedInput<double>* restartScndF_;
      GriddedInput<double>* restartScndN_;
      GriddedInput<double>* restartPrimF_;
      GriddedInput<double>* restartPrimN_;
      GriddedInput<double>* restartUrban_;
      GriddedInput<double>* restartSMA_;
      GriddedInput<double>* restartSMB_;
      GriddedInput<double>* restartC3AnnFert_;
      GriddedInput<double>* restartC4AnnFert_;
      GriddedInput<double>* restartC3PerFert_;
      GriddedInput<double>* restartC4PerFert_;
      GriddedInput<double>* restartC3NfxFert_;
      GriddedInput<double>* restartC3AnnIrrig_;
      GriddedInput<double>* restartC4AnnIrrig_;
      GriddedInput<double>* restartC3PerIrrig_;
      GriddedInput<double>* restartC4PerIrrig_;
      GriddedInput<double>* restartC3NfxIrrig_;
      GriddedInput<double>* restartFlood_;
      GriddedInput<double>* restartFharvC3Per_;
      GriddedInput<double>* restartFharvC4Per_;
      GriddedInput<double>* restartCropbfC3Ann_;
      GriddedInput<double>* restartCropbfC4Ann_;
      GriddedInput<double>* restartCropbfC3Per_;
      GriddedInput<double>* restartCropbfC4Per_;
      GriddedInput<double>* restartCropbfC3Nfx_;
      GriddedInput<double>* restartRndwd_;
      GriddedInput<double>* restartFulwd_;
      GriddedInput<double>* restartCombf_;
      GriddedInput<double>* restartVBH1_;
      GriddedInput<double>* restartSBH1_;
      GriddedInput<double>* restartVBH2_;
      GriddedInput<double>* restartSBH2_;
      GriddedInput<double>* restartSBH3_;
      GriddedInput<double>* restartSC2099_;
   };

}

#endif // GLM_POLITICAL_INPUTS_H