---basic_daisyworld_random.nlogo

来自「NETLOGO」· NLOGO 代码 · 共 885 行 · 第 1/2 页

globals
[ 
;;all of the areas are in percent of screen
 a_b   ;area covered by black daisies
 a_w   ;area covered by white daisies
 a_un  ;area of barren land
 
; Al_b   ;albedo of black daisies
; Al_w   ;albedo of white daisies
; Al_un  ;albedo of barren land
 Al_p   ;albedo of whole planet
 
 g_w   ;growth factor of white daisies
 g_b   ;growth factor of black daisies
 
 T_p   ;planetary temperature
 T_b   ;temperature around black daisies
 T_w   ;temperature around white daisies
 T_dead;temperature of a dead planet
 
 F_HA  ;heat absorption factor
 ;death-rate
 SB_const ;Stefan-Boltzmann constant
 SF_const ;solar flux constant
 ;SL       ;solar luminosity
 worldarea;area of the world
 
 time     ;counter of time
 
 success ;global for communicating with patches 
]

patches-own
[
  status?
]
to init
ca
;;set up areas
set a_un 1
set a_b  0
set a_w  0

;;initialize albedos
set Al_un 0.5
set Al_b  0.25
set Al_w  0.75
set Al_p  0.5 ;no daisies 

;;initalise ideal temperatures for daisies
;set idealblacktemp 22
;set idealwhitetemp 22

;;initialize constants
set F_HA 20
set death-rate 0.3
set SB_const 5.669e-8;W/degK^4
set SF_const 917;W/m^2

;;initalize temperature and solar luminosity
set SL 0.6
set T_dead (((SL * SF_const * (1 - Al_p) / SB_const) ^ 0.25 ) - 273) ;temperature of a dead planet (degC)
ifelse(autotemp)
[
set T_p T_dead
]
[
set T_p 20
set SL 0.6
]


set time 0

set worldarea 4 * ( screen-edge-x + 1 ) * (screen-edge-y + 1)
ask patches[set status? "barren"]
end


to iterate

;;recalculate luminosity and global temp
if(autotemp)
[ set SL (0.6 + (time * (0.001))) ] ;recalculate luminosity
ifelse(a_un < 1)
[
set T_p ((SL * SF_const * (1 - Al_p) / SB_const) ^ 0.25) - 273 ;recalculate global temperature
]
[
set T_p (((SL * SF_const * (1 - 0.5) / SB_const) ^ 0.25 ) - 273) ;temperature of a dead planet (degC)
]
set T_dead (((SL * SF_const * (1 - 0.5) / SB_const) ^ 0.25 ) - 273) ;temperature of a dead planet (degC)
;;recalculate local temperatures
set T_b (F_HA * (Al_p - Al_b) + T_p)
set T_w (F_HA * (Al_p - Al_w) + T_p)

set g_b (1 - 0.003265 * ((idealblacktemp - T_b) ^ 2))
set g_w (1 - 0.003265 * ((idealwhitetemp - T_w) ^ 2))


grow-daisies

let da_b a_b * (a_un * g_b - death-rate) + 0.01
let da_w a_w * (a_un * g_w - death-rate) + 0.01

;;recalculate planetary albedo and areas 
set a_b (count patches with[status? = "black"]) / worldarea
set a_w (count patches with[status? = "white"]) / worldarea
set a_un 1 - a_b - a_w

set Al_p ((a_un) * (Al_un) + (a_b) * (Al_b) + (a_w) * (Al_w))

;;increment time, color daisies
set time (time + 1)
ask patches[color-daisies]


set-current-plot "Daisies"
set-plot-pen-interval 1
set-current-plot-pen "black"
plot a_b
set-current-plot-pen "white"
plot a_w
set-current-plot-pen "barren"
plot a_un

set-current-plot "Temperature"
set-plot-pen-interval 1
set-current-plot-pen "temp"
plot T_p
set-current-plot-pen "temp-without-daisies"
plot T_dead

;file-open "c:\\daisyms.out"
;file-write time
;file-write a_b file-write a_w file-write a_un file-write T_p file-write T_dead
;file-print " "
;file-close

end




to grow-daisies

if(a_b = 0) and (g_b > 0)
[
start-new-black-patch
]
if(a_w = 0) and (g_w > 0)
[
start-new-white-patch
]

