Home

Context Navigation

source: palm/trunk/SOURCE/nudging_mod.f90 @ 2232

Last change on this file since 2232 was 2232, checked in by suehring, 7 years ago
Adjustments according new topography and surface-modelling concept implemented
Property svn:keywords set to `Id`
File size: 21.4 KB

Line
1	!> @file nudging_mod.f90
2	!------------------------------------------------------------------------------!
3	! This file is part of PALM.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the
6	! terms of the GNU General Public License as published by the Free Software
7	! Foundation, either version 3 of the License, or (at your option) any later
8	! version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2017 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! ------------------
22	! Adjustments to new topography concept
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: nudging_mod.f90 2232 2017-05-30 17:47:52Z suehring $
27	!
28	! 2000 2016-08-20 18:09:15Z knoop
29	! Forced header and separation lines into 80 columns
30	!
31	! 1850 2016-04-08 13:29:27Z maronga
32	! Module renamed
33	!
34	!
35	! 1757 2016-02-22 15:49:32Z maronga
36	! Bugfix: allow for using higher vertical resolution in nudging file than grid
37	! spacing in the LES model
38	!
39	! 1682 2015-10-07 23:56:08Z knoop
40	! Code annotations made doxygen readable
41	!
42	! 1398 2014-05-07 11:15:00Z heinze
43	! Subroutine nudge_ref is extended to set u_init and v_init to the current
44	! nudging profiles
45	!
46	! 1382 2014-04-30 12:15:41Z boeske
47	! Changed the weighting factor that is used in the summation of nudging
48	! tendencies for profile data output from weight_pres to weight_substep,
49	! added Neumann boundary conditions for profile data output of nudging terms at
50	! nzt+1
51	!
52	! 1380 2014-04-28 12:40:45Z heinze
53	! Subroutine nudge_ref added to account for proper upper scalar boundary
54	! conditions in case of nudging
55	!
56	! 1365 2014-04-22 15:03:56Z boeske
57	! Variable t renamed nt, variable currtnudge renamed tmp_tnudge,
58	! summation of nudging tendencies for data output added
59	! +sums_ls_l, tmp_tend
60	! Added new subroutine calc_tnudge, which calculates the current nudging time
61	! scale at each time step
62	!
63	! 1355 2014-04-10 10:21:29Z heinze
64	! Error message specified.
65	!
66	! 1353 2014-04-08 15:21:23Z heinze
67	! REAL constants provided with KIND-attribute
68	!
69	! 1320 2014-03-20 08:40:49Z raasch
70	! ONLY-attribute added to USE-statements,
71	! kind-parameters added to all INTEGER and REAL declaration statements,
72	! kinds are defined in new module kinds,
73	! old module precision_kind is removed,
74	! revision history before 2012 removed,
75	! comment fields (!:) to be used for variable explanations added to
76	! all variable declaration statements
77	!
78	! 1318 2014-03-17 13:35:16Z raasch
79	! module interfaces removed
80	!
81	! 1268 2013-12-12 09:47:53Z heinze
82	! bugfix: argument of calc_mean_profile corrected
83	!
84	! 1251 2013-11-07 08:14:30Z heinze
85	! bugfix: calculate dtm and dtp also in vector version
86	!
87	! 1249 2013-11-06 10:45:47Z heinze
88	! remove call of user module
89	! reformatting
90	!
91	! 1241 2013-10-30 11:36:58Z heinze
92	! Initial revision
93	!
94	! Description:
95	! ------------
96	!> Nudges u, v, pt and q to given profiles on a relaxation timescale tnudge.
97	!> Profiles are read in from NUDGIN_DATA. Code is based on Neggers et al. (2012)
98	!> and also part of DALES and UCLA-LES.
99	!--------------------------------------------------------------------------------!
100	MODULE nudge_mod
101
102
103	PRIVATE
104	PUBLIC init_nudge, calc_tnudge, nudge, nudge_ref
105	SAVE
106
107	INTERFACE nudge
108	MODULE PROCEDURE nudge
109	MODULE PROCEDURE nudge_ij
110	END INTERFACE nudge
111
112	CONTAINS
113
114	!------------------------------------------------------------------------------!
115	! Description:
116	! ------------
117	!> @todo Missing subroutine description.
118	!------------------------------------------------------------------------------!
119	SUBROUTINE init_nudge
120
121	USE arrays_3d, &
122	ONLY: ptnudge, qnudge, timenudge, tmp_tnudge, tnudge, unudge, &
123	vnudge, wnudge, zu
124
125	USE control_parameters, &
126	ONLY: dt_3d, lptnudge, lqnudge, lunudge, lvnudge, lwnudge, &
127	message_string, ntnudge
128
129	USE indices, &
130	ONLY: nzb, nzt
131
132	USE kinds
133
134	IMPLICIT NONE
135
136
137	INTEGER(iwp) :: finput = 90 !<
138	INTEGER(iwp) :: ierrn !<
139	INTEGER(iwp) :: k !<
140	INTEGER(iwp) :: nt !<
141
142	CHARACTER(1) :: hash !<
143
144	REAL(wp) :: highheight !<
145	REAL(wp) :: highqnudge !<
146	REAL(wp) :: highptnudge !<
147	REAL(wp) :: highunudge !<
148	REAL(wp) :: highvnudge !<
149	REAL(wp) :: highwnudge !<
150	REAL(wp) :: hightnudge !<
151
152	REAL(wp) :: lowheight !<
153	REAL(wp) :: lowqnudge !<
154	REAL(wp) :: lowptnudge !<
155	REAL(wp) :: lowunudge !<
156	REAL(wp) :: lowvnudge !<
157	REAL(wp) :: lowwnudge !<
158	REAL(wp) :: lowtnudge !<
159
160	REAL(wp) :: fac !<
161
162	ALLOCATE( ptnudge(nzb:nzt+1,1:ntnudge), qnudge(nzb:nzt+1,1:ntnudge), &
163	tnudge(nzb:nzt+1,1:ntnudge), unudge(nzb:nzt+1,1:ntnudge), &
164	vnudge(nzb:nzt+1,1:ntnudge), wnudge(nzb:nzt+1,1:ntnudge) )
165
166	ALLOCATE( tmp_tnudge(nzb:nzt) )
167
168	ALLOCATE( timenudge(0:ntnudge) )
169
170	ptnudge = 0.0_wp; qnudge = 0.0_wp; tnudge = 0.0_wp; unudge = 0.0_wp
171	vnudge = 0.0_wp; wnudge = 0.0_wp; timenudge = 0.0_wp
172	!
173	!-- Initialize array tmp_nudge with a current nudging time scale of 6 hours
174	tmp_tnudge = 21600.0_wp
175
176	nt = 0
177	OPEN ( finput, FILE='NUDGING_DATA', STATUS='OLD', &
178	FORM='FORMATTED', IOSTAT=ierrn )
179
180	IF ( ierrn /= 0 ) THEN
181	message_string = 'file NUDGING_DATA does not exist'
182	CALL message( 'nudging', 'PA0365', 1, 2, 0, 6, 0 )
183	ENDIF
184
185	ierrn = 0
186
187	rloop:DO
188	nt = nt + 1
189	hash = "#"
190	ierrn = 1 ! not zero
191	!
192	!-- Search for the next line consisting of "# time",
193	!-- from there onwards the profiles will be read
194	DO WHILE ( .NOT. ( hash == "#" .AND. ierrn == 0 ) )
195
196	READ ( finput, *, IOSTAT=ierrn ) hash, timenudge(nt)
197	IF ( ierrn < 0 ) EXIT rloop
198
199	ENDDO
200
201	ierrn = 0
202	READ ( finput, *, IOSTAT=ierrn ) lowheight, lowtnudge, lowunudge, &
203	lowvnudge, lowwnudge , lowptnudge, &
204	lowqnudge
205
206	IF ( ierrn /= 0 ) THEN
207	message_string = 'errors in file NUDGING_DATA'
208	CALL message( 'nudging', 'PA0366', 1, 2, 0, 6, 0 )
209	ENDIF
210
211	ierrn = 0
212	READ ( finput, *, IOSTAT=ierrn ) highheight, hightnudge, highunudge, &
213	highvnudge, highwnudge , highptnudge, &
214	highqnudge
215
216	IF ( ierrn /= 0 ) THEN
217	message_string = 'errors in file NUDGING_DATA'
218	CALL message( 'nudging', 'PA0366', 1, 2, 0, 6, 0 )
219	ENDIF
220
221	DO k = nzb, nzt+1
222	DO WHILE ( highheight < zu(k) )
223	lowheight = highheight
224	lowtnudge = hightnudge
225	lowunudge = highunudge
226	lowvnudge = highvnudge
227	lowwnudge = highwnudge
228	lowptnudge = highptnudge
229	lowqnudge = highqnudge
230
231	ierrn = 0
232	READ ( finput, *, IOSTAT=ierrn ) highheight , hightnudge , &
233	highunudge , highvnudge , &
234	highwnudge , highptnudge, &
235	highqnudge
236	IF (ierrn /= 0 ) THEN
237	WRITE( message_string, * ) 'zu(nzt+1) = ', zu(nzt+1), 'm is ',&
238	'higher than the maximum height in NUDING_DATA which ', &
239	'is ', lowheight, 'm. Interpolation on PALM ', &
240	'grid is not possible.'
241	CALL message( 'nudging', 'PA0364', 1, 2, 0, 6, 0 )
242	ENDIF
243	ENDDO
244
245	!
246	!-- Interpolation of prescribed profiles in space
247
248	fac = ( highheight - zu(k) ) / ( highheight - lowheight )
249
250	tnudge(k,nt) = fac * lowtnudge + ( 1.0_wp - fac ) * hightnudge
251	unudge(k,nt) = fac * lowunudge + ( 1.0_wp - fac ) * highunudge
252	vnudge(k,nt) = fac * lowvnudge + ( 1.0_wp - fac ) * highvnudge
253	wnudge(k,nt) = fac * lowwnudge + ( 1.0_wp - fac ) * highwnudge
254	ptnudge(k,nt) = fac * lowptnudge + ( 1.0_wp - fac ) * highptnudge
255	qnudge(k,nt) = fac * lowqnudge + ( 1.0_wp - fac ) * highqnudge
256	ENDDO
257
258	ENDDO rloop
259
260	CLOSE ( finput )
261
262	!
263	!-- Prevent nudging if nudging profiles exhibt too small values
264	!-- not used so far
265	lptnudge = ANY( ABS( ptnudge ) > 1.0e-8_wp )
266	lqnudge = ANY( ABS( qnudge ) > 1.0e-8_wp )
267	lunudge = ANY( ABS( unudge ) > 1.0e-8_wp )
268	lvnudge = ANY( ABS( vnudge ) > 1.0e-8_wp )
269	lwnudge = ANY( ABS( wnudge ) > 1.0e-8_wp )
270
271	END SUBROUTINE init_nudge
272
273
274	!------------------------------------------------------------------------------!
275	! Description:
276	! ------------
277	!> @todo Missing subroutine description.
278	!------------------------------------------------------------------------------!
279	SUBROUTINE calc_tnudge ( time )
280
281	USE arrays_3d, &
282	ONLY: timenudge, tmp_tnudge, tnudge
283
284	USE control_parameters, &
285	ONLY: dt_3d
286
287	USE indices, &
288	ONLY: nzb, nzt
289
290	USE kinds
291
292	IMPLICIT NONE
293
294
295	REAL(wp) :: dtm !<
296	REAL(wp) :: dtp !<
297	REAL(wp) :: time !<
298
299	INTEGER(iwp) :: k !<
300	INTEGER(iwp) :: nt !<
301
302	nt = 1
303	DO WHILE ( time > timenudge(nt) )
304	nt = nt+1
305	ENDDO
306	IF ( time /= timenudge(1) ) THEN
307	nt = nt-1
308	ENDIF
309
310	dtm = ( time - timenudge(nt) ) / ( timenudge(nt+1) - timenudge(nt) )
311	dtp = ( timenudge(nt+1) - time ) / ( timenudge(nt+1) - timenudge(nt) )
312
313	DO k = nzb, nzt
314	tmp_tnudge(k) = MAX( dt_3d, tnudge(k,nt) * dtp + tnudge(k,nt+1) * dtm )
315	ENDDO
316
317	END SUBROUTINE calc_tnudge
318
319	!------------------------------------------------------------------------------!
320	! Description:
321	! ------------
322	!> Call for all grid points
323	!------------------------------------------------------------------------------!
324	SUBROUTINE nudge ( time, prog_var )
325
326	USE arrays_3d, &
327	ONLY: pt, ptnudge, q, qnudge, tend, timenudge, tmp_tnudge, tnudge, &
328	u, unudge, v, vnudge
329
330	USE control_parameters, &
331	ONLY: dt_3d, intermediate_timestep_count, message_string
332
333	USE indices, &
334	ONLY: nxl, nxr, nys, nyn, nzb, nzt, wall_flags_0
335
336	USE kinds
337
338	USE statistics, &
339	ONLY: hom, sums_ls_l, weight_substep
340
341	IMPLICIT NONE
342
343	CHARACTER (LEN=*) :: prog_var !<
344
345	REAL(wp) :: tmp_tend !<
346	REAL(wp) :: dtm !<
347	REAL(wp) :: dtp !<
348	REAL(wp) :: time !<
349
350	INTEGER(iwp) :: i !<
351	INTEGER(iwp) :: j !<
352	INTEGER(iwp) :: k !<
353	INTEGER(iwp) :: nt !<
354
355
356	nt = 1
357	DO WHILE ( time > timenudge(nt) )
358	nt = nt+1
359	ENDDO
360	IF ( time /= timenudge(1) ) THEN
361	nt = nt-1
362	ENDIF
363
364	dtm = ( time - timenudge(nt) ) / ( timenudge(nt+1) - timenudge(nt) )
365	dtp = ( timenudge(nt+1) - time ) / ( timenudge(nt+1) - timenudge(nt) )
366
367	SELECT CASE ( prog_var )
368
369	CASE ( 'u' )
370
371	DO i = nxl, nxr
372	DO j = nys, nyn
373
374	DO k = nzb+1, nzt
375
376	tmp_tend = - ( hom(k,1,1,0) - ( unudge(k,nt) * dtp + &
377	unudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
378
379	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
380	MERGE( 1.0_wp, 0.0_wp, &
381	BTEST( wall_flags_0(k,j,i), 1 ) )
382
383	sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend * &
384	weight_substep(intermediate_timestep_count)
385	ENDDO
386
387	sums_ls_l(nzt+1,6) = sums_ls_l(nzt,6)
388
389	ENDDO
390	ENDDO
391
392	CASE ( 'v' )
393
394	DO i = nxl, nxr
395	DO j = nys, nyn
396
397	DO k = nzb+1, nzt
398
399	tmp_tend = - ( hom(k,1,2,0) - ( vnudge(k,nt) * dtp + &
400	vnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
401
402	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
403	MERGE( 1.0_wp, 0.0_wp, &
404	BTEST( wall_flags_0(k,j,i), 2 ) )
405
406	sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend * &
407	weight_substep(intermediate_timestep_count)
408	ENDDO
409
410	sums_ls_l(nzt+1,7) = sums_ls_l(nzt,7)
411
412	ENDDO
413	ENDDO
414
415	CASE ( 'pt' )
416
417	DO i = nxl, nxr
418	DO j = nys, nyn
419
420	DO k = nzb+1, nzt
421
422	tmp_tend = - ( hom(k,1,4,0) - ( ptnudge(k,nt) * dtp + &
423	ptnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
424
425	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
426	MERGE( 1.0_wp, 0.0_wp, &
427	BTEST( wall_flags_0(k,j,i), 0 ) )
428
429	sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend * &
430	weight_substep(intermediate_timestep_count)
431	ENDDO
432
433	sums_ls_l(nzt+1,4) = sums_ls_l(nzt,4)
434
435	ENDDO
436	ENDDO
437
438	CASE ( 'q' )
439
440	DO i = nxl, nxr
441	DO j = nys, nyn
442
443	DO k = nzb+1, nzt
444
445	tmp_tend = - ( hom(k,1,41,0) - ( qnudge(k,nt) * dtp + &
446	qnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
447
448	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
449	MERGE( 1.0_wp, 0.0_wp, &
450	BTEST( wall_flags_0(k,j,i), 0 ) )
451
452	sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend * &
453	weight_substep(intermediate_timestep_count)
454	ENDDO
455
456	sums_ls_l(nzt+1,5) = sums_ls_l(nzt,5)
457
458	ENDDO
459	ENDDO
460
461	CASE DEFAULT
462	message_string = 'unknown prognostic variable "' // prog_var // '"'
463	CALL message( 'nudge', 'PA0367', 1, 2, 0, 6, 0 )
464
465	END SELECT
466
467	END SUBROUTINE nudge
468
469
470	!------------------------------------------------------------------------------!
471	! Description:
472	! ------------
473	!> Call for grid point i,j
474	!------------------------------------------------------------------------------!
475
476	SUBROUTINE nudge_ij( i, j, time, prog_var )
477
478	USE arrays_3d, &
479	ONLY: pt, ptnudge, q, qnudge, tend, timenudge, tmp_tnudge, tnudge, &
480	u, unudge, v, vnudge
481
482	USE control_parameters, &
483	ONLY: dt_3d, intermediate_timestep_count, message_string
484
485	USE indices, &
486	ONLY: nxl, nxr, nys, nyn, nzb, nzt, wall_flags_0
487
488	USE kinds
489
490	USE statistics, &
491	ONLY: hom, sums_ls_l, weight_substep
492
493	IMPLICIT NONE
494
495
496	CHARACTER (LEN=*) :: prog_var !<
497
498	REAL(wp) :: tmp_tend !<
499	REAL(wp) :: dtm !<
500	REAL(wp) :: dtp !<
501	REAL(wp) :: time !<
502
503	INTEGER(iwp) :: i !<
504	INTEGER(iwp) :: j !<
505	INTEGER(iwp) :: k !<
506	INTEGER(iwp) :: nt !<
507
508
509	nt = 1
510	DO WHILE ( time > timenudge(nt) )
511	nt = nt+1
512	ENDDO
513	IF ( time /= timenudge(1) ) THEN
514	nt = nt-1
515	ENDIF
516
517	dtm = ( time - timenudge(nt) ) / ( timenudge(nt+1) - timenudge(nt) )
518	dtp = ( timenudge(nt+1) - time ) / ( timenudge(nt+1) - timenudge(nt) )
519
520	SELECT CASE ( prog_var )
521
522	CASE ( 'u' )
523
524	DO k = nzb+1, nzt
525
526	tmp_tend = - ( hom(k,1,1,0) - ( unudge(k,nt) * dtp + &
527	unudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
528
529	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
530	MERGE( 1.0_wp, 0.0_wp, &
531	BTEST( wall_flags_0(k,j,i), 1 ) )
532
533	sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend &
534	* weight_substep(intermediate_timestep_count)
535	ENDDO
536
537	sums_ls_l(nzt+1,6) = sums_ls_l(nzt,6)
538
539	CASE ( 'v' )
540
541	DO k = nzb+1, nzt
542
543	tmp_tend = - ( hom(k,1,2,0) - ( vnudge(k,nt) * dtp + &
544	vnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
545
546	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
547	MERGE( 1.0_wp, 0.0_wp, &
548	BTEST( wall_flags_0(k,j,i), 2 ) )
549
550	sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend &
551	* weight_substep(intermediate_timestep_count)
552	ENDDO
553
554	sums_ls_l(nzt+1,7) = sums_ls_l(nzt,7)
555
556	CASE ( 'pt' )
557
558	DO k = nzb+1, nzt
559
560	tmp_tend = - ( hom(k,1,4,0) - ( ptnudge(k,nt) * dtp + &
561	ptnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
562
563	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
564	MERGE( 1.0_wp, 0.0_wp, &
565	BTEST( wall_flags_0(k,j,i), 0 ) )
566
567	sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend &
568	* weight_substep(intermediate_timestep_count)
569	ENDDO
570
571	sums_ls_l(nzt+1,4) = sums_ls_l(nzt,4)
572
573
574	CASE ( 'q' )
575
576	DO k = nzb+1, nzt
577
578	tmp_tend = - ( hom(k,1,41,0) - ( qnudge(k,nt) * dtp + &
579	qnudge(k,nt+1) * dtm ) ) / tmp_tnudge(k)
580
581	tend(k,j,i) = tend(k,j,i) + tmp_tend * &
582	MERGE( 1.0_wp, 0.0_wp, &
583	BTEST( wall_flags_0(k,j,i), 0 ) )
584
585	sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend &
586	* weight_substep(intermediate_timestep_count)
587	ENDDO
588
589	sums_ls_l(nzt+1,5) = sums_ls_l(nzt,5)
590
591	CASE DEFAULT
592	message_string = 'unknown prognostic variable "' // prog_var // '"'
593	CALL message( 'nudge', 'PA0367', 1, 2, 0, 6, 0 )
594
595	END SELECT
596
597
598	END SUBROUTINE nudge_ij
599
600
601	!------------------------------------------------------------------------------!
602	! Description:
603	! ------------
604	!> @todo Missing subroutine description.
605	!------------------------------------------------------------------------------!
606	SUBROUTINE nudge_ref ( time )
607
608	USE arrays_3d, &
609	ONLY: time_vert, ptnudge, pt_init, qnudge, q_init, unudge, u_init, &
610	vnudge, v_init
611
612	USE kinds
613
614
615	IMPLICIT NONE
616
617	INTEGER(iwp) :: nt !<
618
619	REAL(wp) :: fac !<
620	REAL(wp), INTENT(in) :: time !<
621
622	!
623	!-- Interpolation in time of NUDGING_DATA for pt_init and q_init. This is
624	!-- needed for correct upper boundary conditions for pt and q and in case that
625	! large scale subsidence as well as scalar Rayleigh-damping are used
626	nt = 1
627	DO WHILE ( time > time_vert(nt) )
628	nt = nt + 1
629	ENDDO
630	IF ( time /= time_vert(nt) ) THEN
631	nt = nt - 1
632	ENDIF
633
634	fac = ( time-time_vert(nt) ) / ( time_vert(nt+1)-time_vert(nt) )
635
636	pt_init = ptnudge(:,nt) + fac * ( ptnudge(:,nt+1) - ptnudge(:,nt) )
637	q_init = qnudge(:,nt) + fac * ( qnudge(:,nt+1) - qnudge(:,nt) )
638	u_init = unudge(:,nt) + fac * ( unudge(:,nt+1) - unudge(:,nt) )
639	v_init = vnudge(:,nt) + fac * ( vnudge(:,nt+1) - vnudge(:,nt) )
640
641	END SUBROUTINE nudge_ref
642
643	END MODULE nudge_mod

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |