Datum

Datei

Version

Änderung

11.01.2006

bmix.f90

5.0

- neue sub. get_lower_bound_species (Interpolieren von species aus init-File für unterste Schicht)
- neue sub. add_to_species (Werte aus EDGAR Emissionsdateien zu species hinzufügen)
- Hauptprogramm verkürzt: Code in Subroutines ausgelagert

11.01.2006

bmix.inp_tmpl

5.0

zusätzliche Angaben für EDGAR Dateien

11.01.2006

bmix.f90

5.1

Syntaxfehler (auf JUMP) korrigiert

16.01.2006

lib_mix.f90

4.3

- neue Funktion "theta2thetadelta"
- nicht mehr genutzte sub. "vert_interpolation" gelöscht

16.01.2006

bmix.f90

5.2

neue Routinen "zeta2alt" und "alt2zeta"; in sub. "add_to_species" wird "theta2thetadelta" und "zeta2alt" für die Berechnung eines Faktors genutzt

16.01.2006

Makefile

3.19

qhull-Library muss nach utils-Library im Link-Pfad stehen

17.01.2006

bmix.f90

5.3

sub. add_to_species: Felder am Ende der Schleife freigeben (mit pgf90 nicht aufgefallen, aber Programm-Absturz auf dem JUMP!)

07.02.2006

bmix.f90

5.4

Das Bound-File für die oberste Schicht soll entweder 2-dimensionale Variablen (lat,time) oder 3-dimensionale Variablen (lat,lon,time) enthalten können. In sub. get_bound_species wird überprüft, ob "lon" als dritte Dimension vorhanden ist. Ist das der Fall, wird statt der Routine "interpol_time_grid" die neue Subroutine "interpol_time_grid_3d" für eine 3-dimensionale Interpolation aufgerufen.

20.02.2006

bmix.f90

5.5

sub. get_lower_bound_species: Es wird ein missing value gesetzt, wenn die Triangel nicht gefunden wurde.

20.03.2006

mix.f90

4.9

ohne Mischen: nthetas=0 an Ausgaberoutine übergeben!

28.03.2006

bmix.f90

5.6

in sub. replace_lower_layer: Sub. get_lower_bound_species wird nur bei Berücksichtung der Orographie (delta_theta>0) aufgerufen.

02.05.2006

bmix.f90

5.7

sub. replace_lower_layer: logische Var. "use_oro" zur korrekten Unterscheidung (mit/ohne Beruecksichtigung der Orographie); Bisher wurde auf delta_theta>0 abgefragt. Das war aber immer .true., da delta_theta bei nicht Berücksichtigung der Orographie auf theta_in_down-theta_down gesetzt wird!

29.05.2006

age.f90

5.1

neues Programm

29.05.2006

Makefile

3.20

Compilation von "age" ergänzt

21.06.2006

bmix.f90

5.8

sub. replace_lower_layer: Aufruf von get_orography vor der Verzweigung (damit theta0 allociert ist)

22.11.2006

bmix.f90

5.9

Änderungen aus Seitenzweig (Version 5.7.2.5) eingefügt

23.11.2006

bmix.f90

5.10

Aufruf der neuen Wichtungsroutine determ_weight_new (utils:lib_interpol.f90)

23.11.2006

age.f90

5.2

Änderungen aus Seitenzweig (Version 5.1.2.1) eingefügt

01.12.2006

lib_io.f90

4.8

Änderungen aus Seitenzweig (Vers. 4.7.2.4) eingefügt

01.12.2006

dynmod: types_m.f90

4.6

Änderungen aus Seitenzweig (Vers. 4.5.2.4) eingefügt

01.12.2006

dynmod: ap_m_access.f90

4.9

Änderungen aus Seitenzweig (Vers. 4.8.2.3) eingefügt

01.12.2006

mix.f90

4.10

save_tp wird nicht mehr genutzt (s. Seitenzweig Vers. 4.9.2.2), sonst Probleme mit geänderten Ein-und Ausgaberoutinen

08.12.2006

cmdl/make_bound_file.f90

1.2

Hiermit können die NOAA CMDL Files in ein NetCDF Boundary File konvertiert werden: Kopie aus Seitenzweig (1.1.2.2.)

08.12.2006

cmdl/add_to_bound_files.f90

1.2

Hiermit können die nach Netcdf konvertierten NOAA CMDL Files um einen angegebenen Zeitraum erweitert werden: Kopie aus Seitenzweig (1.1.2.3)

08.12.2006

