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

Last change on this file since 4017 was 4017, checked in by schwenkel, 5 years ago

Modularization of all lagrangian particle model code components

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