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

Last change on this file since 4015 was 4008, checked in by moh.hefny, 5 years ago

Bugfix in check radiation related variables

  • Property svn:keywords set to Id
File size: 73.0 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 4008 2019-05-30 09:50:11Z raasch $
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 land_surface_model_mod,                                                &
246        ONLY:  lsm_parin,                                                      &
247               lsm_check_parameters,                                           &
248               lsm_check_data_output_pr,                                       &
249               lsm_check_data_output,                                          &
250               lsm_init_arrays,                                                &
251               lsm_init,                                                       &
252               lsm_header,                                                     &
253               lsm_swap_timelevel,                                             &
254               lsm_3d_data_averaging,                                          &
255               lsm_data_output_2d,                                             &
256               lsm_rrd_local,                                                  &
257               lsm_wrd_local
258
259    USE lsf_nudging_mod,                                                       &
260        ONLY:  lsf_nudging_check_parameters,                                   &
261               lsf_nudging_check_data_output_pr,                               &
262               lsf_init,                                                       &
263               nudge_init,                                                     &
264               lsf_nudging_header
265
266    USE multi_agent_system_mod,                                                &
267        ONLY:  mas_parin
268
269    USE nesting_offl_mod,                                                      &
270        ONLY:  nesting_offl_parin,                                             &
271               nesting_offl_check_parameters,                                  &
272               nesting_offl_header
273
274    USE ocean_mod,                                                             &
275        ONLY:  ocean_parin,                                                    &
276               ocean_check_parameters,                                         &
277               ocean_check_data_output_pr,                                     &
278               ocean_check_data_output,                                        &
279               ocean_init_arrays,                                              &
280               ocean_init,                                                     &
281               ocean_header,                                                   &
282               ocean_actions,                                                  &
283               ocean_prognostic_equations,                                     &
284               ocean_swap_timelevel,                                           &
285               ocean_3d_data_averaging,                                        &
286               ocean_data_output_2d,                                           &
287               ocean_data_output_3d,                                           &
288               ocean_rrd_global,                                               &
289               ocean_wrd_global,                                               &
290               ocean_rrd_local,                                                &
291               ocean_wrd_local
292
293    USE plant_canopy_model_mod,                                                &
294         ONLY: pcm_parin,                                                      &
295               pcm_check_parameters,                                           &
296               pcm_check_data_output,                                          &
297               pcm_init,                                                       &
298               pcm_header,                                                     &
299               pcm_data_output_3d
300
301    USE radiation_model_mod,                                                   &
302        ONLY:  radiation,                                                      &
303               radiation_parin,                                                &
304               radiation_check_parameters,                                     &
305               radiation_check_data_output_ts,                                 &
306               radiation_check_data_output_pr,                                 &
307               radiation_check_data_output,                                    &
308               radiation_init,                                                 &
309               radiation_header,                                               &
310               radiation_3d_data_averaging,                                    &
311               radiation_data_output_2d,                                       &
312               radiation_data_output_3d,                                       &
313               radiation_rrd_local,                                            &
314               radiation_wrd_local
315
316    USE salsa_mod,                                                             &
317        ONLY:  salsa_parin,                                                    &
318               salsa_check_parameters,                                         &
319               salsa_check_data_output,                                        &
320               salsa_init_arrays,                                              &
321               salsa_init,                                                     &
322               salsa_header,                                                   &
323               salsa_actions,                                                  &
324               salsa_non_advective_processes,                                  &
325               salsa_exchange_horiz_bounds,                                    &
326               salsa_prognostic_equations,                                     &
327               salsa_swap_timelevel,                                           &
328               salsa_3d_data_averaging,                                        &
329               salsa_data_output_2d,                                           &
330               salsa_data_output_3d,                                           &
331               salsa_rrd_local,                                                &
332               salsa_wrd_local
333
334    USE spectra_mod,                                                           &
335        ONLY:  calculate_spectra,                                              &
336               spectra_parin,                                                  &
337               spectra_check_parameters,                                       &
338               spectra_header
339
340    USE surface_data_output_mod,                                               &
341        ONLY:  surface_data_output_parin,                                      &
342               surface_data_output_check_parameters,                           &
343               surface_data_output_init_arrays,                                &
344               surface_data_output_rrd_local,                                  &
345               surface_data_output_rrd_global,                                 &
346               surface_data_output_wrd_local,                                  &
347               surface_data_output_wrd_global
348
349    USE synthetic_turbulence_generator_mod,                                    &
350        ONLY:  stg_parin,                                                      &
351               stg_check_parameters,                                           &
352               stg_header,                                                     &
353               stg_rrd_global,                                                 &
354               stg_wrd_global
355
356    USE urban_surface_mod,                                                     &
357        ONLY:  usm_parin,                                                      &
358               usm_check_parameters,                                           &
359               usm_check_data_output,                                          &
360               usm_init_arrays,                                                &
361               usm_init,                                                       &
362               usm_swap_timelevel,                                             &
363               usm_3d_data_averaging,                                          &
364               usm_rrd_local,                                                  &
365               usm_wrd_local
366
367    USE user,                                                                  &
368        ONLY:  user_module_enabled,                                            &
369               user_parin,                                                     &
370               user_check_parameters,                                          &
371               user_check_data_output_ts,                                      &
372               user_check_data_output_pr,                                      &
373               user_check_data_output,                                         &
374               user_init,                                                      &
375               user_init_arrays,                                               &
376               user_header,                                                    &
377               user_actions,                                                   &
378               user_3d_data_averaging,                                         &
379               user_data_output_2d,                                            &
380               user_data_output_3d,                                            &
381               user_statistics,                                                &
382               user_rrd_global,                                                &
383               user_rrd_local,                                                 &
384               user_wrd_global,                                                &
385               user_wrd_local,                                                 &
386               user_last_actions
387
388    USE virtual_measurement_mod,                                               &
389        ONLY:  vm_check_parameters,                                            &
390               vm_init,                                                        &
391               vm_last_actions,                                                &
392               vm_parin
393
394    USE wind_turbine_model_mod,                                                &
395        ONLY:  wtm_parin,                                                      &
396               wtm_check_parameters,                                           &
397               wtm_init_arrays,                                                &
398               wtm_init,                                                       &
399               wtm_actions,                                                    &
400               wtm_rrd_global,                                                 &
401               wtm_wrd_global
402
403    IMPLICIT NONE
404
405    PRIVATE
406
407!
408!-- Public functions
409    PUBLIC                                                                     &
410       module_interface_parin,                                                 &
411       module_interface_check_parameters,                                      &
412       module_interface_check_data_output_ts,                                  &
413       module_interface_check_data_output_pr,                                  &
414       module_interface_check_data_output,                                     &
415       module_interface_init_masks,                                            &
416       module_interface_define_netcdf_grid,                                    &
417       module_interface_init_arrays,                                           &
418       module_interface_init,                                                  &
419       module_interface_init_checks,                                           &
420       module_interface_header,                                                &
421       module_interface_actions,                                               &
422       module_interface_non_advective_processes,                               &
423       module_interface_exchange_horiz,                                        &
424       module_interface_prognostic_equations,                                  &
425       module_interface_swap_timelevel,                                        &
426       module_interface_3d_data_averaging,                                     &
427       module_interface_data_output_2d,                                        &
428       module_interface_data_output_3d,                                        &
429       module_interface_statistics,                                            &
430       module_interface_rrd_global,                                            &
431       module_interface_wrd_global,                                            &
432       module_interface_rrd_local,                                             &
433       module_interface_wrd_local,                                             &
434       module_interface_last_actions
435
436
437    INTERFACE module_interface_parin
438       MODULE PROCEDURE module_interface_parin
439    END INTERFACE module_interface_parin
440
441    INTERFACE module_interface_check_parameters
442       MODULE PROCEDURE module_interface_check_parameters
443    END INTERFACE module_interface_check_parameters
444
445    INTERFACE module_interface_check_data_output_ts
446       MODULE PROCEDURE module_interface_check_data_output_ts
447    END INTERFACE module_interface_check_data_output_ts
448
449    INTERFACE module_interface_check_data_output_pr
450       MODULE PROCEDURE module_interface_check_data_output_pr
451    END INTERFACE module_interface_check_data_output_pr
452
453    INTERFACE module_interface_check_data_output
454       MODULE PROCEDURE module_interface_check_data_output
455    END INTERFACE module_interface_check_data_output
456
457    INTERFACE module_interface_init_masks
458       MODULE PROCEDURE module_interface_init_masks
459    END INTERFACE module_interface_init_masks
460
461    INTERFACE module_interface_define_netcdf_grid
462       MODULE PROCEDURE module_interface_define_netcdf_grid
463    END INTERFACE module_interface_define_netcdf_grid
464
465    INTERFACE module_interface_init_arrays
466       MODULE PROCEDURE module_interface_init_arrays
467    END INTERFACE module_interface_init_arrays
468
469    INTERFACE module_interface_init
470       MODULE PROCEDURE module_interface_init
471    END INTERFACE module_interface_init
472
473    INTERFACE module_interface_init_checks
474       MODULE PROCEDURE module_interface_init_checks
475    END INTERFACE module_interface_init_checks
476
477    INTERFACE module_interface_header
478       MODULE PROCEDURE module_interface_header
479    END INTERFACE module_interface_header
480
481    INTERFACE module_interface_actions
482       MODULE PROCEDURE module_interface_actions
483       MODULE PROCEDURE module_interface_actions_ij
484    END INTERFACE module_interface_actions
485
486    INTERFACE module_interface_non_advective_processes
487       MODULE PROCEDURE module_interface_non_advective_processes
488       MODULE PROCEDURE module_interface_non_advective_processes_ij
489    END INTERFACE module_interface_non_advective_processes
490   
491    INTERFACE module_interface_exchange_horiz
492       MODULE PROCEDURE module_interface_exchange_horiz
493    END INTERFACE module_interface_exchange_horiz
494   
495    INTERFACE module_interface_prognostic_equations
496       MODULE PROCEDURE module_interface_prognostic_equations
497       MODULE PROCEDURE module_interface_prognostic_equations_ij
498    END INTERFACE module_interface_prognostic_equations
499
500    INTERFACE module_interface_swap_timelevel
501       MODULE PROCEDURE module_interface_swap_timelevel
502    END INTERFACE module_interface_swap_timelevel
503
504    INTERFACE module_interface_3d_data_averaging
505       MODULE PROCEDURE module_interface_3d_data_averaging
506    END INTERFACE module_interface_3d_data_averaging
507
508    INTERFACE module_interface_data_output_2d
509       MODULE PROCEDURE module_interface_data_output_2d
510    END INTERFACE module_interface_data_output_2d
511
512    INTERFACE module_interface_data_output_3d
513       MODULE PROCEDURE module_interface_data_output_3d
514    END INTERFACE module_interface_data_output_3d
515
516    INTERFACE module_interface_statistics
517       MODULE PROCEDURE module_interface_statistics
518    END INTERFACE module_interface_statistics
519
520    INTERFACE module_interface_rrd_global
521       MODULE PROCEDURE module_interface_rrd_global
522    END INTERFACE module_interface_rrd_global
523
524    INTERFACE module_interface_wrd_global
525       MODULE PROCEDURE module_interface_wrd_global
526    END INTERFACE module_interface_wrd_global
527
528    INTERFACE module_interface_rrd_local
529       MODULE PROCEDURE module_interface_rrd_local
530    END INTERFACE module_interface_rrd_local
531
532    INTERFACE module_interface_wrd_local
533       MODULE PROCEDURE module_interface_wrd_local
534    END INTERFACE module_interface_wrd_local
535
536    INTERFACE module_interface_last_actions
537       MODULE PROCEDURE module_interface_last_actions
538    END INTERFACE module_interface_last_actions
539
540
541 CONTAINS
542
543
544!------------------------------------------------------------------------------!
545! Description:
546! ------------
547!> Read module-specific parameter namelists
548!------------------------------------------------------------------------------!
549 SUBROUTINE module_interface_parin
550
551
552    IF ( debug_output )  CALL debug_message( 'reading module-specific parameters', 'start' )
553
554    CALL bio_parin
555    CALL bcm_parin
556    CALL chem_parin
557    CALL flight_parin ! ToDo: rename module to match filename
558    CALL gust_parin
559    CALL im_parin
560    CALL lsm_parin
561    ! ToDo: create parin routine for large_scale_forcing and nudging (should be seperate modules or new module switch)
562    CALL mas_parin
563    CALL nesting_offl_parin
564    CALL ocean_parin
565    CALL pcm_parin
566    CALL radiation_parin
567    CALL salsa_parin
568    CALL spectra_parin
569    CALL surface_data_output_parin
570    CALL stg_parin
571    CALL user_parin ! ToDo: make user code a single Fortran module
572    CALL usm_parin
573    CALL vm_parin
574    CALL wtm_parin
575
576    CALL package_parin ! ToDo: deprecated, needs to be dissolved
577
578    IF ( debug_output )  CALL debug_message( 'reading module-specific parameters', 'end' )
579
580
581 END SUBROUTINE module_interface_parin
582
583
584!------------------------------------------------------------------------------!
585! Description:
586! ------------
587!> Perform module-specific initialization checks
588!------------------------------------------------------------------------------!
589 SUBROUTINE module_interface_check_parameters
590
591
592    IF ( debug_output )  CALL debug_message( 'checking module-specific parameters', 'start' )
593
594    IF ( bulk_cloud_model )     CALL bcm_check_parameters
595    IF ( air_chemistry )        CALL chem_check_parameters
596    IF ( gust_module_enabled )  CALL gust_check_parameters
597    IF ( indoor_model )         CALL im_check_parameters
598    IF ( land_surface )         CALL lsm_check_parameters
599    IF ( large_scale_forcing  .OR.  nudging )  CALL lsf_nudging_check_parameters ! ToDo: create single module switch
600    IF ( nesting_offline )      CALL nesting_offl_check_parameters
601    IF ( ocean_mode )           CALL ocean_check_parameters
602    IF ( plant_canopy )         CALL pcm_check_parameters
603    IF ( radiation )            CALL radiation_check_parameters
604    IF ( salsa )                CALL salsa_check_parameters
605    IF ( calculate_spectra )    CALL spectra_check_parameters
606    IF ( surface_output )       CALL surface_data_output_check_parameters
607    IF ( syn_turb_gen )         CALL stg_check_parameters
608    IF ( urban_surface )        CALL usm_check_parameters
609    IF ( virtual_measurement )  CALL vm_check_parameters
610    IF ( wind_turbine )         CALL wtm_check_parameters
611    IF ( user_module_enabled )  CALL user_check_parameters
612
613    IF ( debug_output )  CALL debug_message( 'checking module-specific parameters', 'end' )
614
615
616 END SUBROUTINE module_interface_check_parameters
617
618
619!------------------------------------------------------------------------------!
620! Description:
621! ------------
622!> Check module-specific data output of timeseries
623!------------------------------------------------------------------------------!
624 SUBROUTINE module_interface_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
625
626
627    INTEGER(iwp),      INTENT(IN)    ::  dots_max !< variable output array index
628    INTEGER(iwp),      INTENT(INOUT)    ::  dots_num !< variable output array index
629    CHARACTER (LEN=*), DIMENSION(dots_max), INTENT(INOUT) :: dots_label
630    CHARACTER (LEN=*), DIMENSION(dots_max), INTENT(INOUT) :: dots_unit
631
632
633    IF ( debug_output )  CALL debug_message( 'checking module-specific data output ts', 'start' )
634
635    IF ( radiation )  THEN
636       CALL radiation_check_data_output_ts( dots_max, dots_num )
637    ENDIF
638
639    IF ( user_module_enabled )  THEN
640       CALL user_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
641    ENDIF
642
643    IF ( debug_output )  CALL debug_message( 'checking module-specific data output ts', 'end' )
644
645
646 END SUBROUTINE module_interface_check_data_output_ts
647
648
649!------------------------------------------------------------------------------!
650! Description:
651! ------------
652!> Check module-specific data output of profiles
653!------------------------------------------------------------------------------!
654 SUBROUTINE module_interface_check_data_output_pr( variable, var_count, unit,  &
655                                                   dopr_unit )
656
657
658    CHARACTER (LEN=*), INTENT(IN)    ::  variable  !< variable name
659    INTEGER(iwp),      INTENT(IN)    ::  var_count !< variable output array index
660    CHARACTER (LEN=*), INTENT(INOUT) ::  unit      !< physical unit of variable
661    CHARACTER (LEN=*), INTENT(OUT)   ::  dopr_unit !< local value of dopr_unit
662
663
664    IF ( debug_output )  CALL debug_message( 'checking module-specific data output pr', 'start' )
665
666    IF ( unit == 'illegal' .AND.  bulk_cloud_model )  THEN
667       CALL bcm_check_data_output_pr( variable, var_count, unit, dopr_unit )
668    ENDIF
669
670    IF ( unit == 'illegal' .AND.  air_chemistry )  THEN
671       CALL chem_check_data_output_pr( variable, var_count, unit, dopr_unit )
672    ENDIF
673
674    IF ( unit == 'illegal'  .AND.  gust_module_enabled  )  THEN
675       CALL gust_check_data_output_pr( variable, var_count, unit, dopr_unit )
676    ENDIF
677
678    IF ( unit == 'illegal' )  THEN ! ToDo: add module switch if possible
679       CALL lsm_check_data_output_pr( variable, var_count, unit, dopr_unit )
680    ENDIF
681
682    IF ( unit == 'illegal' )  THEN ! ToDo: add module switch if possible
683       CALL lsf_nudging_check_data_output_pr( variable, var_count, unit, dopr_unit )
684    ENDIF
685
686    IF ( unit == 'illegal'  .AND.  ocean_mode )  THEN
687       CALL ocean_check_data_output_pr( variable, var_count, unit, dopr_unit )
688    ENDIF
689
690    IF ( unit == 'illegal'  .AND.  radiation )  THEN
691       CALL radiation_check_data_output_pr( variable, var_count, unit, dopr_unit )
692    ENDIF
693
694    IF ( unit == 'illegal'  .AND.  user_module_enabled )  THEN
695       unit = '' ! ToDo: Seems like a hack. Find a general soultion!
696       CALL user_check_data_output_pr( variable, var_count, unit, dopr_unit )
697    ENDIF
698
699    IF ( debug_output )  CALL debug_message( 'checking module-specific data output pr', 'end' )
700
701
702 END SUBROUTINE module_interface_check_data_output_pr
703
704!------------------------------------------------------------------------------!
705! Description:
706! ------------
707!> Check module-specific 2D and 3D data output
708!------------------------------------------------------------------------------!
709 SUBROUTINE module_interface_check_data_output( variable, unit, i, j, ilen, k )
710
711
712    CHARACTER (LEN=*), INTENT(IN)    ::  variable !< variable name
713    CHARACTER (LEN=*), INTENT(INOUT) ::  unit     !< physical unit of variable
714
715    INTEGER(iwp),      INTENT(IN)    :: i         !< ToDo: remove dummy argument, instead pass string from data_output
716    INTEGER(iwp),      INTENT(IN)    :: j         !< average quantity? 0 = no, 1 = yes
717    INTEGER(iwp),      INTENT(IN)    :: ilen      !< ToDo: remove dummy argument, instead pass string from data_output
718    INTEGER(iwp),      INTENT(IN)    :: k         !< ToDo: remove dummy argument, instead pass string from data_output
719
720
721    IF ( debug_output )  CALL debug_message( 'checking module-specific data output 2d/3d', 'start' )
722
723    IF ( unit == 'illegal'  .AND.  biometeorology )  THEN
724       CALL bio_check_data_output( variable, unit, i, j, ilen, k )
725    ENDIF
726
727    IF ( unit == 'illegal'  .AND.  bulk_cloud_model  )  THEN
728       CALL bcm_check_data_output( variable, unit )
729    ENDIF
730
731    IF ( unit == 'illegal'  .AND.  air_chemistry                               &
732         .AND.  (variable(1:3) == 'kc_' .OR. variable(1:3) == 'em_') )  THEN  ! ToDo: remove aditional conditions
733       CALL chem_check_data_output( variable, unit, i, ilen, k )
734    ENDIF
735
736    IF ( unit == 'illegal'  .AND.  gust_module_enabled  )  THEN
737       CALL gust_check_data_output( variable, unit )
738    ENDIF
739
740    IF ( unit == 'illegal' )  THEN  ! ToDo: add module switch if possible
741       CALL lsm_check_data_output( variable, unit, i, ilen, k )
742    ENDIF
743
744    IF ( unit == 'illegal'  .AND.  ocean_mode )  THEN
745       CALL ocean_check_data_output( variable, unit )
746    ENDIF
747
748    IF ( unit == 'illegal'  .AND.  plant_canopy                                &
749         .AND.  variable(1:4) == 'pcm_' )  THEN  ! ToDo: remove aditional conditions
750       CALL pcm_check_data_output( variable, unit )
751    ENDIF
752
753    IF ( unit == 'illegal'  .AND.  radiation )  THEN
754       CALL radiation_check_data_output( variable, unit, i, ilen, k )
755    ENDIF
756
757    IF ( unit == 'illegal' .AND. salsa ) THEN
758       CALL salsa_check_data_output( variable, unit )
759    ENDIF
760
761    IF ( unit == 'illegal' .AND. indoor_model ) THEN
762       CALL im_check_data_output( variable, unit )
763    ENDIF
764   
765    IF ( unit == 'illegal'  .AND.  urban_surface                      &
766        .AND.  variable(1:4) == 'usm_' )  THEN  ! ToDo: remove aditional conditions
767       CALL usm_check_data_output( variable, unit )
768    ENDIF
769
770    IF ( unit == 'illegal'  .AND.  user_module_enabled )  THEN
771       unit = ''
772       CALL user_check_data_output( variable, unit )
773    ENDIF
774
775    IF ( debug_output )  CALL debug_message( 'checking module-specific data output 2d/3d', 'end' )
776
777
778 END SUBROUTINE module_interface_check_data_output
779
780
781!------------------------------------------------------------------------------!
782!
783! Description:
784! ------------
785!> Interface for init_masks. ToDo: get rid of these redundant calls!
786!------------------------------------------------------------------------------!
787 SUBROUTINE module_interface_init_masks( variable, unit )
788
789
790    CHARACTER (LEN=*), INTENT(IN)    ::  variable !< variable name
791    CHARACTER (LEN=*), INTENT(INOUT) ::  unit     !< physical unit of variable
792
793
794    IF ( debug_output )  CALL debug_message( 'initializing module-specific masks', 'start' )
795
796    IF ( unit == 'illegal'  .AND.  air_chemistry                               &
797         .AND.  (variable(1:3) == 'kc_' .OR. variable(1:3) == 'em_') )  THEN  ! ToDo: remove aditional conditions
798       CALL chem_check_data_output( variable, unit, 0, 0, 0 )
799    ENDIF
800
801    IF ( unit == 'illegal'  .AND.  radiation )  THEN
802       CALL radiation_check_data_output( variable, unit, 0, 0, 0 )
803    ENDIF
804
805    IF ( unit == 'illegal'  .AND.  salsa )  THEN
806       CALL salsa_check_data_output( variable, unit )
807    ENDIF
808
809    IF ( unit == 'illegal'  .AND.  user_module_enabled )  THEN
810       unit = ''
811       CALL user_check_data_output( variable, unit )
812    ENDIF
813
814    IF ( debug_output )  CALL debug_message( 'initializing module-specific masks', 'end' )
815
816
817 END SUBROUTINE module_interface_init_masks
818
819
820!------------------------------------------------------------------------------!
821!
822! Description:
823! ------------
824!> Define appropriate grid for module-specific netcdf output variables.
825!------------------------------------------------------------------------------!
826 SUBROUTINE module_interface_define_netcdf_grid( var, found,                   &
827                                                 grid_x, grid_y, grid_z )
828
829
830    CHARACTER (LEN=*), INTENT(IN)  ::  var    !< variable name
831    LOGICAL,           INTENT(OUT) ::  found  !< indicates if variable was found
832    CHARACTER (LEN=*), INTENT(OUT) ::  grid_x !< netcdf dimension in x-direction
833    CHARACTER (LEN=*), INTENT(OUT) ::  grid_y !< netcdf dimension in y-direction
834    CHARACTER (LEN=*), INTENT(OUT) ::  grid_z !< netcdf dimension in z-direction
835
836
837    IF ( debug_output )  CALL debug_message( 'defining module-specific netcdf grids', 'start' )
838!
839!-- As long as no action is done in this subroutine, initialize strings with
840!-- intent(out) attribute, in order to avoid compiler warnings.
841    found  = .FALSE.
842    grid_x = 'none'
843    grid_y = 'none'
844    grid_z = 'none'
845!
846!-- Use var to avoid compiler warning about unused variable
847    IF ( var == ' ' )  RETURN
848
849    IF ( debug_output )  CALL debug_message( 'defining module-specific netcdf grids', 'end' )
850
851
852 END SUBROUTINE module_interface_define_netcdf_grid
853
854
855!------------------------------------------------------------------------------!
856! Description:
857! ------------
858!> Allocate module-specific arrays and pointers
859!------------------------------------------------------------------------------!
860 SUBROUTINE module_interface_init_arrays
861
862
863    IF ( debug_output )  CALL debug_message( 'initializing module-specific arrays', 'start' )
864
865    IF ( bulk_cloud_model    )  CALL bcm_init_arrays
866    IF ( air_chemistry       )  CALL chem_init_arrays
867    IF ( gust_module_enabled )  CALL gust_init_arrays
868    IF ( land_surface        )  CALL lsm_init_arrays
869    IF ( ocean_mode          )  CALL ocean_init_arrays
870    IF ( salsa               )  CALL salsa_init_arrays
871    IF ( urban_surface       )  CALL usm_init_arrays
872    IF ( surface_output      )  CALL surface_data_output_init_arrays
873    IF ( wind_turbine        )  CALL wtm_init_arrays
874    IF ( user_module_enabled )  CALL user_init_arrays
875
876    IF ( debug_output )  CALL debug_message( 'initializing module-specific arrays', 'end' )
877
878
879 END SUBROUTINE module_interface_init_arrays
880
881
882!------------------------------------------------------------------------------!
883! Description:
884! ------------
885!> Perform module-specific initialization
886!------------------------------------------------------------------------------!
887 SUBROUTINE module_interface_init
888
889
890    IF ( debug_output )  CALL debug_message( 'module-specific initialization', 'start' )
891
892    IF ( biometeorology      )  CALL bio_init
893    IF ( bulk_cloud_model    )  CALL bcm_init
894    IF ( air_chemistry       )  CALL chem_init
895    IF ( virtual_flight      )  CALL flight_init
896    IF ( gust_module_enabled )  CALL gust_init
897    IF ( indoor_model        )  CALL im_init
898    IF ( large_scale_forcing )  CALL lsf_init
899    IF ( land_surface        )  CALL lsm_init
900    IF ( nudging             )  CALL nudge_init
901    IF ( ocean_mode          )  CALL ocean_init
902    IF ( plant_canopy        )  CALL pcm_init
903    IF ( salsa               )  CALL salsa_init
904    IF ( urban_surface       )  CALL usm_init
905    IF ( virtual_measurement )  CALL vm_init
906    IF ( wind_turbine        )  CALL wtm_init
907    IF ( radiation           )  CALL radiation_init
908    IF ( user_module_enabled )  CALL user_init
909
910    IF ( debug_output )  CALL debug_message( 'module-specific initialization', 'end' )
911
912
913 END SUBROUTINE module_interface_init
914
915
916!------------------------------------------------------------------------------!
917! Description:
918! ------------
919!> Perform module-specific post-initialization checks
920!------------------------------------------------------------------------------!
921 SUBROUTINE module_interface_init_checks
922
923
924    IF ( debug_output )  CALL debug_message( 'module-specific post-initialization checks', 'start' )
925
926    IF ( biometeorology      )  CALL bio_init_checks
927
928    IF ( debug_output )  CALL debug_message( 'module-specific post-initialization checks', 'end' )
929
930
931 END SUBROUTINE module_interface_init_checks
932
933
934!------------------------------------------------------------------------------!
935! Description:
936! ------------
937!> Gather module-specific header output
938!------------------------------------------------------------------------------!
939 SUBROUTINE module_interface_header( io )
940
941
942    INTEGER(iwp), INTENT(IN) ::  io  !< unit of the output file
943
944
945    IF ( debug_output )  CALL debug_message( 'module-specific header output', 'start' )
946
947    IF ( biometeorology      )  CALL bio_header ( io )
948    IF ( bulk_cloud_model    )  CALL bcm_header( io )
949    IF ( air_chemistry       )  CALL chem_header ( io )
950    IF ( virtual_flight      )  CALL flight_header( io )
951    IF ( gust_module_enabled )  CALL gust_header( io )
952    IF ( land_surface        )  CALL lsm_header( io )
953    IF ( large_scale_forcing )  CALL lsf_nudging_header( io )
954    IF ( nesting_offline     )  CALL nesting_offl_header( io )
955    IF ( ocean_mode          )  CALL ocean_header( io )
956    IF ( plant_canopy        )  CALL pcm_header( io )
957    IF ( radiation           )  CALL radiation_header( io )
958    IF ( salsa               )  CALL salsa_header( io )
959    IF ( calculate_spectra   )  CALL spectra_header( io )
960    IF ( syn_turb_gen        )  CALL stg_header( io )
961    IF ( user_module_enabled )  CALL user_header( io )
962
963    IF ( debug_output )  CALL debug_message( 'module-specific header output', 'end' )
964
965
966 END SUBROUTINE module_interface_header
967
968
969!------------------------------------------------------------------------------!
970! Description:
971! ------------
972!> Perform module-specific actions while in time-integration (vector-optimized)
973!------------------------------------------------------------------------------!
974 SUBROUTINE module_interface_actions( location )
975
976
977    CHARACTER (LEN=*), INTENT(IN) ::  location !< call location string
978
979
980    IF ( bulk_cloud_model    )  CALL bcm_actions( location )
981    IF ( air_chemistry       )  CALL chem_actions( location )
982    IF ( gust_module_enabled )  CALL gust_actions( location )
983    IF ( ocean_mode          )  CALL ocean_actions( location )
984    IF ( salsa               )  CALL salsa_actions( location )
985    IF ( wind_turbine        )  CALL wtm_actions( location )
986    IF ( user_module_enabled )  CALL user_actions( location )
987
988
989 END SUBROUTINE module_interface_actions
990
991
992!------------------------------------------------------------------------------!
993! Description:
994! ------------
995!> Perform module-specific actions while in time-integration (cache-optimized)
996!------------------------------------------------------------------------------!
997 SUBROUTINE module_interface_actions_ij( i, j, location )
998
999
1000    INTEGER(iwp),      INTENT(IN) ::  i         !< grid index in x-direction
1001    INTEGER(iwp),      INTENT(IN) ::  j         !< grid index in y-direction
1002    CHARACTER (LEN=*), INTENT(IN) ::  location  !< call location string
1003
1004
1005    IF ( bulk_cloud_model    )  CALL bcm_actions( i, j, location )
1006    IF ( air_chemistry       )  CALL chem_actions( i, j, location )
1007    IF ( gust_module_enabled )  CALL gust_actions( i, j, location )
1008    IF ( ocean_mode          )  CALL ocean_actions( i, j, location )
1009    IF ( salsa               )  CALL salsa_actions( i, j, location )
1010    IF ( wind_turbine        )  CALL wtm_actions( i, j, location )
1011    IF ( user_module_enabled )  CALL user_actions( i, j, location )
1012
1013
1014 END SUBROUTINE module_interface_actions_ij
1015
1016
1017!------------------------------------------------------------------------------!
1018! Description:
1019! ------------
1020!> Compute module-specific non_advective_processes (vector-optimized)
1021!------------------------------------------------------------------------------!
1022 SUBROUTINE module_interface_non_advective_processes()
1023
1024
1025    IF ( bulk_cloud_model    )  CALL bcm_non_advective_processes()
1026    IF ( air_chemistry       )  CALL chem_non_advective_processes()
1027    IF ( salsa               )  CALL salsa_non_advective_processes()
1028
1029
1030 END SUBROUTINE module_interface_non_advective_processes
1031
1032
1033!------------------------------------------------------------------------------!
1034! Description:
1035! ------------
1036!> Compute module-specific non_advective_processes (cache-optimized)
1037!------------------------------------------------------------------------------!
1038 SUBROUTINE module_interface_non_advective_processes_ij( i, j )
1039
1040
1041    INTEGER(iwp), INTENT(IN) ::  i            !< grid index in x-direction
1042    INTEGER(iwp), INTENT(IN) ::  j            !< grid index in y-direction
1043
1044
1045    IF ( bulk_cloud_model    )  CALL bcm_non_advective_processes( i, j )
1046    IF ( air_chemistry       )  CALL chem_non_advective_processes( i, j )
1047    IF ( salsa               )  CALL salsa_non_advective_processes( i, j )
1048
1049
1050 END SUBROUTINE module_interface_non_advective_processes_ij
1051 
1052!------------------------------------------------------------------------------!
1053! Description:
1054! ------------
1055!> Exchange horiz for module-specific quantities
1056!------------------------------------------------------------------------------!
1057 SUBROUTINE module_interface_exchange_horiz()
1058
1059
1060    IF ( debug_output_timestep )  CALL debug_message( 'module-specific exchange_horiz', 'start' )
1061
1062    IF ( bulk_cloud_model    )  CALL bcm_exchange_horiz()
1063    IF ( air_chemistry       )  CALL chem_exchange_horiz_bounds()
1064    IF ( salsa               )  CALL salsa_exchange_horiz_bounds()
1065
1066    IF ( debug_output_timestep )  CALL debug_message( 'module-specific exchange_horiz', 'end' )
1067
1068
1069 END SUBROUTINE module_interface_exchange_horiz
1070
1071
1072!------------------------------------------------------------------------------!
1073! Description:
1074! ------------
1075!> Compute module-specific prognostic_equations (vector-optimized)
1076!------------------------------------------------------------------------------!
1077 SUBROUTINE module_interface_prognostic_equations()
1078
1079
1080    IF ( bulk_cloud_model    )  CALL bcm_prognostic_equations()
1081    IF ( air_chemistry       )  CALL chem_prognostic_equations()
1082    IF ( gust_module_enabled )  CALL gust_prognostic_equations()
1083    IF ( ocean_mode          )  CALL ocean_prognostic_equations()
1084    IF ( salsa               )  CALL salsa_prognostic_equations()
1085
1086
1087 END SUBROUTINE module_interface_prognostic_equations
1088
1089
1090!------------------------------------------------------------------------------!
1091! Description:
1092! ------------
1093!> Compute module-specific prognostic_equations (cache-optimized)
1094!------------------------------------------------------------------------------!
1095 SUBROUTINE module_interface_prognostic_equations_ij( i, j, i_omp_start, tn )
1096
1097
1098    INTEGER(iwp), INTENT(IN) ::  i            !< grid index in x-direction
1099    INTEGER(iwp), INTENT(IN) ::  j            !< grid index in y-direction
1100    INTEGER(iwp), INTENT(IN) ::  i_omp_start  !< first loop index of i-loop in prognostic_equations
1101    INTEGER(iwp), INTENT(IN) ::  tn           !< task number of openmp task
1102
1103
1104    IF ( bulk_cloud_model    )  CALL bcm_prognostic_equations( i, j, i_omp_start, tn )
1105    IF ( air_chemistry       )  CALL chem_prognostic_equations( i, j, i_omp_start, tn )
1106    IF ( gust_module_enabled )  CALL gust_prognostic_equations( i, j, i_omp_start, tn )
1107    IF ( ocean_mode          )  CALL ocean_prognostic_equations( i, j, i_omp_start, tn )
1108    IF ( salsa               )  CALL salsa_prognostic_equations( i, j, i_omp_start, tn )
1109
1110
1111 END SUBROUTINE module_interface_prognostic_equations_ij
1112
1113
1114!------------------------------------------------------------------------------!
1115! Description:
1116! ------------
1117!> Swap the timelevel pointers for module-specific arrays
1118!------------------------------------------------------------------------------!
1119 SUBROUTINE module_interface_swap_timelevel ( swap_mode )
1120
1121
1122    INTEGER(iwp), INTENT(IN) :: swap_mode !< determines procedure of pointer swap
1123
1124
1125    IF ( debug_output_timestep )  CALL debug_message( 'module-specific swap timelevel', 'start' )
1126
1127    IF ( bulk_cloud_model    )  CALL bcm_swap_timelevel( swap_mode )
1128    IF ( air_chemistry       )  CALL chem_swap_timelevel( swap_mode )
1129    IF ( gust_module_enabled )  CALL gust_swap_timelevel( swap_mode )
1130    IF ( land_surface        )  CALL lsm_swap_timelevel( swap_mode )
1131    IF ( ocean_mode          )  CALL ocean_swap_timelevel( swap_mode )
1132    IF ( salsa               )  CALL salsa_swap_timelevel( swap_mode )
1133    IF ( urban_surface       )  CALL usm_swap_timelevel( swap_mode )
1134
1135    IF ( debug_output_timestep )  CALL debug_message( 'module-specific swap timelevel', 'end' )
1136
1137
1138 END SUBROUTINE module_interface_swap_timelevel
1139
1140
1141!------------------------------------------------------------------------------!
1142!
1143! Description:
1144! ------------
1145!> Perform module-specific averaging of 3D data
1146!------------------------------------------------------------------------------!
1147 SUBROUTINE module_interface_3d_data_averaging( mode, variable )
1148
1149
1150    CHARACTER (LEN=*), INTENT(IN) ::  mode     !< averaging interface mode
1151    CHARACTER (LEN=*), INTENT(IN) ::  variable !< variable name
1152
1153
1154    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 3d data averaging', 'start' )
1155
1156    IF ( biometeorology      )  CALL bio_3d_data_averaging( mode, variable )
1157    IF ( bulk_cloud_model    )  CALL bcm_3d_data_averaging( mode, variable )
1158    IF ( air_chemistry       )  CALL chem_3d_data_averaging( mode, variable )
1159    IF ( gust_module_enabled )  CALL gust_3d_data_averaging( mode, variable )
1160    IF ( land_surface        )  CALL lsm_3d_data_averaging( mode, variable )
1161    IF ( ocean_mode          )  CALL ocean_3d_data_averaging( mode, variable )
1162    IF ( radiation           )  CALL radiation_3d_data_averaging( mode, variable )
1163    IF ( salsa               )  CALL salsa_3d_data_averaging( mode, variable )
1164    IF ( urban_surface       )  CALL usm_3d_data_averaging( mode, variable )
1165    IF ( user_module_enabled )  CALL user_3d_data_averaging( mode, variable )
1166
1167    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 3d data averaging', 'end' )
1168
1169
1170 END SUBROUTINE module_interface_3d_data_averaging
1171
1172!------------------------------------------------------------------------------!
1173!
1174! Description:
1175! ------------
1176!> Define module-specific 2D output variables
1177!------------------------------------------------------------------------------!
1178 SUBROUTINE module_interface_data_output_2d( av, variable, found, grid, mode,  &
1179                                             local_pf, two_d, nzb_do, nzt_do,  &
1180                                             fill_value )
1181
1182    INTEGER(iwp),      INTENT(IN)    ::  av         !< flag for (non-)average output
1183    CHARACTER (LEN=*), INTENT(IN)    ::  variable   !< variable name
1184    LOGICAL,           INTENT(INOUT) ::  found      !< flag if output variable is found
1185    CHARACTER (LEN=*), INTENT(INOUT) ::  grid       !< name of vertical grid
1186    CHARACTER (LEN=*), INTENT(IN)    ::  mode       !< either 'xy', 'xz' or 'yz'
1187    LOGICAL,           INTENT(OUT)   ::  two_d      !< flag for 2D variables
1188    INTEGER(iwp),      INTENT(IN)    ::  nzb_do     !< vertical output index (bottom) (usually 0)
1189    INTEGER(iwp),      INTENT(IN)    ::  nzt_do     !< vertical output index (top) (usually nz_do3d)
1190    REAL(wp),          INTENT(IN)    ::  fill_value !< to be removed
1191
1192    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf !< ToDo: can also be kind=sp
1193
1194
1195    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 2d data output', 'start' )
1196
1197    IF ( .NOT. found  .AND.  biometeorology )  THEN
1198       CALL bio_data_output_2d(                                                &
1199               av, variable, found, grid, local_pf, two_d, nzb_do, nzt_do      &
1200            )
1201    ENDIF
1202
1203    IF ( .NOT. found  .AND.  bulk_cloud_model )  THEN
1204       CALL bcm_data_output_2d(                                                &
1205               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1206            )
1207    ENDIF
1208
1209    IF ( .NOT. found  .AND.  air_chemistry )  THEN
1210       CALL chem_data_output_2d(                                               &
1211               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
1212            )
1213    ENDIF
1214
1215    IF ( .NOT. found  .AND.  gust_module_enabled )  THEN
1216       CALL gust_data_output_2d(                                               &
1217               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
1218            )
1219    ENDIF
1220
1221    IF ( .NOT. found  .AND.  land_surface )  THEN
1222       CALL lsm_data_output_2d(                                                &
1223               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1224            )
1225    ENDIF
1226
1227    IF ( .NOT. found  .AND.  ocean_mode )  THEN
1228       CALL ocean_data_output_2d(                                              &
1229               av, variable, found, grid, mode, local_pf, nzb_do, nzt_do       &
1230            )
1231    ENDIF
1232
1233    IF ( .NOT. found  .AND.  radiation )  THEN
1234       CALL radiation_data_output_2d(                                          &
1235               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1236            )
1237    ENDIF
1238
1239    IF ( .NOT. found  .AND.  salsa )  THEN
1240       CALL salsa_data_output_2d(                                              &
1241               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
1242            )
1243    ENDIF
1244
1245    IF ( .NOT. found  .AND.  user_module_enabled )  THEN
1246       CALL user_data_output_2d(                                               &
1247               av, variable, found, grid, local_pf, two_d, nzb_do, nzt_do      &
1248            )
1249    ENDIF
1250
1251    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 2d data output', 'end' )
1252
1253
1254 END SUBROUTINE module_interface_data_output_2d
1255
1256
1257!------------------------------------------------------------------------------!
1258!
1259! Description:
1260! ------------
1261!> Define module-specific 3D output variables
1262!------------------------------------------------------------------------------!
1263 SUBROUTINE module_interface_data_output_3d( av, variable, found, local_pf,    &
1264                                             fill_value, resorted, nzb_do, nzt_do )
1265
1266
1267    INTEGER(iwp),      INTENT(IN)    ::  av         !< flag for (non-)average output
1268    CHARACTER (LEN=*), INTENT(IN)    ::  variable   !< variable name
1269    LOGICAL,           INTENT(INOUT) ::  found      !< flag if output variable is found
1270    REAL(wp),          INTENT(IN)    ::  fill_value !< ToDo: refactor
1271    LOGICAL,           INTENT(OUT)   ::  resorted   !< flag if output has been resorted
1272    INTEGER(iwp),      INTENT(IN)    ::  nzb_do     !< vertical output index (bottom) (usually 0)
1273    INTEGER(iwp),      INTENT(IN)    ::  nzt_do     !< vertical output index (top) (usually nz_do3d)
1274
1275    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf
1276
1277
1278    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 3d data output', 'start' )
1279
1280    IF ( .NOT. found  .AND.  biometeorology )  THEN
1281       CALL bio_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1282       resorted = .FALSE.
1283    ENDIF
1284
1285    IF ( .NOT. found  .AND.  bulk_cloud_model )  THEN
1286       CALL bcm_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1287       resorted = .TRUE.
1288    ENDIF
1289
1290    IF ( .NOT. found  .AND.  air_chemistry )  THEN
1291       CALL chem_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1292       resorted = .TRUE.
1293    ENDIF
1294
1295    IF ( .NOT. found  .AND.  gust_module_enabled )  THEN
1296       CALL gust_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1297       resorted = .TRUE.
1298    ENDIF
1299   
1300    IF ( .NOT. found  .AND.  indoor_model )  THEN
1301       CALL im_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1302       resorted = .TRUE.
1303    ENDIF
1304
1305    IF ( .NOT. found  .AND.  ocean_mode )  THEN
1306       CALL ocean_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1307       resorted = .TRUE.
1308    ENDIF
1309
1310    IF ( .NOT. found  .AND.  plant_canopy )  THEN
1311       CALL pcm_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
1312       resorted = .TRUE.
1313    ENDIF
1314
1315    IF ( .NOT. found  .AND.  radiation )  THEN
1316       CALL radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1317       resorted = .TRUE.
1318    ENDIF
1319
1320    IF ( .NOT. found  .AND.  salsa )  THEN
1321       CALL salsa_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1322       resorted = .TRUE.
1323    ENDIF
1324
1325    IF ( .NOT. found  .AND.  user_module_enabled )  THEN
1326       CALL user_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
1327       resorted = .TRUE.
1328    ENDIF
1329
1330    IF ( debug_output_timestep )  CALL debug_message( 'module-specific 3d data output', 'end' )
1331
1332
1333 END SUBROUTINE module_interface_data_output_3d
1334
1335
1336!------------------------------------------------------------------------------!
1337! Description:
1338! ------------
1339!> Compute module-specific profile and timeseries data
1340!------------------------------------------------------------------------------!
1341 SUBROUTINE module_interface_statistics( mode, sr, tn, dots_max )
1342
1343
1344    CHARACTER (LEN=*), INTENT(IN) ::  mode     !< statistical analysis mode
1345    INTEGER(iwp),      INTENT(IN) ::  sr       !<
1346    INTEGER(iwp),      INTENT(IN) ::  tn       !<
1347    INTEGER(iwp),      INTENT(IN) ::  dots_max !< maximum number of timeseries
1348
1349
1350    IF ( debug_output_timestep )  CALL debug_message( 'module-specific statistics', 'start' )
1351
1352    IF ( gust_module_enabled )  CALL gust_statistics( mode, sr, tn, dots_max )
1353    IF ( air_chemistry       )  CALL chem_statistics( mode, sr, tn )
1354    IF ( user_module_enabled )  CALL user_statistics( mode, sr, tn )
1355
1356    IF ( debug_output_timestep )  CALL debug_message( 'module-specific statistics', 'end' )
1357
1358
1359 END SUBROUTINE module_interface_statistics
1360
1361
1362!------------------------------------------------------------------------------!
1363! Description:
1364! ------------
1365!> Read module-specific restart data globaly shared by all MPI ranks
1366!------------------------------------------------------------------------------!
1367 SUBROUTINE module_interface_rrd_global( found )
1368
1369
1370    LOGICAL, INTENT(INOUT) ::  found    !< flag if variable was found
1371
1372
1373    IF ( debug_output )  CALL debug_message( 'module-specific read global restart data', 'start' )
1374
1375    IF ( .NOT. found )  CALL bio_rrd_global( found ) ! ToDo: change interface to pass variable
1376    IF ( .NOT. found )  CALL bcm_rrd_global( found ) ! ToDo: change interface to pass variable
1377    IF ( .NOT. found )  CALL flight_rrd_global( found ) ! ToDo: change interface to pass variable
1378    IF ( .NOT. found )  CALL gust_rrd_global( found ) ! ToDo: change interface to pass variable
1379    IF ( .NOT. found )  CALL ocean_rrd_global( found ) ! ToDo: change interface to pass variable
1380    IF ( .NOT. found )  CALL stg_rrd_global ( found ) ! ToDo: change interface to pass variable
1381    IF ( .NOT. found )  CALL wtm_rrd_global( found ) ! ToDo: change interface to pass variable
1382    IF ( .NOT. found )  CALL surface_data_output_rrd_global( found )
1383    IF ( .NOT. found )  CALL user_rrd_global( found ) ! ToDo: change interface to pass variable
1384
1385    IF ( debug_output )  CALL debug_message( 'module-specific read global restart data', 'end' )
1386
1387
1388 END SUBROUTINE module_interface_rrd_global
1389
1390
1391!------------------------------------------------------------------------------!
1392! Description:
1393! ------------
1394!> Write module-specific restart data globaly shared by all MPI ranks
1395!------------------------------------------------------------------------------!
1396 SUBROUTINE module_interface_wrd_global
1397
1398
1399    IF ( debug_output )  CALL debug_message( 'module-specific write global restart data', 'start' )
1400
1401    IF ( biometeorology )       CALL bio_wrd_global
1402    IF ( bulk_cloud_model )     CALL bcm_wrd_global
1403    IF ( virtual_flight )       CALL flight_wrd_global
1404    IF ( gust_module_enabled )  CALL gust_wrd_global
1405    IF ( ocean_mode )           CALL ocean_wrd_global
1406    IF ( syn_turb_gen )         CALL stg_wrd_global
1407    IF ( wind_turbine )         CALL wtm_wrd_global
1408    IF ( surface_output )       CALL surface_data_output_wrd_global
1409    IF ( user_module_enabled )  CALL user_wrd_global
1410
1411    IF ( debug_output )  CALL debug_message( 'module-specific write global restart data', 'end' )
1412
1413
1414 END SUBROUTINE module_interface_wrd_global
1415
1416
1417!------------------------------------------------------------------------------!
1418! Description:
1419! ------------
1420!> Read module-specific restart data specific to local MPI ranks
1421!------------------------------------------------------------------------------!
1422 SUBROUTINE module_interface_rrd_local( map_index,                             &
1423                                        nxlf, nxlc, nxl_on_file,               &
1424                                        nxrf, nxrc, nxr_on_file,               &
1425                                        nynf, nync, nyn_on_file,               &
1426                                        nysf, nysc, nys_on_file,               &
1427                                        tmp_2d, tmp_3d, found )
1428
1429
1430    INTEGER(iwp), INTENT(IN)  ::  map_index    !<
1431    INTEGER(iwp), INTENT(IN)  ::  nxlc         !<
1432    INTEGER(iwp), INTENT(IN)  ::  nxlf         !<
1433    INTEGER(iwp), INTENT(IN)  ::  nxl_on_file  !<
1434    INTEGER(iwp), INTENT(IN)  ::  nxrc         !<
1435    INTEGER(iwp), INTENT(IN)  ::  nxrf         !<
1436    INTEGER(iwp), INTENT(IN)  ::  nxr_on_file  !<
1437    INTEGER(iwp), INTENT(IN)  ::  nync         !<
1438    INTEGER(iwp), INTENT(IN)  ::  nynf         !<
1439    INTEGER(iwp), INTENT(IN)  ::  nyn_on_file  !<
1440    INTEGER(iwp), INTENT(IN)  ::  nysc         !<
1441    INTEGER(iwp), INTENT(IN)  ::  nysf         !<
1442    INTEGER(iwp), INTENT(IN)  ::  nys_on_file  !<
1443    LOGICAL,      INTENT(INOUT) ::  found        !< flag if variable was found
1444
1445    REAL(wp), DIMENSION(nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp), INTENT(OUT) :: tmp_2d   !<
1446    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   !<
1447
1448
1449    IF ( debug_output )  CALL debug_message( 'module-specific read local restart data', 'start' )
1450
1451    IF ( .NOT. found ) CALL bio_rrd_local(                                     &
1452                               found                                           &
1453                            )
1454
1455    IF ( .NOT. found ) CALL bcm_rrd_local(                                     &
1456                               map_index,                                      &
1457                               nxlf, nxlc, nxl_on_file,                        &
1458                               nxrf, nxrc, nxr_on_file,                        &
1459                               nynf, nync, nyn_on_file,                        &
1460                               nysf, nysc, nys_on_file,                        &
1461                               tmp_2d, tmp_3d, found                           &
1462                            ) ! ToDo: change interface to pass variable
1463
1464    IF ( .NOT. found ) CALL chem_rrd_local(                                    &
1465                               map_index,                                      &
1466                               nxlf, nxlc, nxl_on_file,                        &
1467                               nxrf, nxrc, nxr_on_file,                        &
1468                               nynf, nync, nyn_on_file,                        &
1469                               nysf, nysc, nys_on_file,                        &
1470                               tmp_3d, found                                   &
1471                            ) ! ToDo: change interface to pass variable
1472
1473    IF ( .NOT. found ) CALL gust_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 lsm_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_2d, found                                   &
1489                            ) ! ToDo: change interface to pass variable
1490
1491    IF ( .NOT. found ) CALL ocean_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_3d, found                                   &
1498                            ) ! ToDo: change interface to pass variable
1499
1500    IF ( .NOT. found ) CALL radiation_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_2d, tmp_3d, found                           &
1507                            ) ! ToDo: change interface to pass variable
1508
1509    IF ( .NOT. found ) CALL salsa_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_3d, found                                   &
1516                            ) ! ToDo: change interface to pass variable
1517
1518    IF ( .NOT. found ) CALL usm_rrd_local(                                     &
1519                               map_index,                                      &
1520                               nxlf, nxlc, nxl_on_file,                        &
1521                               nxrf, nxr_on_file,                              &
1522                               nynf, nyn_on_file,                              &
1523                               nysf, nysc, nys_on_file,                        &
1524                               found                                           &
1525                            ) ! ToDo: change interface to pass variable
1526!
1527!-- Surface data do not need overlap data, so do not pass these information.
1528    IF ( .NOT. found ) CALL surface_data_output_rrd_local( found )
1529                           
1530    IF ( .NOT. found ) CALL user_rrd_local(                                    &
1531                               map_index,                                      &
1532                               nxlf, nxlc, nxl_on_file,                        &
1533                               nxrf, nxrc, nxr_on_file,                        &
1534                               nynf, nync, nyn_on_file,                        &
1535                               nysf, nysc, nys_on_file,                        &
1536                               tmp_3d, found                                   &
1537                            ) ! ToDo: change interface to pass variable
1538
1539    IF ( debug_output )  CALL debug_message( 'module-specific read local restart data', 'end' )
1540
1541
1542 END SUBROUTINE module_interface_rrd_local
1543
1544
1545!------------------------------------------------------------------------------!
1546! Description:
1547! ------------
1548!> Write module-specific restart data specific to local MPI ranks
1549!------------------------------------------------------------------------------!
1550 SUBROUTINE module_interface_wrd_local
1551
1552
1553    IF ( debug_output )  CALL debug_message( 'module-specific write local restart data', 'start' )
1554
1555    IF ( biometeorology )       CALL bio_wrd_local
1556    IF ( bulk_cloud_model )     CALL bcm_wrd_local
1557    IF ( air_chemistry )        CALL chem_wrd_local
1558    IF ( gust_module_enabled )  CALL gust_wrd_local
1559    IF ( land_surface )         CALL lsm_wrd_local
1560    IF ( ocean_mode )           CALL ocean_wrd_local
1561    IF ( radiation )            CALL radiation_wrd_local
1562    IF ( salsa )                CALL salsa_wrd_local
1563    IF ( urban_surface )        CALL usm_wrd_local
1564    IF ( surface_output )       CALL surface_data_output_wrd_local
1565    IF ( user_module_enabled )  CALL user_wrd_local
1566
1567    IF ( debug_output )  CALL debug_message( 'module-specific write local restart data', 'end' )
1568
1569
1570 END SUBROUTINE module_interface_wrd_local
1571
1572
1573!------------------------------------------------------------------------------!
1574! Description:
1575! ------------
1576!> Perform module-specific last actions before the program terminates
1577!------------------------------------------------------------------------------!
1578 SUBROUTINE module_interface_last_actions
1579
1580
1581    IF ( debug_output )  CALL debug_message( 'module-specific last actions', 'start' )
1582
1583    IF ( virtual_measurement )  CALL vm_last_actions
1584    IF ( user_module_enabled )  CALL user_last_actions
1585
1586    IF ( debug_output )  CALL debug_message( 'module-specific last actions', 'end' )
1587
1588
1589 END SUBROUTINE module_interface_last_actions
1590
1591
1592 END MODULE module_interface
Note: See TracBrowser for help on using the repository browser.