Home

Context Navigation

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

Last change on this file since 3937 was 3931, checked in by schwenkel, 6 years ago
renamed module interface from non_transport_physics to non_advective_processes
Property svn:keywords set to `Id`
File size: 67.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |