int MDRainSMoistChgDef () {
int ret = 0;
float par;
const char *optStr;
if (_MDOutSMoistChgID != MFUnset) return (_MDOutSMoistChgID);
const char *soilTemperatureOptions [] = { "none", "calculate", (char *) NULL };
int soilTemperatureID;
if (((optStr = MFOptionGet (MDOptSoilTemperature)) == (char *) NULL) ||
((soilTemperatureID = CMoptLookup (soilTemperatureOptions, optStr, true)) == CMfailed)) {
CMmsgPrint(CMmsgUsrError," Soil TemperatureOption not specifed! Options = 'none' or 'calculate'\n");
return CMfailed;
}
if (((optStr = MFOptionGet (MDParSoilMoistALPHA)) != (char *) NULL) && (sscanf (optStr,"%f",&par) == 1)) _MDSoilMoistALPHA = par;
MFDefEntering ("Rainfed Soil Moisture");
if (soilTemperatureID == MFcalculate ) {
if (((ret = MDPermafrostDef()) == CMfailed) ||
((_MDInRelativeIceContent = MFVarGetID ("SoilIceContent_01", "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutLiquidSoilMoistureID = MFVarGetID (MDVarLiquidSoilMoisture, "-", MFOutput, MFState, MFBoundary)) == CMfailed)) return CMfailed;
}
if ((ret = MDIrrGrossDemandDef ()) == CMfailed) return (CMfailed);
if ((ret != MFUnset) && ((_MDInIrrAreaFracID = MDIrrigatedAreaDef ())== CMfailed) ) return (CMfailed);
if (((_MDInPrecipID = MDPrecipitationDef ()) == CMfailed) ||
((_MDInSPackChgID = MDSPackChgDef ()) == CMfailed) ||
((_MDInPotETID = MDRainPotETDef ()) == CMfailed) ||
((_MDInInterceptID = MDRainInterceptDef ()) == CMfailed) ||
((_MDInSoilAvailWaterCapID = MDSoilAvailWaterCapDef ()) == CMfailed) ||
((_MDInAirTMeanID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutEvaptrsID = MFVarGetID (MDVarRainEvapotranspiration, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDOutSoilMoistCellID = MFVarGetID (MDVarRainSoilMoistureCell, "mm", MFOutput, MFState, MFInitial)) == CMfailed) ||
((_MDOutSoilMoistID = MFVarGetID (MDVarRainSoilMoisture, "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutSMoistChgID = MFVarGetID (MDVarRainSoilMoistChange, "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
(MFModelAddFunction (_MDRainSMoistChg) == CMfailed)) return (CMfailed);
MFDefLeaving ("Rainfed Soil Moisture");
return (_MDOutSMoistChgID);
}
int MDRainSMoistChgDef () {
int ret = 0;
float par;
const char *optStr;
/*
// Test if the percolation pathway is active (if its not - PercolationBETA should not be in the Options)
if ((optStr = MFOptionGet(MDParSoilPercolationBETA)) != (char *) NULL) {
if ((_MDOutPercolationID = MFVarGetID (MDVarSoilPercolation, "mm", MFInput, MFFlux, MFBoundary)) == CMfailed) return (CMfailed); // SZ 10212014
}
*/
if (_MDOutSMoistChgID != MFUnset) return (_MDOutSMoistChgID);
const char *soilTemperatureOptions [] = { "none", "calculate", (char *) NULL };
int soilTemperatureID;
if (((optStr = MFOptionGet (MDOptSoilTemperature)) == (char *) NULL) ||
((soilTemperatureID = CMoptLookup (soilTemperatureOptions, optStr, true)) == CMfailed)) {
CMmsgPrint(CMmsgUsrError," Soil TemperatureOption not specifed! Options = 'none' or 'calculate'\n");
return CMfailed;
}
if (((optStr = MFOptionGet (MDParSoilMoistALPHA)) != (char *) NULL) && (sscanf (optStr,"%f",&par) == 1)) _MDSoilMoistALPHA = par;
if (((optStr = MFOptionGet (MDParSoilPercolationBETA)) != (char *) NULL) && (sscanf (optStr,"%f",&par) ==1)) _MDparPercolation = par;
MFDefEntering ("Rainfed Soil Moisture");
if (soilTemperatureID == MFcalculate ) {
if (((ret = MDPermafrostDef()) == CMfailed) ||
((_MDInRelativeIceContent = MFVarGetID ("SoilIceContent_01", "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutLiquidSoilMoistureID = MFVarGetID (MDVarLiquidSoilMoisture, "-", MFOutput, MFState, MFBoundary)) == CMfailed)) return CMfailed;
}
if ((ret = MDIrrGrossDemandDef ()) == CMfailed) return (CMfailed);
if ((ret != MFUnset) && ((_MDInIrrAreaFracID = MDIrrigatedAreaDef ())== CMfailed) ) return (CMfailed);
if (((_MDInPrecipID = MDPrecipitationDef ()) == CMfailed) ||
((_MDInSPackChgID = MDSPackChgDef ()) == CMfailed) ||
((_MDInPotETID = MDRainPotETDef ()) == CMfailed) ||
((_MDInInterceptID = MDRainInterceptDef ()) == CMfailed) ||
((_MDInSoilAvailWaterCapID = MDSoilAvailWaterCapDef ()) == CMfailed) ||
((_MDInAirTMeanID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInRunoffToPervID = MDStormRunoffDef ()) == CMfailed) || // RJS 082812
// ((_MDInH2OFractionID = MFVarGetID (MDVarH2OFracSpatial, "mm", MFInput, MFState, MFBoundary)) == CMfailed) || // RJS 082812
((_MDInStormRunoffTotalID = MFVarGetID (MDVarStormRunoffTotal, "mm", MFInput, MFFlux, MFBoundary)) == CMfailed) || // RJS 082812
((_MDInImpFractionID = MFVarGetID (MDVarImpFracSpatial, "-", MFInput, MFState, MFBoundary)) == CMfailed) || // RJS 082812
((_MDInH2OFractionID = MFVarGetID (MDVarH2OFracSpatial, "-", MFInput, MFState, MFBoundary)) == CMfailed) || // RJS 082812
((_MDOutEvaptrsID = MFVarGetID (MDVarRainEvapotranspiration, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDOutExcessID = MFVarGetID (MDVarExcess, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) || // RJS 091813
((_MDOutExcessNotScaledID = MFVarGetID (MDVarExcessNotScaled, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) || // RJS 091813
((_MDOutEvaptrsNotScaledID = MFVarGetID (MDVarRainETnotScaled, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) || // RJS 082812
((_MDOutSoilMoistNotScaledID = MFVarGetID (MDVarRainSoilMoistureNotScaled, "mm", MFOutput, MFState, MFInitial)) == CMfailed) ||
((_MDOutSoilMoistID = MFVarGetID (MDVarRainSoilMoisture, "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutSMoistChgNotScaledID = MFVarGetID (MDVarRainSoilMoistureChangeNotScaled, "mm", MFOutput, MFState, MFInitial)) == CMfailed) ||
((_MDOutSMoistChgID = MFVarGetID (MDVarRainSoilMoistChange, "mm", MFOutput, MFState, MFBoundary)) == CMfailed) ||
(MFModelAddFunction (_MDRainSMoistChg) == CMfailed)) return (CMfailed);
const char *snowImpervMeltOptions [] = { MDNoneStr, MDCalculateStr,(char *) NULL };
int impervSnowMeltCalcID;
if ((optStr = MFOptionGet (MDOptImperviousMeltCalc) ) == (char *) NULL) {
optStr = MDNoneStr;
//CMmsgPrint(CMmsgWarning," Impervious Snow Fall runoff method not specified - defaulting to none (RainSMoistChgDef).\n");
}
if ((impervSnowMeltCalcID = CMoptLookup (snowImpervMeltOptions,optStr,true)) == CMfailed) {
CMmsgPrint(CMmsgUsrError," Impervious Snow Fall runoff method incorrectly specified. Options are 'none' or 'calculate'.\n");
return (CMfailed);
}
switch (impervSnowMeltCalcID) {
case MFnone:
// Nothing to be done
break;
case MFcalculate:
if (((_MDInImpSnowFallROID = MFVarGetID ( MDVarImpSnowFallRunoff,"mm", MFInput, MFFlux, MFBoundary)) == CMfailed)
) return (CMfailed);
break;
default: MFOptionMessage (MDOptImperviousMeltCalc, optStr, snowImpervMeltOptions); return (CMfailed);
}
MFDefLeaving ("Rainfed Soil Moisture");
return (_MDOutSMoistChgID);
}