cmdl/make_bound_file.inp_tmpl

1.2

Kopie aus Seitenzweig (1.1.2.1)

08.12.2006

cmdl/add_to_bound_files.inp_tmpl

1.2

Kopie aus Seitenzweig (1.1.2.1)

08.12.2006

cmdl/Makefile

1.2

Makefile aus Seitenzweig (1.1.2.2)

11.12.2006

cmdl/make_bound_file.inp_tmpl

1.3

weitere Angaben (user, output filename)

12.12.2006

mix.f90

4.11

Erweiterungen aus dem Seitenzweig (4.9.2.8) eingebaut:
- 2 Faktoren + 2 Level
- Abbruch der Schleife "implicit mixing" wenn Anz. sich nicht mehr als vorgegeben verändert
- temporäre Ausgabedatei löschen
- GPL Header

12.12.2006

lib_mix.f90

4.4

GPL Header

12.12.2006

lib_io.f90

4.9

GPL Header

12.12.2006

mix.inp_tmpl

3.1

- unnötige Angaben gelöscht (Bound-File, start/end-offset und save_temp)
- Zeile 6: Faktor1, Faktor2, Ober- und Untergrenze fuer ersten Faktor
- Zeile 8: Max. rel. Abweichung zwischen n und n_old fuer Abbruch der Schleife

25.05.2007

mix.f90

4.12

zwei Probleme auf dem JUMP behoben:
- in Zeile 439: i => itheta
- temporäre Datei nur von 1. Prozessor löschen lassen

30.05.2007

mix.f90

4.13

Fehler auf JUMP gehoben:
geshiftete Level (grid_switch=1) wurden bei mehreren CPUs nicht korrekt ermittelt (die Hälfte des untersten Intervalls muss wieder subtrahiert werden!)

05.06.2007

bmix.f90

5.11

sub. get_vert_species: Boundfile fuer Nord-/Süd-Rand kann 3-dimensional sein (lon,theta,time)

02.07.2007

mix.f90

4.14

Benutzung von 30-Tage-Monaten möglich

18.09.2007

mix.f90

4.15

Aufruf von nc_write_ap_s mit "dates30"

18.09.2007

bmix.f90

5.12

Benutzung von 30-Tage-Monaten möglich

21.09.2007

cmdl/add_to_bound_files.f90

1.3

Erweitern der Daten in die Vergangenheit:
Es wird zusätzlich ein Startjahr eingelesen und ab diesem bis zu den vorhandenen Daten die Werte des ersten vorhandenen Jahres wiederholt.

21.09.2007

cmdl/add_to_bound_files.inp_tmpl

1.3

zusätzliche Angabe des Startdatums

24.09.2007

mix.f90

4.16

Warnung bei 30-Tage-Monaten ausgeben

24.09.2007

bmix.f90

5.13

Warnung bei 30-Tage-Monaten ausgeben

02.10.2007

bmix.f90

5.13

Tag "bmix-mit-edgar-files" gesetzt

02.10.2007

bmix.f90

5.14

CO-Update aus Edgar-Datenbank wird nicht mehr ausgeführt

02.10.2007

bmix.inp_tmpl

5.1

Angaben zu Edgar-Files gelöscht;
Interpolation, wenn delta_theta>0 (Orographie), nur wenn "y" in letzter Zeile von bmix.inp angegeben ist

04.10.2007

bmix.f90

5.15

sub. alt2zeta und zeta2alt entfernt (wurde für Edgar-Dateien genutzt)

09.10.2007

Makefile

3.21

für "age" werden nicht alle Libraries gelinkt

23.10.2007

replace_lbound.f90

5.1

neues Programm: Unteren Rand bis zu einem angegebenen Level mit Daten aus der angegebenen Datei ersetzen => z. B. Mopitt-Daten für unteren Rand bis 250K ersetzen (mit replace_lbound), aber ansonsten unteren Rand bis 75K ersetzen (mit bmix) (Aufgebaut wie pos_chem mit der Interpolation am unteren Rand wie in bmix)

23.10.2007

replace_lbound.inp_tmpl

5.1

Template für replace_lbound.inp

23.10.2007

Makefile

3.22

Programm replace_lbound ergänzt

26.11.2007

lib_io.f90

4.10

sub. get_lowest_layer: Feldüberschreitung in theta0 abgefangen (Es wird nur noch theta0(ilat,ilon) genutzt, nicht mehr die umliegenden Feldelemente)

30.11.2007

