示例#1
0
    /// <summary>
    /// Respond to a remove biomass event.
    /// </summary>
    internal void OnRemoveBiomass(double removeBiomPheno)
    {
        string existingStage = CurrentStageName;

        if (RewindDueToBiomassRemoved != null)
        {
            FractionBiomassRemoved = removeBiomPheno; // The RewindDueToBiomassRemoved function will use this.

            double ttCritical       = TTInAboveGroundPhase;
            double removeFractPheno = RewindDueToBiomassRemoved.Value;
            double removeTTPheno    = ttCritical * removeFractPheno;

            string msg;
            msg  = "Phenology change:-\r\n";
            msg += "    Fraction DM removed  = " + removeBiomPheno.ToString() + "\r\n";
            msg += "    Fraction TT removed  = " + removeFractPheno.ToString() + "\r\n";
            msg += "    Critical TT          = " + ttCritical.ToString() + "\r\n";
            msg += "    Remove TT            = " + removeTTPheno.ToString() + "\r\n";

            double ttRemaining = removeTTPheno;
            for (int i = Phases.Count - 1; i >= 0; i--)
            {
                Phase Phase = Phases[i];
                if (Phase.TTinPhase > 0)
                {
                    double ttCurrentPhase = Phase.TTinPhase;
                    if (ttRemaining > ttCurrentPhase)
                    {
                        Phase.ResetPhase();
                        ttRemaining       -= ttCurrentPhase;
                        CurrentPhaseIndex -= 1;
                        if (CurrentPhaseIndex < 4)  //FIXME - hack to stop onEmergence being fired which initialises biomass parts
                        {
                            CurrentPhaseIndex = 4;
                            break;
                        }
                    }
                    else
                    {
                        Phase.Add(-ttRemaining);
                        // Return fraction of thermal time we are through the current
                        // phenological phase (0-1)
                        //double frac = Phase.FractionComplete;
                        //if (frac > 0.0 && frac < 1.0)  // Don't skip out of this stage - some have very low targets, eg 1.0 in "maturity"
                        //    currentStage = frac + floor(currentStage);

                        break;
                    }
                }
                else
                { // phase is empty - not interested in it
                }
            }
            Stage = (CurrentPhaseIndex + 1) + CurrentPhase.FractionComplete;

            if (existingStage != CurrentStageName)
            {
                PhaseChangedType PhaseChangedData = new PhaseChangedType();
                PhaseChangedData.OldPhaseName = existingStage;
                PhaseChangedData.NewPhaseName = CurrentPhase.Name;
                PhaseChanged.Invoke(PhaseChangedData);
                MyPaddock.Publish(CurrentPhase.Start);
            }
        }
    }