let da_b a_b * (a_un * g_b - death-rate) + 0.001
let da_w a_w * (a_un * g_w - death-rate) + 0.001
ifelse(da_b > 0)
[
ask patches with[status? = "black"] [ grow-black-daisies ]
]
[
ask patches with[status? = "black"] [if(random-float 1 < (- da_b))[set status? "barren"]]
]

ifelse(da_w > 0)
[
ask patches with[status? = "white"] [ grow-white-daisies ]
]
[
ask patches with[status? = "white"] [if(random-float 1 < (- da_w))[set status? "barren"]]
]

end



to grow-black-daisies
let numblack (count patches with[status? = "black"]) 
let da_b a_b * (a_un * g_b - death-rate) + 0.001
set da_b da_b *  (100 / 8)

ifelse (count neighbors with [not (status? = "barren")] > spread-thresh )
[
 if(random-float 1 < da_b)
 [start-new-black-patch]
]
[
  ask neighbors 
  [
  if(random-float 1 < da_b and status? = "barren")
  [
  set status? "black"
  ]
  
  ]
]

end

to grow-white-daisies
let numwhite (count patches with[status? = "white"]) 
let da_w a_w * (a_un * g_w - death-rate) + 0.001
set da_w da_w * (100 / (8))
ifelse (count neighbors with [not (status? = "barren")] > spread-thresh )
[
if(random-float 1 < da_w)
 [start-new-white-patch]
]
[
  ask neighbors 
  [
  if(random-float 1 < da_w and status? = "barren")
  [
  set status? "white"
  ]
  
  ]
]

end

to start-new-black-patch
set success false
  ;;pick a random x coordinate
  let currxcor (random (2 * screen-edge-x)) - screen-edge-x
  ;;pick a random y coordinate
  let currycor (random (2 * screen-edge-y)) - screen-edge-y
  ask patches with[ pxcor = currxcor and pycor = currycor and status? = "barren"]
  [
    set status? "black"
    set success true
  ]
  ifelse (success)
  [
    ;;do nothing
  ]
  [
    start-new-black-patch
  ]


end

to start-new-white-patch
set success false
  ;;pick a random x coordinate
  let currxcor (random (2 * screen-edge-x)) - screen-edge-x
  ;;pick a random y coordinate
  let currycor (random (2 * screen-edge-y)) - screen-edge-y
  ask patches with[ pxcor = currxcor and pycor = currycor and status? = "barren"]
  [
    set status? "white"
    set success true
  ]
  ifelse (success)
  [
    ;;do nothing
  ]
  [
    start-new-white-patch
  ]

end

to color-daisies

ifelse (status? = "black")
[
set pcolor black
]
[
  ifelse (status? = "white")
  [
  set pcolor white
  ]
  [
  set pcolor brown
  ]
]

end
@#$#@#$#@
GRAPHICS-WINDOW
207
61
522
397
30
30
5.0
1
10
1
1
1
0
1
1
1

CC-WINDOW
5
485
784
580
Command Center
0

BUTTON
452
16
515
49
NIL
init\n
NIL
1
T
OBSERVER
T
NIL

MONITOR
608
11
682
60
Black Area
a_b\n
3
1

MONITOR
688
10
763
59
White Area
a_w
3
1

MONITOR
434
408
517
457
Planet Temp.
T_p
3
1

PLOT
530
61
775
263
Daisies
NIL
NIL
0.0
10.0
0.0
1.0
true
false
PENS
"black" 1.0 0 -16777216 true
"white" 1.0 0 -5325092 true
"barren" 1.0 0 -6459832 true

PLOT
531
268
775
471
Temperature
time ->
Temperature (C)
0.0
10.0
-20.0
100.0
true
false
PENS
"temp" 1.0 0 -10899396 true
"temp-without-daisies" 1.0 0 -6459832 true

BUTTON
375
16
441
49
NIL
iterate\n
T
1
T
OBSERVER
T
NIL

SLIDER
23
22
195
55
idealblacktemp
idealblacktemp
0
100
22
1
1
NIL

SLIDER
23
59
195
92
idealwhitetemp
idealwhitetemp
0
100
22
1
1
NIL

SLIDER
25
117
197
150
Al_w
Al_w
0
1
0.75
0.01