Absturz von bmix auf dem JUMP: ap_s_1(i)%ind wurde nicht initialisiert (nullify)

30.11.2007

bmix.f90

5.16

ap_s_1 wird jetzt in sub. get_lowest_layer allokiert (in der richtigen Größe)

30.11.2007

lib_io.f90

4.11

Feld ap_s_1 wird in der richtigen Größe allokiert (und nullify für pointer "ind")

04.12.2007

Absturz von bmix auf 64-Bit-Maschinen

04.12.2007

bmix.f90

5.17

Größe von ap_s_1 wird durch Aufruf von sub. get_nparts_lowest_layer ermittelt

04.12.2007

lib_io.f90

4.12

neue sub. get_nparts_lowest_layer, keine Allokierungen mehr in sub. get_lowest_layer (jetzt im HP)

10.12.2007

lib_io.f90

4.13

Absturz, wenn PV_INIT nicht in Bound-File vorhanden ist:
sub. nc_read_ap_s: PV_init immer mit Größe "nparts" allokieren

14.12.2007

bmix.f90

5.18

Absturz, wenn für unterstes Level Interpolation vom Init-File eingeschaltet wird:
sub. get_lower_bound_species: Dim. der Felder ap0_coor und ap1_coor angepaßt

17.12.2007

cmdl/make_bound_file.f90

1.3

Erweiterungen für N2O- und F11-Daten:
- Ermittlung von ntimes geändert
- sub. read_species: umrechnen von ppm,ppb oder ppt

14.02.2008

cmdl/add_to_bound_files.f90

1.4

nyeartimes angeben (NICHT berechnen) => Zeitangaben werden sonst evtl. zu ungenau

14.02.2008

cmdl/make_bound_file.f90

1.4

Eingabedateien können mit Verzeichnis angegeben werden

15.02.2008

bmix.f90

5.19

Verwendung eines HCN Boundary-Files möglich

22.02.2008

Makefile

3.23

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

03.03.2008

cmdl/data/*

Dateien unter cmdl/data eingecheckt

04.03.2008

cmdl/input_data/*

Dateien unter cmdl/input_data eingecheckt

04.03.2008

cmdl/cats_f11.f90

1.1

F11-Dateien (unter cmdl/input_data) werden in das Format der CMDL-Files konvertiert (nach cmdl/data)

04.03.2008

cmdl/cats_n2o.f90

1.1

N2O-Dateien (unter cmdl/input_data) werden in das Format der CMDL-Files konvertiert (nach cmdl/data)

04.03.2008

cmdl/ccval.f90

1.1

ccmval-File (unter cmdl/input_data) wird in das Format der CMDL-Files konvertiert (nach cmdl/data)

07.03.2008

cmdl/readme.txt

1.1

Beschreibung der Boundfiles im Verzeichnis "data" und der Erstellung von Boundfiles für bmix aus diesen Dateien

10.03.2008

bmix.f90

5.20

Die Namen der zusätzlichen Boundfiles und die Spezies, die von diesen am oberen oder unteren Rand ersetzt werden sollen, werden von bmix.inp eingelesen.

10.03.2008

bmix.inp_tmpl

5.2

Namen der Boundfiles können optional angegeben werden

14.03.2008

bmix.f90

5.21, 5.22

Kommentare in bmix.inp in Zeilen mit Boundfiles möglich

31.03.2008

lib_io.f90

4.14

Funktion get_nparts_lowest_layer und sub. get_lowest_layer: modulo(lon,360.) verwenden, damit immer lon<360. (sonst evtl. Absturz in bmix)

07.04.2008

lib_io.f90

4.15

sub. nc_read_ap_s_0_info: spec_arr mit Leerstring vorbelegen

07.04.2008

bmix.f90

5.23

In den eingelesenen Bound-files wurde irrtümlich die einglesenen Werte aller Spezies am unteren Rand auf den unteren Bereich der ersten Spezies in spec_arr geschrieben
=> in sub. get_bound_species: c_spec_arr nutzen (nicht globales Feld (spec_arr) überschreiben und an sub.interpol_time_grid[_3d] weitergeben
=> sub. interpol_time_grid u. interpol_time_grid_3d: Position der aktuellen Spezies (in c_spec_arr) im globalen Feld "spec_arr" suchen

07.04.2008

lib_io.f90

4.16

Änderung aus 4.15 wieder rückgängig gemacht (sonst enthält Species im Ausgabedatensatz Kontrollzeichen ?!)

07.04.2008

bmix.f90

5.24

sub.interpol_time_grid: nur mit den ersten 10 Zeichen von spec_arr vergleichen. Da spec_arr nicht mehr vorbelegt wird (s. lib_io.f90), kann es nach der 10. Stelle Kontrollzeichen enthalten.

08.04.2008

bmix.f90

5.25

sub.interpol_time_grid_3d: wie sub. interpol_time_grid (s.07.04.) geändert

10.04.2008

bmix.f90

5.26

Testausgabe entfernt

06.05.2008

mix.f90

4.17

Aufruf der Funktion nc_set_new_history zum Erzeugen einer neuen History, die nur aus dem aktuellen Datum und dem Namen des Programmes (mix) besteht (Bisher wurde die History vom Init-File des Vortags kopiert und somit mit jedem Tag eines Clams-Laufes länger !)

19.05.2008

*.f90

Tag "bmix-ohne-add" gesetzt (Version vor dem Branch "mix-mit-add" mit einer Testversion für Robert)

17.06.2008

cmdl/Makefile

1.3

an Clams-Makefile-Struktur angepaßt (nutze mkincl)

26.06.2008

age.inp_tmpl

5.1

Template für age.inp (bisher nur im Branch new-mix-test eingecheckt)

13.06.2006

mix.f90

4.9.2.1

Schichtweise Bearbeitung geändert

13.06.2006

lib_mix.f90

4.3.2.1

Schichtweise Bearbeitung geändert

13.06.2006

lib_io.f90

4.7.2.1

Schichtweise Bearbeitung geändert

20.06.2006

mix.inp_tmpl

3.0.2.1

2 Faktoren + Level angeben; start_offset, end_offset, save_tp nicht mehr angeben

20.06.2006

mix.f90

4.9.2.2

- Statt einem Faktor für die maximale horizontale Distanz werden jetzt zwei Faktoren und ein Level eingelesen. Unterhalb des angegebenen Levels wird Faktor1 genutzt und oberhalb Faktor2.
- Start- und Endoffset werden nicht mehr eingelesen
- Temperatur und Druck werden nicht mehr mit ausgegeben

20.06.2006

bmix.f90

5.7.2.1

Temperatur und Druck werden nicht mehr mit ausgegeben

20.06.2006

lib_io.f90

4.7.2.2

Temperatur und Druck werden nicht mehr mit ausgegeben

21.06.2006

lib_io.f90

4.7.2.3

sub. nc_create_outfile: logische Variable pv_defined entfernt (da sie immer den Wert .true. hatte)

04.07.2006

mix.f90

4.9.2.3

Anpassungen für 2d-isentropic-Fall

28.08.2006

Erzeugung doppelter Punkte vermeiden

28.08.2006

mix.f90

4.9.2.4

l_min_act an sub. adapt_grid übergeben

28.08.2006

lib_mix.f90

4.3.2.2

"eliminate" bzw. "new_point" nur ausführen, wenn beide Punkte größer sind als "l_min_act" (untere Grenze des kleinen Intervalls)

28.08.2006

dynmod/source:
ap_m_modify.f90

4.0.2.3

Funktion "theta_gt_min" ergänzt, die überprüft, ob die beiden angegebenen Punkte oberhalb eine angegebenen Theta-Levels liegen

08.09.2006

bmix.f90

5.7.2.2

NOAA CMDL boundary files können genutzt werden

08.09.2006

bmix.f90

5.7.2.3

temporäre Ausgabedatei wieder löschen

08.09.2006

mix.f90

4.9.2.5

temporäre Ausgabedatei wieder löschen

08.09.2006

cmdl/make_bound_file.f90

1.1.2.1

Hiermit können die NOAA CMDL Files in ein NetCDF Boundary File konvertiert werden.

08.09.2006

cmdl/Makefile

1.1.2.1

Makefile für make_bound_file

13.09.2006

Stetig ansteigenden Speicherplatz vermeiden

13.09.2006

mix.f90

4.9.2.6

für ap_s_0: Komponenten "ind" und "c" wieder deallokiert; GPL-Header ergänzt

13.09.2006

lib_io.f90

4.7.2.4

einige Hilfsfelder wieder deallokiert

13.09.2006

dynmod/source:
ap_m_access.f90

4.8.2.3

sub. free: Listen mit Indizes der Nachbarn wieder deallokiert

13.09.2006

dynmod/source:
lib_triang.f90

4.1.2.1

sub. qhull: Feld "points" wieder deallokieren

26.09.2006

lib_mix.f90

4.3.2.3

switch an sub. eliminate übergeben

29.09.2006

bmix.f90

5.7.2.4

sub. replace_lower_layer: übersichtlicher gestaltet

17.10.2006

age.f90

5.1.2.1

Wenn EQLAT vorhanden ist, wird LAT mit EQLAT überschrieben

17.10.2006

age.inp_tmpl

1.1.2.1

template für age.inp angelegt

17.10.2006

age.inp_tmpl

1.1.2.2

'offset in days' ergänzt

02.11.2006

mix.f90

4.9.2.7

für max. Distanz: zweites Level einlesen (fac_limit_level2)
Unterhalb von fac_limit_level und oberhalb von fac_limit_level2 wird fac_limit genutzt, zwischen den beiden Leveln fac_limit2

14.11.2006

mix.f90

4.9.2.8

- Einlesen einer max. rel. Abweichung auf adapt_par%r_dev
- implicit mixing: Abbruch der Schleife, wenn rel. Abweichung zwischen n und n_old kleiner als r_dev

14.11.2006

mix.f90

4.9.2.9

in implicit mixing: Vergleich mit <= statt mit <

14.11.2006

mix.inp_tmpl

3.0.2.3

Zeile 6: Ober- und Untergrenze für ersten Faktor;
Zeile 8: Max. rel. Abweichung zwischen n und n_old für Abbruch der Schleife

22.11.2006

bmix.f90

5.7.2.5

MOPITT boundary files können genutzt werden

23.11.2006

bmix.f90

5.7.2.6

Aufruf der neuen Wichtungsroutine determ_weight_new (utils:lib_interpol.f90)

07.12.2006

cmdl/add_to_bound_files.f90

1.1.2.1

Hiermit können die nach Netcdf konvertierten NOAA CMDL Files um einen angegebenen Zeitraum erweitert werden

07.12.2006

cmdl/add_to_bound_files.f90

1.1.2.2

Einlesen von add_to_bound_files.inp

07.12.2006

cmdl/add_to_bound_files.f90

1.1.2.3

Doppelt genaue Matrix benutzen

07.12.2006

cmdl/Makefile

1.1.2.2

Prg. add_to_bound_files ergänzt

07.12.2006

cmdl/make_bound_file.inp_tmpl

1.1.2.1

template für make_bound_file.inp

08.12.2006

cmdl/add_to_bound_files.inp_tmpl

1.1.2.1

template für add_to_bound_files.inp

08.12.2006

cmdl/make_bound_file.f90

1.1.2.2

- Einlesen von make_bound_file.inp;
- times doppelt genau;
- sub. read_species: Suche nach der aktuellen Zeit geändert

11.12.2006

cmdl/make_bound_file.inp_tmpl

1.1.2.2

weitere Angaben (user, output filename)

15.04.2008

Änderungen/Erweiterungen aus dem Hauptzweig in den Nebenzweig (Tag new-mix-test) von mix und dynmod eingebaut.

15.04.2008

mix.f90

4.9.2.10

Änderungen aus Hauptzweig bis Vers. 4.16 eingebaut

15.04.2008

bmix.f90

5.7.2.7

Änderungen aus Hauptzweig bis Vers. 5.26 eingebaut

15.04.2008

lib_io.f90

4.7.2.5

Änderungen aus Hauptzweig bis Vers. 4.16 eingebaut

15.04.2008

dynmod/types_m.f90

4.5.2.5

Änderungen aus Hauptzweig bis Vers. 4.7 eingebaut

15.04.2008

dynmod/ap_m_access.f90

4.8.2.4

Änderungen aus Hauptzweig bis Vers. 4.11 eingebaut

15.04.2008

dynmod/ap_m_modifiy.f90

4.0.2.7

Änderungen aus Hauptzweig bis Vers. 4.2 eingebaut

15.04.2008

dynmod/lib_triang.f90

4.1.2.2

Änderugnen aus Hauptzweig bis Vers. 4.2 eingebaut

06.05.2008

mix.f90

4.9.2.10

Aufruf der Funktion nc_set_new_history zum Erzeugen einer neuen History, die nur aus dem aktuellen Datum und dem Namen des Programmes (mix) besteht (Bisher wurde die History vom Init-File des Vortags kopiert und somit mit jedem Tag eines Clams-Laufes länger !)

19.05.2008

bmix.f90

5.26.2.1

in sub. interpol_time_grid[3d] können die Werte ersetzt oder aufaddiert werden -> welche Operation durchzuführen ist, wird aus bmix.inp eingelesen