Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |