bool CWaterRetentionCapacity::On_Execute(void){
int i,j;
int x,y;
int iField;
int iShape;
int iRows;
float fValue = 0;
float **pData;
int iX, iY;
float fC;
double dSlope,dAspect;
CSG_Shape* pShape;
CSG_Shapes* pShapes = Parameters("SHAPES")->asShapes();
CSG_Grid* pDEM = Parameters("DEM")->asGrid();
m_pRetention = Parameters("RETENTION")->asGrid();
m_pSlope = SG_Create_Grid(pDEM);
m_pOutput = Parameters("OUTPUT")->asShapes();
m_pOutput->Assign(pShapes);
m_pOutput->Add_Field("CCC", SG_DATATYPE_Double);
m_pOutput->Add_Field("CIL", SG_DATATYPE_Double);
m_pOutput->Add_Field(_TL("Permeability"), SG_DATATYPE_Double);
m_pOutput->Add_Field(_TL("Equivalent Moisture"), SG_DATATYPE_Double);
m_pOutput->Add_Field(_TL("Water Retention Capacity"), SG_DATATYPE_Double);
for(y=0; y<Get_NY() && Set_Progress(y); y++){
for(x=0; x<Get_NX(); x++){
if( pDEM->Get_Gradient(x, y, dSlope, dAspect) ){
m_pSlope->Set_Value(x, y, dSlope);
}
else{
m_pSlope->Set_NoData(x, y);
}
}
}
iRows = pShapes->Get_Field_Count() / 5;
pData = new float*[iRows];
for (iShape = 0; iShape < pShapes->Get_Count(); iShape++){
pShape = pShapes->Get_Shape(iShape);
for (i = 0; i< iRows; i++){
pData[i] = new float[5];
for (j = 0; j < 5; j++){
pData[i][j] = 0;
try{
pData[i][j] = pShape->asFloat(j+i*5);
}//try
catch(...){}
}//for
}//for
iX = (int)((pShape->Get_Point(0).x - pDEM->Get_XMin())/pDEM->Get_Cellsize());
iY = (int)((pShape->Get_Point(0).y - pDEM->Get_YMin())/pDEM->Get_Cellsize());
fC = (float)(1. - tan(m_pSlope->asFloat(iX,iY,false)));
pShape = m_pOutput->Get_Shape(iShape);
CalculateWaterRetention(pData, iRows, fC, pShape);
}//for
iField = m_pOutput->Get_Field_Count()-1;
CIDW IDW;
IDW.setParameters(m_pRetention, m_pOutput, iField);
IDW.Interpolate();
CorrectWithSlope();
return true;
}//method
