CropWater.h

Go to the documentation of this file.

//------------------------------------------------------------------------------
//
//  Thomas Berger (main author), Pepijn Schreinemachers, and Thorsten Arnold
//
//
//  Hohenheim University (490d)
//  Professorship for Land Use Economics in the Tropics and Subtropics
//
//  May not be distributed without permission by main author
//
//  File:       CropWater.h (formerly part of agentF.h)
//
//  Contents:   Implemenation of the water-stress crop growth model CROPWAT
//
//
//
//------------------------------------------------------------------------------

#ifndef _CropWater_h
#define _CropWater_h

#include "BasicData.h"
#include "MatrixDouble.h"

#ifndef LP_OLD
        #include "AgentSolution.h"
#else
        #include "lp9.h"
#endif

class InfoHydro{  
public:

        int getColumnInLP()const;
        int getPriceColumnInLP()const;
        int getYieldRowInLP()const;
        int getSoilTypeID()const;
        int getIrrigationPriorityGroup()const;
        int getCropIrrigationType()const;
        int getIrrigationMethod()const;
        int getFormOfStressCoefficientKR()const;
        double getWaterStressCoefficientKY()const;
        double getCropYieldPotential()const;
        double getCropYieldStartValue()const;
        double getPlantWaterDemandsGrossInLiterPerSecond(int m)const;
        double getWaterDemandsGrossInOriginalUnit(int m)const;
        //double getWaterDemandsExpectedReducedByEffectiveRainfall(int m)const;//## H2O
        

        void copyRhsToThis(const InfoHydro& rhs);
        void printToScreen();
        void writeToFile(ofstream& outStream);
        bool readFromFile(ifstream& fileStreamCrop, int numIrrigationPeriods);
        void allocate();
        void initialize(double initVal);
        void write(FILE* strm);


        InfoHydro();
        InfoHydro(double initVal);
        InfoHydro(const InfoHydro&);
        InfoHydro& operator=(const InfoHydro&);
        ~InfoHydro();

private:
        /* TODO 1 -oBerger -cImplement : LP info related to WASIM */
        int columnInLP;
        int priceColumnInLP;
        int yieldRowInLP;
        
        //int soilClassID; ## soil
        int soilTypeID;

        int irrigationPriorityGroup;
        int cropIrrigationType;
        int irrigationMethod;
        int formOfStressCoefficientKR;
        double waterStressCoefficientKY;
        double cropYieldPotential;
        double cropYieldStartValue;
        MatrixDouble waterDemandsGross;
        //MatrixDouble waterDemandsReducedByExpectedEffectiveRainfall;//## h2o
};



class IrrigationMethodSector;
 
class CropWater
{
public:

        virtual double computeCropEffectivePrecipitation(double actualRainfall, double monthlyWaterDemand);

        //virtual void computeCropEffectiveWaterDemandsInMonth //substitute
        virtual void computeCropEffectivePrecipitationInMonth //Tbe 080220
                         (int m, double actualRainfall, 
                         MatrixDouble& vectorOfMonthlyCropEffectivePrecipitation);
        virtual void computeCropIrrigationWaterDemandsInMonth //Tbe 080220
                        (int m,   int secID_,
                        const IrrigationMethodSector& sectorparamsIrrigationMethod, 
                        const MatrixDouble& effectivePrecipitation,
                        MatrixDouble& vectorOfMonthlyCropIrrigationWaterDemands);

        //virtual double computeIrrigationWaterDemandPerHectareForMonth //substitute
        virtual double computeIrrigatedAreaPerMethodAndTotalIrrigationWaterDemand(int m, //Tbe 080220
                        const MatrixDouble& cropWaterData,
                        const MatrixDouble*& irrigationWaterDemands,
                        MatrixDouble& vectorOfIrrigatedAreaByMethod);

        //virtual double allocateIrrigationWaterToCropsPerMonth //substitute
        virtual double allocateIrrigationWaterToCropsInMonth //Tbe 080220
                        (int m,
                        double*& D_sat,
                        const  MatrixDouble*& irrigationWaterDemands,
                        const MatrixDouble& cropWaterData,
                        MatrixDouble& irrigationWaterDemandsPerPriorityGroup,
                        MatrixDouble& percentageIrrigatedPerPriorityGroup,
                        MatrixDouble& irrigationWaterSuppliedToCrops);

        //virtual void computeWaterDeficitForHydroActivity //substitute
        virtual void computeWaterDeficitsOfAllCrops //Tbe 080220
                        (const MatrixDouble& cropWaterData,
                        const MatrixDouble*& effectivePrecipitation,
                        const MatrixDouble& irrigationWaterSuppliedToCrops,
                        MatrixDouble& waterDeficitsOfCropsAsMonthlyKR);
        virtual void computeYieldReductionUnderWaterStress //Tbe 080220
                        (MatrixDouble& cropWaterData,
                        MatrixDouble& waterDeficitsOfCropsAsMonthlyKR);


        //---------------------------------------------------------------------------



        virtual int getNumWaterActivities();
        virtual int getNumIrrigationPriorityGroups();
        virtual int getNumWaterMonths();
        
        virtual int getLpColumnForActivityID(int _actID);
        virtual int getIrrigationTypeForActivityID(int _actID);
        
        virtual int getActivityIDForLpColumn(int _colLP);
        virtual int getIrrigationTypeForLpColumn(int _colLP);

   virtual double getPlantWaterDemandsGrossInLiterPerSecond(int _hydroActID, int m_);
        virtual int getSoilTypeID(int _hydroActID);
        virtual InfoHydro getInfoOnActivityCropHydrology(int arrayIndex);
        const InfoHydro* getPointerToInfoOnActivityCropHydrology(int arrayIndex);
        
        //---------------------------------------------------------------------------
        


        virtual void readInputFromFile(string fn);
        virtual void writeInputIntoFile();

        

        //---------------------------------------------------------------------------



        CropWater();

        //CropWater(const CropWater& cropWater_);

        //CropWater& operator=(const CropWater& cropWater_);

        virtual ~CropWater()
        {       delete [] waterActivities;
                waterActivityLookup.deallocate();
        }
        
        
protected:


        int numWaterActivities;
        
        int numIrrigationPriorityGroups;

        int numWaterMonths;

        double effRainCoeffConstant;            // 1.2995160
        double effRainCoeffDemand;                      //-0.0348437
        double effRainCoeffDemandSquared;// 0.0001663
        double effRainCoeffRain;                        // 0.6439529
        double effRainCoeffRainSquared; //-0.0004772
        double effRainCoeffCross;                       // 0.0010454    
        
        /* TODO 1 -oBerger -cImplement : LP info related to WASIM */
        InfoHydro* waterActivities;//formerly 'actHydro' or 'actHydroVector'


        MatrixDouble waterActivityLookup;




        //double growthRedFunction(double ky, double kr);//substitute
        double computeCropWaterYieldReductionFactor(double ky, double kr);//Tbe 080220
        
};


#endif

---