Differences between revisions 4 and 7 (spanning 3 versions)
Revision 4 as of 2008-07-15 09:21:00
Size: 7861
Editor: NicoleThomas
Comment:
Revision 7 as of 2008-11-20 08:35:58
Size: 8387
Editor: NicoleThomas
Comment:
Deletions are marked like this. Additions are marked like this.
Line 24: Line 24:
|| 04.09.2007 || pos_chem.f90 || 3.1 || Neues Programm pos_chem eingecheckt, in dem die chem. Spezies in einer Init-Datei unterhalb einer Untergrenze und oberhalb einer Obergrenze, durch Werte aus einer Klimatologie-Datei ersetzt werden. Die Variablen aus der Init-Datei, die als Koordinaten genutzt werden sollen werden eingelesen (lat/eqlat und theta/zeta). Für jede Spezies wird der Name in der Init-Datei, der entspr. Name in der Clim-Datei, eine Unter- und eine Obergrenze und der Name der Clim-Datei eingelesen. Die Interpolation von den Klimatologie-Daten, z.B. O3(eqlat,press,month) auf ein Clams-Feld O3(nparts) erfolgt mit Hilfe der Routine Die Positionen eines pos-Files und die eines Flugzeug-Datensatzes werden zusammen in ein neues pos-File geschrieben. Außerdem wird eine Tracer-Variable (TRACER_NEW) ergänzt, die nur auf dem Flugpfad belegt ist.!
interpol_2d_values (utils-Lib, read2d_utils.f90). ||
|| 04.09.2007 || pos_chem.f90 || 3.1 || Neues Programm pos_chem eingecheckt, in dem die chem. Spezies in einer Init-Datei unterhalb einer Untergrenze und oberhalb einer Obergrenze, durch Werte aus einer Klimatologie-Datei ersetzt werden. Die Variablen aus der Init-Datei, die als Koordinaten genutzt werden sollen werden eingelesen (lat/eqlat und theta/zeta). Für jede Spezies wird der Name in der Init-Datei, der entspr. Name in der Clim-Datei, eine Unter- und eine Obergrenze und der Name der Clim-Datei eingelesen. Die Interpolation von den Klimatologie-Daten, z.B. O3(eqlat,press,month) auf ein Clams-Feld O3(nparts) erfolgt mit Hilfe der Routine interpol_2d_values (utils-Lib, read2d_utils.f90). ||
Line 54: Line 53:
|| 28.10.2008 || pos_dyn.f90 || 3.17 || - Aufruf von interpolate_data_met_press durch Aufruf von interpolate_data_met (mit presslev=.true.) ersetzt <<BR>> - use_log jetzt integer:<<BR>> 0= lineare Interpolation (bisher use_log=.false.) <<BR>> 1=log.-lin. Interpolation (bisher use_log=.true,.) <<BR>> 2=doppelt log.-lin. Interpolation (wird im Moment nur für SH genutzt) ||
|| 28.10.2008 || lib_pos.f90 || 3.9 || - sub. interpolate_data_met_press entfernt <<BR>> - sub. interpolate_data_met: Interpolation auf Druckniveaus integriert , uselog jetz integer ||
|| 17.11.2008 || pos_dyn.f90 || 3.18 || ||
|| 17.11.2008 || pos_dyn.f90 || 3.19 || ||
|| 17.11.2008 || pos_dyn.f90 || 3.20 || ||
|| 06.11.2008 || pos_dyn.f90 || 3.21 || ||

Änderungen im Paket pos

13.07.2004

lib_pos.f90

3.6

in function "interpolate_data_met" kann jetzt optional auch log.-linear interpoliert werden.

13.07.2004

pos_dyn.f90

3.5

SH wird log.-linear interpoliert

10.06.2005

pos_3d.f90

3.1

Beliebige Vertikalkoordinate

10.06.2005

pos_area.f90

3.8

Beliebige Vertikalkoordinate

10.06.2005

pos_dyn.f90

3.9

Beliebige Vertikalkoordinate

10.06.2005

pos_voronoi.f90

3.1

Beliebige Vertikalkoordinate

10.06.2005

lib_pos.f90

3.7

Beliebige Vertikalkoordinate

10.06.2005

pos_3d.inp_tmpl

3.1

Angabe der Vertikalkoordianten

10.06.2005

pos_area.inp_tmpl

3.3

Angabe der Vertikalkoordianten

29.06.2005

pos_area.f90

3.9

