Bose-Einstein-Kondensation in magnetischen und optischen Fallen

102 Anhang B. Beschreibung der verwendeten Programme
box_levels.f erstellt eine Datei mit den Eigenwerten im Kastenpotential.
kugel_levels.f ist für das Kugelpotential geschrieben worden. Die Nullstellen der benötigten
halbzahligen Besselfunktionen wurden vorher mit Mathematica berechnet und
in eine Datei bessel_half gespeichert.
hohlkugel_levels.f berechnet die Eigenwerte im Potential der Hohlkugel.
zylinder_levels.f benötigt die Nullstellen der ganzzahligen Besselfunktionen, die sich in
der Datei bessel_full befinden.
zylinder_levels_period.f berechnet die Energieeigenwerte, wenn der Zylinder periodische
Randbedingungen besitzt.
hohlzylinder_levels.f ermittelt die Werte für den Hohlzylinder.
ho_levels_3d.bin.f ist entscheidend für die Simulation der TOP-Falle, da es hier möglich
ist ein anisotropes harmonisches Oszillatorpotential anzugeben.
Die Sortierung der Ergebnisse erfolgt mit sort_zeros.f, in dem die Funktion dsortx der
ESSL Bibliothek von IBM entstammt. Siehe dazu [71].
B.1.3
howave.f und howavet.f
Sowohl die Wellenfunktion des harmonischen Oszillators, als auch die Teilchendichte in
der TOP-Falle berechnet das Programm howave.f. Erzeugt wird eine Ausgabedatei, die
die Dichte und die zugehörigen Raumkoordinaten enthält. Es unterscheidet sich von howavet.f
nur dadurch, daß dieses zusätzlich eine Zeitentwicklung wie sie in Abschnitt 5.4
benötigt wird, erlaubt. Die hierzu notwendige Fouriertransformation entstammt der ESSL-
Bibliothek von IBM [71].
B.2 Hilfsskripte und -programme
B.2.1
Sortierprogramme
Die Ausgabe von howave.f und howavet.f wurde mit der Visualisierungssoftware Gsharp
in anschauliche Bilder umgewandelt. Da dieses Softwarepaket ein ganz bestimmtes Format
der Datenfiles voraussetzt, existieren einige Shellskripte und Fortran-Programme, um
die Daten zu sortieren und zu konvertieren. denssort.f liest eine Datei mit Daten für einen
Raumquadranten und kopiert die entsprechenden Einträge in die anderen Quadranten. Dieses
ist möglich, da die betrachteten Potentiale symmetrisch sind. Anschließend werden
die ausgegebenen Daten mit dem UNIX sort Kommando aufsteigend sortiert. denszaxis.f
extrahiert die letzte Spalte der nun vorhandenen Ausgabedatei, deren Werte der Teilchendichte
an einem Ort entsprechen.