tonyday567 / daisy

daisy

Daisy World

Usage

import Daisy

Development

:r
:set -Wno-type-defaults
:set -Wno-name-shadowing
:set -XOverloadedStrings
:set -XOverloadedLabels
import Optics.Core
import Chart
import Daisy
import Chart.Examples
import Prettychart
(display, quit) <- startChartServer (Just "daisy")

display unitExample

check

cfg = defaultConfig
as = take 1000 $ alphas defaultConfig
b = tempB cfg . resTime <$> as
p = tempP cfg <$> as
a = (1000*) . albedoP cfg . resAlpha <$> as
sendChart $ mempty & set #hudOptions defaultHudOptions & set #chartTree (unnamed [simpleLineChart 0.01 (palette 0) a])
-- tempP cfg <$> (take 2 $ drop 100 $ alphas defaultConfig)
-- tempB cfg . resTime <$> (take 2 $ drop 100 $ alphas defaultConfig)

cfg = defaultConfig
as = take 1000 $ alphas defaultConfig
b = tempB cfg . resTime <$> as
p = tempP cfg <$> as
a = (1000*) . albedoP cfg . resAlpha <$> as
aW = (1000*) . w . resAlpha <$> as
aB = (1000*) . b . resAlpha <$> as
sendChart $ mempty & set #hudOptions defaultHudOptions & set #chartTree (unnamed [simpleLineChart 0.01 (palette 0) aW, simpleLineChart 0.01 (palette 1) aB])
-- tempP cfg <$> (take 2 $ drop 100 $ alphas defaultConfig)
-- tempB cfg . resTime <$> (take 2 $ drop 100 $ alphas defaultConfig)

a

cfg = defaultConfig
albedoP cfg (x0 cfg)
tP cfg 2 (x0 cfg)

sendChart unitExample

reference

https://en.wikipedia.org/wiki/Runge%E2%80%93Kutta_methods