globales Attr. exp_VERTCOOR_ref_level in Ausgabedatei kopieren

26.09.2005

pos_area.f90

3.10

Umstellung auf Parallelverarbeitung (mit nur einem Prozessor)

26.09.2005

pos_dyn.f90

3.10

Umstellung auf Parallelverarbeitung (mit nur einem Prozessor)

19.10.2005

pos_area.f90

3.11

Entropie-Koordinaten für Zeta möglich:
neue Subroutines:
- make_theta_grid2
- interpol (lineare Interpolation eines Parameters)
- int_tabulated (Integration mit Newton-Codes-Formel, s. Num.Rec. Buch S. 126 "Bode's rule")
- calc_zeta (Berechne Zeta aus TEMP und PRESS)
- calc_ind
- find_nb
- find_next_index

19.10.2005

pos_area.inp_tmpl

3.4

zusätzliche Angaben für neue Version

19.10.2005

Makefile

3.9

numerical recipes linken ("spline" und "splint" werden in int_tabulated genutzt)

17.11.2005

pos_area.f90

3.12

Verbesserung der Übersichtlichkeit:
- globale Var. i nach itheta umbenannt
- in Sub. i lokal deklariert
- Übergabeparameter und lokale Variablen in Subroutines getrennt etc.

20.04.2007

pos_dyn.f90

3.11

- Code in Unterprogramme ausgelagert
- Einlesen einer zweiten Parameterliste mit 2 Levelangaben: Die Parameter müssen bereits in der Datei vorhanden sein. Die Werte werden nur unterhalb des ersten und oberhalb des zweiten Levels aus ECMWF-Daten interpoliert, dazwischen bleiben sie unverändert.
- GPL-Header ergänzt

14.06.2007

pos_dyn.inp_tmpl

3.3

Angabe einer zweiten Parameterliste (mit Unter- und Obergrenzen für das Ersetzen aus Isentropic-Dateien)

21.06.2007

pos_3d.f90

3.2

30-Tage-Monate möglich (wenn "dates30" im HP auf true gesetzt wird)

04.09.2007

pos_dyn.f90

3.12

- 1. Parameterliste: Wenn 03 oder SH hinzugefügt oder überschrieben werden sollen, werden sie immer unter dem Namen O3_ECMWF bzw. SH_ECMWF hinzugefügt (bzw. überschrieben)
- 2. Parameterliste: Parameter werden unterhalb des ersten und oberhalb des zweiten angegeben Levels überschrieben (SH, O3, H20). Bei SH kann wie bisher mit SH aus ECMWF überschrieben werden. Bei O3 und H2O wird O3 bzw. SH aus ECMWF genommen und zunächst umgerechnet (von Massenmischungsverhältnis nach Volumenmischungsverhältnis):
SH = SH_ECMWF
H2O = (28.9644/18.015) * SH_ECMWF
O3 = (28.9644/48.0) * O3_ECMWF

04.09.2007

pos_chem.f90

3.1

Neues Programm pos_chem eingecheckt, in dem die chem. Spezies in einer Init-Datei unterhalb einer Untergrenze und oberhalb einer Obergrenze, durch Werte aus einer Klimatologie-Datei ersetzt werden. Die Variablen aus der Init-Datei, die als Koordinaten genutzt werden sollen werden eingelesen (lat/eqlat und theta/zeta). Für jede Spezies wird der Name in der Init-Datei, der entspr. Name in der Clim-Datei, eine Unter- und eine Obergrenze und der Name der Clim-Datei eingelesen. Die Interpolation von den Klimatologie-Daten, z.B. O3(eqlat,press,month) auf ein Clams-Feld O3(nparts) erfolgt mit Hilfe der Routine interpol_2d_values (utils-Lib, read2d_utils.f90).

04.09.2007

pos_chem.inp_tmpl

3.1

Template für pos_chem.inp

04.09.2007

Makefile

3.10

Programm pos_chem ergänzt

07.09.2007

Makefile

3.11

pos_chem zu prog_posTools hinzugefügt, damit es bei "gmake all" compiliert wird

07.09.2007

Makefile

3.12

Compile-Befehl für pos_chem geändert (dynmod wird nicht benötigt)

07.09.2007

pos_chem.f90

3.2

Unterhalb der angegebenen Untergrenze treten falsche Werte auf:
=> missing values auf THETA abfangen
=> Fehler in utils/src/read2d_utils.f90 korrigiert

10.09.2007

pos_dyn.f90

3.13

- Attribute für geänderte Parameter (2. Parameterliste) nicht überschreiben, sonst sind evtl. die Einheiten falsch (für H2O, O3 wird umgerechnet!)
- Attribute für neu hinzugefügtes O3_ECMWF oder SH_ECMWF von O3 bzw. SH auf Met-Datei übernehmen

11.09.2007

pos_dyn.inp_tmpl

3.4

Kommentare ergänzt

18.09.2007

pos_chem.f90

3.3

Benutzung von 30-Tage-Monaten möglich (sub. read_init_file: ymds2js_interface)

18.09.2007

Makefile

3.13

- Compile- und Link-Befehl getrennt
- nur die benötigten Files linken, Compile-Befehle angepasst
- GPL-Header ergänzt

24.09.2007

pos_chem.f90

3.4

Warnung bei 30-Tage-Monaten ausgeben

04.10.2007

pos.tex

3.2

Beschreibungen für pos_dyn, pos_voronoi und pos_chem ergänzt

27.11.2007

Makefile

3.14

F90-Optionen für Link-Step ergänzt (sonst läßt es sich auf JUMP nicht compilieren)

27.11.2007

Makefile

3.15

Compiler-Optionen für Jump (-O3) angegeben

03.12.2007

pos_dyn.f90

3.14

Parallelverarbeitung auf dem JUMP ist jetzt möglich (bisher wurde immer nur eine CPU genutzt)

10.12.2007

pos_chem.f90

3.5

Benutzung von MPI möglich (läuft aber nur auf einer CPU)

28.01.2008

lib_pos.f90

3.8

sub. interpolate_data_met:
- gotos entfernt
- neue sub. interpolate_data_met_press (wenn Druck die Vertikalkoordinate ist)

28.01.2008

pos_dyn.f90

3.15

Vertikalkoordinate "press" möglich: sub. interpolate_parameters und sub. change_parameters angepaßt (und übersichtlicher gemacht)

28.01.2008

pos_dyn.f90

3.16

Testausgaben entfernt

22.02.2008

Makefile

3.16

lokale Testversion (unter $HOME/clams/utils) von utils-Library kann genutzt werden

22.04.2008

add_flightpath.f90

3.1

neues Programm: Die Positionen eines pos-Files und die eines Flugzeug-Datensatzes werden zusammen in ein neues pos-File geschrieben. Außerdem wird eine Tracer-Variable (TRACER_NEW) ergänzt, die nur auf dem Flugpfad belegt ist.

22.04.2008

add_flightpath.inp_tmpl

3.1

Beispiel-Konfigurationsdatei für add_flightpath

22.04.2008

Makefile

3.17

Programm add_flightpath ergänzt

28.04.2008

add_flightpath.f90

3.2

Es können durch mehrfaches Aufrufen dieses Programms mehrere Flugpfade ergänzt werden:
-> schreibe auf TRACER_nn
-> Felder für TRACER_01-TRACER_n-1 müssen entsprechend (mit 0) aufgefüllt werden

26.05.2008

pos_flightpath

3.1

neues Programm: Die Positionen (lat,lon,time) eines Flugpfades werden eingelesen und für vorgegebene Theta-Level (pos_flightpath.inp) in ein pos-File ausgegeben.

26.05.2008

pos_flightpath.inp_tmpl

3.1

Beispiel-Konfigurationsdatei für pos_flightpath

26.05.2008

Makefile

3.18

Programm pos_flightpath ergänzt

15.07.2008

pos_flightpath.f90

3.2

Header ergänzt

15.07.2008

add_flightpath.f90

3.3

Header ergänzt

28.10.2008

pos_dyn.f90

3.17

- Aufruf von interpolate_data_met_press durch Aufruf von interpolate_data_met (mit presslev=.true.) ersetzt
- use_log jetzt integer:
0= lineare Interpolation (bisher use_log=.false.)
1=log.-lin. Interpolation (bisher use_log=.true,.)
2=doppelt log.-lin. Interpolation (wird im Moment nur für SH genutzt)

28.10.2008

lib_pos.f90

3.9

- sub. interpolate_data_met_press entfernt
- sub. interpolate_data_met: Interpolation auf Druckniveaus integriert , uselog jetz integer

17.11.2008

pos_dyn.f90

3.18

17.11.2008

pos_dyn.f90

3.19

17.11.2008

pos_dyn.f90

3.20

06.11.2008

pos_dyn.f90

3.21

pos/Changes (last edited 2013-07-15 07:25:42 by NicoleThomas)