int MDRainPotETSWGdayDef () {
if (_MDOutPetID != MFUnset) return (_MDOutPetID);
MFDefEntering ("Rainfed Potential Evapotranspiration (Shuttleworth - Wallace [day])");
if (((_MDInDayLengthID = MDSRadDayLengthDef ()) == CMfailed) ||
((_MDInI0HDayID = MDSRadI0HDayDef ()) == CMfailed) ||
((_MDInCParamAlbedoID = MDCParamAlbedoDef ()) == CMfailed) ||
((_MDInCParamCHeightID = MDCParamCHeightDef ()) == CMfailed) ||
((_MDInCParamLWidthID = MDCParamLWidthDef ()) == CMfailed) ||
((_MDInCParamRSSID = MDCParamRSSDef ()) == CMfailed) ||
((_MDInCParamR5ID = MDCParamR5Def ()) == CMfailed) ||
((_MDInCParamCDID = MDCParamCDDef ()) == CMfailed) ||
((_MDInCParamCRID = MDCParamCRDef ()) == CMfailed) ||
((_MDInCParamGLMaxID = MDCParamGLMaxDef ()) == CMfailed) ||
((_MDInCParamZ0gID = MDCParamZ0gDef ()) == CMfailed) ||
((_MDInLeafAreaIndexID = MDLeafAreaIndexDef ()) == CMfailed) ||
((_MDInStemAreaIndexID = MDStemAreaIndexDef ()) == CMfailed) ||
((_MDInSolRadID = MDSolarRadDef ()) == CMfailed) ||
((_MDInAtMeanID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInAtMinID = MFVarGetID (MDVarAirTempMinimum, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInVPressID = MFVarGetID (MDVarVaporPressure, "kPa", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInWSpeedID = MFVarGetID (MDVarWindSpeed, "m/s", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutPetID = MFVarGetID (MDVarRainPotEvapotrans, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
(MFModelAddFunction (_MDRainPotETSWGday) == CMfailed)) return (CMfailed);
MFDefLeaving ("Rainfed Potential Evapotranspiration (Shuttleworth - Wallace [day])");
return(_MDOutPetID);
}
int MDWTempRiverRouteDef () {
int optID = MFUnset, waterBalanceID;
const char *optStr;
const char *options [] = { "none", "calculate", (char *) NULL };
if (_MDWTempRiverRouteID != MFUnset) return (_MDWTempRiverRouteID);
MFDefEntering ("Route river temperature");
if (((optStr = MFOptionGet (MDOptReservoirs)) == (char *) NULL) ||
((optID = CMoptLookup ( options, optStr, true)) == CMfailed)) {
CMmsgPrint(CMmsgUsrError,"Reservoir Option not specified! Option none or calculate");
return CMfailed;
}
if (optID==1) {
if ((( waterBalanceID = MDWaterBalanceDef ()) == CMfailed) ||
((_MDInResStorageChangeID = MFVarGetID (MDVarReservoirStorageChange, "km3", MFInput, MFState, MFBoundary)) == CMfailed) || //RJS 071511
((_MDInResStorageID = MFVarGetID (MDVarReservoirStorage, "km3", MFInput, MFState, MFInitial)) == CMfailed) || //RJS 121311 changed from MFBoundary to MFInitial
((_MDInResCapacityID = MFVarGetID (MDVarReservoirCapacity, "km3", MFInput, MFState, MFBoundary)) == CMfailed))
return (CMfailed);
}
//input
if (((_MDInDischargeID = MDDischargeDef ()) == CMfailed) ||
(( waterBalanceID = MDWaterBalanceDef ()) == CMfailed) ||
((_MDInSolarRadID = MDSolarRadDef ()) == CMfailed) ||
((_MDInWTempRiverID = MDWTempRiverDef ()) == CMfailed) ||
((_MDInRiverWidthID = MDRiverWidthDef ()) == CMfailed) ||
((_MDInRunoffVolumeID = MDRunoffVolumeDef ()) == CMfailed) ||
((_MDInWdl_QxTID = MDThermalInputsDef ()) == CMfailed) || //RJS 072011
((_MDInDischargeIncomingID = MFVarGetID (MDVarDischarge0, "m3/s", MFInput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDInWindSpeedID = MFVarGetID (MDVarWindSpeed, "m/s", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInAirTemperatureID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInCloudCoverID = MFVarGetID (MDVarCloudCover, "%", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInRiverStorageChgID = MFVarGetID (MDVarRiverStorageChg, "m3/s", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInRiverStorageID = MFVarGetID (MDVarRiverStorage, "m3", MFInput, MFState, MFInitial)) == CMfailed) ||
((_MDInSnowPackID = MFVarGetID (MDVarSnowPack, "mm", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInWarmingTempID = MFVarGetID (MDVarWarmingTemp, "degC", MFInput, MFState, MFBoundary)) == CMfailed) || //RJS 072011
((_MDInThermalWdlID = MFVarGetID (MDVarThermalWdl, "-", MFInput, MFState, MFBoundary)) == CMfailed) || //RJS 072011
// output
((_MDLocalIn_QxTID = MFVarGetID (MDVarBgcLocalIn_QxT, "m3*degC/d", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDRemoval_QxTID = MFVarGetID (MDVarRemoval_QxT, "m3*degC/d", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDFlux_QxTID = MFVarGetID (MDVarFlux_QxT, "m3*degC/d", MFRoute, MFFlux, MFBoundary)) == CMfailed) ||
((_MDStorage_QxTID = MFVarGetID (MDVarStorage_QxT, "m3*degC", MFOutput, MFState, MFInitial)) == CMfailed) ||
((_MDDeltaStorage_QxTID = MFVarGetID (MDVarDeltaStorage_QxT, "m3*degC/d", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDWTemp_QxTID = MFVarGetID (MDVarWTemp_QxT, "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||
((_MDWTempDeltaT_QxTID = MFVarGetID (MDVarWTempDeltaT_QxT, "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||
//output mixing
((_MDFluxMixing_QxTID = MFVarGetID (MDVarFluxMixing_QxT, "m3*degC/d", MFRoute, MFFlux, MFBoundary)) == CMfailed) ||
((_MDStorageMixing_QxTID = MFVarGetID (MDVarStorageMixing_QxT, "m3*degC", MFOutput, MFState, MFInitial)) == CMfailed) ||
((_MDDeltaStorageMixing_QxTID = MFVarGetID (MDVarDeltaStorageMixing_QxT, "m3*degC/d", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
((_MDWTempMixing_QxTID = MFVarGetID (MDVarWTempMixing_QxT, "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||
(MFModelAddFunction (_MDWTempRiverRoute) == CMfailed)) return (CMfailed);
MFDefLeaving ("Route river temperature");
return (_MDWTemp_QxTID);
}
int MDRainPotETPsTaylorDef () {
if (_MDOutPetID != MFUnset) return (_MDOutPetID);
MFDefEntering ("Rainfed Potential Evapotranspiration (Priestley - Taylor)");
if (((_MDInDayLengthID = MDSRadDayLengthDef ()) == CMfailed) ||
((_MDInI0HDayID = MDSRadI0HDayDef ()) == CMfailed) ||
((_MDInCParamAlbedoID = MDCParamAlbedoDef ()) == CMfailed) ||
((_MDInSolRadID = MDSolarRadDef ()) == CMfailed) ||
((_MDInAtMeanID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInVPressID = MFVarGetID (MDVarVaporPressure, "kPa", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutPetID = MFVarGetID (MDVarRainPotEvapotrans, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed) ||
(MFModelAddFunction (_MDRainPotETPsTaylor) == CMfailed)) return (CMfailed);
MFDefLeaving ("Rainfed Potential Evapotranspiration (Priestley - Taylor)");
return (_MDOutPetID);
}
int MDIrrRefEvapotransFAODef () {
if (_MDOutIrrRefEvapotransID != MFUnset) return (_MDOutIrrRefEvapotransID);
MFDefEntering ("Irrigation Reference ETP (FAO)");
if (((_MDInDayLengthID = MDSRadDayLengthDef ()) == CMfailed) ||
((_MDInI0HDayID = MDSRadI0HDayDef ()) == CMfailed) ||
((_MDInSolRadID = MDSolarRadDef ()) == CMfailed) ||
((_MDInElevationID = MFVarGetID (MDVarMeanElevation, "m", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInAtMeanID = MFVarGetID (MDVarAirTemperature, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInAtMinID = MFVarGetID (MDVarAirTempMinimum, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInAtMaxID = MFVarGetID (MDVarAirTempMaximum, "degC", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInVPressID = MFVarGetID (MDVarVaporPressure, "kPa", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDInWSpeedID = MFVarGetID (MDVarWindSpeed, "m/s", MFInput, MFState, MFBoundary)) == CMfailed) ||
((_MDOutIrrRefEvapotransID = MFVarGetID (MDVarIrrRefEvapotrans, "mm", MFOutput, MFFlux, MFBoundary)) == CMfailed)) return (CMfailed);
if (MFModelAddFunction (_MDIrrRefEvapotransFAO)== CMfailed) return (CMfailed);
MFDefLeaving ("Irrigation Reference ETP (FAO)");
return(_MDOutIrrRefEvapotransID);
}
