Simple empirical fire model - Revision history

Fei: /* Succession.lse */

2006-11-14T15:53:41Z

Succession.lse

← Older revision		Revision as of 15:53, 14 November 2006
Line 126:		Line 126:
	TRANSITIONS = TRUE		TRANSITIONS = TRUE
	TimeSinceFire = MIN(TimeSinceFire + 1, MaxTSF)		TimeSinceFire = MIN(TimeSinceFire + 1, MaxTSF)
	ENDTR		ENDTR

	==Suggested Experiments==		==Suggested Experiments==

	To explore this cellular automata model further, try the following:		To explore this cellular automata model further, try the following:

Fei: /* Succession.lse */

2006-11-14T15:53:05Z

Succession.lse

← Older revision		Revision as of 15:53, 14 November 2006
Line 115:		Line 115:
	ENDNR		ENDNR
	===Succession.lse===		===Succession.lse===
	LSEVENT: SimpleAging		LSEVENT: SimpleAging
			DEFINITIONS
	DEFINITIONS
	GLOBAL CONSTANT: MaxTSF		GLOBAL CONSTANT: MaxTSF

	LAYER: TimeSinceFire		LAYER: TimeSinceFire
	ENDDEF		ENDDEF
			// Succession is modelled every year
	// Succession is modelled every year		RETURNTIME = 365.25
	RETURNTIME = 365.25		// Aging simply increments TSF by one year!
			TRANSITIONS
	// Aging simply increments TSF by one year!
	TRANSITIONS
	TRANSITIONS = TRUE		TRANSITIONS = TRUE
	TimeSinceFire = MIN(TimeSinceFire + 1, MaxTSF)		TimeSinceFire = MIN(TimeSinceFire + 1, MaxTSF)

Fei at 15:52, 14 November 2006

2006-11-14T15:52:28Z

@@ Line 114: / Line 114: @@
    NUMSPREADRECIPIENTS = CLAMP(NORMAL(0.5,1), 1,3)
   ENDNR
 ==Suggested Experiments==
 To explore this cellular automata model further, try the following:

Fei: /* Fire.lse */

2006-11-14T15:50:23Z

Fire.lse

← Older revision		Revision as of 15:50, 14 November 2006
Line 46:		Line 46:
	===Fire.lse===		===Fire.lse===
	SEVENT: SimpleFire		SEVENT: SimpleFire

	DEFINITIONS		DEFINITIONS
	LAYER: TimeSinceFire		LAYER: TimeSinceFire
	GLOBAL VARIABLE: ForestSize, MeanFireSize, MeanFiresPerYear, FireCycle, oldFireCycle		GLOBAL VARIABLE: ForestSize, MeanFireSize, MeanFiresPerYear, FireCycle, oldFireCycle

	CLUSTER VARIABLE: Extent		CLUSTER VARIABLE: Extent
	ENDDEF		ENDDEF

	RETURNTIME		RETURNTIME
	RETURNTIME = 365.25		RETURNTIME = 365.25

	// In case the user tries to change this...		// In case the user tries to change this...
	ForestSize = NUMCELLS		ForestSize = NUMCELLS

	impliedForestSize = MeanFireSize * MeanFiresPerYear * FireCycle		impliedForestSize = MeanFireSize * MeanFiresPerYear * FireCycle

	IF (impliedForestSize NEQ ForestSize)		IF (impliedForestSize NEQ ForestSize)
	// If the user didn't change fire cycle: adjust it.		// If the user didn't change fire cycle: adjust it.
Line 76:		Line 70:
	ENDFN		ENDFN
	oldFireCycle = FireCycle		oldFireCycle = FireCycle

	// Otherwise, the fire cycle was changed so adjust num fires instead		// Otherwise, the fire cycle was changed so adjust num fires instead
	ELSE		ELSE
Line 92:		Line 85:
	ENDFN		ENDFN
	ENDRT		ENDRT


	// Pick number of fires		// Pick number of fires
	// with random ignition		// with random ignition
	NUMCLUSTERS		NUMCLUSTERS
	NUMCLUSTERS = ROUND(NEGEXP(MeanFiresPerYear))		NUMCLUSTERS = ROUND(NEGEXP(MeanFiresPerYear))

	// For each pick a size		// For each pick a size
	Extent = ROUND(NEGEXP(MeanFireSize))		Extent = ROUND(NEGEXP(MeanFireSize))
	ENDNC		ENDNC


	TRANSITIONS		TRANSITIONS
	// Make a transition only if there is still extent to be burned		// Make a transition only if there is still extent to be burned
	// AND if the stand didn't just burn during this event		// AND if the stand didn't just burn during this event
	TRANSITIONS = (Extent >= 1) AND (TimeSinceFire > 0)		TRANSITIONS = (Extent >= 1) AND (TimeSinceFire > 0)

	// Set the time since fire to 0		// Set the time since fire to 0
	TimeSinceFire = 0		TimeSinceFire = 0

	// Decrement the number of cells remaining to burn for this opening		// Decrement the number of cells remaining to burn for this opening
	Extent = Extent - 1		Extent = Extent - 1
	ENDTR		ENDTR


	// Spread time not relevant		// Spread time not relevant
	SPREADTIMESTEP = 0		SPREADTIMESTEP = 0


	/* Allow spread to the 4 cardinal neighbours: up, down, left and right */		/* Allow spread to the 4 cardinal neighbours: up, down, left and right */
	SPREADLOCATION		SPREADLOCATION
Line 126:		Line 108:
	DECISION TimeSinceFire > 0		DECISION TimeSinceFire > 0
	ENDSL		ENDSL

	NUMSPREADRECIPIENTS		NUMSPREADRECIPIENTS
	/* Mid-Complex shapes. This is the number of neighbours to */		/* Mid-Complex shapes. This is the number of neighbours to */

