The program pos_3d

Das Programm pos_3d bietet die Möglichkeit, verschiedenen Standard-Initialisierungsdateien zu erzeugen. In der Konfigurationsdatei kann zwischen vier unterschiedlichen Verteilungen der Startpositionen gewählt werden:

  1. Es werden für mehrere Theta-Level Startpunkte gewählt, die auf einem polar-stereografischen Plot für jedes Level gleichverteilt sind. Hierzu muss der Benutzer die Hemisphäre angeben und den Breitengrad bis zu dem Startpositionen erzeugt werden sollen. Bei der grafischen Ausgabe wird später der entsprechende Pol im Zentrum der Karte liegen und der angegebene Breitengrad am Rand der Grafik. Der Benutzer gibt die Anzahl der gewünschten Startpositionen entlang eines Längengrades (=nx) an. Es werden für jedes angegebene Theta-Level

    $$nx*nx$$

    Punkte untersucht, die im quadratischen Plot-Ausschnitt gleichverteilt sind. Es werden jedoch nur die Punkte als Startpositionen genommen, die innerhalb des durch den Breitengrad vorgegebenen Kreises (bei polar-stereographischer Darstellung) liegen. Die maximale Anzahl ermittelter Startpositionen ist daher

    $$(nlev) * nx * nx * \pi * 0.5^2$$, wobei nlev der Anzahl der Theta-Level entspricht.

  2. Es werden Startpositionen erzeugt, die auf einem Längen- und Breitengrad-Gitter gleichverteilt sind. Derart verteilte Startpositionen sind für eine spätere zylindisch-äquidistante grafische Ausgabe geeignet. Der Benutzer gibt zwei Breitengrade und zwei Längengrade, die den zu betrachtenden Bereich definieren, und die gewünschte Anzahl der Punkte entlang eines Längengrades (=nx) und der Punkte entlang eines Breitengrades (=ny) an. Es werden

    $$nx*ny$$ Startpositionen erzeugt, die in dem durch die Längen- und Breitengrade vorgegebenen Rechteck liegen.

  3. Es werden Startpositionen entlang eines Längengrades erzeugt. Es wird die vorgegebene Anzahl von Startpositionen erzeugt, wobei diese auf dem angegebenen Längengrad zwischen den beiden angegebenen Breitengraden gleichmäßig verteilt werden. Für eine spätere grafische Ausgabe mit traj_quer.pro (Querschnittplot entlang eines Längengrades) ist ein solcher Eingabedatensatz für das Trajektorienprogramm erforderlich.

  4. Es werden Startpositionen entlang eines Breitengrades erzeugt. Es wird eine vorgegebene Anzahl von Startpositionen erzeugt, die auf dem angegebenen Breitengrad zwischen den beiden angegebenen Längengraden gleichmäßig verteilt sind.

The options for the execution of pos_3d are set in the configuration file pos_3d.inp:

ABC                     ! Prefix for output filename
1995, 3, 16, 12, 0      ! time (year, month, day, hour, min)
theta                   ! Name of vertical coordinate
3                       ! number of theta-levels
320, 340, 360           ! theta-levels
1                       ! 1=polar-stereogr.,2=cyl.-equidistant,3=cross-section(longitude),4=cross-section (latitude)
1                       ! Centre for polar-stereogr. Plots, 1=North Pole, -1=South Pole
30                      ! edge latitude
10                      ! number of particles across latitudes
icg112                  ! userid
n.thomas                ! PI

ABC                     ! Prefix for output filename
1995, 3, 16, 12, 0      ! time (year, month, day, hour, min)
theta                   ! Name of vertical coordinate
3                       ! number of theta-levels
320, 340, 360           ! theta-levels
2                       ! 1=polar-stereogr.,2=cyl.-equidistant,3=cross-section(longitude),4=cross-section (latitude)
30 60                   ! edge latitudes
-10 20                  ! edge longitudes
10 10                   ! number of latitudes, number of longitudes
icg112                  ! userid
n.thomas                ! PI

ABC                     ! Prefix for output filename
1995, 3, 16, 12, 0      ! time (year, month, day, hour, min)
theta                   ! Name of vertical coordinate
3                       ! number of theta-levels
320, 340, 360           ! theta-levels
3                       ! 1=polar-stereogr.,2=cyl.-equidistant,3=cross-section(longitude),4=cross-section (latitude)

30 60                   ! edge latitudes
-10                     ! longitude
10                      ! number of particles across latitudes
icg112                  ! userid
n.thomas                ! PI

ABC                     ! Prefix for output filename
1995, 3, 16, 12, 0      ! time (year, month, day, hour, min)
theta                   ! Name of vertical coordinate
3                       ! number of theta-levels
320, 340, 360           ! theta-levels
4                       ! 1=polar-stereogr.,2=cyl.-equidistant,3=cross-section(longitude),4=cross-section (latitude)

-10 10                  ! edge longitudes
60                      ! latitude
10                      ! number of particles across longitudes
icg112                  ! userid
n.thomas                ! PI

pos/pos_3d.f90 (last edited 2010-04-29 08:25:17 by NicoleThomas)