#ifndef GLM_ZONE_H_
#define GLM_ZONE_H_

#include "grid_cell.h"

#define NZ 61 //31 // Number of z dims tracked (z is distance from managed cells)
#define NL 61  //31 // Number of L dims tracked (L is distance from Landsat cells)

class Config;

namespace glm {

   enum priority_t { VRGN_PRIORITY, SCND_PRIORITY };

    class Zone {
   public:
      void   resetTotals ();
      void   precomputeHarvestRatios (Config& config);
      void   calcConvertedForest ();
      void   calcVrgnForestHarvestRatios (int maxZ, double& remaining);
      void   calcScndForestHarvestRatios (int maxZ, double& remaining);


      char   countryName[128];
      int    countryCode;
      int    continentCode;
      int    regionCode;

      int harvestCase;
      priority_t priority;
      priority_t force_priority;
      bool   useZDistance;
      double   includeClearingInHarvest;
 
      // units MgC
      double harvestDemand;
      double fuelwoodFraction;
      double floodedChange;
      double harvestRemaining;
      double totalAvailBiomass;
      double totalHarvested;
      double newHarvestDemand;

      double vrgnForestBiomass[NZ];
      double vrgnForestHarvestRatio[NZ];
      double totalVrgnForestBiomass;
      double vrgnNonForestBiomass;
      double vrgnNonForestHarvestRatio;
      double scndZdisForestBiomass[NZ];
      double totalScndForestBiomass;
      double scndMatureZdisForestBiomass[NZ];
      double scndZdisForestHarvestRatio[NZ];      
      double totalScndMatureForestBiomass;
      double harvestedScndMatureForestBiomass;
      double harvestedScndYoungForestBiomass;
      double targetScndMatureForestBiomassHarvest;
      double targetScndYoungForestBiomassHarvest;
      double targetVrgnForestBiomass;
      double harvestedVrgnForestBiomass;
      double scndYoungForestBiomass;
      double scndYoungForestHarvestRatio;
      double scndNonForestBiomass;
      double scndNonForestHarvestRatio;


      double fh_sbh2;
      double fh_sbh3;
      double fh_vbh2;
      double predict_b;

      double converted_forest_land;
      double converted_forest_land_total;

      double fertilizerInput[N_SUB_LU_TYPES];

      double tillage;
      
      double c4ann_to_c3nfx;
      double c4ann_to_c3ann;
      double c3ann_to_c3nfx;

      double fertDataC3ANN;
      double fertDataC4ANN;
      double fertDataC3PER;
      double fertDataC4PER;
      double fertDataC3NFX;
      
      double fertDataC3ANNPrev;
      double fertDataC4ANNPrev;
      double fertDataC3PERPrev;
      double fertDataC4PERPrev;
      double fertDataC3NFXPrev;
      
      double irrigData;
      double irrigDataPrev;

      double cropBiofuelsC3ANN;
      double cropBiofuelsC4ANN;
      double cropBiofuelsC3PER;
      double cropBiofuelsC4PER;
      double cropBiofuelsC3NFX;

      double flowvc_prime;
      double flowsc_prime;
      double flowcs_prime;
      double flowvp_prime;
      double flowsp_prime;
      double flowps_prime;
      double flowvu_prime;
      double flowsu_prime;
      double flowus_prime;
      double reg_avail;
      double reg_unmet;
      double sum;
      
      // Crop functional type (CFT) fraction of grid cell ag area
      double c3AnnualFraction;
      double c3PerennialFraction;
      double c4AnnualFraction;
      double c4PerennialFraction;
      double NfixingFraction;
      
      double c3annTotal;
      double c4annTotal;
      double c3perTotal;
      double c4perTotal;
      double c3nfxTotal;
      double cropTotal;
      
      double c3annFertCurrent;
      double c4annFertCurrent;
      double c3perFertCurrent;
      double c4perFertCurrent;
      double c3nfxFertCurrent;
      
      double irrigCurrent;
      double irrigAvail;
      double irrigUnmet;
      double irrigChange;
      double irrigIncFrac1;
      double irrigIncFrac2;
      double irrigDecFrac;
      double irrigNew;
      double irrigCurrentNew;
      double irrigDefaultChange;
      double irrigNeededChange;
      
   };

} // end glm namespace

#endif // GLM_ZONE_H_