First plot of a basic daisyworld simulation

This is the first plot of a result of the daisyworld model in Fortran. It can be compared to one of the original article (need to cite properly later):

daisyworld_1

The wiggles in the plot are most probably related to the stability of the integration in time. We need to check the algorithm implemented in the program. This work has been done with Halime, Liv, Livia, and Kira.

Experiment zur Diffusion von Wärme

Entwurf:

  1. Die IR Kamera eichen:
  2. +-----------------------------+
    |  Kalt (Eis/Wasser Mischung) |
    +-----------------------------+
           t_1   |   | 
           t_2   | S |
           t_3   | t |                  |
           t_4   | a | (Metall)         | IR Kamera   
           t_... | b |                  | 
           t_n   |   |
    +-----------------------------+
    |  Warm (Geheiztes Wasser)    |
    +-----------------------------+
    
    • Den Stab mit einer Kamera aufnehmen, bis sich die Temperaturverteilung nicht mehr ändert (Bild ändert sich nicht mehr).
    • Temperaturprofil des Stabes messen und Graustufen oder Falschfarben des Bildes den Temperaturen zuordnen.
  3. Zeitliches Verhalten bestimmen:
  4. Einen Stab mit Raumtemperatur in das geheizte Wasser stellen und die Erwärmung filmen.

           t_1   |   | 
           t_2   | S |
           t_3   | t |                  |
           t_4   | a | (Metall)         | IR Kamera   
           t_... | b |                  | 
           t_n   |   |
    +-----------------------------+
    |  Warm (Geheiztes Wasser)    |
    +-----------------------------+
    

    Daraus sollte man eine Funktion dT/dt = f(x, T_u, dT/dx, t)
    ableiten können.

Daisyworld status

After implementing the Daisyworld model in a first version, we (Halime, Liv, and me) had to spend the last week for debugging.
We encountered a first problem with the Euler time integration scheme used to solve the differential equation. I setup an issue on Gitlab so that people can work on that.
Next Thursday we will start to work on upgrading the model to one dimension. For that I need the equation that describes the temperature diffusion between the 0d Daisyworlds. To get the equation, we will run an experiment: we heat an metal stick on one end and film it with an infrared camera.
The other aspect to look into is to use more then two colors of daisies.
On Gitlab are the corresponding issues.
http://klostercluster.dkrz.de/gitlab/s013014/daisyworld