LandUseTable.h

Go to the documentation of this file.

#ifndef _LANDUSETABLE_ 

#define _LANDUSETABLE_ 
class TimeHandler;
#include "CommonTable.h"

// Extracted from first part of WASIM-file modef.h

#include <iostream>
#include <iomanip>
#include <cstdio>
#include <vector>
#include <string.h>


#include"MatrixDouble.h"        

using namespace std; 

#define _MAXCHAR_ 200

typedef void (*functionpointer) (int);



struct KcMonthlyValuesStruct
{
  MatrixDouble kcForAllMonth; 
  int code;
  char nameAct[100];        
  float cutoffFactor;
};



class LUSEtabletype:
public CommonTableType
{
        
private:
  int method;
  
        int ID_cropwatTable;
        bool flag_useCropwatTable;
        KcMonthlyValuesStruct kcMonthly;
        
  int iNumFixPoints;
  int *days;          /* St?tzstellen der Jahresgang-Kurve */
  int *daystemp;        /* nach H?he korrigierten aktuellen St?tzstellen */
  float *rsc,     /* minimal canopy-resistance der Landnutzung als Monatswerte */
            *rsi,     /* minimal interception-resistance der Landnutzung als Monatswerte */
                *rse,     /* minimal soil-resistance der Landnutzung als Monatswerte */
                *albedo,  /* albedo if no snow is covering the crops                     */
        *LAI,     /* Blattflaechenindizes an den 4 Knickpunkten des Jahresganges */
        *z0,      /* effektive Bewuchshoehe der Landnutzung zu den Knickpunkten  */
        *Vbg,     /* Vegetationsbedeckung der Landnutzung an den Knickpunkten    */
        *Wurz,    /* Wurzeltiefe der Landnutzung                                 */
                *AltDep;  /* Altitude dependency in days per m                           */
  float form_par, /* Wurzeldichte fuer Entnahmefunktion 0...                     */
        redu_parDry, /* ab welcher Bodenfeuchte im Trockenbereich beginnt Verdunstungsreduktion einlesen als suction in m, Umrechnen auf absolten Wassergehalt Theta */
                redu_parWet, /* ab welcher Bodenfeuchte im Nassbereich beginnt Verdunstungsreduktion, einlesen als relativer Wassergehalt (0..1), Umrechnen auf absoluten Wassergehalt Theta */
                LimitReduWet; /* minimal relative reduction of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0) */
  float m_fIntercepStor; // specific interception storage capacity in mm -> will be multiplied by LAI to get the effektive storage capacity
  int mon[14];
  
public:
  LUSEtabletype(int index);
  LUSEtabletype(); // defaultconstructor
  ~LUSEtabletype();
private:
  float interpolate(double x1,double y1,double x2,double y2,double x);
  int GetJulDay(int iMonthDay,int iMonth,int iYear);
  void AltitudeCorrection(float dhm);
  int GetFixPoint(int iJulDay,float height);
  void stripParameters(string &row,string sWhatParam,float *param,int maxIndex);
  void stripParametersJulDays(string &row);

public:

  void initialize(ostream *WSstream,string row,int version);
  void initialize(ostream *WSstream,string row,string sVersion);

  void initializeKcEmpty();
  void deallocateKcStruct();
  
  float getFormPar();
  float getReduParDry();
  float getReduParWet();
  float getLimitReduWet();
  float getIntercepStor();
  char* getName();
  int   getCode();
  void  setCode(int code);
  void  printContent();

  bool  getIsActivity();
  void printKcMonthlyValuesStruct();
  int  getActivityNumberInMPMAS();
  void extractActivityNumberFromName();
  void readKcMonthlyValuesStruct(FILE* steuer);
  double getMonthlyKcValue(int m);
  double getCutoffFactor() {return (double) kcMonthly.cutoffFactor;};
  
  
  float GetRsc(int iJulDay, float height=0.0);
  float GetRsc(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetRsi(int iJulDay, float height=0.0);
  float GetRsi(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetRse(int iJulDay, float height=0.0);
  float GetRse(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetZ0(int iJulDay, float height=0.0);
  float GetZ0(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetLAI(int iJulDay, float height=0.0);
  float GetLAI(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetVbg(int iJulDay, float height=0.0);
  float GetVbg(int iMonthDay,int iMonth,int iYear, float height=0.0);

  float GetWurz(int iJulDay, float height=0.0);
  float GetWurz(int iMonthDay,int iMonth,int iYear, float height=0.0);
  
  float GetAlbedo(int iJulDay, float height);
  float GetAlbedo(int iMonthDay,int iMonth, int iYear, float height=0.0);

};


void GetLuseTable(int version,  string , LUSEtabletype*& LUtable);
void GetLuseTable(int version,  LUSEtabletype*& LUtable);
//void GetReferenceTable(int version, string fn,        LUSEtabletype *&LUtable);

#endif

---