Home

Context Navigation

source: palm/trunk/SOURCE/module_interface.f90 @ 4384

Last change on this file since 4384 was 4361, checked in by suehring, 5 years ago
Plant-canopy: Avoid an unnecessary exchange of ghost points; remove unused arrays in pcm_rrd_local
Property svn:keywords set to `Id`
File size: 85.4 KB

Line
1	!> @file module_interface.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-2020 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! -----------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: module_interface.f90 4361 2020-01-07 12:22:38Z Giersch $
27	! Remove unused arrays in pmc_rrd_local
28	!
29	! 4360 2020-01-07 11:25:50Z suehring
30	! Add pcm_rrd_local and pcm_wrd_local
31	!
32	! 4331 2019-12-10 18:25:02Z suehring
33	! Change interface for doq_check_data_output, in order to perform further
34	! output checks.
35	!
36	! 4281 2019-10-29 15:15:39Z schwenkel
37	! Added dynamics boundary conditions
38	!
39	! 4272 2019-10-23 15:18:57Z schwenkel
40	! Further modularization of boundary conditions: moved boundary conditions to
41	! respective modules
42	!
43	! 4268 2019-10-17 11:29:38Z schwenkel
44	! Introduction of module_interface_boundary_conditions
45	!
46	! 4182 2019-08-22 15:20:23Z scharf
47	! Corrected "Former revisions" section
48	!
49	! 4173 2019-08-20 12:04:06Z gronemeier
50	! add vdi_internal_controls
51	!
52	! 4157 2019-08-14 09:19:12Z suehring
53	! Call doq_init from module interface
54	!
55	! 4132 2019-08-02 12:34:17Z suehring
56	! Bugfix in masked data output for diagnostic quantities
57	!
58	! 4131 2019-08-02 11:06:18Z monakurppa
59	! Add output of 3D plant-canopy outputs (merge from branch resler)
60	!
61	! 4048 2019-06-21 21:00:21Z knoop
62	! Moved turbulence_closure_mod calls into this module_interface
63	!
64	! 4047 2019-06-21 18:58:09Z knoop
65	! Introduction of the dynamics module
66	!
67	! 4039 2019-06-18 10:32:41Z suehring
68	! Introduce diagnostic output
69	!
70	! 4028 2019-06-13 12:21:37Z schwenkel
71	! Further modularization of particle code components
72	!
73	! 4017 2019-06-06 12:16:46Z schwenkel
74	! local_pf need INTENT(INOUT) attribute rather than INTENT(OUT). This is
75	! because INTENT(OUT) sets the array to not-defined. Especially for outputs that
76	! are not defined everywhere, e.g. land-surface outputs, this will be
77	! problematic as NaN will be output.
78	!
79	! 3987 2019-05-22 09:52:13Z kanani
80	! Introduce switchable DEBUG file output via debug_message routine
81	!
82	! 3956 2019-05-07 12:32:52Z monakurppa
83	! - Added calls for salsa_non_advective_processes and
84	! salsa_exchange_horiz_bounds
85	! - Moved the call for salsa_data_output_2d/3d before that of
86	! radiation_data_output_2d/3d. radiation_data_output_2d/3d tries to read a
87	! salsa output variable and encounters a segmentation fault for "Ntot" due
88	! to the shortoutput name
89	!
90	! 3931 2019-04-24 16:34:28Z schwenkel
91	! Changed non_transport_physics to non_advective_processes
92	!
93	! 3930 2019-04-24 14:57:18Z forkel
94	! Correct/complete module_interface introduction for chemistry model
95	!
96	! 3887 2019 -04-12 08:47:41Z schwenkel
97	! Changes related to global restructuring of location messages and introduction
98	! of additional debug messages
99	!
100	! 3880 2019 -04-08 21:43:02Z knoop
101	! Add a call for salsa_prognostic_equations
102	!
103	! 3840 2019-03-29 10:35:52Z knoop
104	! bugfix: intent of dummy arguments changed to inout
105	!
106	! 3770 2019-02-28 11:22:32Z moh.hefny
107	! removed unused variables in module_interface_check_data_output_ts
108	!
109	! 3767 08:18:02Z raasch
110	! unused variable file_index removed from subroutine parameter list
111	!
112	! 3766 2019-02-26 16:23:41Z raasch
113	! first argument removed from module_interface_rrd_*, statement added to avoid
114	! compiler warning about unused variable, file reformatted with respect to coding
115	! standards
116	!
117	! 3762 2019-02-25 16:54:16Z suehring
118	! only pass required arguments to surface_data_output_rrd_local
119	!
120	! 3747 2019-02-16 15:15:23Z gronemeier
121	! Call user_init_arrays
122	!
123	! 3745 2019-02-15 18:57:56Z suehring
124	! Add indoor model
125	!
126	! 3744 2019-02-15 18:38:58Z suehring
127	! Removed bio_check_parameters as the method is empty.
128	!
129	! 3735 2019-02-12 09:52:40Z dom_dwd_user
130	! Accepting variable j from check_parameters and passing it to
131	! bio_check_data_output
132	! Add required restart data for surface output module
133	!
134	! 3731 2019-02-11 13:06:27Z suehring
135	! Add check_parameters routine for virtual measurements
136	!
137	! 3711 2019-01-31 13:44:26Z knoop
138	! Introduced module_interface_init_checks for post-init checks
139	!
140	! 3705 2019-01-29 19:56:39Z suehring
141	! Add last_actions for virtual measurements
142	!
143	! 3704 2019-01-29 19:51:41Z suehring
144	! Some interface calls moved to module_interface + cleanup
145	!
146	! 3684 2019-01-20 20:20:58Z knoop
147	! Bugfix: made unit intend INOUT
148	!
149	! 3650 2019-01-04 13:01:33Z kanani
150	! Add restart routines for biometeorology
151	!
152	! 3649 2019-01-02 16:52:21Z suehring
153	! Initialize strings, in order to avoid compiler warnings for non-initialized
154	! characters with intent(out) attribute
155	!
156	! 3648 2019-01-02 16:35:46Z suehring
157	! 3641 2018-12-23 22:10:01Z knoop
158	! Initial implementation of the PALM module interface
159	!
160	!
161	! Description:
162	! ------------
163	!> This is the interface between the PALM model core and all its modules.
164	!>
165	!> @todo Re-format module to be consistent with coding standard
166	!------------------------------------------------------------------------------!
167	MODULE module_interface
168
169	USE indices, &
170	ONLY: nbgp, nxl, nxlg, nxr, nxrg, nys, nysg, nyn, nyng, nzb, nzt
171
172	USE kinds
173
174	!
175	!-- load module-specific control parameters.
176	!-- ToDo: move all of them to respective module or a dedicated central module
177
178	USE dynamics_mod, &
179	ONLY: dynamics_parin, &
180	dynamics_check_parameters, &
181	dynamics_check_data_output_ts, &
182	dynamics_check_data_output_pr, &
183	dynamics_check_data_output, &
184	dynamics_init_masks, &
185	dynamics_define_netcdf_grid, &
186	dynamics_init_arrays, &
187	dynamics_init, &
188	dynamics_init_checks, &
189	dynamics_header, &
190	dynamics_actions, &
191	dynamics_non_advective_processes, &
192	dynamics_exchange_horiz, &
193	dynamics_prognostic_equations, &
194	dynamics_boundary_conditions, &
195	dynamics_swap_timelevel, &
196	dynamics_3d_data_averaging, &
197	dynamics_data_output_2d, &
198	dynamics_data_output_3d, &
199	dynamics_statistics, &
200	dynamics_rrd_global, &
201	dynamics_rrd_local, &
202	dynamics_wrd_global, &
203	dynamics_wrd_local, &
204	dynamics_last_actions
205
206	USE turbulence_closure_mod, &
207	ONLY: tcm_check_parameters, &
208	tcm_check_data_output, &
209	tcm_init_arrays, &
210	tcm_init, &
211	tcm_actions, &
212	tcm_prognostic_equations, &
213	tcm_boundary_conds, &
214	tcm_swap_timelevel, &
215	tcm_3d_data_averaging, &
216	tcm_data_output_2d, &
217	tcm_data_output_3d
218
219	USE control_parameters, &
220	ONLY: air_chemistry, &
221	biometeorology, &
222	debug_output, &
223	debug_output_timestep, &
224	indoor_model, &
225	land_surface, &
226	large_scale_forcing, &
227	nesting_offline, &
228	nudging, &
229	ocean_mode, &
230	plant_canopy, &
231	salsa, &
232	surface_output, &
233	syn_turb_gen, &
234	urban_surface, &
235	vdi_checks, &
236	virtual_flight, &
237	virtual_measurement, &
238	wind_turbine
239
240	!
241	!-- load interface routines of all PALM modules
242	USE biometeorology_mod, &
243	ONLY: bio_parin, &
244	bio_check_data_output, &
245	bio_init, &
246	bio_init_checks, &
247	bio_header, &
248	bio_3d_data_averaging, &
249	bio_data_output_2d, &
250	bio_data_output_3d, &
251	bio_rrd_global, &
252	bio_rrd_local, &
253	bio_wrd_global, &
254	bio_wrd_local
255
256	USE bulk_cloud_model_mod, &
257	ONLY: bulk_cloud_model, &
258	bcm_parin, &
259	bcm_check_parameters, &
260	bcm_check_data_output_pr, &
261	bcm_check_data_output, &
262	bcm_init_arrays, &
263	bcm_init, &
264	bcm_header, &
265	bcm_actions, &
266	bcm_non_advective_processes, &
267	bcm_exchange_horiz, &
268	bcm_prognostic_equations, &
269	bcm_boundary_conditions, &
270	bcm_swap_timelevel, &
271	bcm_3d_data_averaging, &
272	bcm_data_output_2d, &
273	bcm_data_output_3d, &
274	bcm_rrd_global, &
275	bcm_wrd_global, &
276	bcm_rrd_local, &
277	bcm_wrd_local
278
279	USE chemistry_model_mod, &
280	ONLY: chem_parin, &
281	chem_check_parameters, &
282	chem_check_data_output_pr, &
283	chem_check_data_output, &
284	chem_exchange_horiz_bounds, &
285	chem_init_arrays, &
286	chem_init, &
287	chem_header, &
288	chem_actions, &
289	chem_non_advective_processes, &
290	chem_prognostic_equations, &
291	chem_boundary_conditions, &
292	chem_swap_timelevel, &
293	chem_3d_data_averaging, &
294	chem_data_output_2d, &
295	chem_data_output_3d, &
296	chem_statistics, &
297	chem_rrd_local, &
298	chem_wrd_local
299
300	USE diagnostic_output_quantities_mod, &
301	ONLY: doq_3d_data_averaging, &
302	doq_check_data_output, &
303	doq_define_netcdf_grid, &
304	doq_init, &
305	doq_output_2d, &
306	doq_output_3d, &
307	doq_wrd_local
308	! doq_rrd_local, &
309
310	USE flight_mod, &
311	ONLY: flight_parin, &
312	flight_header, &
313	flight_init, &
314	flight_rrd_global, &
315	flight_wrd_global
316
317	USE gust_mod, &
318	ONLY: gust_module_enabled, &
319	gust_parin, &
320	gust_check_parameters, &
321	gust_check_data_output_pr, &
322	gust_check_data_output, &
323	gust_init_arrays, &
324	gust_init, &
325	gust_header, &
326	gust_actions, &
327	gust_prognostic_equations, &
328	gust_swap_timelevel, &
329	gust_3d_data_averaging, &
330	gust_data_output_2d, &
331	gust_data_output_3d, &
332	gust_statistics, &
333	gust_rrd_global, &
334	gust_wrd_global, &
335	gust_rrd_local, &
336	gust_wrd_local
337
338	USE indoor_model_mod, &
339	ONLY: im_parin, &
340	im_check_data_output, &
341	im_check_parameters, &
342	im_data_output_3d, &
343	im_init
344
345	USE lagrangian_particle_model_mod, &
346	ONLY: lpm_parin, &
347	lpm_header, &
348	lpm_check_parameters, &
349	lpm_init_arrays, &
350	lpm_init, &
351	lpm_actions, &
352	lpm_rrd_global, &
353	lpm_rrd_local, &
354	lpm_wrd_local, &
355	lpm_wrd_global
356
357	USE land_surface_model_mod, &
358	ONLY: lsm_parin, &
359	lsm_check_parameters, &
360	lsm_check_data_output_pr, &
361	lsm_check_data_output, &
362	lsm_init_arrays, &
363	lsm_init, &
364	lsm_header, &
365	lsm_swap_timelevel, &
366	lsm_3d_data_averaging, &
367	lsm_data_output_2d, &
368	lsm_rrd_local, &
369	lsm_wrd_local
370
371	USE lsf_nudging_mod, &
372	ONLY: lsf_nudging_check_parameters, &
373	lsf_nudging_check_data_output_pr, &
374	lsf_init, &
375	nudge_init, &
376	lsf_nudging_header
377
378	USE multi_agent_system_mod, &
379	ONLY: mas_parin
380
381	USE nesting_offl_mod, &
382	ONLY: nesting_offl_parin, &
383	nesting_offl_check_parameters, &
384	nesting_offl_header
385
386	USE ocean_mod, &
387	ONLY: ocean_parin, &
388	ocean_check_parameters, &
389	ocean_check_data_output_pr, &
390	ocean_check_data_output, &
391	ocean_init_arrays, &
392	ocean_init, &
393	ocean_header, &
394	ocean_actions, &
395	ocean_prognostic_equations, &
396	ocean_boundary_conditions, &
397	ocean_swap_timelevel, &
398	ocean_3d_data_averaging, &
399	ocean_data_output_2d, &
400	ocean_data_output_3d, &
401	ocean_rrd_global, &
402	ocean_wrd_global, &
403	ocean_rrd_local, &
404	ocean_wrd_local
405
406	USE particle_attributes, &
407	ONLY: particle_advection
408
409	USE plant_canopy_model_mod, &
410	ONLY: pcm_parin, &
411	pcm_check_parameters, &
412	pcm_check_data_output, &
413	pcm_init, &
414	pcm_header, &
415	pcm_3d_data_averaging, &
416	pcm_data_output_3d, &
417	pcm_rrd_local, &
418	pcm_wrd_local
419
420	USE radiation_model_mod, &
421	ONLY: radiation, &
422	radiation_parin, &
423	radiation_check_parameters, &
424	radiation_check_data_output_ts, &
425	radiation_check_data_output_pr, &
426	radiation_check_data_output, &
427	radiation_init, &
428	radiation_header, &
429	radiation_3d_data_averaging, &
430	radiation_data_output_2d, &
431	radiation_data_output_3d, &
432	radiation_rrd_local, &
433	radiation_wrd_local
434
435	USE salsa_mod, &
436	ONLY: salsa_parin, &
437	salsa_check_parameters, &
438	salsa_check_data_output_pr, &
439	salsa_check_data_output, &
440	salsa_init_arrays, &
441	salsa_init, &
442	salsa_header, &
443	salsa_actions, &
444	salsa_non_advective_processes, &
445	salsa_exchange_horiz_bounds, &
446	salsa_prognostic_equations, &
447	salsa_boundary_conditions, &
448	salsa_swap_timelevel, &
449	salsa_3d_data_averaging, &
450	salsa_data_output_2d, &
451	salsa_data_output_3d, &
452	salsa_statistics, &
453	salsa_rrd_local, &
454	salsa_wrd_local
455
456	USE spectra_mod, &
457	ONLY: calculate_spectra, &
458	spectra_parin, &
459	spectra_check_parameters, &
460	spectra_header
461
462	USE surface_data_output_mod, &
463	ONLY: surface_data_output_parin, &
464	surface_data_output_check_parameters, &
465	surface_data_output_init_arrays, &
466	surface_data_output_rrd_local, &
467	surface_data_output_rrd_global, &
468	surface_data_output_wrd_local, &
469	surface_data_output_wrd_global
470
471	USE synthetic_turbulence_generator_mod, &
472	ONLY: stg_parin, &
473	stg_check_parameters, &
474	stg_header, &
475	stg_rrd_global, &
476	stg_wrd_global
477
478	USE urban_surface_mod, &
479	ONLY: usm_parin, &
480	usm_check_parameters, &
481	usm_check_data_output, &
482	usm_init_arrays, &
483	usm_init, &
484	usm_swap_timelevel, &
485	usm_3d_data_averaging, &
486	usm_rrd_local, &
487	usm_wrd_local
488
489	USE vdi_internal_controls, &
490	ONLY: vdi_actions
491
492	USE virtual_measurement_mod, &
493	ONLY: vm_check_parameters, &
494	vm_init, &
495	vm_last_actions, &
496	vm_parin
497
498	USE wind_turbine_model_mod, &
499	ONLY: wtm_parin, &
500	wtm_check_parameters, &
501	wtm_init_arrays, &
502	wtm_init, &
503	wtm_actions, &
504	wtm_rrd_global, &
505	wtm_wrd_global
506
507	USE user, &
508	ONLY: user_module_enabled, &
509	user_parin, &
510	user_check_parameters, &
511	user_check_data_output_ts, &
512	user_check_data_output_pr, &
513	user_check_data_output, &
514	user_init, &
515	user_init_arrays, &
516	user_header, &
517	user_actions, &
518	user_3d_data_averaging, &
519	user_data_output_2d, &
520	user_data_output_3d, &
521	user_statistics, &
522	user_rrd_global, &
523	user_rrd_local, &
524	user_wrd_global, &
525	user_wrd_local, &
526	user_last_actions
527
528	IMPLICIT NONE
529
530	PRIVATE
531
532	!
533	!-- Public functions
534	PUBLIC &
535	module_interface_parin, &
536	module_interface_check_parameters, &
537	module_interface_check_data_output_ts, &
538	module_interface_check_data_output_pr, &
539	module_interface_check_data_output, &
540	module_interface_init_masks, &
541	module_interface_define_netcdf_grid, &
542	module_interface_init_arrays, &
543	module_interface_init, &
544	module_interface_init_checks, &
545	module_interface_header, &
546	module_interface_actions, &
547	module_interface_non_advective_processes, &
548	module_interface_exchange_horiz, &
549	module_interface_prognostic_equations, &
550	module_interface_boundary_conditions, &
551	module_interface_swap_timelevel, &
552	module_interface_3d_data_averaging, &
553	module_interface_data_output_2d, &
554	module_interface_data_output_3d, &
555	module_interface_statistics, &
556	module_interface_rrd_global, &
557	module_interface_wrd_global, &
558	module_interface_rrd_local, &
559	module_interface_wrd_local, &
560	module_interface_last_actions
561
562
563	INTERFACE module_interface_parin
564	MODULE PROCEDURE module_interface_parin
565	END INTERFACE module_interface_parin
566
567	INTERFACE module_interface_check_parameters
568	MODULE PROCEDURE module_interface_check_parameters
569	END INTERFACE module_interface_check_parameters
570
571	INTERFACE module_interface_check_data_output_ts
572	MODULE PROCEDURE module_interface_check_data_output_ts
573	END INTERFACE module_interface_check_data_output_ts
574
575	INTERFACE module_interface_check_data_output_pr
576	MODULE PROCEDURE module_interface_check_data_output_pr
577	END INTERFACE module_interface_check_data_output_pr
578
579	INTERFACE module_interface_check_data_output
580	MODULE PROCEDURE module_interface_check_data_output
581	END INTERFACE module_interface_check_data_output
582
583	INTERFACE module_interface_init_masks
584	MODULE PROCEDURE module_interface_init_masks
585	END INTERFACE module_interface_init_masks
586
587	INTERFACE module_interface_define_netcdf_grid
588	MODULE PROCEDURE module_interface_define_netcdf_grid
589	END INTERFACE module_interface_define_netcdf_grid
590
591	INTERFACE module_interface_init_arrays
592	MODULE PROCEDURE module_interface_init_arrays
593	END INTERFACE module_interface_init_arrays
594
595	INTERFACE module_interface_init
596	MODULE PROCEDURE module_interface_init
597	END INTERFACE module_interface_init
598
599	INTERFACE module_interface_init_checks
600	MODULE PROCEDURE module_interface_init_checks
601	END INTERFACE module_interface_init_checks
602
603	INTERFACE module_interface_header
604	MODULE PROCEDURE module_interface_header
605	END INTERFACE module_interface_header
606
607	INTERFACE module_interface_actions
608	MODULE PROCEDURE module_interface_actions
609	MODULE PROCEDURE module_interface_actions_ij
610	END INTERFACE module_interface_actions
611
612	INTERFACE module_interface_non_advective_processes
613	MODULE PROCEDURE module_interface_non_advective_processes
614	MODULE PROCEDURE module_interface_non_advective_processes_ij
615	END INTERFACE module_interface_non_advective_processes
616
617	INTERFACE module_interface_exchange_horiz
618	MODULE PROCEDURE module_interface_exchange_horiz
619	END INTERFACE module_interface_exchange_horiz
620
621	INTERFACE module_interface_prognostic_equations
622	MODULE PROCEDURE module_interface_prognostic_equations
623	MODULE PROCEDURE module_interface_prognostic_equations_ij
624	END INTERFACE module_interface_prognostic_equations
625
626	INTERFACE module_interface_swap_timelevel
627	MODULE PROCEDURE module_interface_swap_timelevel
628	END INTERFACE module_interface_swap_timelevel
629
630	INTERFACE module_interface_boundary_conditions
631	MODULE PROCEDURE module_interface_boundary_conditions
632	END INTERFACE module_interface_boundary_conditions
633
634	INTERFACE module_interface_3d_data_averaging
635	MODULE PROCEDURE module_interface_3d_data_averaging
636	END INTERFACE module_interface_3d_data_averaging
637
638	INTERFACE module_interface_data_output_2d
639	MODULE PROCEDURE module_interface_data_output_2d
640	END INTERFACE module_interface_data_output_2d
641
642	INTERFACE module_interface_data_output_3d
643	MODULE PROCEDURE module_interface_data_output_3d
644	END INTERFACE module_interface_data_output_3d
645
646	INTERFACE module_interface_statistics
647	MODULE PROCEDURE module_interface_statistics
648	END INTERFACE module_interface_statistics
649
650	INTERFACE module_interface_rrd_global
651	MODULE PROCEDURE module_interface_rrd_global
652	END INTERFACE module_interface_rrd_global
653
654	INTERFACE module_interface_wrd_global
655	MODULE PROCEDURE module_interface_wrd_global
656	END INTERFACE module_interface_wrd_global
657
658	INTERFACE module_interface_rrd_local
659	MODULE PROCEDURE module_interface_rrd_local
660	END INTERFACE module_interface_rrd_local
661
662	INTERFACE module_interface_wrd_local
663	MODULE PROCEDURE module_interface_wrd_local
664	END INTERFACE module_interface_wrd_local
665
666	INTERFACE module_interface_last_actions
667	MODULE PROCEDURE module_interface_last_actions
668	END INTERFACE module_interface_last_actions
669
670
671	CONTAINS
672
673
674	!------------------------------------------------------------------------------!
675	! Description:
676	! ------------
677	!> Read module-specific parameter namelists
678	!------------------------------------------------------------------------------!
679	SUBROUTINE module_interface_parin
680
681
682	IF ( debug_output ) CALL debug_message( 'reading module-specific parameters', 'start' )
683
684	CALL dynamics_parin
685
686	CALL bio_parin
687	CALL bcm_parin
688	CALL chem_parin
689	CALL flight_parin ! ToDo: rename module to match filename
690	CALL gust_parin
691	CALL im_parin
692	CALL lpm_parin
693	CALL lsm_parin
694	! ToDo: create parin routine for large_scale_forcing and nudging (should be seperate modules or new module switch)
695	CALL mas_parin
696	CALL nesting_offl_parin
697	CALL ocean_parin
698	CALL pcm_parin
699	CALL radiation_parin
700	CALL salsa_parin
701	CALL spectra_parin
702	CALL surface_data_output_parin
703	CALL stg_parin
704	CALL usm_parin
705	CALL vm_parin
706	CALL wtm_parin
707
708	CALL user_parin
709
710	IF ( debug_output ) CALL debug_message( 'reading module-specific parameters', 'end' )
711
712
713	END SUBROUTINE module_interface_parin
714
715
716	!------------------------------------------------------------------------------!
717	! Description:
718	! ------------
719	!> Perform module-specific initialization checks
720	!------------------------------------------------------------------------------!
721	SUBROUTINE module_interface_check_parameters
722
723
724	IF ( debug_output ) CALL debug_message( 'checking module-specific parameters', 'start' )
725
726	CALL dynamics_check_parameters
727	CALL tcm_check_parameters
728
729	IF ( bulk_cloud_model ) CALL bcm_check_parameters
730	IF ( air_chemistry ) CALL chem_check_parameters
731	IF ( gust_module_enabled ) CALL gust_check_parameters
732	IF ( indoor_model ) CALL im_check_parameters
733	IF ( particle_advection ) CALL lpm_check_parameters
734	IF ( land_surface ) CALL lsm_check_parameters
735	IF ( large_scale_forcing .OR. nudging ) CALL lsf_nudging_check_parameters ! ToDo: create single module switch
736	IF ( nesting_offline ) CALL nesting_offl_check_parameters
737	IF ( ocean_mode ) CALL ocean_check_parameters
738	IF ( plant_canopy ) CALL pcm_check_parameters
739	IF ( radiation ) CALL radiation_check_parameters
740	IF ( salsa ) CALL salsa_check_parameters
741	IF ( calculate_spectra ) CALL spectra_check_parameters
742	IF ( surface_output ) CALL surface_data_output_check_parameters
743	IF ( syn_turb_gen ) CALL stg_check_parameters
744	IF ( urban_surface ) CALL usm_check_parameters
745	IF ( virtual_measurement ) CALL vm_check_parameters
746	IF ( wind_turbine ) CALL wtm_check_parameters
747
748	IF ( user_module_enabled ) CALL user_check_parameters
749
750	IF ( debug_output ) CALL debug_message( 'checking module-specific parameters', 'end' )
751
752
753	END SUBROUTINE module_interface_check_parameters
754
755
756	!------------------------------------------------------------------------------!
757	! Description:
758	! ------------
759	!> Check module-specific data output of timeseries
760	!------------------------------------------------------------------------------!
761	SUBROUTINE module_interface_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
762
763
764	INTEGER(iwp), INTENT(IN) :: dots_max !< variable output array index
765	INTEGER(iwp), INTENT(INOUT) :: dots_num !< variable output array index
766	CHARACTER (LEN=*), DIMENSION(dots_max), INTENT(INOUT) :: dots_label
767	CHARACTER (LEN=*), DIMENSION(dots_max), INTENT(INOUT) :: dots_unit
768
769
770	IF ( debug_output ) CALL debug_message( 'checking module-specific data output ts', 'start' )
771
772	CALL dynamics_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
773
774	IF ( radiation ) THEN
775	CALL radiation_check_data_output_ts( dots_max, dots_num )
776	ENDIF
777
778	IF ( user_module_enabled ) THEN
779	CALL user_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
780	ENDIF
781
782	IF ( debug_output ) CALL debug_message( 'checking module-specific data output ts', 'end' )
783
784
785	END SUBROUTINE module_interface_check_data_output_ts
786
787
788	!------------------------------------------------------------------------------!
789	! Description:
790	! ------------
791	!> Check module-specific data output of profiles
792	!------------------------------------------------------------------------------!
793	SUBROUTINE module_interface_check_data_output_pr( variable, var_count, unit, &
794	dopr_unit )
795
796
797	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
798	INTEGER(iwp), INTENT(IN) :: var_count !< variable output array index
799	CHARACTER (LEN=*), INTENT(INOUT) :: unit !< physical unit of variable
800	CHARACTER (LEN=*), INTENT(OUT) :: dopr_unit !< local value of dopr_unit
801
802
803	IF ( debug_output ) CALL debug_message( 'checking module-specific data output pr', 'start' )
804
805	CALL dynamics_check_data_output_pr( variable, var_count, unit, dopr_unit )
806
807	IF ( unit == 'illegal' .AND. bulk_cloud_model ) THEN
808	CALL bcm_check_data_output_pr( variable, var_count, unit, dopr_unit )
809	ENDIF
810
811	IF ( unit == 'illegal' .AND. air_chemistry ) THEN
812	CALL chem_check_data_output_pr( variable, var_count, unit, dopr_unit )
813	ENDIF
814
815	IF ( unit == 'illegal' .AND. gust_module_enabled ) THEN
816	CALL gust_check_data_output_pr( variable, var_count, unit, dopr_unit )
817	ENDIF
818
819	IF ( unit == 'illegal' ) THEN ! ToDo: add module switch if possible
820	CALL lsm_check_data_output_pr( variable, var_count, unit, dopr_unit )
821	ENDIF
822
823	IF ( unit == 'illegal' ) THEN ! ToDo: add module switch if possible
824	CALL lsf_nudging_check_data_output_pr( variable, var_count, unit, dopr_unit )
825	ENDIF
826
827	IF ( unit == 'illegal' .AND. ocean_mode ) THEN
828	CALL ocean_check_data_output_pr( variable, var_count, unit, dopr_unit )
829	ENDIF
830
831	IF ( unit == 'illegal' .AND. radiation ) THEN
832	CALL radiation_check_data_output_pr( variable, var_count, unit, dopr_unit )
833	ENDIF
834
835	IF ( unit == 'illegal' .AND. salsa ) THEN
836	CALL salsa_check_data_output_pr( variable, var_count, unit, dopr_unit )
837	ENDIF
838
839	IF ( unit == 'illegal' .AND. user_module_enabled ) THEN
840	unit = '' ! ToDo: Seems like a hack. Find a general soultion!
841	CALL user_check_data_output_pr( variable, var_count, unit, dopr_unit )
842	ENDIF
843
844	IF ( debug_output ) CALL debug_message( 'checking module-specific data output pr', 'end' )
845
846
847	END SUBROUTINE module_interface_check_data_output_pr
848
849	!------------------------------------------------------------------------------!
850	! Description:
851	! ------------
852	!> Check module-specific 2D and 3D data output
853	!------------------------------------------------------------------------------!
854	SUBROUTINE module_interface_check_data_output( variable, unit, i, j, ilen, k )
855
856
857	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
858	CHARACTER (LEN=*), INTENT(INOUT) :: unit !< physical unit of variable
859
860	INTEGER(iwp), INTENT(IN) :: i !< ToDo: remove dummy argument, instead pass string from data_output
861	INTEGER(iwp), INTENT(IN) :: j !< average quantity? 0 = no, 1 = yes
862	INTEGER(iwp), INTENT(IN) :: ilen !< ToDo: remove dummy argument, instead pass string from data_output
863	INTEGER(iwp), INTENT(IN) :: k !< ToDo: remove dummy argument, instead pass string from data_output
864
865
866	IF ( debug_output ) CALL debug_message( 'checking module-specific data output 2d/3d', 'start' )
867
868	CALL dynamics_check_data_output( variable, unit )
869
870	CALL tcm_check_data_output( variable, unit )
871
872	IF ( unit == 'illegal' .AND. biometeorology ) THEN
873	CALL bio_check_data_output( variable, unit, i, j, ilen, k )
874	ENDIF
875
876	IF ( unit == 'illegal' .AND. bulk_cloud_model ) THEN
877	CALL bcm_check_data_output( variable, unit )
878	ENDIF
879
880	IF ( unit == 'illegal' .AND. air_chemistry &
881	.AND. (variable(1:3) == 'kc_' .OR. variable(1:3) == 'em_') ) THEN ! ToDo: remove aditional conditions
882	CALL chem_check_data_output( variable, unit, i, ilen, k )
883	ENDIF
884
885	IF ( unit == 'illegal' ) THEN
886	CALL doq_check_data_output( variable, unit, i, ilen, k )
887	ENDIF
888
889	IF ( unit == 'illegal' .AND. gust_module_enabled ) THEN
890	CALL gust_check_data_output( variable, unit )
891	ENDIF
892
893	IF ( unit == 'illegal' ) THEN ! ToDo: add module switch if possible
894	CALL lsm_check_data_output( variable, unit, i, ilen, k )
895	ENDIF
896
897	IF ( unit == 'illegal' .AND. ocean_mode ) THEN
898	CALL ocean_check_data_output( variable, unit )
899	ENDIF
900
901	IF ( unit == 'illegal' .AND. plant_canopy &
902	.AND. variable(1:4) == 'pcm_' ) THEN ! ToDo: remove aditional conditions
903	CALL pcm_check_data_output( variable, unit )
904	ENDIF
905
906	IF ( unit == 'illegal' .AND. radiation ) THEN
907	CALL radiation_check_data_output( variable, unit, i, ilen, k )
908	ENDIF
909
910	IF ( unit == 'illegal' .AND. salsa ) THEN
911	CALL salsa_check_data_output( variable, unit )
912	ENDIF
913
914	IF ( unit == 'illegal' .AND. indoor_model ) THEN
915	CALL im_check_data_output( variable, unit )
916	ENDIF
917
918	IF ( unit == 'illegal' .AND. urban_surface &
919	.AND. variable(1:4) == 'usm_' ) THEN ! ToDo: remove aditional conditions
920	CALL usm_check_data_output( variable, unit )
921	ENDIF
922
923	IF ( unit == 'illegal' .AND. user_module_enabled ) THEN
924	unit = ''
925	CALL user_check_data_output( variable, unit )
926	ENDIF
927
928	IF ( debug_output ) CALL debug_message( 'checking module-specific data output 2d/3d', 'end' )
929
930
931	END SUBROUTINE module_interface_check_data_output
932
933
934	!------------------------------------------------------------------------------!
935	!
936	! Description:
937	! ------------
938	!> Interface for init_masks. ToDo: get rid of these redundant calls!
939	!------------------------------------------------------------------------------!
940	SUBROUTINE module_interface_init_masks( variable, unit )
941
942
943	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
944	CHARACTER (LEN=*), INTENT(INOUT) :: unit !< physical unit of variable
945
946
947	IF ( debug_output ) CALL debug_message( 'initializing module-specific masks', 'start' )
948
949	CALL dynamics_init_masks( variable, unit )
950
951	IF ( unit == 'illegal' .AND. air_chemistry &
952	.AND. (variable(1:3) == 'kc_' .OR. variable(1:3) == 'em_') ) THEN ! ToDo: remove aditional conditions
953	CALL chem_check_data_output( variable, unit, 0, 0, 0 )
954	ENDIF
955
956	IF ( unit == 'illegal' ) THEN
957	CALL doq_check_data_output( variable, unit )
958	ENDIF
959
960	IF ( unit == 'illegal' .AND. radiation ) THEN
961	CALL radiation_check_data_output( variable, unit, 0, 0, 0 )
962	ENDIF
963
964	IF ( unit == 'illegal' .AND. salsa ) THEN
965	CALL salsa_check_data_output( variable, unit )
966	ENDIF
967
968	IF ( unit == 'illegal' .AND. user_module_enabled ) THEN
969	unit = ''
970	CALL user_check_data_output( variable, unit )
971	ENDIF
972
973	IF ( debug_output ) CALL debug_message( 'initializing module-specific masks', 'end' )
974
975
976	END SUBROUTINE module_interface_init_masks
977
978
979	!------------------------------------------------------------------------------!
980	!
981	! Description:
982	! ------------
983	!> Define appropriate grid for module-specific netcdf output variables.
984	!------------------------------------------------------------------------------!
985	SUBROUTINE module_interface_define_netcdf_grid( var, found, &
986	grid_x, grid_y, grid_z )
987
988
989	CHARACTER (LEN=*), INTENT(IN) :: var !< variable name
990	LOGICAL, INTENT(OUT) :: found !< indicates if variable was found
991	CHARACTER (LEN=*), INTENT(OUT) :: grid_x !< netcdf dimension in x-direction
992	CHARACTER (LEN=*), INTENT(OUT) :: grid_y !< netcdf dimension in y-direction
993	CHARACTER (LEN=*), INTENT(OUT) :: grid_z !< netcdf dimension in z-direction
994
995
996	IF ( debug_output ) CALL debug_message( 'defining module-specific netcdf grids', 'start' )
997	!
998	!-- As long as no action is done in this subroutine, initialize strings with
999	!-- intent(out) attribute, in order to avoid compiler warnings.
1000	found = .FALSE.
1001	grid_x = 'none'
1002	grid_y = 'none'
1003	grid_z = 'none'
1004	!
1005	!-- Use var to avoid compiler warning about unused variable
1006	IF ( var == ' ' ) RETURN
1007
1008	IF ( debug_output ) CALL debug_message( 'defining module-specific netcdf grids', 'end' )
1009
1010
1011	END SUBROUTINE module_interface_define_netcdf_grid
1012
1013
1014	!------------------------------------------------------------------------------!
1015	! Description:
1016	! ------------
1017	!> Allocate module-specific arrays and pointers
1018	!------------------------------------------------------------------------------!
1019	SUBROUTINE module_interface_init_arrays
1020
1021
1022	IF ( debug_output ) CALL debug_message( 'initializing module-specific arrays', 'start' )
1023
1024	CALL dynamics_init_arrays
1025	CALL tcm_init_arrays
1026
1027	IF ( bulk_cloud_model ) CALL bcm_init_arrays
1028	IF ( air_chemistry ) CALL chem_init_arrays
1029	IF ( gust_module_enabled ) CALL gust_init_arrays
1030	IF ( particle_advection ) CALL lpm_init_arrays
1031	IF ( land_surface ) CALL lsm_init_arrays
1032	IF ( ocean_mode ) CALL ocean_init_arrays
1033	IF ( salsa ) CALL salsa_init_arrays
1034	IF ( urban_surface ) CALL usm_init_arrays
1035	IF ( surface_output ) CALL surface_data_output_init_arrays
1036	IF ( wind_turbine ) CALL wtm_init_arrays
1037
1038	IF ( user_module_enabled ) CALL user_init_arrays
1039
1040	IF ( debug_output ) CALL debug_message( 'initializing module-specific arrays', 'end' )
1041
1042
1043	END SUBROUTINE module_interface_init_arrays
1044
1045
1046	!------------------------------------------------------------------------------!
1047	! Description:
1048	! ------------
1049	!> Perform module-specific initialization
1050	!------------------------------------------------------------------------------!
1051	SUBROUTINE module_interface_init
1052
1053
1054	IF ( debug_output ) CALL debug_message( 'module-specific initialization', 'start' )
1055
1056	CALL dynamics_init
1057	CALL tcm_init
1058
1059	IF ( biometeorology ) CALL bio_init
1060	IF ( bulk_cloud_model ) CALL bcm_init
1061	IF ( air_chemistry ) CALL chem_init
1062	IF ( virtual_flight ) CALL flight_init
1063	IF ( gust_module_enabled ) CALL gust_init
1064	IF ( indoor_model ) CALL im_init
1065	IF ( particle_advection ) CALL lpm_init
1066	IF ( large_scale_forcing ) CALL lsf_init
1067	IF ( land_surface ) CALL lsm_init
1068	IF ( nudging ) CALL nudge_init
1069	IF ( ocean_mode ) CALL ocean_init
1070	IF ( plant_canopy ) CALL pcm_init
1071	IF ( salsa ) CALL salsa_init
1072	IF ( urban_surface ) CALL usm_init
1073	IF ( virtual_measurement ) CALL vm_init
1074	IF ( wind_turbine ) CALL wtm_init
1075	IF ( radiation ) CALL radiation_init
1076
1077	CALL doq_init
1078
1079	IF ( user_module_enabled ) CALL user_init
1080
1081	IF ( debug_output ) CALL debug_message( 'module-specific initialization', 'end' )
1082
1083
1084	END SUBROUTINE module_interface_init
1085
1086
1087	!------------------------------------------------------------------------------!
1088	! Description:
1089	! ------------
1090	!> Perform module-specific post-initialization checks
1091	!------------------------------------------------------------------------------!
1092	SUBROUTINE module_interface_init_checks
1093
1094
1095	IF ( debug_output ) CALL debug_message( 'module-specific post-initialization checks', 'start' )
1096
1097	CALL dynamics_init_checks
1098
1099	IF ( biometeorology ) CALL bio_init_checks
1100
1101	IF ( debug_output ) CALL debug_message( 'module-specific post-initialization checks', 'end' )
1102
1103
1104	END SUBROUTINE module_interface_init_checks
1105
1106
1107	!------------------------------------------------------------------------------!
1108	! Description:
1109	! ------------
1110	!> Gather module-specific header output
1111	!------------------------------------------------------------------------------!
1112	SUBROUTINE module_interface_header( io )
1113
1114
1115	INTEGER(iwp), INTENT(IN) :: io !< unit of the output file
1116
1117
1118	IF ( debug_output ) CALL debug_message( 'module-specific header output', 'start' )
1119
1120	CALL dynamics_header( io )
1121
1122	IF ( biometeorology ) CALL bio_header ( io )
1123	IF ( bulk_cloud_model ) CALL bcm_header( io )
1124	IF ( air_chemistry ) CALL chem_header ( io )
1125	IF ( virtual_flight ) CALL flight_header( io )
1126	IF ( gust_module_enabled ) CALL gust_header( io )
1127	IF ( particle_advection ) CALL lpm_header( io )
1128	IF ( land_surface ) CALL lsm_header( io )
1129	IF ( large_scale_forcing ) CALL lsf_nudging_header( io )
1130	IF ( nesting_offline ) CALL nesting_offl_header( io )
1131	IF ( ocean_mode ) CALL ocean_header( io )
1132	IF ( plant_canopy ) CALL pcm_header( io )
1133	IF ( radiation ) CALL radiation_header( io )
1134	IF ( salsa ) CALL salsa_header( io )
1135	IF ( calculate_spectra ) CALL spectra_header( io )
1136	IF ( syn_turb_gen ) CALL stg_header( io )
1137
1138	IF ( user_module_enabled ) CALL user_header( io )
1139
1140	IF ( debug_output ) CALL debug_message( 'module-specific header output', 'end' )
1141
1142
1143	END SUBROUTINE module_interface_header
1144
1145
1146	!------------------------------------------------------------------------------!
1147	! Description:
1148	! ------------
1149	!> Perform module-specific actions while in time-integration (vector-optimized)
1150	!------------------------------------------------------------------------------!
1151	SUBROUTINE module_interface_actions( location )
1152
1153
1154	CHARACTER (LEN=*), INTENT(IN) :: location !< call location string
1155
1156	CALL dynamics_actions( location )
1157	CALL tcm_actions( location )
1158
1159	IF ( bulk_cloud_model ) CALL bcm_actions( location )
1160	IF ( air_chemistry ) CALL chem_actions( location )
1161	IF ( gust_module_enabled ) CALL gust_actions( location )
1162	IF ( particle_advection ) CALL lpm_actions( location )
1163	IF ( ocean_mode ) CALL ocean_actions( location )
1164	IF ( salsa ) CALL salsa_actions( location )
1165	IF ( wind_turbine ) CALL wtm_actions( location )
1166
1167	IF ( user_module_enabled ) CALL user_actions( location )
1168	IF ( vdi_checks ) CALL vdi_actions( location )
1169
1170
1171	END SUBROUTINE module_interface_actions
1172
1173
1174	!------------------------------------------------------------------------------!
1175	! Description:
1176	! ------------
1177	!> Perform module-specific actions while in time-integration (cache-optimized)
1178	!------------------------------------------------------------------------------!
1179	SUBROUTINE module_interface_actions_ij( i, j, location )
1180
1181
1182	INTEGER(iwp), INTENT(IN) :: i !< grid index in x-direction
1183	INTEGER(iwp), INTENT(IN) :: j !< grid index in y-direction
1184	CHARACTER (LEN=*), INTENT(IN) :: location !< call location string
1185
1186	CALL dynamics_actions( i, j, location )
1187	CALL tcm_actions( i, j, location )
1188
1189	IF ( bulk_cloud_model ) CALL bcm_actions( i, j, location )
1190	IF ( air_chemistry ) CALL chem_actions( i, j, location )
1191	IF ( gust_module_enabled ) CALL gust_actions( i, j, location )
1192	IF ( ocean_mode ) CALL ocean_actions( i, j, location )
1193	IF ( salsa ) CALL salsa_actions( i, j, location )
1194	IF ( wind_turbine ) CALL wtm_actions( i, j, location )
1195
1196	IF ( user_module_enabled ) CALL user_actions( i, j, location )
1197
1198
1199	END SUBROUTINE module_interface_actions_ij
1200
1201
1202	!------------------------------------------------------------------------------!
1203	! Description:
1204	! ------------
1205	!> Compute module-specific non_advective_processes (vector-optimized)
1206	!------------------------------------------------------------------------------!
1207	SUBROUTINE module_interface_non_advective_processes
1208
1209
1210	CALL dynamics_non_advective_processes
1211
1212	IF ( bulk_cloud_model ) CALL bcm_non_advective_processes
1213	IF ( air_chemistry ) CALL chem_non_advective_processes
1214	IF ( salsa ) CALL salsa_non_advective_processes
1215
1216
1217	END SUBROUTINE module_interface_non_advective_processes
1218
1219
1220	!------------------------------------------------------------------------------!
1221	! Description:
1222	! ------------
1223	!> Compute module-specific non_advective_processes (cache-optimized)
1224	!------------------------------------------------------------------------------!
1225	SUBROUTINE module_interface_non_advective_processes_ij( i, j )
1226
1227
1228	INTEGER(iwp), INTENT(IN) :: i !< grid index in x-direction
1229	INTEGER(iwp), INTENT(IN) :: j !< grid index in y-direction
1230
1231	CALL dynamics_non_advective_processes( i, j )
1232
1233	IF ( bulk_cloud_model ) CALL bcm_non_advective_processes( i, j )
1234	IF ( air_chemistry ) CALL chem_non_advective_processes( i, j )
1235	IF ( salsa ) CALL salsa_non_advective_processes( i, j )
1236
1237
1238	END SUBROUTINE module_interface_non_advective_processes_ij
1239
1240	!------------------------------------------------------------------------------!
1241	! Description:
1242	! ------------
1243	!> Exchange horiz for module-specific quantities
1244	!------------------------------------------------------------------------------!
1245	SUBROUTINE module_interface_exchange_horiz
1246
1247
1248	IF ( debug_output_timestep ) CALL debug_message( 'module-specific exchange_horiz', 'start' )
1249
1250	CALL dynamics_exchange_horiz
1251
1252	IF ( bulk_cloud_model ) CALL bcm_exchange_horiz
1253	IF ( air_chemistry ) CALL chem_exchange_horiz_bounds
1254	IF ( salsa ) CALL salsa_exchange_horiz_bounds
1255
1256	IF ( debug_output_timestep ) CALL debug_message( 'module-specific exchange_horiz', 'end' )
1257
1258
1259	END SUBROUTINE module_interface_exchange_horiz
1260
1261
1262	!------------------------------------------------------------------------------!
1263	! Description:
1264	! ------------
1265	!> Compute module-specific prognostic_equations (vector-optimized)
1266	!------------------------------------------------------------------------------!
1267	SUBROUTINE module_interface_prognostic_equations
1268
1269
1270	CALL dynamics_prognostic_equations
1271	CALL tcm_prognostic_equations
1272
1273	IF ( bulk_cloud_model ) CALL bcm_prognostic_equations
1274	IF ( air_chemistry ) CALL chem_prognostic_equations
1275	IF ( gust_module_enabled ) CALL gust_prognostic_equations
1276	IF ( ocean_mode ) CALL ocean_prognostic_equations
1277	IF ( salsa ) CALL salsa_prognostic_equations
1278
1279
1280	END SUBROUTINE module_interface_prognostic_equations
1281
1282
1283	!------------------------------------------------------------------------------!
1284	! Description:
1285	! ------------
1286	!> Compute module-specific prognostic_equations (cache-optimized)
1287	!------------------------------------------------------------------------------!
1288	SUBROUTINE module_interface_prognostic_equations_ij( i, j, i_omp_start, tn )
1289
1290
1291	INTEGER(iwp), INTENT(IN) :: i !< grid index in x-direction
1292	INTEGER(iwp), INTENT(IN) :: j !< grid index in y-direction
1293	INTEGER(iwp), INTENT(IN) :: i_omp_start !< first loop index of i-loop in prognostic_equations
1294	INTEGER(iwp), INTENT(IN) :: tn !< task number of openmp task
1295
1296	CALL dynamics_prognostic_equations( i, j, i_omp_start, tn )
1297	CALL tcm_prognostic_equations( i, j, i_omp_start, tn )
1298
1299	IF ( bulk_cloud_model ) CALL bcm_prognostic_equations( i, j, i_omp_start, tn )
1300	IF ( air_chemistry ) CALL chem_prognostic_equations( i, j, i_omp_start, tn )
1301	IF ( gust_module_enabled ) CALL gust_prognostic_equations( i, j, i_omp_start, tn )
1302	IF ( ocean_mode ) CALL ocean_prognostic_equations( i, j, i_omp_start, tn )
1303	IF ( salsa ) CALL salsa_prognostic_equations( i, j, i_omp_start, tn )
1304
1305
1306	END SUBROUTINE module_interface_prognostic_equations_ij
1307
1308	!------------------------------------------------------------------------------!
1309	! Description:
1310	! ------------
1311	!> Compute module-specific boundary conditions
1312	!------------------------------------------------------------------------------!
1313	SUBROUTINE module_interface_boundary_conditions
1314
1315
1316	IF ( debug_output_timestep ) CALL debug_message( 'module-specific boundary_conditions', 'start' )
1317
1318	CALL dynamics_boundary_conditions
1319	CALL tcm_boundary_conds
1320
1321	IF ( bulk_cloud_model ) CALL bcm_boundary_conditions
1322	IF ( air_chemistry ) CALL chem_boundary_conditions
1323	IF ( ocean_mode ) CALL ocean_boundary_conditions
1324	IF ( salsa ) CALL salsa_boundary_conditions
1325
1326	IF ( debug_output_timestep ) CALL debug_message( 'module-specific boundary_conditions', 'end' )
1327
1328
1329	END SUBROUTINE module_interface_boundary_conditions
1330
1331	!------------------------------------------------------------------------------!
1332	! Description:
1333	! ------------
1334	!> Swap the timelevel pointers for module-specific arrays
1335	!------------------------------------------------------------------------------!
1336	SUBROUTINE module_interface_swap_timelevel ( swap_mode )
1337
1338
1339	INTEGER(iwp), INTENT(IN) :: swap_mode !< determines procedure of pointer swap
1340
1341
1342	IF ( debug_output_timestep ) CALL debug_message( 'module-specific swap timelevel', 'start' )
1343
1344	CALL dynamics_swap_timelevel( swap_mode )
1345	CALL tcm_swap_timelevel( swap_mode )
1346
1347	IF ( bulk_cloud_model ) CALL bcm_swap_timelevel( swap_mode )
1348	IF ( air_chemistry ) CALL chem_swap_timelevel( swap_mode )
1349	IF ( gust_module_enabled ) CALL gust_swap_timelevel( swap_mode )
1350	IF ( land_surface ) CALL lsm_swap_timelevel( swap_mode )
1351	IF ( ocean_mode ) CALL ocean_swap_timelevel( swap_mode )
1352	IF ( salsa ) CALL salsa_swap_timelevel( swap_mode )
1353	IF ( urban_surface ) CALL usm_swap_timelevel( swap_mode )
1354
1355	IF ( debug_output_timestep ) CALL debug_message( 'module-specific swap timelevel', 'end' )
1356
1357
1358	END SUBROUTINE module_interface_swap_timelevel
1359
1360
1361	!------------------------------------------------------------------------------!
1362	!
1363	! Description:
1364	! ------------
1365	!> Perform module-specific averaging of 3D data
1366	!------------------------------------------------------------------------------!
1367	SUBROUTINE module_interface_3d_data_averaging( mode, variable )
1368
1369
1370	CHARACTER (LEN=*), INTENT(IN) :: mode !< averaging interface mode
1371	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
1372
1373
1374	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 3d data averaging', 'start' )
1375
1376	CALL dynamics_3d_data_averaging( mode, variable )
1377	CALL tcm_3d_data_averaging( mode, variable )
1378
1379	IF ( biometeorology ) CALL bio_3d_data_averaging( mode, variable )
1380	IF ( bulk_cloud_model ) CALL bcm_3d_data_averaging( mode, variable )
1381	IF ( air_chemistry ) CALL chem_3d_data_averaging( mode, variable )
1382	CALL doq_3d_data_averaging( mode, variable ) ! ToDo: this seems to be not according to the design
1383	IF ( gust_module_enabled ) CALL gust_3d_data_averaging( mode, variable )
1384	IF ( land_surface ) CALL lsm_3d_data_averaging( mode, variable )
1385	IF ( ocean_mode ) CALL ocean_3d_data_averaging( mode, variable )
1386	IF ( plant_canopy ) CALL pcm_3d_data_averaging( mode, variable )
1387	IF ( radiation ) CALL radiation_3d_data_averaging( mode, variable )
1388	IF ( salsa ) CALL salsa_3d_data_averaging( mode, variable )
1389	IF ( urban_surface ) CALL usm_3d_data_averaging( mode, variable )
1390
1391	IF ( user_module_enabled ) CALL user_3d_data_averaging( mode, variable )
1392
1393	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 3d data averaging', 'end' )
1394
1395
1396	END SUBROUTINE module_interface_3d_data_averaging
1397
1398	!------------------------------------------------------------------------------!
1399	!
1400	! Description:
1401	! ------------
1402	!> Define module-specific 2D output variables
1403	!------------------------------------------------------------------------------!
1404	SUBROUTINE module_interface_data_output_2d( av, variable, found, grid, mode, &
1405	local_pf, two_d, nzb_do, nzt_do, &
1406	fill_value )
1407
1408	INTEGER(iwp), INTENT(IN) :: av !< flag for (non-)average output
1409	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
1410	LOGICAL, INTENT(INOUT) :: found !< flag if output variable is found
1411	CHARACTER (LEN=*), INTENT(INOUT) :: grid !< name of vertical grid
1412	CHARACTER (LEN=*), INTENT(IN) :: mode !< either 'xy', 'xz' or 'yz'
1413	LOGICAL, INTENT(OUT) :: two_d !< flag for 2D variables
1414	INTEGER(iwp), INTENT(IN) :: nzb_do !< vertical output index (bottom) (usually 0)
1415	INTEGER(iwp), INTENT(IN) :: nzt_do !< vertical output index (top) (usually nz_do3d)
1416	REAL(wp), INTENT(IN) :: fill_value !< to be removed
1417
1418	REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) :: local_pf !< ToDo: can also be kind=sp
1419
1420
1421	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 2d data output', 'start' )
1422
1423	CALL dynamics_data_output_2d( &
1424	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
1425	)
1426
1427	IF ( .NOT. found ) THEN
1428	CALL tcm_data_output_2d( &
1429	av, variable, found, grid, mode, local_pf, nzb_do, nzt_do &
1430	)
1431	ENDIF
1432
1433	IF ( .NOT. found .AND. biometeorology ) THEN
1434	CALL bio_data_output_2d( &
1435	av, variable, found, grid, local_pf, two_d, nzb_do, nzt_do &
1436	)
1437	ENDIF
1438
1439	IF ( .NOT. found .AND. bulk_cloud_model ) THEN
1440	CALL bcm_data_output_2d( &
1441	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1442	)
1443	ENDIF
1444
1445	IF ( .NOT. found .AND. air_chemistry ) THEN
1446	CALL chem_data_output_2d( &
1447	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
1448	)
1449	ENDIF
1450
1451	IF ( .NOT. found ) THEN
1452	CALL doq_output_2d( &
1453	av, variable, found, grid, mode, local_pf, two_d, &
1454	nzb_do, nzt_do, fill_value )
1455	ENDIF
1456
1457	IF ( .NOT. found .AND. gust_module_enabled ) THEN
1458	CALL gust_data_output_2d( &
1459	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
1460	)
1461	ENDIF
1462
1463	IF ( .NOT. found .AND. land_surface ) THEN
1464	CALL lsm_data_output_2d( &
1465	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1466	)
1467	ENDIF
1468
1469	IF ( .NOT. found .AND. ocean_mode ) THEN
1470	CALL ocean_data_output_2d( &
1471	av, variable, found, grid, mode, local_pf, nzb_do, nzt_do &
1472	)
1473	ENDIF
1474
1475	IF ( .NOT. found .AND. radiation ) THEN
1476	CALL radiation_data_output_2d( &
1477	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1478	)
1479	ENDIF
1480
1481	IF ( .NOT. found .AND. salsa ) THEN
1482	CALL salsa_data_output_2d( &
1483	av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1484	)
1485	ENDIF
1486
1487	IF ( .NOT. found .AND. user_module_enabled ) THEN
1488	CALL user_data_output_2d( &
1489	av, variable, found, grid, local_pf, two_d, nzb_do, nzt_do &
1490	)
1491	ENDIF
1492
1493	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 2d data output', 'end' )
1494
1495
1496	END SUBROUTINE module_interface_data_output_2d
1497
1498
1499	!------------------------------------------------------------------------------!
1500	!
1501	! Description:
1502	! ------------
1503	!> Define module-specific 3D output variables
1504	!------------------------------------------------------------------------------!
1505	SUBROUTINE module_interface_data_output_3d( av, variable, found, local_pf, &
1506	fill_value, resorted, nzb_do, nzt_do )
1507
1508
1509	INTEGER(iwp), INTENT(IN) :: av !< flag for (non-)average output
1510	CHARACTER (LEN=*), INTENT(IN) :: variable !< variable name
1511	LOGICAL, INTENT(INOUT) :: found !< flag if output variable is found
1512	REAL(wp), INTENT(IN) :: fill_value !< ToDo: refactor
1513	LOGICAL, INTENT(OUT) :: resorted !< flag if output has been resorted
1514	INTEGER(iwp), INTENT(IN) :: nzb_do !< vertical output index (bottom) (usually 0)
1515	INTEGER(iwp), INTENT(IN) :: nzt_do !< vertical output index (top) (usually nz_do3d)
1516
1517	REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) :: local_pf
1518
1519
1520	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 3d data output', 'start' )
1521
1522	CALL dynamics_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1523	resorted = .FALSE.
1524
1525	IF ( .NOT. found ) THEN
1526	CALL tcm_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1527	resorted = .TRUE.
1528	ENDIF
1529
1530	IF ( .NOT. found .AND. biometeorology ) THEN
1531	CALL bio_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1532	resorted = .FALSE.
1533	ENDIF
1534
1535	IF ( .NOT. found .AND. bulk_cloud_model ) THEN
1536	CALL bcm_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1537	resorted = .TRUE.
1538	ENDIF
1539
1540	IF ( .NOT. found .AND. air_chemistry ) THEN
1541	CALL chem_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1542	resorted = .TRUE.
1543	ENDIF
1544
1545	IF ( .NOT. found ) THEN
1546	CALL doq_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1547	resorted = .TRUE.
1548	ENDIF
1549
1550	IF ( .NOT. found .AND. gust_module_enabled ) THEN
1551	CALL gust_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1552	resorted = .TRUE.
1553	ENDIF
1554
1555	IF ( .NOT. found .AND. indoor_model ) THEN
1556	CALL im_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1557	resorted = .TRUE.
1558	ENDIF
1559
1560	IF ( .NOT. found .AND. ocean_mode ) THEN
1561	CALL ocean_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1562	resorted = .TRUE.
1563	ENDIF
1564
1565	IF ( .NOT. found .AND. plant_canopy ) THEN
1566	CALL pcm_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1567	resorted = .TRUE.
1568	ENDIF
1569
1570	IF ( .NOT. found .AND. radiation ) THEN
1571	CALL radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1572	resorted = .TRUE.
1573	ENDIF
1574
1575	IF ( .NOT. found .AND. salsa ) THEN
1576	CALL salsa_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1577	resorted = .TRUE.
1578	ENDIF
1579
1580	IF ( .NOT. found .AND. user_module_enabled ) THEN
1581	CALL user_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1582	resorted = .TRUE.
1583	ENDIF
1584
1585	IF ( debug_output_timestep ) CALL debug_message( 'module-specific 3d data output', 'end' )
1586
1587
1588	END SUBROUTINE module_interface_data_output_3d
1589
1590
1591	!------------------------------------------------------------------------------!
1592	! Description:
1593	! ------------
1594	!> Compute module-specific profile and timeseries data
1595	!------------------------------------------------------------------------------!
1596	SUBROUTINE module_interface_statistics( mode, sr, tn, dots_max )
1597
1598
1599	CHARACTER (LEN=*), INTENT(IN) :: mode !< statistical analysis mode
1600	INTEGER(iwp), INTENT(IN) :: sr !<
1601	INTEGER(iwp), INTENT(IN) :: tn !<
1602	INTEGER(iwp), INTENT(IN) :: dots_max !< maximum number of timeseries
1603
1604
1605	IF ( debug_output_timestep ) CALL debug_message( 'module-specific statistics', 'start' )
1606
1607	CALL dynamics_statistics( mode, sr, tn )
1608
1609	IF ( gust_module_enabled ) CALL gust_statistics( mode, sr, tn, dots_max )
1610	IF ( air_chemistry ) CALL chem_statistics( mode, sr, tn )
1611	IF ( salsa ) CALL salsa_statistics( mode, sr, tn )
1612
1613	IF ( user_module_enabled ) CALL user_statistics( mode, sr, tn )
1614
1615	IF ( debug_output_timestep ) CALL debug_message( 'module-specific statistics', 'end' )
1616
1617
1618	END SUBROUTINE module_interface_statistics
1619
1620
1621	!------------------------------------------------------------------------------!
1622	! Description:
1623	! ------------
1624	!> Read module-specific restart data globaly shared by all MPI ranks
1625	!------------------------------------------------------------------------------!
1626	SUBROUTINE module_interface_rrd_global( found )
1627
1628
1629	LOGICAL, INTENT(INOUT) :: found !< flag if variable was found
1630
1631
1632	IF ( debug_output ) CALL debug_message( 'module-specific read global restart data', 'start' )
1633
1634	CALL dynamics_rrd_global( found ) ! ToDo: change interface to pass variable
1635
1636	IF ( .NOT. found ) CALL bio_rrd_global( found ) ! ToDo: change interface to pass variable
1637	IF ( .NOT. found ) CALL bcm_rrd_global( found ) ! ToDo: change interface to pass variable
1638	IF ( .NOT. found ) CALL flight_rrd_global( found ) ! ToDo: change interface to pass variable
1639	IF ( .NOT. found ) CALL gust_rrd_global( found ) ! ToDo: change interface to pass variable
1640	IF ( .NOT. found ) CALL lpm_rrd_global( found ) ! ToDo: change interface to pass variable
1641	IF ( .NOT. found ) CALL ocean_rrd_global( found ) ! ToDo: change interface to pass variable
1642	IF ( .NOT. found ) CALL stg_rrd_global ( found ) ! ToDo: change interface to pass variable
1643	IF ( .NOT. found ) CALL wtm_rrd_global( found ) ! ToDo: change interface to pass variable
1644	IF ( .NOT. found ) CALL surface_data_output_rrd_global( found )
1645
1646	IF ( .NOT. found ) CALL user_rrd_global( found ) ! ToDo: change interface to pass variable
1647
1648	IF ( debug_output ) CALL debug_message( 'module-specific read global restart data', 'end' )
1649
1650
1651	END SUBROUTINE module_interface_rrd_global
1652
1653
1654	!------------------------------------------------------------------------------!
1655	! Description:
1656	! ------------
1657	!> Write module-specific restart data globaly shared by all MPI ranks
1658	!------------------------------------------------------------------------------!
1659	SUBROUTINE module_interface_wrd_global
1660
1661
1662	IF ( debug_output ) CALL debug_message( 'module-specific write global restart data', 'start' )
1663
1664	CALL dynamics_wrd_global
1665
1666	IF ( biometeorology ) CALL bio_wrd_global
1667	IF ( bulk_cloud_model ) CALL bcm_wrd_global
1668	IF ( virtual_flight ) CALL flight_wrd_global
1669	IF ( gust_module_enabled ) CALL gust_wrd_global
1670	IF ( ocean_mode ) CALL ocean_wrd_global
1671	IF ( syn_turb_gen ) CALL stg_wrd_global
1672	IF ( wind_turbine ) CALL wtm_wrd_global
1673	IF ( surface_output ) CALL surface_data_output_wrd_global
1674
1675	IF ( user_module_enabled ) CALL user_wrd_global
1676
1677	IF ( debug_output ) CALL debug_message( 'module-specific write global restart data', 'end' )
1678
1679
1680	END SUBROUTINE module_interface_wrd_global
1681
1682
1683	!------------------------------------------------------------------------------!
1684	! Description:
1685	! ------------
1686	!> Read module-specific restart data specific to local MPI ranks
1687	!------------------------------------------------------------------------------!
1688	SUBROUTINE module_interface_rrd_local( map_index, &
1689	nxlf, nxlc, nxl_on_file, &
1690	nxrf, nxrc, nxr_on_file, &
1691	nynf, nync, nyn_on_file, &
1692	nysf, nysc, nys_on_file, &
1693	tmp_2d, tmp_3d, found )
1694
1695
1696	INTEGER(iwp), INTENT(IN) :: map_index !<
1697	INTEGER(iwp), INTENT(IN) :: nxlc !<
1698	INTEGER(iwp), INTENT(IN) :: nxlf !<
1699	INTEGER(iwp), INTENT(IN) :: nxl_on_file !<
1700	INTEGER(iwp), INTENT(IN) :: nxrc !<
1701	INTEGER(iwp), INTENT(IN) :: nxrf !<
1702	INTEGER(iwp), INTENT(IN) :: nxr_on_file !<
1703	INTEGER(iwp), INTENT(IN) :: nync !<
1704	INTEGER(iwp), INTENT(IN) :: nynf !<
1705	INTEGER(iwp), INTENT(IN) :: nyn_on_file !<
1706	INTEGER(iwp), INTENT(IN) :: nysc !<
1707	INTEGER(iwp), INTENT(IN) :: nysf !<
1708	INTEGER(iwp), INTENT(IN) :: nys_on_file !<
1709	LOGICAL, INTENT(INOUT) :: found !< flag if variable was found
1710
1711	REAL(wp), &
1712	DIMENSION(nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp), &
1713	INTENT(OUT) :: tmp_2d !<
1714	REAL(wp), &
1715	DIMENSION(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp), &
1716	INTENT(OUT) :: tmp_3d !<
1717
1718
1719	IF ( debug_output ) CALL debug_message( 'module-specific read local restart data', 'start' )
1720
1721	CALL dynamics_rrd_local( &
1722	map_index, &
1723	nxlf, nxlc, nxl_on_file, &
1724	nxrf, nxrc, nxr_on_file, &
1725	nynf, nync, nyn_on_file, &
1726	nysf, nysc, nys_on_file, &
1727	tmp_2d, tmp_3d, found &
1728	) ! ToDo: change interface to pass variable
1729
1730	IF ( .NOT. found ) CALL bio_rrd_local( &
1731	found &
1732	)
1733
1734	IF ( .NOT. found ) CALL bcm_rrd_local( &
1735	map_index, &
1736	nxlf, nxlc, nxl_on_file, &
1737	nxrf, nxrc, nxr_on_file, &
1738	nynf, nync, nyn_on_file, &
1739	nysf, nysc, nys_on_file, &
1740	tmp_2d, tmp_3d, found &
1741	) ! ToDo: change interface to pass variable
1742
1743	IF ( .NOT. found ) CALL chem_rrd_local( &
1744	map_index, &
1745	nxlf, nxlc, nxl_on_file, &
1746	nxrf, nxrc, nxr_on_file, &
1747	nynf, nync, nyn_on_file, &
1748	nysf, nysc, nys_on_file, &
1749	tmp_3d, found &
1750	) ! ToDo: change interface to pass variable
1751
1752	! IF ( .NOT. found ) CALL doq_rrd_local( &
1753	! map_index, &
1754	! nxlf, nxlc, nxl_on_file, &
1755	! nxrf, nxrc, nxr_on_file, &
1756	! nynf, nync, nyn_on_file, &
1757	! nysf, nysc, nys_on_file, &
1758	! tmp_3d_non_standard, found &
1759	! ) ! ToDo: change interface to pass variable CALL doq_wrd_local
1760
1761	IF ( .NOT. found ) CALL gust_rrd_local( &
1762	map_index, &
1763	nxlf, nxlc, nxl_on_file, &
1764	nxrf, nxrc, nxr_on_file, &
1765	nynf, nync, nyn_on_file, &
1766	nysf, nysc, nys_on_file, &
1767	tmp_2d, tmp_3d, found &
1768	) ! ToDo: change interface to pass variable
1769
1770	IF ( .NOT. found ) CALL lpm_rrd_local( &
1771	map_index, &
1772	nxlf, nxlc, nxl_on_file, &
1773	nxrf, nxrc, nxr_on_file, &
1774	nynf, nync, nyn_on_file, &
1775	nysf, nysc, nys_on_file, &
1776	tmp_3d, found &
1777	) ! ToDo: change interface to pass variable
1778
1779	IF ( .NOT. found ) CALL lsm_rrd_local( &
1780	map_index, &
1781	nxlf, nxlc, nxl_on_file, &
1782	nxrf, nxrc, nxr_on_file, &
1783	nynf, nync, nyn_on_file, &
1784	nysf, nysc, nys_on_file, &
1785	tmp_2d, found &
1786	) ! ToDo: change interface to pass variable
1787
1788	IF ( .NOT. found ) CALL pcm_rrd_local( &
1789	map_index, &
1790	nxlf, nxlc, nxl_on_file, &
1791	nxrf, nxrc, nxr_on_file, &
1792	nynf, nync, nyn_on_file, &
1793	nysf, nysc, nys_on_file, &
1794	found &
1795	)
1796
1797	IF ( .NOT. found ) CALL ocean_rrd_local( &
1798	map_index, &
1799	nxlf, nxlc, nxl_on_file, &
1800	nxrf, nxrc, nxr_on_file, &
1801	nynf, nync, nyn_on_file, &
1802	nysf, nysc, nys_on_file, &
1803	tmp_3d, found &
1804	) ! ToDo: change interface to pass variable
1805
1806	IF ( .NOT. found ) CALL radiation_rrd_local( &
1807	map_index, &
1808	nxlf, nxlc, nxl_on_file, &
1809	nxrf, nxrc, nxr_on_file, &
1810	nynf, nync, nyn_on_file, &
1811	nysf, nysc, nys_on_file, &
1812	tmp_2d, tmp_3d, found &
1813	) ! ToDo: change interface to pass variable
1814
1815	IF ( .NOT. found ) CALL salsa_rrd_local( &
1816	map_index, &
1817	nxlf, nxlc, nxl_on_file, &
1818	nxrf, nxrc, nxr_on_file, &
1819	nynf, nync, nyn_on_file, &
1820	nysf, nysc, nys_on_file, &
1821	tmp_3d, found &
1822	) ! ToDo: change interface to pass variable
1823
1824	IF ( .NOT. found ) CALL usm_rrd_local( &
1825	map_index, &
1826	nxlf, nxlc, nxl_on_file, &
1827	nxrf, nxr_on_file, &
1828	nynf, nyn_on_file, &
1829	nysf, nysc, nys_on_file, &
1830	found &
1831	) ! ToDo: change interface to pass variable
1832	!
1833	!-- Surface data do not need overlap data, so do not pass these information.
1834	IF ( .NOT. found ) CALL surface_data_output_rrd_local( found )
1835
1836	IF ( .NOT. found ) CALL user_rrd_local( &
1837	map_index, &
1838	nxlf, nxlc, nxl_on_file, &
1839	nxrf, nxrc, nxr_on_file, &
1840	nynf, nync, nyn_on_file, &
1841	nysf, nysc, nys_on_file, &
1842	tmp_3d, found &
1843	) ! ToDo: change interface to pass variable
1844
1845	IF ( debug_output ) CALL debug_message( 'module-specific read local restart data', 'end' )
1846
1847
1848	END SUBROUTINE module_interface_rrd_local
1849
1850
1851	!------------------------------------------------------------------------------!
1852	! Description:
1853	! ------------
1854	!> Write module-specific restart data specific to local MPI ranks
1855	!------------------------------------------------------------------------------!
1856	SUBROUTINE module_interface_wrd_local
1857
1858
1859	IF ( debug_output ) CALL debug_message( 'module-specific write local restart data', 'start' )
1860
1861	CALL dynamics_wrd_local
1862
1863	IF ( biometeorology ) CALL bio_wrd_local
1864	IF ( bulk_cloud_model ) CALL bcm_wrd_local
1865	IF ( air_chemistry ) CALL chem_wrd_local
1866	CALL doq_wrd_local
1867	IF ( gust_module_enabled ) CALL gust_wrd_local
1868	IF ( particle_advection ) CALL lpm_wrd_local
1869	IF ( land_surface ) CALL lsm_wrd_local
1870	IF ( plant_canopy ) CALL pcm_wrd_local
1871	IF ( ocean_mode ) CALL ocean_wrd_local
1872	IF ( radiation ) CALL radiation_wrd_local
1873	IF ( salsa ) CALL salsa_wrd_local
1874	IF ( urban_surface ) CALL usm_wrd_local
1875	IF ( surface_output ) CALL surface_data_output_wrd_local
1876
1877	IF ( user_module_enabled ) CALL user_wrd_local
1878
1879	IF ( debug_output ) CALL debug_message( 'module-specific write local restart data', 'end' )
1880
1881
1882	END SUBROUTINE module_interface_wrd_local
1883
1884
1885	!------------------------------------------------------------------------------!
1886	! Description:
1887	! ------------
1888	!> Perform module-specific last actions before the program terminates
1889	!------------------------------------------------------------------------------!
1890	SUBROUTINE module_interface_last_actions
1891
1892
1893	IF ( debug_output ) CALL debug_message( 'module-specific last actions', 'start' )
1894
1895	CALL dynamics_last_actions
1896
1897	IF ( virtual_measurement ) CALL vm_last_actions
1898
1899	IF ( user_module_enabled ) CALL user_last_actions
1900
1901	IF ( debug_output ) CALL debug_message( 'module-specific last actions', 'end' )
1902
1903
1904	END SUBROUTINE module_interface_last_actions
1905
1906
1907	END MODULE module_interface

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |