Example #1
0
 public FireEffects(int numberOfSeverities)
 {
     ReductionsTable = new FireReductions[numberOfSeverities+1];  //will ignore zero
     
     for(int i=0; i <= numberOfSeverities; i++)
     {
         ReductionsTable[i] = new FireReductions();
     }
 }
Example #2
0
        //---------------------------------------------------------------------

        protected override IInputParameters Parse()
        {
            ReadLandisDataVar();

            int numLitterTypes     = 4;
            int numFunctionalTypes = 25;

            Parameters parameters = new Parameters(ecoregionDataset, speciesDataset, numLitterTypes, numFunctionalTypes);

            InputVar <int> timestep = new InputVar <int>(Names.Timestep);

            ReadVar(timestep);
            parameters.Timestep = timestep.Value;

            InputVar <SeedingAlgorithms> seedAlg = new InputVar <SeedingAlgorithms>(Names.SeedingAlgorithm);

            ReadVar(seedAlg);
            parameters.SeedAlgorithm = seedAlg.Value;

            //---------------------------------------------------------------------------------

            InputVar <string> initCommunities = new InputVar <string>("InitialCommunities");

            ReadVar(initCommunities);
            parameters.InitialCommunities = initCommunities.Value;

            InputVar <string> communitiesMap = new InputVar <string>("InitialCommunitiesMap");

            ReadVar(communitiesMap);
            parameters.InitialCommunitiesMap = communitiesMap.Value;


            InputVar <string> climateConfigFile = new InputVar <string>(Names.ClimateConfigFile);

            ReadVar(climateConfigFile);
            parameters.ClimateConfigFile = climateConfigFile.Value;

            InputVar <bool> calimode = new InputVar <bool>(Names.CalibrateMode);

            if (ReadOptionalVar(calimode))
            {
                parameters.CalibrateMode = calimode.Value;
            }
            else
            {
                parameters.CalibrateMode = false;
            }

            InputVar <double> spinMort = new InputVar <double>("SpinupMortalityFraction");

            if (ReadOptionalVar(spinMort))
            {
                parameters.SpinupMortalityFraction = spinMort.Value;
            }
            else
            {
                parameters.SpinupMortalityFraction = 0.0;
            }

            InputVar <string> wt = new InputVar <string>("WaterDecayFunction");

            ReadVar(wt);
            parameters.WType = WParse(wt.Value);

            InputVar <double> pea = new InputVar <double>("ProbEstablishAdjust");

            if (ReadOptionalVar(pea))
            {
                parameters.ProbEstablishAdjustment = pea.Value;
            }
            else
            {
                parameters.ProbEstablishAdjustment = 1.0;
            }

            //InputVar<string> soilCarbonMaps = new InputVar<string>("SoilCarbonMapNames");
            //if (ReadOptionalVar(soilCarbonMaps))
            //{
            //    PlugIn.SoilCarbonMapNames = soilCarbonMaps.Value;

            //    InputVar<int> soilCarbonMapFreq = new InputVar<int>("SoilCarbonMapFrequency");
            //    ReadVar(soilCarbonMapFreq);
            //    PlugIn.SoilCarbonMapFrequency = soilCarbonMapFreq.Value;

            //}

            //InputVar<string> soilNitrogenMaps = new InputVar<string>("SoilNitrogenMapNames");
            //if (ReadOptionalVar(soilNitrogenMaps))
            //{
            //    PlugIn.SoilNitrogenMapNames = soilNitrogenMaps.Value;

            //    InputVar<int> soilNitrogenMapFreq = new InputVar<int>("SoilNitrogenMapFrequency");
            //    ReadVar(soilNitrogenMapFreq);
            //    PlugIn.SoilNitrogenMapFrequency = soilNitrogenMapFreq.Value;

            //}

            InputVar <string> anppMaps = new InputVar <string>("ANPPMapNames");

            if (ReadOptionalVar(anppMaps))
            {
                PlugIn.ANPPMapNames = anppMaps.Value;

                InputVar <int> anppMapFreq = new InputVar <int>("ANPPMapFrequency");
                ReadVar(anppMapFreq);
                PlugIn.ANPPMapFrequency = anppMapFreq.Value;
            }

            InputVar <string> aneeMaps = new InputVar <string>("ANEEMapNames");

            if (ReadOptionalVar(aneeMaps))
            {
                PlugIn.ANEEMapNames = aneeMaps.Value;

                InputVar <int> aneeMapFreq = new InputVar <int>("ANEEMapFrequency");
                ReadVar(aneeMapFreq);
                PlugIn.ANEEMapFrequency = aneeMapFreq.Value;
            }

            InputVar <string> soilCarbonMaps = new InputVar <string>("SoilCarbonMapNames");

            if (ReadOptionalVar(soilCarbonMaps))
            {
                PlugIn.SoilCarbonMapNames = soilCarbonMaps.Value;

                InputVar <int> soilCarbonMapFreq = new InputVar <int>("SoilCarbonMapFrequency");
                ReadVar(soilCarbonMapFreq);
                PlugIn.SoilCarbonMapFrequency = soilCarbonMapFreq.Value;
            }

            InputVar <string> soilNitrogenMaps = new InputVar <string>("SoilNitrogenMapNames");

            if (ReadOptionalVar(soilNitrogenMaps))
            {
                PlugIn.SoilNitrogenMapNames = soilNitrogenMaps.Value;

                InputVar <int> soilNitrogenMapFreq = new InputVar <int>("SoilNitrogenMapFrequency");
                ReadVar(soilNitrogenMapFreq);
                PlugIn.SoilNitrogenMapFrequency = soilNitrogenMapFreq.Value;
            }

            InputVar <string> totalCMaps = new InputVar <string>("TotalCMapNames");

            if (ReadOptionalVar(totalCMaps))
            {
                PlugIn.TotalCMapNames = totalCMaps.Value;

                InputVar <int> totalCMapFreq = new InputVar <int>("TotalCMapFrequency");
                ReadVar(totalCMapFreq);
                PlugIn.TotalCMapFrequency = totalCMapFreq.Value;
            }
            //--------------------------
            //  MinRelativeBiomass table

            ReadName("AvailableLightBiomass");

            List <IEcoregion> ecoregions    = ReadEcoregions();
            string            lastEcoregion = ecoregions[ecoregions.Count - 1].Name;

            InputVar <byte> shadeClassVar = new InputVar <byte>("Shade Class");

            for (byte shadeClass = 1; shadeClass <= 5; shadeClass++)
            {
                if (AtEndOfInput)
                {
                    throw NewParseException("Expected a line with available light class {0}", shadeClass);
                }

                StringReader currentLine = new StringReader(CurrentLine);
                ReadValue(shadeClassVar, currentLine);
                if (shadeClassVar.Value.Actual != shadeClass)
                {
                    throw new InputValueException(shadeClassVar.Value.String,
                                                  "Expected the available light class {0}", shadeClass);
                }

                foreach (IEcoregion ecoregion in ecoregions)
                {
                    InputVar <Percentage> MinRelativeBiomass = new InputVar <Percentage>("Ecoregion " + ecoregion.Name);
                    ReadValue(MinRelativeBiomass, currentLine);
                    parameters.SetMinRelativeBiomass(shadeClass, ecoregion, MinRelativeBiomass.Value);
                }

                CheckNoDataAfter("the Ecoregion " + lastEcoregion + " column",
                                 currentLine);
                GetNextLine();
            }

            //----------------------------------------------------------
            //  Read table of sufficient light probabilities.
            //  Available light classes are in increasing order.
            ReadName("LightEstablishmentTable");

            InputVar <byte>   sc  = new InputVar <byte>("Available Light Class");
            InputVar <double> pl0 = new InputVar <double>("Probability of Germination - Light Level 0");
            InputVar <double> pl1 = new InputVar <double>("Probability of Germination - Light Level 1");
            InputVar <double> pl2 = new InputVar <double>("Probability of Germination - Light Level 2");
            InputVar <double> pl3 = new InputVar <double>("Probability of Germination - Light Level 3");
            InputVar <double> pl4 = new InputVar <double>("Probability of Germination - Light Level 4");
            InputVar <double> pl5 = new InputVar <double>("Probability of Germination - Light Level 5");

            int previousNumber = 0;

            while (!AtEndOfInput && CurrentName != Names.SpeciesParameters &&
                   previousNumber != 6)
            {
                StringReader currentLine = new StringReader(CurrentLine);

                ISufficientLight suffLight = new SufficientLight();

                ReadValue(sc, currentLine);
                suffLight.ShadeClass = sc.Value;

                //  Check that the current shade class is 1 more than
                //  the previous number (numbers are must be in increasing order).
                if (sc.Value.Actual != (byte)previousNumber + 1)
                {
                    throw new InputValueException(sc.Value.String,
                                                  "Expected the severity number {0}",
                                                  previousNumber + 1);
                }
                previousNumber = (int)sc.Value.Actual;

                ReadValue(pl0, currentLine);
                suffLight.ProbabilityLight0 = pl0.Value;

                ReadValue(pl1, currentLine);
                suffLight.ProbabilityLight1 = pl1.Value;

                ReadValue(pl2, currentLine);
                suffLight.ProbabilityLight2 = pl2.Value;

                ReadValue(pl3, currentLine);
                suffLight.ProbabilityLight3 = pl3.Value;

                ReadValue(pl4, currentLine);
                suffLight.ProbabilityLight4 = pl4.Value;

                ReadValue(pl5, currentLine);
                suffLight.ProbabilityLight5 = pl5.Value;

                parameters.LightClassProbabilities.Add(suffLight);

                CheckNoDataAfter("the " + pl5.Name + " column",
                                 currentLine);
                GetNextLine();
            }
            if (parameters.LightClassProbabilities.Count == 0)
            {
                throw NewParseException("No sufficient light probabilities defined.");
            }
            if (previousNumber != 5)
            {
                throw NewParseException("Expected shade class {0}", previousNumber + 1);
            }

            //-------------------------
            //  Species Parameters table

            ReadName("SpeciesParameters");

            speciesLineNums.Clear();  //  If parser re-used (i.e., for testing purposes)

            InputVar <int>    ft            = new InputVar <int>("Functional Type");
            InputVar <bool>   nt            = new InputVar <bool>("Nitrogen Fixer");
            InputVar <int>    gddmn         = new InputVar <int>("Growing Degree Day Minimum");
            InputVar <int>    gddmx         = new InputVar <int>("Growing Degree Day Maximum");
            InputVar <int>    mjt           = new InputVar <int>("Minimum January Temperature");
            InputVar <double> maxd          = new InputVar <double>("Maximum Allowable Drought");
            InputVar <double> leafLongevity = new InputVar <double>("Leaf Longevity");
            InputVar <bool>   epicorm       = new InputVar <bool>("Epicormic:  Y/N");
            InputVar <double> leafLignin    = new InputVar <double>("Leaf Percent Lignin");
            InputVar <double> wLignin       = new InputVar <double>("Wood Percent Lignin");
            InputVar <double> crLignin      = new InputVar <double>("Coarse Root Percent Lignin");
            InputVar <double> frLignin      = new InputVar <double>("Fine Root Percent Lignin");
            InputVar <double> leafCN        = new InputVar <double>("Leaf CN Ratio");
            InputVar <double> woodCN        = new InputVar <double>("Wood CN Ratio");
            InputVar <double> cRootCN       = new InputVar <double>("Coarse Root CN Ratio");
            InputVar <double> foliarCN      = new InputVar <double>("Foliage CN Ratio");
            InputVar <double> fRootCN       = new InputVar <double>("Fine Root CN Ratio");
            string            lastColumn    = "the " + fRootCN.Name + " column";

            while (!AtEndOfInput && CurrentName != Names.FunctionalGroupParameters)
            {
                StringReader currentLine = new StringReader(CurrentLine);
                ISpecies     species     = ReadSpecies(currentLine);

                ReadValue(ft, currentLine);
                parameters.SetFunctionalType(species, ft.Value);

                ReadValue(nt, currentLine);
                parameters.NFixer[species] = nt.Value;

                ReadValue(gddmn, currentLine);
                parameters.SetGDDmin(species, gddmn.Value);

                ReadValue(gddmx, currentLine);
                parameters.SetGDDmax(species, gddmx.Value);

                ReadValue(mjt, currentLine);
                parameters.SetMinJanTemp(species, mjt.Value);

                ReadValue(maxd, currentLine);
                parameters.SetMaxDrought(species, maxd.Value);

                ReadValue(leafLongevity, currentLine);
                parameters.SetLeafLongevity(species, leafLongevity.Value);

                ReadValue(epicorm, currentLine);
                parameters.Epicormic[species] = epicorm.Value;

                ReadValue(leafLignin, currentLine);
                parameters.SetLeafLignin(species, leafLignin.Value);

                ReadValue(frLignin, currentLine);
                parameters.SetFineRootLignin(species, frLignin.Value);

                ReadValue(wLignin, currentLine);
                parameters.SetWoodLignin(species, wLignin.Value);

                ReadValue(crLignin, currentLine);
                parameters.SetCoarseRootLignin(species, crLignin.Value);

                ReadValue(leafCN, currentLine);
                parameters.SetLeafCN(species, leafCN.Value);

                ReadValue(fRootCN, currentLine);
                parameters.SetFineRootCN(species, fRootCN.Value);

                ReadValue(woodCN, currentLine);
                parameters.SetWoodCN(species, woodCN.Value);

                ReadValue(cRootCN, currentLine);
                parameters.SetCoarseRootCN(species, cRootCN.Value);

                ReadValue(foliarCN, currentLine);
                parameters.SetFoliageLitterCN(species, foliarCN.Value);

                CheckNoDataAfter(lastColumn, currentLine);
                GetNextLine();
            }

            //--------- Read In Functional Group Table -------------------------------
            PlugIn.ModelCore.UI.WriteLine("   Begin parsing FUNCTIONAL GROUP table.");

            ReadName(Names.FunctionalGroupParameters);
            string InitialEcoregionParameters = "InitialEcoregionParameters";

            InputVar <string> ftname             = new InputVar <string>("Name");
            InputVar <int>    ftindex            = new InputVar <int>("Index (< 25)");
            InputVar <double> ppdf1              = new InputVar <double>("PPDF(1)");
            InputVar <double> ppdf2              = new InputVar <double>("PPDF(2)");
            InputVar <double> ppdf3              = new InputVar <double>("PPDF(3)");
            InputVar <double> ppdf4              = new InputVar <double>("PPDF(4)");
            InputVar <double> fcfleaf            = new InputVar <double>("FCFRAC: Leaf");
            InputVar <double> btolai             = new InputVar <double>("BTOLAI");
            InputVar <double> klai               = new InputVar <double>("KLAI");
            InputVar <double> maxlai             = new InputVar <double>("MAXLAI");
            InputVar <double> mwm                = new InputVar <double>("Monthly Wood Mortality");
            InputVar <double> wdr                = new InputVar <double>("Wood Decay Rate");
            InputVar <double> mortCurveShapeParm = new InputVar <double>("Mortality Curve Shape Parameter");
            InputVar <int>    leafNeedleDrop     = new InputVar <int>("Leaf or Needle Drop Month");

            InputVar <double> ppr2 = new InputVar <double>("PPRPTS2");
            InputVar <double> ppr3 = new InputVar <double>("PPRPTS3");
            InputVar <double> coarseRootFraction = new InputVar <double>("CRootFrac");
            InputVar <double> fineRootFraction   = new InputVar <double>("FRootFrac");

            while (!AtEndOfInput && CurrentName != InitialEcoregionParameters)
            {
                StringReader currentLine = new StringReader(CurrentLine);

                ReadValue(ftname, currentLine);

                ReadValue(ftindex, currentLine);
                int ln = (int)ftindex.Value.Actual;

                if (ln >= numFunctionalTypes)
                {
                    throw new InputValueException(ftindex.Value.String,
                                                  "The index:  {0} exceeds the allowable number of functional groups, {1}",
                                                  ftindex.Value.String, numFunctionalTypes - 1);
                }


                //IEditableFunctionalType funcTParms = new EditableFunctionalType();
                FunctionalType funcTParms = new FunctionalType();
                parameters.FunctionalTypes[ln] = funcTParms;

                ReadValue(ppdf1, currentLine);
                funcTParms.PPDF1 = ppdf1.Value;

                ReadValue(ppdf2, currentLine);
                funcTParms.PPDF2 = ppdf2.Value;

                ReadValue(ppdf3, currentLine);
                funcTParms.PPDF3 = ppdf3.Value;

                ReadValue(ppdf4, currentLine);
                funcTParms.PPDF4 = ppdf4.Value;

                ReadValue(fcfleaf, currentLine);
                funcTParms.FCFRACleaf = fcfleaf.Value;

                ReadValue(btolai, currentLine);
                funcTParms.BTOLAI = btolai.Value;

                ReadValue(klai, currentLine);
                funcTParms.KLAI = klai.Value;

                ReadValue(maxlai, currentLine);
                funcTParms.MAXLAI = maxlai.Value;

                ReadValue(ppr2, currentLine);
                funcTParms.PPRPTS2 = ppr2.Value;

                ReadValue(ppr3, currentLine);
                funcTParms.PPRPTS3 = ppr3.Value;

                ReadValue(wdr, currentLine);
                funcTParms.WoodDecayRate = wdr.Value;

                ReadValue(mwm, currentLine);
                funcTParms.MonthlyWoodMortality = mwm.Value;

                ReadValue(mortCurveShapeParm, currentLine);
                funcTParms.MortCurveShape = mortCurveShapeParm.Value;

                ReadValue(leafNeedleDrop, currentLine);
                funcTParms.LeafNeedleDrop = leafNeedleDrop.Value;

                ReadValue(coarseRootFraction, currentLine);
                funcTParms.CoarseRootFraction = coarseRootFraction.Value;

                ReadValue(fineRootFraction, currentLine);
                funcTParms.FineRootFraction = fineRootFraction.Value;

                //PlugIn.ModelCore.UI.WriteLine("PPRPTS2={0}.", parameters.FunctionalTypeTable[ln].PPRPTS2);

                CheckNoDataAfter("the " + fineRootFraction.Name + " column", currentLine);
                GetNextLine();
            }


            //--------- Read In FIRST Ecoregion Table ---------------------------
            PlugIn.ModelCore.UI.WriteLine("   Begin reading INITIAL ECOREGION parameters.");
            ReadName(InitialEcoregionParameters);

            InputVar <string>        ecoregionName = new InputVar <string>("Ecoregion");
            InputVar <double>        iS1surfC      = new InputVar <double>("Initial SOM1 surface C");
            InputVar <double>        iS1surfN      = new InputVar <double>("Initial SOM1 surface N");
            InputVar <double>        iS1soilC      = new InputVar <double>("Initial SOM1 soil C");
            InputVar <double>        iS1soilN      = new InputVar <double>("Initial SOM1 soil N");
            InputVar <double>        iS2C          = new InputVar <double>("Initial SOM2 (intermediate turnover) C");
            InputVar <double>        iS2N          = new InputVar <double>("Initial SOM2 (intermediate turnover) N");
            InputVar <double>        iS3C          = new InputVar <double>("Initial SOM3 (slow turnover) C");
            InputVar <double>        iS3N          = new InputVar <double>("Initial SOM3 (slow turnover) N");
            InputVar <double>        iMN           = new InputVar <double>("Initial Mineral (available) N");
            Dictionary <string, int> lineNumbers2  = new Dictionary <string, int>();

            while (!AtEndOfInput && CurrentName != Names.EcoregionParameters)
            {
                StringReader currentLine = new StringReader(CurrentLine);

                ReadValue(ecoregionName, currentLine);

                IEcoregion ecoregion = GetEcoregion(ecoregionName.Value,
                                                    lineNumbers2);

                ReadValue(iS1surfC, currentLine);
                parameters.SetInitSOM1surfC(ecoregion, iS1surfC.Value);

                ReadValue(iS1surfN, currentLine);
                parameters.SetInitSOM1surfN(ecoregion, iS1surfN.Value);

                ReadValue(iS1soilC, currentLine);
                parameters.SetInitSOM1soilC(ecoregion, iS1soilC.Value);

                ReadValue(iS1soilN, currentLine);
                parameters.SetInitSOM1soilN(ecoregion, iS1soilN.Value);

                ReadValue(iS2C, currentLine);
                parameters.SetInitSOM2C(ecoregion, iS2C.Value);

                ReadValue(iS2N, currentLine);
                parameters.SetInitSOM2N(ecoregion, iS2N.Value);

                ReadValue(iS3C, currentLine);
                parameters.SetInitSOM3C(ecoregion, iS3C.Value);

                ReadValue(iS3N, currentLine);
                parameters.SetInitSOM3N(ecoregion, iS3N.Value);

                ReadValue(iMN, currentLine);
                parameters.SetInitMineralN(ecoregion, iMN.Value);

                CheckNoDataAfter("the " + iMN.Name + " column", currentLine);

                GetNextLine();
            }

            //--------- Read In SECOND Ecoregion Table ---------------------------
            // First, read table of additional parameters for ecoregions
            PlugIn.ModelCore.UI.WriteLine("   Begin reading FIXED ECOREGION parameters.");
            ReadName(Names.EcoregionParameters);

            InputVar <double> pclay  = new InputVar <double>("Percent Clay");
            InputVar <double> psand  = new InputVar <double>("Percent Sand");
            InputVar <int>    sd     = new InputVar <int>("Soil Depth");
            InputVar <double> fc     = new InputVar <double>("Field Capacity");
            InputVar <double> wp     = new InputVar <double>("Wilting Point");
            InputVar <double> sff    = new InputVar <double>("Storm Flow Fraction");
            InputVar <double> bff    = new InputVar <double>("Base Flow Fraction");
            InputVar <double> drain  = new InputVar <double>("Drain Fraction");
            InputVar <double> ans    = new InputVar <double>("Atmospheric N Slope");
            InputVar <double> ani    = new InputVar <double>("Atmospheric N Intercept");
            InputVar <double> lat    = new InputVar <double>("Latitude");
            InputVar <double> drsoms = new InputVar <double>("Decay Rate Surf");
            InputVar <double> drsom1 = new InputVar <double>("Decay Rate SOM1");
            InputVar <double> drsom2 = new InputVar <double>("Decay Rate SOM2");
            InputVar <double> drsom3 = new InputVar <double>("Decay Rate SOM3");
            InputVar <double> denits = new InputVar <double>("Denitrification");

            Dictionary <string, int> lineNumbers = new Dictionary <string, int>();

            while (!AtEndOfInput && CurrentName != Names.FireReductionParameters)
            {
                StringReader currentLine = new StringReader(CurrentLine);

                ReadValue(ecoregionName, currentLine);

                IEcoregion ecoregion = GetEcoregion(ecoregionName.Value,
                                                    lineNumbers);

                ReadValue(sd, currentLine);
                parameters.SetSoilDepth(ecoregion, sd.Value);

                ReadValue(pclay, currentLine);
                parameters.SetPercentClay(ecoregion, pclay.Value);

                ReadValue(psand, currentLine);
                parameters.SetPercentSand(ecoregion, psand.Value);

                ReadValue(fc, currentLine);
                parameters.SetFieldCapacity(ecoregion, fc.Value);

                ReadValue(wp, currentLine);
                parameters.SetWiltingPoint(ecoregion, wp.Value);

                ReadValue(sff, currentLine);
                parameters.SetStormFlowFraction(ecoregion, sff.Value);

                ReadValue(bff, currentLine);
                parameters.SetBaseFlowFraction(ecoregion, bff.Value);

                ReadValue(drain, currentLine);
                parameters.SetDrain(ecoregion, drain.Value);

                ReadValue(ans, currentLine);
                parameters.SetAtmosNslope(ecoregion, ans.Value);

                ReadValue(ani, currentLine);
                parameters.SetAtmosNintercept(ecoregion, ani.Value);

                ReadValue(lat, currentLine);
                parameters.SetLatitude(ecoregion, lat.Value);

                ReadValue(drsoms, currentLine);
                parameters.SetDecayRateSurf(ecoregion, drsoms.Value);

                ReadValue(drsom1, currentLine);
                parameters.SetDecayRateSOM1(ecoregion, drsom1.Value);

                ReadValue(drsom2, currentLine);
                parameters.SetDecayRateSOM2(ecoregion, drsom2.Value);

                ReadValue(drsom3, currentLine);
                parameters.SetDecayRateSOM3(ecoregion, drsom3.Value);

                ReadValue(denits, currentLine);
                parameters.SetDenitrif(ecoregion, denits.Value);


                CheckNoDataAfter("the " + drsom3.Name + " column", currentLine);

                GetNextLine();
            }
            //--------- Read In Fire Reductions Table ---------------------------
            PlugIn.ModelCore.UI.WriteLine("   Begin reading FIRE REDUCTION parameters.");
            ReadName(Names.FireReductionParameters);

            InputVar <int>    frindex = new InputVar <int>("Fire Severity Index MUST = 1-5");
            InputVar <double> wred    = new InputVar <double>("Wood Reduction");
            InputVar <double> lred    = new InputVar <double>("Litter Reduction");

            while (!AtEndOfInput && CurrentName != Names.MonthlyMaxNPP &&
                   CurrentName != Names.HarvestReductionParameters)
            {
                StringReader currentLine = new StringReader(CurrentLine);

                ReadValue(frindex, currentLine);
                int ln = (int)frindex.Value.Actual;

                if (ln < 1 || ln > 5)
                {
                    throw new InputValueException(ftindex.Value.String,
                                                  "The fire severity index:  {0} must be 1-5,",
                                                  frindex.Value.String);
                }


                FireReductions inputFireReduction = new FireReductions();  // ignoring severity = zero
                parameters.FireReductionsTable[ln] = inputFireReduction;

                ReadValue(wred, currentLine);
                inputFireReduction.WoodReduction = wred.Value;

                ReadValue(lred, currentLine);
                inputFireReduction.LitterReduction = lred.Value;

                CheckNoDataAfter("the " + lred.Name + " column", currentLine);

                GetNextLine();
            }
            //--------- Read In Harvest Reductions Table ---------------------------
            InputVar <string> hreds = new InputVar <string>("HarvestReductions");

            ReadOptionalName(Names.HarvestReductionParameters);
            {
                PlugIn.ModelCore.UI.WriteLine("   Begin reading HARVEST REDUCTION parameters.");
                InputVar <string> prescriptionName = new InputVar <string>("Prescription");
                InputVar <double> wred_pr          = new InputVar <double>("Wood Reduction");
                InputVar <double> lred_pr          = new InputVar <double>("Litter Reduction");

                //lineNumbers.Clear();
                List <string> prescriptionNames = new List <string>();
                //Dictionary<int, int> DisturbanceTypeLineNumbers = new Dictionary<int, int>();

                while (!AtEndOfInput && CurrentName != Names.MonthlyMaxNPP)
                {
                    HarvestReductions harvReduction = new HarvestReductions();
                    parameters.HarvestReductionsTable.Add(harvReduction);

                    StringReader currentLine = new StringReader(CurrentLine);

                    ReadValue(prescriptionName, currentLine);
                    harvReduction.PrescriptionName = prescriptionName.Value;

                    ReadValue(wred_pr, currentLine);
                    harvReduction.WoodReduction = wred.Value;

                    ReadValue(lred_pr, currentLine);
                    harvReduction.LitterReduction = lred.Value;

                    GetNextLine();
                }
            }

            //---------------------------------------------------------------------

            ParseBiomassParameters(parameters, Names.AgeOnlyDisturbanceParms,
                                   Names.DynamicChange);

            string            lastParameter           = null;
            InputVar <string> ageOnlyDisturbanceParms = new InputVar <string>(Names.AgeOnlyDisturbanceParms);

            ReadVar(ageOnlyDisturbanceParms);
            parameters.AgeOnlyDisturbanceParms = ageOnlyDisturbanceParms.Value;


            //  Climate Change table (optional)
            if (ReadOptionalName(Names.DynamicChange))
            {
                ReadDynamicTable(parameters.DynamicUpdates);
            }
            else if (lastParameter != null)
            {
                CheckNoDataAfter(lastParameter);
            }

            return(parameters); //.GetComplete();
        }