Fei: /* Model.sel */

2006-11-14T15:49:44Z

Model.sel

← Older revision		Revision as of 15:49, 14 November 2006
Line 34:		Line 34:
	Global Constants:		Global Constants:
	MaxTSF = 200		MaxTSF = 200

	Output Model Bounds:		Output Model Bounds:
	TimeSinceFire: MaxTSF		TimeSinceFire: MaxTSF
Line 44:		Line 43:
	oldFireCycle = FireCycle OFF		oldFireCycle = FireCycle OFF
	Output Frequency: 365.250000		Output Frequency: 365.250000



	===Fire.lse===		===Fire.lse===

Fei: /* Fire.scn */

2006-11-14T15:49:34Z

Fire.scn

← Older revision		Revision as of 15:49, 14 November 2006
Line 19:		Line 19:
	Model Dimensions: initialTimeSinceFire		Model Dimensions: initialTimeSinceFire
	FireModel.sel		FireModel.sel

	Minimize Initial State		Minimize Initial State
	Tile		Tile

Fei at 15:49, 14 November 2006

2006-11-14T15:49:18Z

Show changes

Fei at 15:44, 14 November 2006

2006-11-14T15:44:53Z

← Older revision		Revision as of 15:44, 14 November 2006
Line 15:		Line 15:

	===Fire.scn===		===Fire.scn===
	~~Scenario Information~~		initialTimeSinceFire = ..\gisData\cell\initialTSF1
	initialTimeSinceFire = ..\gisData\cell\initialTSF1		Model Dimensions: initialTimeSinceFire
	Model Dimensions: initialTimeSinceFire		FireModel.sel
	FireModel.sel

	Minimize Initial State		Minimize Initial State
	Tile		Tile

	===Model.sel===		===Model.sel===
	Seles Model		Seles Model Information
	Time Units: ~~Step kiloStep 1000 1000~~		Time Units: Day Year 365.25
	Landscape Events:		Landscape Events:
	~~Producer~~.lse ~~DEBUG~~		Succession.lse
	~~Consumer~~.lse ~~DEBUG~~		Fire.lse
	~~Update~~.lse ~~DEBUG~~		ReportResults.lse
	~~Spatial Variables~~:		Variable-Input View Maps:
	~~LSState[~~3]		Species = initialSpecies
	~~PrevLSState[~~3]		TimeSinceFire = initialTimeSinceFire
	~~Legends~~:		Variable-Output View Maps:
	~~LSStateTypes = {0~~:~~Uncolonised, 1:Producer, 2~~:~~Consumer}~~		Species
			TimeSinceFire
			Global Constants:
			MaxTSF = 200

			rAspen = 1
			rPine = 2
			rSpruce = 3
			rBalsam = 4
			Aspen = 1
			Pine = 2
			Spruce = 3
			Balsam = 4

			Output Model Bounds:
			TimeSinceFire: MaxTSF
			Species: 4
	Global Variables:		Global Variables:
	~~pInitialProducer~~ = ~~0.01~~		ForestSize = 250000
	~~pInitialConsumer~~ = ~~0.01~~		MeanFireSize = 100
	~~pSpreadProducer~~ = ~~0.1~~		MeanFiresPerYear = 10
	~~pSpreadConsumer~~ = 0.2		FireCycle = ForestSize /(MeanFireSize * MeanFiresPerYear)
	~~pDeathProducer~~ = 0.01		oldFireCycle = FireCycle OFF
	~~pDeathConsumer~~ = 0.05
	~~// SpreadRateProducer = 1~~		PatchBased = TRUE
	~~// SpreadRateConsumer~~ = 1		SProbAspen2Spruce = 0.1
	Output Frequency: 1		SProbPine = 0.01
			SProbPine2Spruce = 0.5

			RecordResults = TRUE
			Output Frequency: 365.250000


	===consumer.lse===		===consumer.lse===

Fei: /* Consume.scn */

2006-11-14T15:41:44Z

Consume.scn

← Older revision		Revision as of 15:41, 14 November 2006
Line 14:		Line 14:
	In the following sections we will examine all of the model files for this model. Note that instead of downloading the zip file above, you could just copy the text in the boxes below into a text editor and save it with the appropriate name (Section title). Opening the resulting .scn file in the SELES simulator would run this model.		In the following sections we will examine all of the model files for this model. Note that instead of downloading the zip file above, you could just copy the text in the boxes below into a text editor and save it with the appropriate name (Section title). Opening the resulting .scn file in the SELES simulator would run this model.

	===~~Consume~~.scn===		===Fire.scn===
	~~Seles~~ Scenario		Scenario Information
	Model.sel		initialTimeSinceFire = ..\gisData\cell\initialTSF1
	~~SimPriority Low Priority~~		Model Dimensions: initialTimeSinceFire
			FireModel.sel

			Minimize Initial State
			Tile

	===Model.sel===		===Model.sel===

Fei at 15:38, 14 November 2006

2006-11-14T15:38:21Z

← Older revision		Revision as of 15:38, 14 November 2006
Line 4:		Line 4:
	==Screenshot==		==Screenshot==

	[[Image:~~Consumer~~.gif]]		[[Image:TimeSinceFire.gif]]

	==Download Model==		==Download Model==

	Download the .scn, .sel and .lse files by clicking on the following link: [[Image:~~ProducerConsumer~~.zip]]		Download the .scn, .sel and .lse files by clicking on the following link: [[Image:SimpleFireModel.zip]]

	==Model Code Exploration==		==Model Code Exploration==