自学教程：C++ CMoptLookup函数代码示例

int MDWetlandAreaDef (){		int  optID = MFUnset;	const char *optStr, *optName = MDVarWetlandAreaFraction;	const char *options [] = { MDNoneStr, MDInputStr, MDCalculateStr, (char *) NULL };    if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	if ((optID == MDnone) || (_MDWetlandAreaFracID != MFUnset)) return (_MDWetlandAreaFracID);	MFDefEntering ("WetlandArea");//	printf ("Wetland Area optID  =%i  /n",optID);	if (_MDWetlandAreaFracID != MFUnset) return (_MDWetlandAreaFracID);	switch (optID) {				case MDinput:				    if ((_MDWetlandAreaFracID       = MFVarGetID (MDVarWetlandAreaFraction,          "[-]",   MFInput, MFState, MFBoundary))  == CMfailed)				    	return (CMfailed);				    break;				case MDcalculate:										if ((MFModelAddFunction (_MDWetlandAreaFraction) == CMfailed)) return (CMfailed);					break;				default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}		MFDefLeaving("WetlandArea");	return (_MDWetlandAreaFracID);	}

int MDRiverbedShapeExponentDef () {	int  optID = MFUnset;	const char *optStr, *optName = MDOptRiverbed;	const char *options [] = { MDInputStr, "slope-independent", "slope-dependent", (char *) NULL };	if (_MDOutRiverbedShapeExponentID != MFUnset) return (_MDOutRiverbedShapeExponentID);	MFDefEntering ("Riverbed Geometry");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options,optStr,true);	switch (optID) {		case MDinput:			if (((_MDOutRiverbedAvgDepthMeanID  = MFVarGetID (MDVarRiverbedAvgDepthMean,  "m",      MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedWidthMeanID     = MFVarGetID (MDVarRiverbedWidthMean,     "m",      MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedVelocityMeanID  = MFVarGetID (MDVarRiverbedVelocityMean,  "m/s",    MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedShapeExponentID = MFVarGetID (MDVarRiverbedShapeExponent, MFNoUnit, MFInput,  MFState, MFBoundary)) == CMfailed))				return (CMfailed);			break;		case MDdependent:			if ((_MDInRiverbedSlopeID           = MFVarGetID (MDVarRiverbedSlope,         "m/km",   MFInput,  MFState, MFBoundary)) == CMfailed)				return (CMfailed);		case MDindependent:			if (((_MDInDischMeanID = MDDischMeanDef ()) == CMfailed) ||			    ((_MDOutRiverbedAvgDepthMeanID  = MFVarGetID (MDVarRiverbedAvgDepthMean,  "m",      MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedWidthMeanID     = MFVarGetID (MDVarRiverbedWidthMean,     "m",      MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedVelocityMeanID  = MFVarGetID (MDVarRiverbedVelocityMean,  "m/s",    MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutRiverbedShapeExponentID = MFVarGetID (MDVarRiverbedShapeExponent, MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDRiverbedShapeExponent) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Riverbed Geometry");	return (_MDOutRiverbedShapeExponentID);}

int MDSmallReservoirCapacityDef () {	int  optID = MFUnset;	const char *optStr, *optName = MDVarSmallResCapacity;	const char *options [] = { MDNoneStr, MDInputStr, MDCalculateStr, (char *) NULL };	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);			if ((optID == MDnone) || (_MDOutSmallResCapacityID != MFUnset)) return (_MDOutSmallResCapacityID);	MFDefEntering("SmallReservoirCapacity");	if ((_MDInIrrGrossDemandID = MDIrrGrossDemandDef  ()) != MFUnset) {		switch (optID) {			case MDinput:			    if ((_MDOutSmallResCapacityID       = MFVarGetID (MDVarSmallResCapacity,          "mm",   MFInput, MFState, MFBoundary))  == CMfailed)			    	return (CMfailed);			    break;			case MDcalculate:				if ((_MDInIrrGrossDemandID == CMfailed) ||				    ((_MDInIrrAreaID                    = MDIrrigatedAreaDef ())==  CMfailed) ||				    ((_MDInRainSurfRunoffID             = MFVarGetID (MDVarRainSurfRunoff,            "mm",  MFInput,  MFFlux,  MFBoundary)) == CMfailed) ||				    ((_MDOutRainSurfRunoffAccumulatedID = MFVarGetID ("__SurfaceROAccumulated",       "mm",  MFOutput, MFState, MFInitial))  == CMfailed) ||				    ((_MDOutIrrGrossDemandAccumulatedID = MFVarGetID ("__GrossDemandAccumulated",     "mm",  MFOutput, MFState, MFInitial))  == CMfailed) ||				    ((_MDInSmallResStorageFractionID    = MFVarGetID (MDVarSmallReservoirStorageFrac, "-",   MFInput,MFState,MFBoundary))  == CMfailed) ||				    ((_MDOutSmallResCapacityID          = MFVarGetID (MDVarSmallResCapacity,          "mm",  MFOutput, MFState, MFInitial))  == CMfailed) ||				    (MFModelAddFunction (_MDSmallReservoirCapacity) == CMfailed)) return (CMfailed);				break;			default: MFOptionMessage (optName, optStr, options); return (CMfailed);		}	}	MFDefLeaving("SmallReservoirCapacity");	return (_MDOutSmallResCapacityID);}

int MDWTempGrdWaterDef () {	int  optID = MDinput;																												// RJS 061312	const char *optStr, *optName = MDOptGrdWaterTemp;																					// RJS 061312	const char *options [] = { MDCalculateStr, MDInputStr,  (char *) NULL };															// RJS 061312	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);									// RJS 061312	if (_MDOutWTempGrdWaterID != MFUnset) return (_MDOutWTempGrdWaterID);	MFDefEntering ("Groundwater temperature");	switch (optID) {																													// RJS 061312	case MDcalculate:																													// RJS 061312	if (((_MDInWTempSurfRunoffID = MDWTempSurfRunoffDef  ()) == CMfailed) ||	    ((_MDInRainRechargeID    = MDRainInfiltrationDef ()) == CMfailed) ||	    ((_MDInIrrReturnFlowID   = MDIrrReturnFlowDef    ()) == CMfailed) ||	    ((_MDOutGrdWaterID       = MDBaseFlowDef         ()) == CMfailed) ||	    ((_MDInAirTempID         = MFVarGetID (MDVarAirTemperature,  "degC", MFInput,  MFState, MFBoundary))  == CMfailed) ||	    ((_MDOutWTempGrdWaterID  = MFVarGetID (MDVarWTempGrdWater, "degC", MFOutput, MFState,  MFInitial)) == CMfailed) ||	    (MFModelAddFunction (_MDWTempGrdWater) == CMfailed)) return (CMfailed);	break;																																// RJS 061312	case MDinput:	_MDOutWTempGrdWaterID = MFVarGetID (MDVarWTempGrdWater, "degC", MFInput,  MFState, MFBoundary); break;				// RJS 061312	MFInitial changed to MFBoundary in order to read in																													// RJS 061312	default: MFOptionMessage (optName, optStr, options); return (CMfailed);																// RJS 061312	}																																	// RJS 061312	MFDefLeaving ("Groundwater temperature");	return (_MDOutWTempGrdWaterID);}

int MDRainInterceptDef () {	int optID = MFUnset;	const char *optStr, *optName = MDVarRainInterception;	const char *options [] = { MDNoneStr, MDInputStr, MDCalculateStr, (char *) NULL };	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options,optStr,true);	if ((optID == MDnone) || (_MDOutInterceptID != MFUnset)) return (_MDOutInterceptID);	MFDefEntering ("Rainfed Intercept");	switch (optID) {		case MDinput: _MDOutInterceptID = MFVarGetID (MDVarRainInterception,  "mm", MFInput, MFFlux, false); break;		case MDcalc:			if (((_MDInPrecipID        = MDPrecipitationDef ()) == CMfailed) ||			    	    	    ((_MDInSPackChgID      = MDSPackChgDef      ()) == CMfailed) ||			    ((_MDInLeafAreaIndexID = MDLeafAreaIndexDef ()) == CMfailed) ||			    ((_MDInStemAreaIndexID = MDStemAreaIndexDef ()) == CMfailed) ||			    ((_MDInPetID           = MDRainPotETDef     ()) == CMfailed) ||			    ((_MDOutInterceptID    = MFVarGetID (MDVarRainInterception,   "mm",     MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDRainIntercept) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Rainfed Intercept");	return (_MDOutInterceptID); }

int MDReservoirDef () {	int optID = MFUnset;	const char *optStr, *optName = MDOptReservoirs;	const char *options [] = { MDNoneStr, MDCalculateStr, (char *) NULL }; 	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true); 	if ((optID == MDnone) || (_MDOutResReleaseID != MFUnset)) return (_MDOutResReleaseID);	MFDefEntering ("Reservoirs");	switch (optID) {		case MDcalculate:			if (((_MDInDischMeanID      = MDDischMeanDef ())   == CMfailed) ||			    ((_MDInDischargeID      = MDDischLevel2Def ()) == CMfailed) ||			    ((_MDInResCapacityID    = MFVarGetID (MDVarReservoirCapacity,      "km3",  MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutResStorageID    = MFVarGetID (MDVarReservoirStorage,       "km3",  MFOutput, MFState, MFInitial))  == CMfailed) ||			    ((_MDOutResStorageChgID = MFVarGetID (MDVarReservoirStorageChange, "km3",  MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutResReleaseID    = MFVarGetID (MDVarReservoirRelease,       "m3/s", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDReservoir) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Reservoirs");	return (_MDOutResReleaseID); }

int MDRainPotETDef () {	int optID = MFUnset;	const char *optStr, *optName = MDVarRainPotEvapotrans;	const char *options [] = { MDInputStr, "Hamon", "Jensen", "PsTaylor", "Pstd", "PMday", "PMdn", "SWGday", "SWGdn", "Turc", (char *) NULL };	if (_MDPotETID != MFUnset) return (_MDPotETID);	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options,optStr,true);	MFDefEntering ("Rainfed Potential Evapotranspiration");	switch (optID) {		case MDinput:    _MDPotETID = MFVarGetID (MDVarRainPotEvapotrans, "mm", MFInput, MFFlux, false); break;		case MDHamon:    _MDPotETID = MDRainPotETHamonDef    (); break;		case MDJensen:   _MDPotETID = MDRainPotETJensenDef   (); break;		case MDPsTaylor: _MDPotETID = MDRainPotETPsTaylorDef (); break;		case MDPstd:     _MDPotETID = MDRainPotETPstdDef     (); break;		case MDPMday:    _MDPotETID = MDRainPotETPMdayDef    (); break;		case MDPMdn:     _MDPotETID = MDRainPotETPMdnDef     (); break;		case MDSWGday:   _MDPotETID = MDRainPotETSWGdayDef   (); break;		case MDSWGdn:    _MDPotETID = MDRainPotETSWGdnDef    (); break;		case MDTurc:     _MDPotETID = MDRainPotETTurcDef     (); break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Rainfed Potential Evapotranspiration");	return (_MDPotETID);}

int MDDischLevel3MuskingumCoeffDef () {	int  optID = MFUnset;	const char *optStr, *optName = MDOptMuskingum;	const char *options [] = { MDInputStr, "static", (char *) NULL };	if (_MDOutMuskingumC0ID != MFUnset) return (_MDOutMuskingumC0ID);	MFDefEntering ("Muskingum Coefficients");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	switch (optID) {		case MDinput:			if (((_MDOutMuskingumC0ID = MFVarGetID (MDVarMuskingumC0, MFNoUnit,   MFInput,  MFState,  MFBoundary)) == CMfailed) ||			    ((_MDOutMuskingumC1ID = MFVarGetID (MDVarMuskingumC1, MFNoUnit,   MFInput,  MFState,  MFBoundary)) == CMfailed) ||			    ((_MDOutMuskingumC2ID = MFVarGetID (MDVarMuskingumC2, MFNoUnit,   MFInput,  MFState,  MFBoundary)) == CMfailed))				return (CMfailed);			break;		case MDstatic:			if (((_MDInRiverbedShapeExponentID  = MDRiverbedShapeExponentDef ()) == CMfailed) ||			    ((_MDInRiverbedWidthMeanID      = MFVarGetID (MDVarRiverbedWidthMean,    "m",      MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDInRiverbedAvgDepthMeanID   = MFVarGetID (MDVarRiverbedAvgDepthMean, "m",      MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDInRiverbedVelocityMeanID   = MFVarGetID (MDVarRiverbedVelocityMean, "m/s",    MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDInRiverbedSlopeID          = MFVarGetID (MDVarRiverbedSlope,        "m/km",   MFInput,  MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutMuskingumC0ID           = MFVarGetID (MDVarMuskingumC0,          MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutMuskingumC1ID           = MFVarGetID (MDVarMuskingumC1,          MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutMuskingumC2ID           = MFVarGetID (MDVarMuskingumC2,          MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    ((_MDOutCourantID               = MFVarGetID ("Courant",                 MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDDischRouteMuskingumCoeff) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Muskingum Coefficients");	return (_MDOutMuskingumC0ID);}

int MDGrossRadDef () {    int optID = MFUnset;    const char *optStr, *optName = MDVarGrossRadiance;    const char *options [] = { MDInputStr, "standard", "Otto", (char *) NULL };    float par;    if (_MDOutGrossRadID != MFUnset) return (_MDOutGrossRadID);    MFDefEntering ("Gross Radiance");    if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);    switch (optID) {    case MDinput:        _MDOutGrossRadID = MFVarGetID (MDVarGrossRadiance, "MJ/m^2", MFInput,  MFFlux,  MFBoundary);        break;    case MDstandard:        if (((optStr = MFOptionGet (MDParGrossRadTAU)) != (char *) NULL) && (sscanf (optStr,"%f",&par) == 1))            _MDGrossRadStdTAU = par;        if (((_MDOutGrossRadID    = MFVarGetID (MDVarGrossRadiance, "MJ/m^2", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                (MFModelAddFunction (_MDGrossRadianceStd)  == CMfailed)) return (CMfailed);        break;    case MDOtto:        if (((_MDOutGrossRadID    = MFVarGetID (MDVarGrossRadiance, "MJ/m^2", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                (MFModelAddFunction (_MDGrossRadianceOtto) == CMfailed)) return (CMfailed);        break;    default:        MFOptionMessage (optName, optStr, options);        return (CMfailed);    }    MFDefLeaving ("Gross Radiance");    return (_MDOutGrossRadID);}

int MDBaseFlowDef () {	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 ((_MDInSoilPercolationID = MFVarGetID (MDVarSoilPercolation,     "mm", MFInput, MFFlux,   MFBoundary)) == CMfailed) return (CMfailed); // SZ 10212014        }                int  optID = MFUnset;                                                                                   // RJS 060214	const char *optName = MDVarRunoff;                                                                      // RJS 060214	const char *options [] = { "spatially", MDCalculateStr, MDPnETStr, (char *) NULL };                                // RJS 060214	if (_MDOutBaseFlowID != MFUnset) return (_MDOutBaseFlowID);	MFDefEntering ("Base flow");	if ((optStr  = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);    // RJS 060214        switch (optID) {            case MDspatial:                if ((_MDInGroundWatBETAID     = MFVarGetID(MDParGroundWatBETA,       "1/d", MFInput, MFState, MFBoundary)) == CMfailed) return (CMfailed);             case MDcalculate: 	                 if (_MDInGroundWatBETAID == MFUnset) {                    if (((optStr = MFOptionGet (MDParGroundWatBETA))  != (char *) NULL) && (sscanf (optStr,"%f",&par) == 1)) _MDGroundWatBETA = par;                 }                if (((_MDInRechargeID       = MDRainInfiltrationDef ()) == CMfailed) ||                    ((_MDInIrrGrossDemandID = MDIrrGrossDemandDef   ()) == CMfailed)) return (CMfailed);                if ( _MDInIrrGrossDemandID != MFUnset) {                        if (((_MDInSmallResReleaseID    = MDSmallReservoirReleaseDef ()) == CMfailed) ||                            ((_MDInIrrAreaFracID        = MDIrrigatedAreaDef         ()) ==  CMfailed) ||                            ((_MDInIrrReturnFlowID      = MFVarGetID (MDVarIrrReturnFlow,     "mm", MFInput,  MFFlux,  MFBoundary)) == CMfailed) ||                            ((_MDOutIrrUptakeExternalID = MFVarGetID (MDVarIrrUptakeExternal, "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                            ((_MDOutIrrUptakeGrdWaterID = MDIrrUptakeGrdWaterDef     ()) == CMfailed))                            return CMfailed;                 }                if (((_MDOutGrdWatID                = MFVarGetID (MDVarGroundWater,         "mm", MFOutput, MFState, MFInitial))  == CMfailed) ||                    ((_MDOutGrdWatChgID             = MFVarGetID (MDVarGroundWaterChange,   "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                    ((_MDOutGrdWatRechargeID        = MFVarGetID (MDVarGroundWaterRecharge, "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                    ((_MDOutGrdWatUptakeID          = MFVarGetID (MDVarGroundWaterUptake,   "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                    ((_MDOutBaseFlowID              = MFVarGetID (MDVarBaseFlow,            "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                     (MFModelAddFunction (_MDBaseFlow) == CMfailed)) return (CMfailed);                break;            case MDPnET:                if (_MDInGroundWatBETAID == MFUnset) {                    if (((optStr = MFOptionGet (MDParGroundWatBETA))  != (char *) NULL) && (sscanf (optStr,"%f",&par) == 1)) _MDGroundWatBETA = par;                 }                 if (((_MDInRechargeID               = MDRainInfiltrationDef ()) == CMfailed) ||                     ((_MDOutGrdWatID                = MFVarGetID (MDVarGroundWater,         "mm", MFOutput, MFState, MFInitial))  == CMfailed) ||                     ((_MDOutGrdWatChgID             = MFVarGetID (MDVarGroundWaterChange,   "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                     ((_MDOutGrdWatRechargeID        = MFVarGetID (MDVarGroundWaterRecharge, "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                     ((_MDOutGrdWatUptakeID          = MFVarGetID (MDVarGroundWaterUptake,   "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                     ((_MDOutBaseFlowID              = MFVarGetID (MDVarBaseFlow,            "mm", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                     (MFModelAddFunction (_MDBaseFlow2) == CMfailed)) return (CMfailed);														// RJS 061312                break;	              }	MFDefLeaving ("Base flow ");	return (_MDOutBaseFlowID);}

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 MDRunoffDef () {	int  optID = MFUnset;	const char *optStr, *optName = MDVarRunoff;	const char *options [] = { MDInputStr, MDCalculateStr, "corrected", MDInput2Str, (char *) NULL };	if (_MDOutRunoffID != MFUnset) return (_MDOutRunoffID);	MFDefEntering ("Runoff");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	switch (optID) {//		case MDinput: _MDOutRunoffID = MFVarGetID (MDVarRunoff,         "mm",     MFInput,  MFFlux, MFBoundary); break;						// RJS commented out 061312		case MDinput: 																										       			// RJS 061312			if (((_MDInBaseFlowID          = MFVarGetID (MDVarBaseFlow,          "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	// RJS 061312//				((_MDInRunoffPoolReleaseID = MFVarGetID (MDVarRunoffPoolRelease, "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	// RJS 061312, commented out 082812				((_MDInTotalSurfRunoffID   = MFVarGetID (MDVarTotalSurfRunoff,   "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	// RJS 082812		                                ((_MDOutPropROStormWaterID   = MFVarGetID (MDVarPropROStormWater,    "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                                ((_MDOutPropROSurfaceWaterID = MFVarGetID (MDVarPropROSurfaceWater,  "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                                 ((_MDOutPropROGroundWaterID  = MFVarGetID (MDVarPropROGroundWater,   "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||  				((_MDOutTotalSurfRunoffID  = MFVarGetID (MDVarTotalSurfRunoff, "mm",     MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||				((_MDOutRunoffID           = MFVarGetID (MDVarRunoff,            "mm",    MFOutput, MFFlux, MFBoundary)) == CMfailed) ||	// RJS 061312				(MFModelAddFunction (_MDRunoffInput) == CMfailed)) return (CMfailed);														// RJS 061312			break;																															// RJS 061312                case MDinput2:                        if (((_MDInBaseFlowID   = MDBaseFlowDef   ()) == CMfailed) ||                            ((_MDInRunoffPoolReleaseID = MDSurfRunoffPoolDef ()) == CMfailed) ||		// RJS 042612                            ((_MDOutPropROStormWaterID   = MFVarGetID (MDVarPropROStormWater,    "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                            ((_MDOutPropROSurfaceWaterID = MFVarGetID (MDVarPropROSurfaceWater,  "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                             ((_MDOutPropROGroundWaterID  = MFVarGetID (MDVarPropROGroundWater,   "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||  			    ((_MDOutTotalSurfRunoffID  = MFVarGetID (MDVarTotalSurfRunoff, "mm",     MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||                            ((_MDOutRunoffID           = MFVarGetID (MDVarRunoff,            "mm",    MFOutput, MFFlux, MFBoundary)) == CMfailed) ||	// RJS 061312                            (MFModelAddFunction (_MDRunoffInput2) == CMfailed)) return (CMfailed);														// RJS 061312			break;	                case MDcorrected:			if ((_MDInRunoffCorrID  = MFVarGetID (MDVarRunoffCorretion, MFNoUnit, MFInput,  MFState, MFBoundary)) == CMfailed)				return (CMfailed);			break;	// RJS 082812		case MDcalculate:					if (((_MDInBaseFlowID   = MDBaseFlowDef   ()) == CMfailed) ||//			    ((_MDInSurfRunoffID = MDSurfRunoffDef ()) == CMfailed) ||				//commented out RJS 042612				((_MDInRunoffPoolReleaseID = MDSurfRunoffPoolDef ()) == CMfailed) ||		// RJS 042612//				((_MDInStormRunoffTotalID  = MDStormRunoffDef ()) == CMfailed) ||		// RJS 082812				((_MDInStormRunoffTotalID  = MFVarGetID (MDVarStormRunoffTotal, "mm",    MFInput,  MFFlux,  MFBoundary)) == CMfailed) ||                                ((_MDOutPropROStormWaterID   = MFVarGetID (MDVarPropROStormWater,    "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                                ((_MDOutPropROSurfaceWaterID = MFVarGetID (MDVarPropROSurfaceWater,  "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||       // RJS 100313                                 ((_MDOutPropROGroundWaterID  = MFVarGetID (MDVarPropROGroundWater,   "-",     MFOutput, MFState, MFBoundary)) == CMfailed) ||  				((_MDOutRunoffID           = MFVarGetID (MDVarRunoff,          "mm",     MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||				((_MDOutTotalSurfRunoffID  = MFVarGetID (MDVarTotalSurfRunoff, "mm",     MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||				(MFModelAddFunction (_MDRunoff) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving  ("Runoff");	return (_MDOutRunoffID);}

NCurlType NCurlGetType (const char *url) {	int ret;	NCurlType  codes [] = { NCurlFile, NCurlHttp, NCurlFtp, NCurlInvalid };	const char *opts [] = { "file",  "http",  "ftp", NULL };	if ((ret = CMoptLookup (opts, url, false)) == CMfailed) {		CMmsgPrint (CMmsgUsrError,"Invalid url: /"%s/"!/n", url);		return (NCurlInvalid);	}	return (codes [ret]);}

int MDWTempRiverDef () {	const char *optStr;        	if (_MDOutWTempRiverID != MFUnset) return (_MDOutWTempRiverID);	MFDefEntering ("River temperature");                int  optID = MFUnset;                                                                                   // RJS 060214	const char *optName = MDVarRunoff;                                                                      // RJS 060214	const char *options [] = { MDCalculateStr, MDInputStr, MDInput2Str, (char *) NULL };                    // RJS 060214            if ((optStr  = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);    // RJS 060214                switch (optID) {        case MDcalculate: 	            	if (            ((_MDInSurfRunoffID      = MDRainSurfRunoffDef ()) == CMfailed) ||//	    ((_MDInBaseFlowID        = MDBaseFlowDef       ()) == CMfailed) ||          // commented out 051614 so that input BaseFlow isn't overwritten            ((_MDInBaseFlowID        = MFVarGetID (MDVarBaseFlow,          "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||     //	    ((_MDInWTempSurfRunoffID = MDWTempSurfRunoffDef ()) == CMfailed) ||		// commented out RJS 060512	    ((_MDInWTempSurfRunoffPoolID = MDWTempSurfRunoffPoolDef ()) == CMfailed) ||		// RJS 060512	    ((_MDInWTempGrdWaterID   = MDWTempGrdWaterDef   ()) == CMfailed) ||	    ((_MDInTotalSurfRunoffID = MFVarGetID (MDVarTotalSurfRunoff, "mm",  MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	//RJS 082812	    ((_MDOutWTempRiverID     = MFVarGetID (MDVarWTempRiver,    "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||		    (MFModelAddFunction (_MDWTempRiver) == CMfailed)) return (CMfailed);                break;       case MDinput:        if (            ((_MDInTotalSurfRunoffID   = MFVarGetID (MDVarTotalSurfRunoff, "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||            ((_MDInBaseFlowID        = MFVarGetID (MDVarBaseFlow,       "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||     	    ((_MDInWTempSurfRunoffPoolID = MDWTempSurfRunoffPoolDef ()) == CMfailed) ||		// RJS 060512	    ((_MDInWTempGrdWaterID   = MDWTempGrdWaterDef   ()) == CMfailed) ||//	    ((_MDInTotalSurfRunoffID = MFVarGetID (MDVarTotalSurfRunoff, "mm",  MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	//RJS 082812	    ((_MDOutWTempRiverID     = MFVarGetID (MDVarWTempRiver,    "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||		    (MFModelAddFunction (_MDWTempRiver) == CMfailed)) return (CMfailed);                break;         case MDinput2:        if (            ((_MDInSurfRunoffID      = MFVarGetID (MDVarRainSurfRunoff, "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||            ((_MDInBaseFlowID        = MFVarGetID (MDVarBaseFlow,       "mm",    MFInput,  MFFlux, MFBoundary)) == CMfailed) ||     	    ((_MDInWTempSurfRunoffPoolID = MDWTempSurfRunoffPoolDef ()) == CMfailed) ||		// RJS 060512	    ((_MDInWTempGrdWaterID   = MDWTempGrdWaterDef   ()) == CMfailed) ||	    ((_MDInTotalSurfRunoffID = MFVarGetID (MDVarTotalSurfRunoff, "mm",  MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	//RJS 082812	    ((_MDOutWTempRiverID     = MFVarGetID (MDVarWTempRiver,    "degC", MFOutput, MFState, MFBoundary)) == CMfailed) ||		    (MFModelAddFunction (_MDWTempRiver) == CMfailed)) return (CMfailed);                break;           }	MFDefLeaving ("River temperature");	return (_MDOutWTempRiverID);}

int MDIrrSoilMoistChgDef() {	int optID = MFUnset, ret;	const char *optStr, *optName = MDOptIrrigation;	const char *options [] = { MDNoneStr, MDInputStr, MDCalculateStr, (char *) NULL };	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);			if ((optID == MDnone) || (optID == MDinput) || (_MDOutIrrSoilMoistChgID != MFUnset)) return (_MDOutIrrSoilMoistChgID);	if ((ret = MDIrrGrossDemandDef ()) == CMfailed) return (CMfailed);	if (ret == MFUnset) return (MFUnset);	_MDOutIrrSoilMoistChgID = MFVarGetID (MDVarIrrSoilMoistChange,     "mm",   MFInput, MFFlux,  MFBoundary);    return (_MDOutIrrSoilMoistChgID);}

int MDBgcDINRoutingDef () {	int optID = MFUnset;	const char *optStr;	const char *options [] = { "none", "calculate", (char *) NULL };     //const char *options [] = { MDNoneStr, (char *) NULL };	if (_MDFlux_DINID != MFUnset)	return (_MDFlux_DINID);	MFDefEntering ("DIN Routing Calculation");	if (((optStr = MFOptionGet (MDOptReservoirs))  == (char *) NULL) || ((optID = CMoptLookup ( options, optStr, true)) == CMfailed)) {						CMmsgPrint(CMmsgUsrError,"Reservoir Option not specified! Option none or calculate /n");						return CMfailed;					}		    if (optID==1){		  //  	printf ("Resoption=%i/n",optID);		    if (((_MDInResStorageID       = MFVarGetID (MDVarReservoirStorage,       "km3",   MFInput, MFState, MFInitial))  == CMfailed) ||		        ((_MDInResStorageChangeID = MFVarGetID (MDVarReservoirStorageChange, "km3/s", MFInput, MFState, MFBoundary)) == CMfailed))		    	return CMfailed;		    }			       // Input	if (((_MDInDischargeID            = MDDischargeDef    ()) == CMfailed) ||        ((_MDInRiverWidthID           = MDRiverWidthDef   ()) == CMfailed) ||	    ((_MDInRunoffVolumeID         = MDRunoffVolumeDef ()) == CMfailed) ||	    ((_MDInRunoffID               = MDRunoffDef ()) == CMfailed) ||	    ((_MDInWTempRiverRouteID      = MDWTempRiverRouteDef ()) == CMfailed) ||	    ((_MDInRiverStorageChgID      = MFVarGetID (MDVarRiverStorageChg,         "m3/s",   MFInput,  MFState, MFBoundary)) == CMfailed) ||	    ((_MDInRiverStorageID         = MFVarGetID (MDVarRiverStorage,            "m3",     MFInput,  MFState, MFInitial))  == CMfailed) ||	    ((_MDInPointScenarioID        = MFVarGetID (MDVarPointScenario,           MFNoUnit, MFInput,  MFState, MFBoundary)) == CMfailed) ||        ((_MDNonPoint_DINID           = MFVarGetID (MDVarBGCNonPoint_DIN,         "kg/m3",  MFInput,  MFState, MFBoundary)) == CMfailed) ||	    ((_MDPointSources_DINID       = MFVarGetID (MDVarBGCPointSources_DIN,     "kg/day", MFInput,  MFState, MFBoundary)) == CMfailed) ||	    	    	    // Output	    ((_MDLocalIn_DINID            = MFVarGetID (MDVarDINLocalIn,              "kg/d",   MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||	    ((_MDRemoval_DINID            = MFVarGetID (MDVarDINRemoval,              "kg/d",   MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||        ((_MDConc_DINID               = MFVarGetID (MDVarDINConcentration,        "kg/m3",  MFOutput, MFState, MFBoundary)) == CMfailed) ||	    ((_MDStorage_DINID            = MFVarGetID (MDVarDINStorage,              "kg",     MFOutput, MFState, MFInitial))  == CMfailed) ||	    ((_MDDeltaStorage_DINID       = MFVarGetID (MDVarDINDeltaStorage,         "kg/day", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||        ((_MDFlux_DINID               = MFVarGetID (MDVarDINFlux  ,               "kg/day", MFRoute,  MFFlux,  MFBoundary)) == CMfailed) ||        ((_MDConcMixing_DINID         = MFVarGetID (MDVarDINConcentration_Mixing, "kg/m3",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDStorageMixing_DINID      = MFVarGetID (MDVarDINStorage_Mixing,       "kg",     MFOutput, MFState, MFInitial))  == CMfailed) ||        ((_MDDeltaStorageMixing_DINID = MFVarGetID (MDVarDINDeltaStorage_Mixing,  "kg/day", MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||        ((_MDFluxMixing_DINID         = MFVarGetID (MDVarDINFlux_Mixing ,         "kg/day", MFRoute,  MFFlux,  MFBoundary)) == CMfailed) ||        (MFModelAddFunction (_MDDINRouting) == CMfailed)) return (CMfailed); 	MFDefLeaving ("DIN Routing Calculation");	return (_MDFlux_DINID); }

int MDStormRunoffDef () {	if (_MDOutStormRunoffTotalID != MFUnset) return (_MDOutStormRunoffTotalID);	MFDefEntering ("Storm Runoff");	if (((_MDInPrecipID             = MDPrecipitationDef ()) == CMfailed) ||            ((_MDInH2OFractionID        = MFVarGetID (MDVarH2OFracSpatial,  "-",   MFInput,  MFState, MFBoundary)) == CMfailed) ||		//commented out 082812	    ((_MDInImpFractionID        = MFVarGetID (MDVarImpFracSpatial,  "-",    MFInput,  MFState, MFBoundary)) == CMfailed) ||	    ((_MDInSPackChgID           = MFVarGetID (MDVarSnowPackChange,  "mm",   MFInput,  MFFlux, MFBoundary)) == CMfailed) ||	    ((_MDInHCIAID               = MFVarGetID (MDVarHCIA,            "-",    MFInput,  MFState, MFBoundary)) == CMfailed) ||	    ((_MDOutStormRunoffImpID    = MFVarGetID (MDVarStormRunoffImp,  "mm",   MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||	    ((_MDOutStormRunoffH2OID    = MFVarGetID (MDVarStormRunoffH2O,  "mm",   MFOutput, MFFlux,  MFBoundary)) == CMfailed) ||		//commented out 082812	    ((_MDOutPrecipPervID        = MFVarGetID (MDVarPrecipPerv,      "mm",   MFOutput, MFFlux, MFBoundary)) == CMfailed) ||	    ((_MDOutStormRunoffTotalID  = MFVarGetID (MDVarStormRunoffTotal,"mm",   MFOutput, MFFlux, MFBoundary))  == CMfailed) ||            ((_MDOutH2OFractionID       = MFVarGetID (MDVarH2OFracSpatial,  "-",    MFOutput,  MFState, MFBoundary)) == CMfailed) ||		// RJS 100113	    ((_MDOutImpFractionID       = MFVarGetID (MDVarImpFracSpatial,  "-",    MFOutput,  MFState, MFBoundary)) == CMfailed) ||           // RJS 100113            ((_MDOutRunofftoPervID      = MFVarGetID (MDVarRunofftoPerv,    "mm",   MFOutput, MFFlux, MFBoundary))  == CMfailed)            ) return (CMfailed);        // Add in check for development ... then read-in Impervious and HCIA ... then add model function        enum { MFnone, MFcalculate };        const char *optStr;        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 (StormRunoff)./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);        }	// Include StormRunoff        if ((MFModelAddFunction (_MDStormRunoff) == CMfailed)) return (CMfailed);	MFDefLeaving ("Storm Runoff");	return (_MDOutRunofftoPervID);}

int MDBgcRiverGPPDef () {	int optID = MFUnset;    const char *optStr;		const char *options [] = { "none", "calculate", (char *) NULL };	if (_MDRiverGPPID != MFUnset) return (_MDRiverGPPID);	MFDefEntering ("Calculate river GPP");		if (((optStr = MFOptionGet (MDOptReservoirs))  == (char *) NULL) || ((optID = CMoptLookup ( options, optStr, true)) == CMfailed)) {				CMmsgPrint(CMmsgUsrError,"Reservoir Option not specified! Option none or calculate /n");				return CMfailed;			}    if (optID==1){//  	printf ("Resoption=%i/n",optID);    if (((_MDInResStorageID           = MFVarGetID (MDVarReservoirStorage,       "km3",     MFInput, MFState, MFInitial))   == CMfailed) ||        ((_MDInResStorageChangeID     = MFVarGetID (MDVarReservoirStorageChange, "km3/s",   MFInput, MFState, MFBoundary))  == CMfailed))    	return CMfailed;    }	//input	if (//((_MDInDischargeID            = MDDischargeDef    ()) == CMfailed) ||		 ((_MDInRiverLightID          = MDBgcRiverLightDef ()) == CMfailed) ||         ((_MDInDischargeID           = MFVarGetID (MDVarDischarge,              "m3/s",    MFInput, MFState,  MFBoundary)) == CMfailed) ||         ((_MDInRiverWidthID          = MFVarGetID (MDVarRiverWidth,             "m",       MFInput, MFState,  MFBoundary)) == CMfailed) ||         ((_MDInPAR2BottomID          = MFVarGetID (MDVarPAR2Bottom,             "MJ/m2/d", MFInput,  MFState, MFBoundary)) == CMfailed) ||// output     	((_MDBenthicAlgaeCHLID        = MFVarGetID (MDVarBenthicAlgaeCHL,        "g/m2",    MFOutput, MFState, MFInitial))  == CMfailed) ||        ((_MDBenthicAlgaeCID          = MFVarGetID (MDVarBenthicAlgaeC,          "g/m2",    MFOutput, MFState, MFInitial))  == CMfailed) ||        ((_MDBenthicGPPID             = MFVarGetID (MDVarBenthicGPP,             "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicRaID              = MFVarGetID (MDVarBenthicRa,              "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicNPPID             = MFVarGetID (MDVarBenthicNPP,             "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicMortalityID       = MFVarGetID (MDVarBenthicMortality,       "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicAlgaeC_REACHID    = MFVarGetID (MDVarBenthicAlgaeC_REACH,    "g/m2",    MFOutput, MFState, MFInitial))  == CMfailed) ||        ((_MDBenthicGPP_REACHID       = MFVarGetID (MDVarBenthicGPP_REACH,       "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicRa_REACHID        = MFVarGetID (MDVarBenthicRa_REACH,        "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicNPP_REACHID       = MFVarGetID (MDVarBenthicNPP_REACH,       "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||        ((_MDBenthicMortality_REACHID = MFVarGetID (MDVarBenthicMortality_REACH, "g/m2/d",  MFOutput, MFState, MFBoundary)) == CMfailed) ||       (MFModelAddFunction (_MDRiverGPP) == CMfailed)) return (CMfailed);      	   MFDefLeaving ("Calculate river GPP");	   return (_MDRiverGPPID);}

int MDDischLevel1Def() {   char *optStr;   const char *options [] = { MDNoneStr, (char *) NULL };   	if (_MDOutDischLevel1ID != MFUnset) return (_MDOutDischLevel1ID);	MFDefEntering ("Discharge Level 1");	if ((_MDInDischLevel2ID = MDDischLevel2Def ()) == CMfailed) return (CMfailed);	if (((optStr = MFOptionGet (MDOptReservoirs)) != (char *) NULL) && (CMoptLookup (options,optStr,true) == CMfailed)) {		if ((_MDInDischReleasedID = MDReservoirDef ()) == CMfailed) return (CMfailed);	}	if (((_MDOutDischLevel1ID = MFVarGetID ("__DischLevel1", "m3/s",  MFOutput,  MFState, MFBoundary)) == CMfailed) ||	    (MFModelAddFunction (_MDDischLevel1) == CMfailed)) return (CMfailed);	MFDefLeaving ("Discharge Level 1");	return (_MDOutDischLevel1ID);}

int MDSpecCondDef () {    int  optID = MFUnset;											//SZ 08212014    const char *optStr, *optName = MDOptSpecConductance;								//SZ 08212014    const char *options [] = { MDCalculateStr, MDInputStr, MDInput2Str, MDNoneStr, (char *) NULL };		//SZ 08212014	MFDefEntering ("Specific Conductance Routing");	    if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);  //SZ 08212014    switch (optID) {	//SZ 08212014        case MDcalculate:        if (//            ((_MDInDischargeID                  = MDDischargeDef()) == CMfailed ) ||            ((_MDInDINFluxID                    = MDDINDef ()) == CMfailed) ||     // Needed for merging with upstream	    ((_MDInWTempRiverID                 = MFVarGetID (MDVarWTemp_QxT,              "degC",      MFInput, MFState, MFBoundary)) == CMfailed)   ||            ((_MDInRiverWidthID                 = MDRiverWidthDef ())     == CMfailed) ||//          ((_MDInLitterFall_POCID             = MDLitterFallDef ()) == CMfailed) ||//          ((_MDInLocalLoad_DOCID              = MFVarGetID (MDVarLocalLoadDOC,                             "kg/d",    MFInput,  MFFlux,  MFBoundary))   == CMfailed) ||            ((_MDInDischarge0ID                 = MFVarGetID (MDVarDischarge0,                               "m3/s",    MFInput,  MFState, MFBoundary))   == CMfailed) ||		    ((_MDInDischargeID                  = MFVarGetID (MDVarDischarge,                                "m3/s",    MFInput,  MFState, MFBoundary))   == CMfailed) ||            ((_MDInRiverStorageID               = MFVarGetID (MDVarRiverStorage,                           "m3/day",    MFInput,  MFState, MFInitial))    == CMfailed) ||	            ((_MDInBaseFlowID                   = MFVarGetID (MDVarBaseFlow,                                  "mm",    MFInput,  MFFlux,MFBoundary))    == CMfailed) ||            ((_MDInStormRunoffTotalID           = MFVarGetID (MDVarStormRunoffTotal,                          "mm",   MFInput, MFFlux, MFBoundary))  == CMfailed) ||            ((_MDInRunoffPoolReleaseID          = MFVarGetID (MDVarRunoffPoolRelease,                        "mm",  MFInput, MFFlux, MFBoundary))  == CMfailed) ||            ((_MDInAirTemperatureID             = MFVarGetID (MDVarAirTemperature,         "degC",       MFInput,  MFState, MFBoundary)) == CMfailed) ||            ((_MDInSubFractionID                = MFVarGetID (MDVarLandUseSpatialSub,                         "-",    MFInput,  MFState, MFBoundary))   == CMfailed) ||            ((_MDOutLocalLoadSCID               = MFVarGetID (MDVarLocalLoadSC,               "ic/d",   MFOutput,  MFFlux,  MFBoundary))   == CMfailed) ||             ((_MDOutPostConc_SCID               = MFVarGetID (MDVarPostSpecCond,    	      "uS/cm",   MFOutput,  MFFlux, MFBoundary))   == CMfailed) ||	                           ((_MDOutStoreWater_SCID             = MFVarGetID (MDVarStoreWaterSC,              "ic/day",   MFOutput,  MFState, MFInitial))    == CMfailed) ||	            ((_MDOutFlux_SCID                   = MFVarGetID (MDVarFluxSC,                    "ic/day",   MFRoute,   MFFlux,  MFBoundary))   == CMfailed) ||	              (MFModelAddFunction (_MDSpecCond) == CMfailed)) return (CMfailed);            break;        default:             MFOptionMessage (optName, optStr, options); return (CMfailed);    }    	MFDefLeaving ("Specific Conductance Routing");    return (_MDOutFlux_SCID);}

int MDBgcRiverLightDef () {	int optID = MFUnset;    const char *optStr;		const char *options [] = { "none", "calculate", (char *) NULL };	if (_MDRiverLightID != MFUnset) return (_MDRiverLightID);	MFDefEntering ("Calculate river light");		if (((optStr = MFOptionGet (MDOptReservoirs))  == (char *) NULL) || ((optID = CMoptLookup ( options, optStr, true)) == CMfailed)) {				CMmsgPrint(CMmsgUsrError,"Reservoir Option not specified! Option none or calculate /n");				return CMfailed;			}    if (optID==1){  //  	printf ("Resoption=%i/n",optID);    if (((_MDInResStorageID        = MFVarGetID (MDVarReservoirStorage,       "km3",   MFInput, MFState, MFInitial))  == CMfailed) ||        ((_MDInResStorageChangeID  = MFVarGetID (MDVarReservoirStorageChange, "km3/s", MFInput, MFState, MFBoundary)) == CMfailed))    	return CMfailed;    }			//input	if (//((_MDInDischargeID          = MDDischargeDef    ()) == CMfailed) ||         // ((_MDInSolarRadID         = MDSolarRadDef     ()) == CMfailed) ||         ((_MDInFluxDOCID           = MDBgcDOCRoutingDef   ()) == CMfailed) ||         ((_MDInSolarRadID          = MFVarGetID (MDVarSolarRadiation,     "MJ/m2/d", MFInput,  MFFlux,  MFBoundary)) == CMfailed) ||         ((_MDInDischargeID         = MFVarGetID (MDVarDischarge,          "m3/s",    MFInput,  MFState, MFBoundary)) == CMfailed) ||         ((_MDInConcDOCID           = MFVarGetID (MDVarDOCConcentration,   "kg/m3",   MFInput,  MFState, MFBoundary)) == CMfailed) ||         ((_MDInRiverWidthID        = MFVarGetID (MDVarRiverWidth,         "m",       MFInput,  MFState, MFBoundary)) == CMfailed) ||         ((_MDInRiverbedWidthMeanID = MFVarGetID (MDVarRiverbedWidthMean,  "m",       MFInput,  MFState, MFBoundary)) == CMfailed) || 	     ((_MDInRiverDepthID        = MFVarGetID (MDVarRiverDepth,         "m",       MFInput,  MFState, MFBoundary)) == CMfailed) || 	     ((_MDInKoppenID            = MFVarGetID (MDVarKoppen,             MFNoUnit,  MFInput,  MFState, MFBoundary)) == CMfailed) ||        // output        ((_MDPAR2BottomID           = MFVarGetID (MDVarPAR2Bottom,         "MJ/m2/d", MFOutput, MFState, MFBoundary))  == CMfailed) ||        ((_MDPAR2ReachID            = MFVarGetID (MDVarPAR2Reach,          "MJ/d",    MFOutput, MFState, MFBoundary))  == CMfailed) ||       (MFModelAddFunction (_MDRiverLight) == CMfailed)) return (CMfailed);      	   MFDefLeaving ("Calculate river light");	   return (_MDRiverLightID);}

int MDCParamZ0gDef () {	int optID = MFUnset;	const char *optStr, *optName = MDVarCParamZ0g;	const char *options [] = { MDInputStr, MDLookupStr, (char *) NULL };	if (_MDOutCParamZ0gID != MFUnset) return (_MDOutCParamZ0gID);	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options,optStr,true);	MFDefEntering ("Z0g");	switch (optID) {		case MDinput: _MDOutCParamZ0gID = MFVarGetID (MDVarCParamZ0g,  "m", MFInput, MFState, false); break;		case MDlookup:			if (((_MDInCoverID = MDLandCoverDef ()) == CMfailed) ||	    		 ((_MDOutCParamZ0gID = MFVarGetID (MDVarCParamZ0g, "m", MFOutput, MFState, false)) == CMfailed) ||	    		 (MFModelAddFunction (_MDCParamZ0g) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Z0g");	return (_MDOutCParamZ0gID); }

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 MDDischMeanDef () {	int  optID = MFUnset;	const char *optStr, *optName = MDVarDischMean;	const char *options [] = { MDInputStr, MDCalculateStr, (char *) NULL };	if (_MDOutDischMeanID != MFUnset) return (_MDOutDischMeanID);	MFDefEntering ("Discharge Mean");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	switch (optID) {		case MDinput: _MDOutDischMeanID  = MFVarGetID (MDVarDischMean, "m3/s", MFInput,  MFState, MFBoundary); break;		case MDcalculate:			if (((_MDInAvgNStepsID       = MDAvgNStepsDef ())   == CMfailed) ||			    ((_MDInAccumDischargeID  = MDAccumRunoffDef ()) == CMfailed) ||			    ((_MDOutDischMeanID      = MFVarGetID (MDVarDischMean, "m3/s", MFOutput, MFState, MFInitial))  == CMfailed) ||			    (MFModelAddFunction (_MDDischMean) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Discharge Mean");	return (_MDOutDischMeanID);}

int MDRunoffVolumeDef () {	int optID = MFUnset;	const char *optStr, *optName = MDVarRunoffVolume;	const char *options [] = { MDInputStr, MDCalculateStr, (char *) NULL };	if (_MDOutRunoffVolumeID != MFUnset) return (_MDOutRunoffVolumeID);	MFDefEntering ("Runoff Volume");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	switch (optID) {		case MDinput:			_MDOutRunoffVolumeID = MFVarGetID (MDVarRunoffVolume, "m3/s", MFInput, MFState, MFBoundary);			break;		case MDcalculate:			if (((_MDInRunoffID        = MDRunoffDef ()) == CMfailed) ||			    ((_MDOutRunoffVolumeID = MFVarGetID (MDVarRunoffVolume, "m3/s", MFOutput, MFState, MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDRunoffVolume) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving  ("Runoff Volume");	return (_MDOutRunoffVolumeID);}

int MDCParamAlbedoDef () {	int optID = MFUnset;	const char *optStr, *optName = MDVarCParamAlbedo;	const char *options [] = { MDInputStr, MDLookupStr, (char *) NULL };	if (_MDOutCParamAlbedoID != MFUnset) return (_MDOutCParamAlbedoID);	MFDefEntering ("Albedo");	if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options,optStr,true);	switch (optID) {		case MDinput:  _MDOutCParamAlbedoID = MFVarGetID (MDVarCParamAlbedo,  MFNoUnit, MFInput, MFState, MFBoundary); break;		case MDlookup:			if (((_MDInCoverID    = MDLandCoverDef ()) == CMfailed) ||			    ((_MDInSnowPackID = MDSPackChgDef  ()) == CMfailed) ||			    ((_MDOutCParamAlbedoID = MFVarGetID (MDVarCParamAlbedo, MFNoUnit, MFOutput, MFState, MFBoundary)) == CMfailed) ||			    (MFModelAddFunction (_MDCParamAlbedo) == CMfailed)) return (CMfailed);			break;		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	MFDefLeaving ("Albedo");	return (_MDOutCParamAlbedoID); }

int main (int argc,char *argv []) {	//printf("at start of main/n");	int argNum;	//int DELAY=10;	//int i;	int  optID = MDbankfullQcalc; //CHANGED may2011	const char *optStr, *optName = MDOptModel;	const char *options [] = { "pet", "surplus", "infiltration", "runoff", "discharge",  "watertemp", "balance", "geometry", "bgc", "bgc_DIN","bgc_DINPLUSBIOMASS", "bgc_DOC", "fecal", "sedimentflux","bedloadflux", "BQARTinputs" , "BQARTpreprocess","bankfullQcalc","riverbedshape",(char *) NULL };	argNum = MFOptionParse (argc,argv);			if ((optStr = MFOptionGet (optName)) != (char *) NULL) optID = CMoptLookup (options, optStr, true);	//printf("before switch/n");	switch (optID) {		case MDpet:          return (MFModelRun (argc,argv,argNum,MDRainPotETDef));		case MDsurplus:      return (MFModelRun (argc,argv,argNum,MDRainWaterSurplusDef));		case MDinfiltration: return (MFModelRun (argc,argv,argNum,MDRainInfiltrationDef));		case MDrunoff:       return (MFModelRun (argc,argv,argNum,MDRunoffDef));		case MDdischarge:    return (MFModelRun (argc,argv,argNum,MDDischargeDef));		case MDwatertemp:    return (MFModelRun (argc,argv,argNum,MDWTempRiverRouteDef));		case MDbalance:      return (MFModelRun (argc,argv,argNum,MDWaterBalanceDef));		case MDgeometry:     return (MFModelRun (argc,argv,argNum,MDRiverWidthDef));		case MDbgc:          return (MFModelRun (argc,argv,argNum,MDBgcRoutingDef));		case MDbgc_DOC:      return (MFModelRun (argc,argv,argNum,MDBgcDOCRoutingDef));		case MDbgc_DIN:      return (MFModelRun (argc,argv,argNum,MDBgcDINRoutingDef));		case MDbgc_DINPLUSBIOMASS: return (MFModelRun (argc,argv,argNum,MDBgcDINPlusBiomassRoutingDef));		case MDsedimentflux:    return (MFModelRun (argc,argv,argNum,MDSedimentFluxDef));		case MDbedloadflux:    return (MFModelRun (argc,argv,argNum,MDBedloadFluxDef));			case MDBQARTpreprocess:    return (MFModelRun (argc,argv,argNum,MDBQARTpreprocessDef));			case MDbankfullQcalc:    return (MFModelRun (argc,argv,argNum,MDBankfullQcalcDef));			case MDRiverbedShapeExponent:    return (MFModelRun (argc,argv,argNum,MDRiverbedShapeExponentDef));		default: MFOptionMessage (optName, optStr, options); return (CMfailed);	}	return (CMfailed);}