geodynamics/aspect — pull request #108 by jaustermann (2014-05-30 17:54:45+00:00)
Title: S40RTS initial temperature field

Hi Wolfgang,

Ich schick dir hier einen pull request fuer das S40RTS Modell, welches man nutzen kann um das Anfangstemperaturfeld zu bestimmen. Angehaengt hab ich auch einen screenshot der den Vergeich zwischen dem originalen S40RTS (der matlab plot) und Aspect S40RTS (der paraview plot) zeigt. Du siehst es unterscheidet sich ein bisschen wegen der geringeren Aufloesung in Aspect und weil in Aspect eine Randbedingung fuer die Oberflaeche und core-mantle-boundary festgelegt ist (dort ist die perturbation 0). Ansonsten ist es gleich. Ich hab auch noch S20RTS mit reingenommen, d.h. man kann zwischen den zwei Tomographie Modellen waehlen.
Das berechnen dieser initial condition dauert relativ lange, v.a. wenn man zu hoeheren Aufloesungen geht. Das liegt daran, dass man fuer jeden grid point 18081 Komponenten aufsummieren muss. Fuer die cubic spline interpolation kann man nicht nur die umliegenden 4 Punkte benutzen, da man die Randbedingungen mit einbeziehen muss. Je nachdem was fuer eine Anwendung man hat gibt es Moeglichkeiten den Prozess zu vereinfachen / beschleunigen, aber das kommt eben auf die Anwendung an deswegen hab ich es jetzt mal beim allgemeinsten Fall gelassen.
Soll ich dir hierfuer noch einen test schreiben? Selbst bei der groebsten Aufloesung wuerde das wohl ein bisschen dauern, aber ich kann dir was schicken wenn du magst.
Ich lasse grad einen run mit hoher Aufloesung laufen und schick dir dann nochmal einen screenshot wenn das fertig ist. Sag mir bescheid falls ihr etwas davon fuer ein cookbook benutzen wollt, bzw. was ich dafuer machen soll.

Zu guter letzt ... Ich bin mir noch nicht ganz sicher ob ich das alles so richtig mach mit Git. Ich hoffe es passt. Du wirst sehen, dass ich am Anfang zweit test commits gemacht haben, die keine Veraenderung bringen sollten (ich konnte sie leider nicht aus meiner history loeschen). Ausserdem zeigt es an wenn ich den master branch merge, was passiert wenn ich mit 'git push origin master' meinen master branch auf github aktualisiere. Ist das so okay? Letztendlich geht es hier nur um den letzten commit.

Liebe Gruesse,
Jacky

![screen shot 2014-05-30 at 1 32 31 pm](https://cloud.githubusercontent.com/assets/7633875/3133991/761265c4-e820-11e3-924a-f305d28c471b.png)