Home

Context Navigation

source: palm/trunk/SOURCE/surface_output_mod.f90 @ 3614

Last change on this file since 3614 was 3614, checked in by raasch, 6 years ago
unused variables removed, abort renamed inifor_abort to avoid intrinsic problem in Fortran
Property svn:keywords set to `Id`
File size: 204.7 KB

Rev	Line
[3421]	1	!> @file surface_output_mod.f90
	2	!------------------------------------------------------------------------------!
	3	! This file is part of the PALM model system.
	4	!
	5	! PALM is free software: you can redistribute it and/or modify it under the
	6	! terms of the GNU General Public License as published by the Free Software
	7	! Foundation, either version 3 of the License, or (at your option) any later
	8	! version.
	9	!
	10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
	11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
	12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
	13	!
	14	! You should have received a copy of the GNU General Public License along with
	15	! PALM. If not, see <http://www.gnu.org/licenses/>.
	16	!
	17	! Copyright 1997-2018 Leibniz Universitaet Hannover
	18	!------------------------------------------------------------------------------!
	19	!
	20	! Current revisions:
	21	! ------------------
	22	!
	23	!
	24	! Former revisions:
	25	! -----------------
	26	! $Id: surface_output_mod.f90 3614 2018-12-10 07:05:46Z raasch $
[3614]	27	! unused variables removed
	28	!
	29	! 3572 2018-11-28 11:40:28Z suehring
[3572]	30	! Added short- and longwave radiation flux arrays (e.g. diffuse, direct,
	31	! reflected, resedual) for all surfaces (M. Salim)
	32	!
	33	! 3494 2018-11-06 14:51:27Z suehring
[3494]	34	! Bugfix in gathering surface data from different types and orientation.
	35	! Output of total number of surfaces and vertices added.
	36	! String length is output, for more convinient post-processing.
	37	! Last actions added.
	38	!
	39	! 3483 2018-11-02 14:19:26Z raasch
[3483]	40	! bugfix: missed directives for MPI added
	41	!
	42	! 3421 2018-10-24 18:39:32Z gronemeier
[3421]	43	! Add output variables
	44	! Write data into binary file
	45	!
	46	! 3420 2018-10-24 17:30:08Z gronemeier
[3494]	47	! Initial implementation from Klaus Ketelsen and Matthias Suehring
[3421]	48	!
[3494]	49	!
	50	! Authors:
	51	! --------
	52	! @author Klaus Ketelsen, Matthias Suehring
	53	!
[3421]	54	! Description:
	55	! ------------
	56	!> Generate output for surface data.
[3494]	57	!>
	58	!> @todo Create namelist file for post-processing tool.
[3421]	59	!------------------------------------------------------------------------------!
	60
	61	MODULE surface_output_mod
	62
	63	USE kinds
	64
	65	USE arrays_3d, &
	66	ONLY: zu, zw
[3494]	67
[3421]	68	USE control_parameters, &
	69	ONLY: surface_data_output
[3494]	70
[3421]	71	USE grid_variables, &
	72	ONLY: dx,dy
[3494]	73
[3421]	74	USE indices, &
	75	ONLY: nxl, nxr, nys, nyn, nzb, nzt
[3494]	76
	77	USE pegrid
	78
[3421]	79	USE surface_mod, &
	80	ONLY: surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, &
	81	surf_usm_h, surf_usm_v
	82
	83	IMPLICIT NONE
	84
[3494]	85	TYPE surf_out !< data structure which contains all surfaces elements of all types on subdomain
[3421]	86
[3494]	87	INTEGER(iwp) :: ns !< number of surface elements on subdomain
	88	INTEGER(iwp) :: ns_total !< total number of surface elements
	89	INTEGER(iwp) :: npoints !< number of points / vertices which define a surface element (on subdomain)
	90	INTEGER(iwp) :: npoints_total !< total number of points / vertices which define a surface element
[3421]	91
	92	REAL(wp) :: fillvalue = -9999.9_wp !< fillvalue for surface elements which are not defined
	93
	94	REAL(wp), DIMENSION(:), ALLOCATABLE :: var_out !< output variables
	95	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: var_av !< variables used for averaging
	96	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: points !< points / vertices of a surface element
	97	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: polygons !< polygon data of a surface element
	98	END TYPE surf_out
	99
	100	CHARACTER(LEN=100), DIMENSION(300) :: data_output_surf = ' ' !< namelist variable which describes the output variables
	101	CHARACTER(LEN=100), DIMENSION(0:1,300) :: dosurf = ' ' !< internal variable which describes the output variables and separates
	102	!< averaged from non-averaged output
	103
	104	INTEGER(iwp) :: average_count_surf = 0 !< number of ensemble members used for averaging
	105	INTEGER(iwp) :: dosurf_no(0:1) = 0 !< number of surface output quantities
	106
	107	LOGICAL :: first_output(0:1) = .FALSE. !< true if first output was already called
	108
	109	REAL(wp) :: averaging_interval_surf = 9999999.9_wp !< averaging interval
	110	REAL(wp) :: dt_dosurf = 9999999.9_wp !< time interval for instantaneous data output
	111	REAL(wp) :: dt_dosurf_av = 9999999.9_wp !< time interval for averaged data output
	112	REAL(wp) :: skip_time_dosurf = 0.0_wp !< skip time for instantaneous data output
	113	REAL(wp) :: skip_time_dosurf_av = 0.0_wp !< skip time for averaged data output
	114	REAL(wp) :: time_dosurf = 0.0_wp !< internal counter variable to check for instantaneous data output
	115	REAL(wp) :: time_dosurf_av = 0.0_wp !< internal counter variable to check for averaged data output
	116
	117
	118	TYPE(surf_out) :: surfaces !< variable which contains all required output information
	119
	120	SAVE
	121
	122	PRIVATE
	123
	124	INTERFACE surface_output
	125	MODULE PROCEDURE surface_output
	126	END INTERFACE surface_output
	127
	128	INTERFACE surface_output_averaging
	129	MODULE PROCEDURE surface_output_averaging
	130	END INTERFACE surface_output_averaging
	131
	132	INTERFACE surface_output_check_parameters
	133	MODULE PROCEDURE surface_output_check_parameters
	134	END INTERFACE surface_output_check_parameters
	135
	136	INTERFACE surface_output_init
	137	MODULE PROCEDURE surface_output_init
	138	END INTERFACE surface_output_init
	139
[3494]	140	INTERFACE surface_output_last_action
	141	MODULE PROCEDURE surface_output_last_action
	142	END INTERFACE surface_output_last_action
	143
[3421]	144	INTERFACE surface_output_parin
	145	MODULE PROCEDURE surface_output_parin
	146	END INTERFACE surface_output_parin
	147
	148	!
	149	!--Public subroutines
	150	PUBLIC surface_output, surface_output_averaging, &
	151	surface_output_check_parameters, surface_output_init, &
[3494]	152	surface_output_last_action, surface_output_parin
[3421]	153	!
	154	!--Public variables
	155	PUBLIC average_count_surf, averaging_interval_surf, dt_dosurf, dt_dosurf_av,&
	156	skip_time_dosurf, skip_time_dosurf_av, time_dosurf, time_dosurf_av
	157
	158	CONTAINS
	159
	160	!------------------------------------------------------------------------------!
	161	! Description:
	162	! ------------
	163	!> Initialization surface-data output data structure: calculation of vertices
	164	!> and polygon data for the surface elements, allocation of required arrays.
	165	!------------------------------------------------------------------------------!
	166	SUBROUTINE surface_output_init
[3494]	167
[3421]	168	IMPLICIT NONE
	169
[3494]	170	INTEGER(iwp) :: i !< grid index in x-direction, also running variable for counting non-average data output
	171	INTEGER(iwp) :: j !< grid index in y-direction, also running variable for counting average data output
	172	INTEGER(iwp) :: k !< grid index in z-direction
	173	INTEGER(iwp) :: l !< running index for surface-element orientation
	174	INTEGER(iwp) :: m !< running index for surface elements
	175	INTEGER(iwp) :: npg !< counter variable for all surface elements ( or polygons )
	176	INTEGER(iwp) :: point_index_count !< local counter variable for point index
[3421]	177
[3494]	178	INTEGER(iwp), DIMENSION(0:numprocs-1) :: num_points_on_pe !< array which contains the number of points on all mpi ranks
[3421]	179	INTEGER(iwp), ALLOCATABLE, DIMENSION(:,:,:) :: point_index !< dummy array used to check where the reference points for surface polygons are located
	180
	181	!
	182	!-- Determine the number of surface elements on subdomain
	183	surfaces%ns = surf_def_h(0)%ns + surf_lsm_h%ns + surf_usm_h%ns & !horizontal upward-facing
	184	+ surf_def_h(1)%ns & !horizontal downard-facing
	185	+ surf_def_v(0)%ns + surf_lsm_v(0)%ns + surf_usm_v(0)%ns & !northward-facing
	186	+ surf_def_v(1)%ns + surf_lsm_v(1)%ns + surf_usm_v(1)%ns & !southward-facing
	187	+ surf_def_v(2)%ns + surf_lsm_v(2)%ns + surf_usm_v(2)%ns & !westward-facing
	188	+ surf_def_v(3)%ns + surf_lsm_v(3)%ns + surf_usm_v(3)%ns !eastward-facing
	189	!
	190	!-- Allocate output variable and set to _FillValue attribute
	191	ALLOCATE ( surfaces%var_out(1:surfaces%ns) )
	192	surfaces%var_out = surfaces%fillvalue
	193	!
	194	!-- If there is an output of time average output variables, allocate the
	195	!-- required array.
	196	IF ( dosurf_no(1) > 0 ) THEN
	197	ALLOCATE ( surfaces%var_av(1:surfaces%ns,1:dosurf_no(1)) )
	198	surfaces%var_av = 0.0_wp
	199	ENDIF
	200	!
	201	!-- Initialize points and polygon data.
	202	!-- In a first step, count the number of points which are defining
	203	!-- the surfaces and the polygons.
	204	ALLOCATE( point_index(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	205	point_index = -1
	206	!
	207	!-- Horizontal default surfaces
	208	surfaces%npoints = 0
	209	DO l = 0, 1
	210	DO m = 1, surf_def_h(0)%ns
	211	!
	212	!-- Determine the indices of the respective grid cell inside the topography
	213	i = surf_def_h(0)%i(m) + surf_def_h(0)%ioff
	214	j = surf_def_h(0)%j(m) + surf_def_h(0)%joff
	215	k = surf_def_h(0)%k(m) + surf_def_h(0)%koff
	216	!
	217	!-- Check if the vertices that define the surface element are already
	218	!-- defined, if not, increment the counter.
[3494]	219	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	220	surfaces%npoints = surfaces%npoints + 1
[3494]	221	point_index(k,j,i) = surfaces%npoints - 1
[3421]	222	ENDIF
[3494]	223	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	224	surfaces%npoints = surfaces%npoints + 1
[3494]	225	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	226	ENDIF
[3494]	227	IF ( point_index(k,j+1,i+1) < 0 ) THEN
[3421]	228	surfaces%npoints = surfaces%npoints + 1
[3494]	229	point_index(k,j+1,i+1) = surfaces%npoints - 1
[3421]	230	ENDIF
[3494]	231	IF ( point_index(k,j+1,i) < 0 ) THEN
	232	surfaces%npoints = surfaces%npoints + 1
	233	point_index(k,j+1,i) = surfaces%npoints - 1
	234	ENDIF
[3421]	235	ENDDO
	236	ENDDO
	237	DO m = 1, surf_lsm_h%ns
	238	i = surf_lsm_h%i(m) + surf_lsm_h%ioff
	239	j = surf_lsm_h%j(m) + surf_lsm_h%joff
	240	k = surf_lsm_h%k(m) + surf_lsm_h%koff
	241
[3494]	242	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	243	surfaces%npoints = surfaces%npoints + 1
[3494]	244	point_index(k,j,i) = surfaces%npoints - 1
[3421]	245	ENDIF
[3494]	246	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	247	surfaces%npoints = surfaces%npoints + 1
[3494]	248	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	249	ENDIF
[3494]	250	IF ( point_index(k,j+1,i+1) < 0 ) THEN
	251	surfaces%npoints = surfaces%npoints + 1
	252	point_index(k,j+1,i+1) = surfaces%npoints - 1
	253	ENDIF
	254	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	255	surfaces%npoints = surfaces%npoints + 1
[3494]	256	point_index(k,j+1,i) = surfaces%npoints - 1
	257	ENDIF
[3421]	258	ENDDO
	259	DO m = 1, surf_usm_h%ns
	260	i = surf_usm_h%i(m) + surf_usm_h%ioff
	261	j = surf_usm_h%j(m) + surf_usm_h%joff
	262	k = surf_usm_h%k(m) + surf_usm_h%koff
	263
[3494]	264	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	265	surfaces%npoints = surfaces%npoints + 1
[3494]	266	point_index(k,j,i) = surfaces%npoints - 1
[3421]	267	ENDIF
[3494]	268	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	269	surfaces%npoints = surfaces%npoints + 1
[3494]	270	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	271	ENDIF
[3494]	272	IF ( point_index(k,j+1,i+1) < 0 ) THEN
	273	surfaces%npoints = surfaces%npoints + 1
	274	point_index(k,j+1,i+1) = surfaces%npoints - 1
	275	ENDIF
	276	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	277	surfaces%npoints = surfaces%npoints + 1
[3494]	278	point_index(k,j+1,i) = surfaces%npoints - 1
[3421]	279	ENDIF
	280	ENDDO
	281	!
	282	!-- Vertical surfaces
	283	DO l = 0, 3
	284	DO m = 1, surf_def_v(l)%ns
	285	!
	286	!-- Determine the indices of the respective grid cell inside the topography.
	287	!-- Please note, j-index for north-facing surfaces ( l==0 ) is
	288	!-- identical to the reference j-index outside the grid box.
	289	!-- Equivalent for east-facing surfaces and i-index.
[3494]	290	i = surf_def_v(l)%i(m) + MERGE( surf_def_v(l)%ioff, 0, l == 3 )
	291	j = surf_def_v(l)%j(m) + MERGE( surf_def_v(l)%joff, 0, l == 1 )
[3421]	292	k = surf_def_v(l)%k(m) + surf_def_v(l)%koff
	293	!
	294	!-- Lower left /front vertex
[3494]	295	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	296	surfaces%npoints = surfaces%npoints + 1
[3494]	297	point_index(k,j,i) = surfaces%npoints - 1
[3421]	298	ENDIF
	299	!
	300	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	301	IF ( l == 0 .OR. l == 1 ) THEN
	302	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	303	surfaces%npoints = surfaces%npoints + 1
[3494]	304	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	305	ENDIF
[3494]	306	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	307	surfaces%npoints = surfaces%npoints + 1
[3494]	308	point_index(k+1,j,i+1) = surfaces%npoints - 1
[3421]	309	ENDIF
	310	!
	311	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	312	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	313	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	314	surfaces%npoints = surfaces%npoints + 1
[3494]	315	point_index(k,j+1,i) = surfaces%npoints - 1
[3421]	316	ENDIF
[3494]	317	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	318	surfaces%npoints = surfaces%npoints + 1
[3494]	319	point_index(k+1,j+1,i) = surfaces%npoints - 1
[3421]	320	ENDIF
[3494]	321	ENDIF
	322	!
	323	!-- Upper left / front vertex
	324	IF ( point_index(k+1,j,i) < 0 ) THEN
	325	surfaces%npoints = surfaces%npoints + 1
	326	point_index(k+1,j,i) = surfaces%npoints - 1
	327	ENDIF
[3421]	328	ENDDO
	329	DO m = 1, surf_lsm_v(l)%ns
[3494]	330	i = surf_lsm_v(l)%i(m) + MERGE( surf_lsm_v(l)%ioff, 0, l == 3 )
	331	j = surf_lsm_v(l)%j(m) + MERGE( surf_lsm_v(l)%joff, 0, l == 1 )
[3421]	332	k = surf_lsm_v(l)%k(m) + surf_lsm_v(l)%koff
	333	!
	334	!-- Lower left /front vertex
[3494]	335	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	336	surfaces%npoints = surfaces%npoints + 1
[3494]	337	point_index(k,j,i) = surfaces%npoints - 1
[3421]	338	ENDIF
	339	!
	340	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	341	IF ( l == 0 .OR. l == 1 ) THEN
	342	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	343	surfaces%npoints = surfaces%npoints + 1
[3494]	344	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	345	ENDIF
[3494]	346	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	347	surfaces%npoints = surfaces%npoints + 1
[3494]	348	point_index(k+1,j,i+1) = surfaces%npoints - 1
[3421]	349	ENDIF
	350	!
	351	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	352	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	353	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	354	surfaces%npoints = surfaces%npoints + 1
[3494]	355	point_index(k,j+1,i) = surfaces%npoints - 1
[3421]	356	ENDIF
[3494]	357	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	358	surfaces%npoints = surfaces%npoints + 1
[3494]	359	point_index(k+1,j+1,i) = surfaces%npoints - 1
[3421]	360	ENDIF
[3494]	361	ENDIF
	362	!
	363	!-- Upper left / front vertex
	364	IF ( point_index(k+1,j,i) < 0 ) THEN
	365	surfaces%npoints = surfaces%npoints + 1
	366	point_index(k+1,j,i) = surfaces%npoints - 1
	367	ENDIF
[3421]	368	ENDDO
[3494]	369
[3421]	370	DO m = 1, surf_usm_v(l)%ns
[3494]	371	i = surf_usm_v(l)%i(m) + MERGE( surf_usm_v(l)%ioff, 0, l == 3 )
	372	j = surf_usm_v(l)%j(m) + MERGE( surf_usm_v(l)%joff, 0, l == 1 )
[3421]	373	k = surf_usm_v(l)%k(m) + surf_usm_v(l)%koff
	374	!
	375	!-- Lower left /front vertex
[3494]	376	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	377	surfaces%npoints = surfaces%npoints + 1
[3494]	378	point_index(k,j,i) = surfaces%npoints - 1
[3421]	379	ENDIF
	380	!
	381	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	382	IF ( l == 0 .OR. l == 1 ) THEN
	383	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	384	surfaces%npoints = surfaces%npoints + 1
[3494]	385	point_index(k,j,i+1) = surfaces%npoints - 1
[3421]	386	ENDIF
[3494]	387	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	388	surfaces%npoints = surfaces%npoints + 1
[3494]	389	point_index(k+1,j,i+1) = surfaces%npoints - 1
[3421]	390	ENDIF
	391	!
	392	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	393	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	394	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	395	surfaces%npoints = surfaces%npoints + 1
[3494]	396	point_index(k,j+1,i) = surfaces%npoints - 1
[3421]	397	ENDIF
[3494]	398	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	399	surfaces%npoints = surfaces%npoints + 1
[3494]	400	point_index(k+1,j+1,i) = surfaces%npoints - 1
[3421]	401	ENDIF
[3494]	402	ENDIF
	403	!
	404	!-- Upper left / front vertex
	405	IF ( point_index(k+1,j,i) < 0 ) THEN
	406	surfaces%npoints = surfaces%npoints + 1
	407	point_index(k+1,j,i) = surfaces%npoints - 1
	408	ENDIF
[3421]	409	ENDDO
[3494]	410
[3421]	411	ENDDO
	412	!
	413	!-- Allocate the number of points and polygons. Note, the number of polygons
	414	!-- is identical to the number of surfaces elements, whereas the number
	415	!-- of points (vertices) which define the polygons can be larger.
	416	ALLOCATE( surfaces%points(3,1:surfaces%npoints) )
[3494]	417	ALLOCATE( surfaces%polygons(5,1:surfaces%ns) )
[3421]	418	!
[3494]	419	!-- Note, PARAVIEW expects consecutively ordered points which can be
	420	!-- unambiguously identified.
	421	!-- Hence, all PEs know where they should start counting, depending on the
	422	!-- number of points on the other PE's with lower MPI rank.
	423	#if defined( __parallel )
	424	CALL MPI_ALLGATHER( surfaces%npoints, 1, MPI_INTEGER, &
	425	num_points_on_pe, 1, MPI_INTEGER, comm2d, ierr )
	426	#else
	427	num_points_on_pe = surfaces%npoints
	428	#endif
	429
	430	!
[3421]	431	!-- After the number of vertices is counted, repeat the loops and define the
[3494]	432	!-- vertices. Start with the horizontal default surfaces
	433	!-- First, however, determine the offset where couting of points should be
	434	!-- started, which is the sum of points of all PE's with lower MPI rank.
	435	i = 0
	436	point_index_count = 0
	437	DO WHILE ( i < myid .AND. i <= SIZE( num_points_on_pe ) )
	438	point_index_count = point_index_count + num_points_on_pe(i)
	439	i = i + 1
	440	ENDDO
	441
	442
[3421]	443	surfaces%npoints = 0
	444	point_index = -1
	445	npg = 0
[3494]	446
[3421]	447	DO l = 0, 1
	448	DO m = 1, surf_def_h(0)%ns
	449	!
	450	!-- Determine the indices of the respective grid cell inside the topography
	451	i = surf_def_h(0)%i(m) + surf_def_h(0)%ioff
	452	j = surf_def_h(0)%j(m) + surf_def_h(0)%joff
	453	k = surf_def_h(0)%k(m) + surf_def_h(0)%koff
	454	!
	455	!-- Check if the vertices that define the surface element are already
	456	!-- defined, if not, increment the counter.
[3494]	457	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	458	surfaces%npoints = surfaces%npoints + 1
[3494]	459	point_index(k,j,i) = point_index_count
	460	point_index_count = point_index_count + 1
[3421]	461	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	462	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	463	surfaces%points(3,surfaces%npoints) = zw(k)
	464	ENDIF
[3494]	465	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	466	surfaces%npoints = surfaces%npoints + 1
[3494]	467	point_index(k,j,i+1) = point_index_count
	468	point_index_count = point_index_count + 1
[3421]	469	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	470	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	471	surfaces%points(3,surfaces%npoints) = zw(k)
	472	ENDIF
[3494]	473	IF ( point_index(k,j+1,i+1) < 0 ) THEN
[3421]	474	surfaces%npoints = surfaces%npoints + 1
[3494]	475	point_index(k,j+1,i+1) = point_index_count
	476	point_index_count = point_index_count + 1
[3421]	477	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	478	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	479	surfaces%points(3,surfaces%npoints) = zw(k)
	480	ENDIF
[3494]	481	IF ( point_index(k,j+1,i) < 0 ) THEN
	482	surfaces%npoints = surfaces%npoints + 1
	483	point_index(k,j+1,i) = point_index_count
	484	point_index_count = point_index_count + 1
	485	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	486	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	487	surfaces%points(3,surfaces%npoints) = zw(k)
	488	ENDIF
[3421]	489
[3494]	490	npg = npg + 1
[3421]	491	surfaces%polygons(1,npg) = 4
	492	surfaces%polygons(2,npg) = point_index(k,j,i)
[3494]	493	surfaces%polygons(3,npg) = point_index(k,j,i+1)
[3421]	494	surfaces%polygons(4,npg) = point_index(k,j+1,i+1)
[3494]	495	surfaces%polygons(5,npg) = point_index(k,j+1,i)
[3421]	496	ENDDO
	497	ENDDO
	498	DO m = 1, surf_lsm_h%ns
	499	i = surf_lsm_h%i(m) + surf_lsm_h%ioff
	500	j = surf_lsm_h%j(m) + surf_lsm_h%joff
	501	k = surf_lsm_h%k(m) + surf_lsm_h%koff
[3494]	502	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	503	surfaces%npoints = surfaces%npoints + 1
[3494]	504	point_index(k,j,i) = point_index_count
	505	point_index_count = point_index_count + 1
[3421]	506	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	507	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	508	surfaces%points(3,surfaces%npoints) = zw(k)
	509	ENDIF
[3494]	510	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	511	surfaces%npoints = surfaces%npoints + 1
[3494]	512	point_index(k,j,i+1) = point_index_count
	513	point_index_count = point_index_count + 1
[3421]	514	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	515	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	516	surfaces%points(3,surfaces%npoints) = zw(k)
	517	ENDIF
[3494]	518	IF ( point_index(k,j+1,i+1) < 0 ) THEN
[3421]	519	surfaces%npoints = surfaces%npoints + 1
[3494]	520	point_index(k,j+1,i+1) = point_index_count
	521	point_index_count = point_index_count + 1
[3421]	522	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	523	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	524	surfaces%points(3,surfaces%npoints) = zw(k)
	525	ENDIF
[3494]	526	IF ( point_index(k,j+1,i) < 0 ) THEN
	527	surfaces%npoints = surfaces%npoints + 1
	528	point_index(k,j+1,i) = point_index_count
	529	point_index_count = point_index_count + 1
	530	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	531	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	532	surfaces%points(3,surfaces%npoints) = zw(k)
	533	ENDIF
[3421]	534
[3494]	535	npg = npg + 1
[3421]	536	surfaces%polygons(1,npg) = 4
	537	surfaces%polygons(2,npg) = point_index(k,j,i)
[3494]	538	surfaces%polygons(3,npg) = point_index(k,j,i+1)
[3421]	539	surfaces%polygons(4,npg) = point_index(k,j+1,i+1)
[3494]	540	surfaces%polygons(5,npg) = point_index(k,j+1,i)
[3421]	541	ENDDO
[3494]	542
[3421]	543	DO m = 1, surf_usm_h%ns
	544	i = surf_usm_h%i(m) + surf_usm_h%ioff
	545	j = surf_usm_h%j(m) + surf_usm_h%joff
	546	k = surf_usm_h%k(m) + surf_usm_h%koff
	547
[3494]	548	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	549	surfaces%npoints = surfaces%npoints + 1
[3494]	550	point_index(k,j,i) = point_index_count
	551	point_index_count = point_index_count + 1
[3421]	552	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	553	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	554	surfaces%points(3,surfaces%npoints) = zw(k)
	555	ENDIF
[3494]	556	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	557	surfaces%npoints = surfaces%npoints + 1
[3494]	558	point_index(k,j,i+1) = point_index_count
	559	point_index_count = point_index_count + 1
[3421]	560	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	561	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	562	surfaces%points(3,surfaces%npoints) = zw(k)
	563	ENDIF
[3494]	564	IF ( point_index(k,j+1,i+1) < 0 ) THEN
[3421]	565	surfaces%npoints = surfaces%npoints + 1
[3494]	566	point_index(k,j+1,i+1) = point_index_count
	567	point_index_count = point_index_count + 1
[3421]	568	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	569	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	570	surfaces%points(3,surfaces%npoints) = zw(k)
	571	ENDIF
[3494]	572	IF ( point_index(k,j+1,i) < 0 ) THEN
	573	surfaces%npoints = surfaces%npoints + 1
	574	point_index(k,j+1,i) = point_index_count
	575	point_index_count = point_index_count + 1
	576	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	577	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
	578	surfaces%points(3,surfaces%npoints) = zw(k)
	579	ENDIF
[3421]	580
[3494]	581	npg = npg + 1
[3421]	582	surfaces%polygons(1,npg) = 4
	583	surfaces%polygons(2,npg) = point_index(k,j,i)
[3494]	584	surfaces%polygons(3,npg) = point_index(k,j,i+1)
[3421]	585	surfaces%polygons(4,npg) = point_index(k,j+1,i+1)
[3494]	586	surfaces%polygons(5,npg) = point_index(k,j+1,i)
[3421]	587	ENDDO
[3494]	588
[3421]	589	DO l = 0, 3
	590	DO m = 1, surf_def_v(l)%ns
	591	!
	592	!-- Determine the indices of the respective grid cell inside the topography.
	593	!-- Please note, j-index for north-facing surfaces ( l==0 ) is
	594	!-- identical to the reference j-index outside the grid box.
	595	!-- Equivalent for east-facing surfaces and i-index.
[3494]	596	i = surf_def_v(l)%i(m) + MERGE( surf_def_v(l)%ioff, 0, l == 3 )
	597	j = surf_def_v(l)%j(m) + MERGE( surf_def_v(l)%joff, 0, l == 1 )
[3421]	598	k = surf_def_v(l)%k(m) + surf_def_v(l)%koff
	599	!
	600	!-- Lower left /front vertex
[3494]	601	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	602	surfaces%npoints = surfaces%npoints + 1
[3494]	603	point_index(k,j,i) = point_index_count
	604	point_index_count = point_index_count + 1
[3421]	605	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	606	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	607	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	608	ENDIF
	609	!
	610	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	611	IF ( l == 0 .OR. l == 1 ) THEN
	612	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	613	surfaces%npoints = surfaces%npoints + 1
[3494]	614	point_index(k,j,i+1) = point_index_count
	615	point_index_count = point_index_count + 1
[3421]	616	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	617	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	618	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	619	ENDIF
[3494]	620	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	621	surfaces%npoints = surfaces%npoints + 1
[3494]	622	point_index(k+1,j,i+1) = point_index_count
	623	point_index_count = point_index_count + 1
[3421]	624	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	625	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	626	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	627	ENDIF
	628	!
	629	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	630	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	631	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	632	surfaces%npoints = surfaces%npoints + 1
[3494]	633	point_index(k,j+1,i) = point_index_count
	634	point_index_count = point_index_count + 1
[3421]	635	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	636	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	637	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	638	ENDIF
[3494]	639	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	640	surfaces%npoints = surfaces%npoints + 1
[3494]	641	point_index(k+1,j+1,i) = point_index_count
	642	point_index_count = point_index_count + 1
[3421]	643	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	644	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	645	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	646	ENDIF
[3494]	647	ENDIF
	648	!
	649	!-- Upper left / front vertex
	650	IF ( point_index(k+1,j,i) < 0 ) THEN
	651	surfaces%npoints = surfaces%npoints + 1
	652	point_index(k+1,j,i) = point_index_count
	653	point_index_count = point_index_count + 1
	654	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	655	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	656	surfaces%points(3,surfaces%npoints) = zw(k)
	657	ENDIF
[3421]	658
	659	npg = npg + 1
[3494]	660	IF ( l == 0 .OR. l == 1 ) THEN
[3421]	661	surfaces%polygons(1,npg) = 4
	662	surfaces%polygons(2,npg) = point_index(k,j,i)
	663	surfaces%polygons(3,npg) = point_index(k,j,i+1)
	664	surfaces%polygons(4,npg) = point_index(k+1,j,i+1)
	665	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	666	ELSE
	667	surfaces%polygons(1,npg) = 4
	668	surfaces%polygons(2,npg) = point_index(k,j,i)
	669	surfaces%polygons(3,npg) = point_index(k,j+1,i)
	670	surfaces%polygons(4,npg) = point_index(k+1,j+1,i)
	671	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	672	ENDIF
	673
	674	ENDDO
[3494]	675
[3421]	676	DO m = 1, surf_lsm_v(l)%ns
[3494]	677	i = surf_lsm_v(l)%i(m) + MERGE( surf_lsm_v(l)%ioff, 0, l == 3 )
	678	j = surf_lsm_v(l)%j(m) + MERGE( surf_lsm_v(l)%joff, 0, l == 1 )
[3421]	679	k = surf_lsm_v(l)%k(m) + surf_lsm_v(l)%koff
	680	!
	681	!-- Lower left /front vertex
[3494]	682	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	683	surfaces%npoints = surfaces%npoints + 1
[3494]	684	point_index(k,j,i) = point_index_count
	685	point_index_count = point_index_count + 1
[3421]	686	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	687	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	688	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	689	ENDIF
	690	!
	691	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	692	IF ( l == 0 .OR. l == 1 ) THEN
	693	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	694	surfaces%npoints = surfaces%npoints + 1
[3494]	695	point_index(k,j,i+1) = point_index_count
	696	point_index_count = point_index_count + 1
[3421]	697	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	698	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	699	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	700	ENDIF
[3494]	701	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	702	surfaces%npoints = surfaces%npoints + 1
[3494]	703	point_index(k+1,j,i+1) = point_index_count
	704	point_index_count = point_index_count + 1
[3421]	705	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	706	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	707	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	708	ENDIF
	709	!
	710	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	711	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	712	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	713	surfaces%npoints = surfaces%npoints + 1
[3494]	714	point_index(k,j+1,i) = point_index_count
	715	point_index_count = point_index_count + 1
[3421]	716	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	717	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	718	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	719	ENDIF
[3494]	720	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	721	surfaces%npoints = surfaces%npoints + 1
[3494]	722	point_index(k+1,j+1,i) = point_index_count
	723	point_index_count = point_index_count + 1
[3421]	724	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	725	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	726	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	727	ENDIF
	728	ENDIF
[3494]	729	!
	730	!-- Upper left / front vertex
	731	IF ( point_index(k+1,j,i) < 0 ) THEN
	732	surfaces%npoints = surfaces%npoints + 1
	733	point_index(k+1,j,i) = point_index_count
	734	point_index_count = point_index_count + 1
	735	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	736	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	737	surfaces%points(3,surfaces%npoints) = zw(k)
	738	ENDIF
[3421]	739
	740	npg = npg + 1
[3494]	741	IF ( l == 0 .OR. l == 1 ) THEN
[3421]	742	surfaces%polygons(1,npg) = 4
	743	surfaces%polygons(2,npg) = point_index(k,j,i)
	744	surfaces%polygons(3,npg) = point_index(k,j,i+1)
	745	surfaces%polygons(4,npg) = point_index(k+1,j,i+1)
	746	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	747	ELSE
	748	surfaces%polygons(1,npg) = 4
	749	surfaces%polygons(2,npg) = point_index(k,j,i)
	750	surfaces%polygons(3,npg) = point_index(k,j+1,i)
	751	surfaces%polygons(4,npg) = point_index(k+1,j+1,i)
	752	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	753	ENDIF
	754	ENDDO
	755	DO m = 1, surf_usm_v(l)%ns
[3494]	756	i = surf_usm_v(l)%i(m) + MERGE( surf_usm_v(l)%ioff, 0, l == 3 )
	757	j = surf_usm_v(l)%j(m) + MERGE( surf_usm_v(l)%joff, 0, l == 1 )
[3421]	758	k = surf_usm_v(l)%k(m) + surf_usm_v(l)%koff
	759	!
	760	!-- Lower left /front vertex
[3494]	761	IF ( point_index(k,j,i) < 0 ) THEN
[3421]	762	surfaces%npoints = surfaces%npoints + 1
[3494]	763	point_index(k,j,i) = point_index_count
	764	point_index_count = point_index_count + 1
[3421]	765	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	766	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	767	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	768	ENDIF
	769	!
	770	!-- Upper / lower right index for north- and south-facing surfaces
[3494]	771	IF ( l == 0 .OR. l == 1 ) THEN
	772	IF ( point_index(k,j,i+1) < 0 ) THEN
[3421]	773	surfaces%npoints = surfaces%npoints + 1
[3494]	774	point_index(k,j,i+1) = point_index_count
	775	point_index_count = point_index_count + 1
[3421]	776	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	777	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	778	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	779	ENDIF
[3494]	780	IF ( point_index(k+1,j,i+1) < 0 ) THEN
[3421]	781	surfaces%npoints = surfaces%npoints + 1
[3494]	782	point_index(k+1,j,i+1) = point_index_count
	783	point_index_count = point_index_count + 1
[3421]	784	surfaces%points(1,surfaces%npoints) = ( i + 1 - 0.5_wp ) * dx
	785	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
[3494]	786	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	787	ENDIF
	788	!
	789	!-- Upper / lower front index for east- and west-facing surfaces
[3494]	790	ELSEIF ( l == 2 .OR. l == 3 ) THEN
	791	IF ( point_index(k,j+1,i) < 0 ) THEN
[3421]	792	surfaces%npoints = surfaces%npoints + 1
[3494]	793	point_index(k,j+1,i) = point_index_count
	794	point_index_count = point_index_count + 1
[3421]	795	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	796	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	797	surfaces%points(3,surfaces%npoints) = zw(k-1)
[3421]	798	ENDIF
[3494]	799	IF ( point_index(k+1,j+1,i) < 0 ) THEN
[3421]	800	surfaces%npoints = surfaces%npoints + 1
[3494]	801	point_index(k+1,j+1,i) = point_index_count
	802	point_index_count = point_index_count + 1
[3421]	803	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	804	surfaces%points(2,surfaces%npoints) = ( j + 1 - 0.5_wp ) * dy
[3494]	805	surfaces%points(3,surfaces%npoints) = zw(k)
[3421]	806	ENDIF
	807	ENDIF
[3494]	808	!
	809	!-- Upper left / front vertex
	810	IF ( point_index(k+1,j,i) < 0 ) THEN
	811	surfaces%npoints = surfaces%npoints + 1
	812	point_index(k+1,j,i) = point_index_count
	813	point_index_count = point_index_count + 1
	814	surfaces%points(1,surfaces%npoints) = ( i - 0.5_wp ) * dx
	815	surfaces%points(2,surfaces%npoints) = ( j - 0.5_wp ) * dy
	816	surfaces%points(3,surfaces%npoints) = zw(k)
	817	ENDIF
[3421]	818
	819	npg = npg + 1
[3494]	820	IF ( l == 0 .OR. l == 1 ) THEN
[3421]	821	surfaces%polygons(1,npg) = 4
	822	surfaces%polygons(2,npg) = point_index(k,j,i)
	823	surfaces%polygons(3,npg) = point_index(k,j,i+1)
	824	surfaces%polygons(4,npg) = point_index(k+1,j,i+1)
	825	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	826	ELSE
	827	surfaces%polygons(1,npg) = 4
	828	surfaces%polygons(2,npg) = point_index(k,j,i)
	829	surfaces%polygons(3,npg) = point_index(k,j+1,i)
	830	surfaces%polygons(4,npg) = point_index(k+1,j+1,i)
	831	surfaces%polygons(5,npg) = point_index(k+1,j,i)
	832	ENDIF
	833	ENDDO
[3494]	834
[3421]	835	ENDDO
	836	!
	837	!-- Deallocate temporary dummy variable
	838	DEALLOCATE ( point_index )
	839	!
[3494]	840	!-- Sum-up total number of surface elements and vertices on domain. This
	841	!-- will be needed for post-processing.
	842	surfaces%npoints_total = 0
	843	#if defined( __parallel )
	844	CALL MPI_ALLREDUCE( surfaces%npoints, surfaces%npoints_total, 1, &
	845	MPI_INTEGER, MPI_SUM, comm2d, ierr )
	846	CALL MPI_ALLREDUCE( surfaces%ns, surfaces%ns_total, 1, &
	847	MPI_INTEGER, MPI_SUM, comm2d, ierr )
	848	#else
	849	surfaces%npoints_total = surfaces%npoints
	850	surfaces%ns_total = surfaces%ns
	851	#endif
	852
[3421]	853	END SUBROUTINE surface_output_init
	854
	855	!------------------------------------------------------------------------------!
	856	! Description:
	857	! ------------
	858	!> Routine for controlling the data output. Surface data is collected from
	859	!> different types of surfaces (default, natural, urban) and different
	860	!> orientation and written to one 1D-output array. Further, NetCDF routines
	861	!> are called to write the surface data in the respective NetCDF files.
	862	!------------------------------------------------------------------------------!
	863	SUBROUTINE surface_output( av )
	864
	865	USE control_parameters, &
	866	ONLY: io_blocks, io_group, simulated_time
	867
	868	USE pegrid, &
	869	ONLY: comm2d, ierr
	870
	871
	872	IMPLICIT NONE
	873
	874	CHARACTER(LEN=100) :: trimvar = ' ' !< dummy for single output variable
	875
	876	INTEGER(iwp) :: av !< id indicating average or non-average data output
	877	INTEGER(iwp) :: i !< loop index
	878	INTEGER(iwp) :: n_out !< counter variables for surface output
	879
	880	!
	881	!-- Return, if nothing to output
	882	IF ( dosurf_no(av) == 0 ) RETURN
	883	!
	884	!-- Open files
	885	CALL check_open( 25+av )
	886	!
	887	!-- Write coordinates
	888	IF ( .NOT. first_output(av) ) THEN
	889	DO i = 0, io_blocks-1
	890	IF ( i == io_group ) THEN
	891	WRITE ( 25+av ) surfaces%npoints
[3494]	892	WRITE ( 25+av ) surfaces%npoints_total
[3421]	893	WRITE ( 25+av ) surfaces%ns
[3494]	894	WRITE ( 25+av ) surfaces%ns_total
[3421]	895	WRITE ( 25+av ) surfaces%points
	896	WRITE ( 25+av ) surfaces%polygons
	897	ENDIF
[3483]	898	#if defined( __parallel )
[3421]	899	CALL MPI_BARRIER( comm2d, ierr )
[3483]	900	#endif
[3421]	901	first_output(av) = .TRUE.
	902	ENDDO
	903	ENDIF
	904
	905	n_out = 0
	906	DO WHILE ( dosurf(av,n_out+1)(1:1) /= ' ' )
	907
	908	n_out = n_out + 1
	909	trimvar = TRIM( dosurf(av,n_out) )
	910	!
	911	!-- Set the output array to the _FillValue in case it is not
	912	!-- defined for each type of surface.
	913	surfaces%var_out = surfaces%fillvalue
	914	SELECT CASE ( trimvar )
	915
	916	CASE ( 'us' )
	917	!
	918	!-- Output of instantaneous data
	919	IF ( av == 0 ) THEN
	920	CALL surface_output_collect( surf_def_h(0)%us, &
[3494]	921	surf_def_h(1)%us, &
	922	surf_lsm_h%us, &
	923	surf_usm_h%us, &
	924	surf_def_v(0)%us, &
	925	surf_lsm_v(0)%us, &
	926	surf_usm_v(0)%us, &
	927	surf_def_v(1)%us, &
	928	surf_lsm_v(1)%us, &
	929	surf_usm_v(1)%us, &
	930	surf_def_v(2)%us, &
	931	surf_lsm_v(2)%us, &
	932	surf_usm_v(2)%us, &
	933	surf_def_v(3)%us, &
	934	surf_lsm_v(3)%us, &
	935	surf_usm_v(3)%us )
[3421]	936	ELSE
	937	!
	938	!-- Output of averaged data
	939	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	940	REAL( average_count_surf, KIND=wp )
	941	surfaces%var_av(:,n_out) = 0.0_wp
	942
	943	ENDIF
	944
	945	CASE ( 'ts' )
	946	!
	947	!-- Output of instantaneous data
	948	IF ( av == 0 ) THEN
	949	CALL surface_output_collect( surf_def_h(0)%ts, &
	950	surf_def_h(1)%ts, &
	951	surf_lsm_h%ts, &
	952	surf_usm_h%ts, &
	953	surf_def_v(0)%ts, &
	954	surf_lsm_v(0)%ts, &
	955	surf_usm_v(0)%ts, &
	956	surf_def_v(1)%ts, &
	957	surf_lsm_v(1)%ts, &
	958	surf_usm_v(1)%ts, &
	959	surf_def_v(2)%ts, &
	960	surf_lsm_v(2)%ts, &
	961	surf_usm_v(2)%ts, &
	962	surf_def_v(3)%ts, &
	963	surf_lsm_v(3)%ts, &
	964	surf_usm_v(3)%ts )
	965	ELSE
	966	!
	967	!-- Output of averaged data
	968	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	969	REAL( average_count_surf, KIND=wp )
	970	surfaces%var_av(:,n_out) = 0.0_wp
	971
	972	ENDIF
	973
	974	CASE ( 'qs' )
	975	!
	976	!-- Output of instantaneous data
	977	IF ( av == 0 ) THEN
	978	CALL surface_output_collect( surf_def_h(0)%qs, &
	979	surf_def_h(1)%qs, &
	980	surf_lsm_h%qs, &
	981	surf_usm_h%qs, &
	982	surf_def_v(0)%qs, &
	983	surf_lsm_v(0)%qs, &
	984	surf_usm_v(0)%qs, &
	985	surf_def_v(1)%qs, &
	986	surf_lsm_v(1)%qs, &
	987	surf_usm_v(1)%qs, &
	988	surf_def_v(2)%qs, &
	989	surf_lsm_v(2)%qs, &
	990	surf_usm_v(2)%qs, &
	991	surf_def_v(3)%qs, &
	992	surf_lsm_v(3)%qs, &
	993	surf_usm_v(3)%qs )
	994	ELSE
	995	!
	996	!-- Output of averaged data
	997	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	998	REAL( average_count_surf, KIND=wp )
	999	surfaces%var_av(:,n_out) = 0.0_wp
	1000
	1001	ENDIF
	1002
	1003	CASE ( 'ss' )
	1004	!
	1005	!-- Output of instantaneous data
	1006	IF ( av == 0 ) THEN
	1007	CALL surface_output_collect( surf_def_h(0)%ss, &
	1008	surf_def_h(1)%ss, &
	1009	surf_lsm_h%ss, &
	1010	surf_usm_h%ss, &
	1011	surf_def_v(0)%ss, &
	1012	surf_lsm_v(0)%ss, &
	1013	surf_usm_v(0)%ss, &
	1014	surf_def_v(1)%ss, &
	1015	surf_lsm_v(1)%ss, &
	1016	surf_usm_v(1)%ss, &
	1017	surf_def_v(2)%ss, &
	1018	surf_lsm_v(2)%ss, &
	1019	surf_usm_v(2)%ss, &
	1020	surf_def_v(3)%ss, &
	1021	surf_lsm_v(3)%ss, &
	1022	surf_usm_v(3)%ss )
	1023	ELSE
	1024	!
	1025	!-- Output of averaged data
	1026	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1027	REAL( average_count_surf, KIND=wp )
	1028	surfaces%var_av(:,n_out) = 0.0_wp
	1029
	1030	ENDIF
	1031
	1032	CASE ( 'qcs' )
	1033	!
	1034	!-- Output of instantaneous data
	1035	IF ( av == 0 ) THEN
	1036	CALL surface_output_collect( surf_def_h(0)%qcs, &
	1037	surf_def_h(1)%qcs, &
	1038	surf_lsm_h%qcs, &
	1039	surf_usm_h%qcs, &
	1040	surf_def_v(0)%qcs, &
	1041	surf_lsm_v(0)%qcs, &
	1042	surf_usm_v(0)%qcs, &
	1043	surf_def_v(1)%qcs, &
	1044	surf_lsm_v(1)%qcs, &
	1045	surf_usm_v(1)%qcs, &
	1046	surf_def_v(2)%qcs, &
	1047	surf_lsm_v(2)%qcs, &
	1048	surf_usm_v(2)%qcs, &
	1049	surf_def_v(3)%qcs, &
	1050	surf_lsm_v(3)%qcs, &
	1051	surf_usm_v(3)%qcs )
	1052	ELSE
	1053	!
	1054	!-- Output of averaged data
	1055	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1056	REAL( average_count_surf, KIND=wp )
	1057	surfaces%var_av(:,n_out) = 0.0_wp
	1058
	1059	ENDIF
	1060
	1061	CASE ( 'ncs' )
	1062	!
	1063	!-- Output of instantaneous data
	1064	IF ( av == 0 ) THEN
	1065	CALL surface_output_collect( surf_def_h(0)%ncs, &
	1066	surf_def_h(1)%ncs, &
	1067	surf_lsm_h%ncs, &
	1068	surf_usm_h%ncs, &
	1069	surf_def_v(0)%ncs, &
	1070	surf_lsm_v(0)%ncs, &
	1071	surf_usm_v(0)%ncs, &
	1072	surf_def_v(1)%ncs, &
	1073	surf_lsm_v(1)%ncs, &
	1074	surf_usm_v(1)%ncs, &
	1075	surf_def_v(2)%ncs, &
	1076	surf_lsm_v(2)%ncs, &
	1077	surf_usm_v(2)%ncs, &
	1078	surf_def_v(3)%ncs, &
	1079	surf_lsm_v(3)%ncs, &
	1080	surf_usm_v(3)%ncs )
	1081	ELSE
	1082	!
	1083	!-- Output of averaged data
	1084	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1085	REAL( average_count_surf, KIND=wp )
	1086	surfaces%var_av(:,n_out) = 0.0_wp
	1087
	1088	ENDIF
	1089
	1090	CASE ( 'qrs' )
	1091	!
	1092	!-- Output of instantaneous data
	1093	IF ( av == 0 ) THEN
	1094	CALL surface_output_collect( surf_def_h(0)%qrs, &
	1095	surf_def_h(1)%qrs, &
	1096	surf_lsm_h%qrs, &
	1097	surf_usm_h%qrs, &
	1098	surf_def_v(0)%qrs, &
	1099	surf_lsm_v(0)%qrs, &
	1100	surf_usm_v(0)%qrs, &
	1101	surf_def_v(1)%qrs, &
	1102	surf_lsm_v(1)%qrs, &
	1103	surf_usm_v(1)%qrs, &
	1104	surf_def_v(2)%qrs, &
	1105	surf_lsm_v(2)%qrs, &
	1106	surf_usm_v(2)%qrs, &
	1107	surf_def_v(3)%qrs, &
	1108	surf_lsm_v(3)%qrs, &
	1109	surf_usm_v(3)%qrs )
	1110	ELSE
	1111	!
	1112	!-- Output of averaged data
	1113	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1114	REAL( average_count_surf, KIND=wp )
	1115	surfaces%var_av(:,n_out) = 0.0_wp
	1116
	1117	ENDIF
	1118
	1119	CASE ( 'nrs' )
	1120	!
	1121	!-- Output of instantaneous data
	1122	IF ( av == 0 ) THEN
	1123	CALL surface_output_collect( surf_def_h(0)%nrs, &
	1124	surf_def_h(1)%nrs, &
	1125	surf_lsm_h%nrs, &
	1126	surf_usm_h%nrs, &
	1127	surf_def_v(0)%nrs, &
	1128	surf_lsm_v(0)%nrs, &
	1129	surf_usm_v(0)%nrs, &
	1130	surf_def_v(1)%nrs, &
	1131	surf_lsm_v(1)%nrs, &
	1132	surf_usm_v(1)%nrs, &
	1133	surf_def_v(2)%nrs, &
	1134	surf_lsm_v(2)%nrs, &
	1135	surf_usm_v(2)%nrs, &
	1136	surf_def_v(3)%nrs, &
	1137	surf_lsm_v(3)%nrs, &
	1138	surf_usm_v(3)%nrs )
	1139	ELSE
	1140	!
	1141	!-- Output of averaged data
	1142	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1143	REAL( average_count_surf, KIND=wp )
	1144	surfaces%var_av(:,n_out) = 0.0_wp
	1145
	1146	ENDIF
	1147
	1148	CASE ( 'ol' )
	1149	!
	1150	!-- Output of instantaneous data
	1151	IF ( av == 0 ) THEN
	1152	CALL surface_output_collect( surf_def_h(0)%ol, &
	1153	surf_def_h(1)%ol, &
	1154	surf_lsm_h%ol, &
	1155	surf_usm_h%ol, &
	1156	surf_def_v(0)%ol, &
	1157	surf_lsm_v(0)%ol, &
	1158	surf_usm_v(0)%ol, &
	1159	surf_def_v(1)%ol, &
	1160	surf_lsm_v(1)%ol, &
	1161	surf_usm_v(1)%ol, &
	1162	surf_def_v(2)%ol, &
	1163	surf_lsm_v(2)%ol, &
	1164	surf_usm_v(2)%ol, &
	1165	surf_def_v(3)%ol, &
	1166	surf_lsm_v(3)%ol, &
	1167	surf_usm_v(3)%ol )
	1168	ELSE
	1169	!
	1170	!-- Output of averaged data
	1171	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1172	REAL( average_count_surf, KIND=wp )
	1173	surfaces%var_av(:,n_out) = 0.0_wp
	1174
	1175	ENDIF
	1176
	1177	CASE ( 'z0' )
	1178	!
	1179	!-- Output of instantaneous data
	1180	IF ( av == 0 ) THEN
	1181	CALL surface_output_collect( surf_def_h(0)%z0, &
	1182	surf_def_h(1)%z0, &
	1183	surf_lsm_h%z0, &
	1184	surf_usm_h%z0, &
	1185	surf_def_v(0)%z0, &
	1186	surf_lsm_v(0)%z0, &
	1187	surf_usm_v(0)%z0, &
	1188	surf_def_v(1)%z0, &
	1189	surf_lsm_v(1)%z0, &
	1190	surf_usm_v(1)%z0, &
	1191	surf_def_v(2)%z0, &
	1192	surf_lsm_v(2)%z0, &
	1193	surf_usm_v(2)%z0, &
	1194	surf_def_v(3)%z0, &
	1195	surf_lsm_v(3)%z0, &
	1196	surf_usm_v(3)%z0 )
	1197	ELSE
	1198	!
	1199	!-- Output of averaged data
	1200	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1201	REAL( average_count_surf, KIND=wp )
	1202	surfaces%var_av(:,n_out) = 0.0_wp
	1203
	1204	ENDIF
	1205
	1206	CASE ( 'z0h' )
	1207	!
	1208	!-- Output of instantaneous data
	1209	IF ( av == 0 ) THEN
	1210	CALL surface_output_collect( surf_def_h(0)%z0h, &
	1211	surf_def_h(1)%z0h, &
	1212	surf_lsm_h%z0h, &
	1213	surf_usm_h%z0h, &
	1214	surf_def_v(0)%z0h, &
	1215	surf_lsm_v(0)%z0h, &
	1216	surf_usm_v(0)%z0h, &
	1217	surf_def_v(1)%z0h, &
	1218	surf_lsm_v(1)%z0h, &
	1219	surf_usm_v(1)%z0h, &
	1220	surf_def_v(2)%z0h, &
	1221	surf_lsm_v(2)%z0h, &
	1222	surf_usm_v(2)%z0h, &
	1223	surf_def_v(3)%z0h, &
	1224	surf_lsm_v(3)%z0h, &
	1225	surf_usm_v(3)%z0h )
	1226	ELSE
	1227	!
	1228	!-- Output of averaged data
	1229	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1230	REAL( average_count_surf, KIND=wp )
	1231	surfaces%var_av(:,n_out) = 0.0_wp
	1232
	1233	ENDIF
	1234
	1235	CASE ( 'z0q' )
	1236	!
	1237	!-- Output of instantaneous data
	1238	IF ( av == 0 ) THEN
	1239	CALL surface_output_collect( surf_def_h(0)%z0q, &
	1240	surf_def_h(1)%z0q, &
	1241	surf_lsm_h%z0q, &
	1242	surf_usm_h%z0q, &
	1243	surf_def_v(0)%z0q, &
	1244	surf_lsm_v(0)%z0q, &
	1245	surf_usm_v(0)%z0q, &
	1246	surf_def_v(1)%z0q, &
	1247	surf_lsm_v(1)%z0q, &
	1248	surf_usm_v(1)%z0q, &
	1249	surf_def_v(2)%z0q, &
	1250	surf_lsm_v(2)%z0q, &
	1251	surf_usm_v(2)%z0q, &
	1252	surf_def_v(3)%z0q, &
	1253	surf_lsm_v(3)%z0q, &
	1254	surf_usm_v(3)%z0q )
	1255	ELSE
	1256	!
	1257	!-- Output of averaged data
	1258	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1259	REAL( average_count_surf, KIND=wp )
	1260	surfaces%var_av(:,n_out) = 0.0_wp
	1261
	1262	ENDIF
	1263
	1264	CASE ( 'theta1' )
	1265	!
	1266	!-- Output of instantaneous data
	1267	IF ( av == 0 ) THEN
	1268	CALL surface_output_collect( surf_def_h(0)%pt1, &
	1269	surf_def_h(1)%pt1, &
	1270	surf_lsm_h%pt1, &
	1271	surf_usm_h%pt1, &
	1272	surf_def_v(0)%pt1, &
	1273	surf_lsm_v(0)%pt1, &
	1274	surf_usm_v(0)%pt1, &
	1275	surf_def_v(1)%pt1, &
	1276	surf_lsm_v(1)%pt1, &
	1277	surf_usm_v(1)%pt1, &
	1278	surf_def_v(2)%pt1, &
	1279	surf_lsm_v(2)%pt1, &
	1280	surf_usm_v(2)%pt1, &
	1281	surf_def_v(3)%pt1, &
	1282	surf_lsm_v(3)%pt1, &
	1283	surf_usm_v(3)%pt1 )
	1284	ELSE
	1285	!
	1286	!-- Output of averaged data
	1287	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1288	REAL( average_count_surf, KIND=wp )
	1289	surfaces%var_av(:,n_out) = 0.0_wp
	1290
	1291	ENDIF
	1292
	1293	CASE ( 'qv1' )
	1294	!
	1295	!-- Output of instantaneous data
	1296	IF ( av == 0 ) THEN
	1297	CALL surface_output_collect( surf_def_h(0)%qv1, &
	1298	surf_def_h(1)%qv1, &
	1299	surf_lsm_h%qv1, &
	1300	surf_usm_h%qv1, &
	1301	surf_def_v(0)%qv1, &
	1302	surf_lsm_v(0)%qv1, &
	1303	surf_usm_v(0)%qv1, &
	1304	surf_def_v(1)%qv1, &
	1305	surf_lsm_v(1)%qv1, &
	1306	surf_usm_v(1)%qv1, &
	1307	surf_def_v(2)%qv1, &
	1308	surf_lsm_v(2)%qv1, &
	1309	surf_usm_v(2)%qv1, &
	1310	surf_def_v(3)%qv1, &
	1311	surf_lsm_v(3)%qv1, &
	1312	surf_usm_v(3)%qv1 )
	1313	ELSE
	1314	!
	1315	!-- Output of averaged data
	1316	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1317	REAL( average_count_surf, KIND=wp )
	1318	surfaces%var_av(:,n_out) = 0.0_wp
	1319
	1320	ENDIF
	1321
	1322	CASE ( 'thetav1' )
	1323	!
	1324	!-- Output of instantaneous data
	1325	IF ( av == 0 ) THEN
	1326	CALL surface_output_collect( surf_def_h(0)%vpt1, &
	1327	surf_def_h(1)%vpt1, &
	1328	surf_lsm_h%vpt1, &
	1329	surf_usm_h%vpt1, &
	1330	surf_def_v(0)%vpt1, &
	1331	surf_lsm_v(0)%vpt1, &
	1332	surf_usm_v(0)%vpt1, &
	1333	surf_def_v(1)%vpt1, &
	1334	surf_lsm_v(1)%vpt1, &
	1335	surf_usm_v(1)%vpt1, &
	1336	surf_def_v(2)%vpt1, &
	1337	surf_lsm_v(2)%vpt1, &
	1338	surf_usm_v(2)%vpt1, &
	1339	surf_def_v(3)%vpt1, &
	1340	surf_lsm_v(3)%vpt1, &
	1341	surf_usm_v(3)%vpt1 )
	1342	ELSE
	1343	!
	1344	!-- Output of averaged data
	1345	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1346	REAL( average_count_surf, KIND=wp )
	1347	surfaces%var_av(:,n_out) = 0.0_wp
	1348
	1349	ENDIF
	1350
	1351	CASE ( 'usws' )
	1352	!
	1353	!-- Output of instantaneous data
	1354	IF ( av == 0 ) THEN
	1355	CALL surface_output_collect( surf_def_h(0)%usws, &
	1356	surf_def_h(1)%usws, &
	1357	surf_lsm_h%usws, &
	1358	surf_usm_h%usws, &
	1359	surf_def_v(0)%usws, &
	1360	surf_lsm_v(0)%usws, &
	1361	surf_usm_v(0)%usws, &
	1362	surf_def_v(1)%usws, &
	1363	surf_lsm_v(1)%usws, &
	1364	surf_usm_v(1)%usws, &
	1365	surf_def_v(2)%usws, &
	1366	surf_lsm_v(2)%usws, &
	1367	surf_usm_v(2)%usws, &
	1368	surf_def_v(3)%usws, &
	1369	surf_lsm_v(3)%usws, &
	1370	surf_usm_v(3)%usws )
	1371	ELSE
	1372	!
	1373	!-- Output of averaged data
	1374	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1375	REAL( average_count_surf, KIND=wp )
	1376	surfaces%var_av(:,n_out) = 0.0_wp
	1377
	1378	ENDIF
	1379
	1380	CASE ( 'vsws' )
	1381	!
	1382	!-- Output of instantaneous data
	1383	IF ( av == 0 ) THEN
	1384	CALL surface_output_collect( surf_def_h(0)%vsws, &
	1385	surf_def_h(1)%vsws, &
	1386	surf_lsm_h%vsws, &
	1387	surf_usm_h%vsws, &
	1388	surf_def_v(0)%vsws, &
	1389	surf_lsm_v(0)%vsws, &
	1390	surf_usm_v(0)%vsws, &
	1391	surf_def_v(1)%vsws, &
	1392	surf_lsm_v(1)%vsws, &
	1393	surf_usm_v(1)%vsws, &
	1394	surf_def_v(2)%vsws, &
	1395	surf_lsm_v(2)%vsws, &
	1396	surf_usm_v(2)%vsws, &
	1397	surf_def_v(3)%vsws, &
	1398	surf_lsm_v(3)%vsws, &
	1399	surf_usm_v(3)%vsws )
	1400	ELSE
	1401	!
	1402	!-- Output of averaged data
	1403	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1404	REAL( average_count_surf, KIND=wp )
	1405	surfaces%var_av(:,n_out) = 0.0_wp
	1406
	1407	ENDIF
	1408
	1409	CASE ( 'shf' )
	1410	!
	1411	!-- Output of instantaneous data
	1412	IF ( av == 0 ) THEN
	1413	CALL surface_output_collect( surf_def_h(0)%shf, &
	1414	surf_def_h(1)%shf, &
	1415	surf_lsm_h%shf, &
	1416	surf_usm_h%shf, &
	1417	surf_def_v(0)%shf, &
	1418	surf_lsm_v(0)%shf, &
	1419	surf_usm_v(0)%shf, &
	1420	surf_def_v(1)%shf, &
	1421	surf_lsm_v(1)%shf, &
	1422	surf_usm_v(1)%shf, &
	1423	surf_def_v(2)%shf, &
	1424	surf_lsm_v(2)%shf, &
	1425	surf_usm_v(2)%shf, &
	1426	surf_def_v(3)%shf, &
	1427	surf_lsm_v(3)%shf, &
	1428	surf_usm_v(3)%shf )
	1429	ELSE
	1430	!
	1431	!-- Output of averaged data
	1432	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1433	REAL( average_count_surf, KIND=wp )
	1434	surfaces%var_av(:,n_out) = 0.0_wp
	1435	ENDIF
	1436
	1437	CASE ( 'qsws' )
	1438	!
	1439	!-- Output of instantaneous data
	1440	IF ( av == 0 ) THEN
	1441	CALL surface_output_collect( surf_def_h(0)%qsws, &
	1442	surf_def_h(1)%qsws, &
	1443	surf_lsm_h%qsws, &
	1444	surf_usm_h%qsws, &
	1445	surf_def_v(0)%qsws, &
	1446	surf_lsm_v(0)%qsws, &
	1447	surf_usm_v(0)%qsws, &
	1448	surf_def_v(1)%qsws, &
	1449	surf_lsm_v(1)%qsws, &
	1450	surf_usm_v(1)%qsws, &
	1451	surf_def_v(2)%qsws, &
	1452	surf_lsm_v(2)%qsws, &
	1453	surf_usm_v(2)%qsws, &
	1454	surf_def_v(3)%qsws, &
	1455	surf_lsm_v(3)%qsws, &
	1456	surf_usm_v(3)%qsws )
	1457	ELSE
	1458	!
	1459	!-- Output of averaged data
	1460	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1461	REAL( average_count_surf, KIND=wp )
	1462	surfaces%var_av(:,n_out) = 0.0_wp
	1463
	1464	ENDIF
	1465
	1466	CASE ( 'ssws' )
	1467	!
	1468	!-- Output of instantaneous data
	1469	IF ( av == 0 ) THEN
	1470	CALL surface_output_collect( surf_def_h(0)%ssws, &
	1471	surf_def_h(1)%ssws, &
	1472	surf_lsm_h%ssws, &
	1473	surf_usm_h%ssws, &
	1474	surf_def_v(0)%ssws, &
	1475	surf_lsm_v(0)%ssws, &
	1476	surf_usm_v(0)%ssws, &
	1477	surf_def_v(1)%ssws, &
	1478	surf_lsm_v(1)%ssws, &
	1479	surf_usm_v(1)%ssws, &
	1480	surf_def_v(2)%ssws, &
	1481	surf_lsm_v(2)%ssws, &
	1482	surf_usm_v(2)%ssws, &
	1483	surf_def_v(3)%ssws, &
	1484	surf_lsm_v(3)%ssws, &
	1485	surf_usm_v(3)%ssws )
	1486	ELSE
	1487	!
	1488	!-- Output of averaged data
	1489	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1490	REAL( average_count_surf, KIND=wp )
	1491	surfaces%var_av(:,n_out) = 0.0_wp
	1492
	1493	ENDIF
	1494
	1495	CASE ( 'qcsws' )
	1496	!
	1497	!-- Output of instantaneous data
	1498	IF ( av == 0 ) THEN
	1499	CALL surface_output_collect( surf_def_h(0)%qcsws, &
	1500	surf_def_h(1)%qcsws, &
	1501	surf_lsm_h%qcsws, &
	1502	surf_usm_h%qcsws, &
	1503	surf_def_v(0)%qcsws, &
	1504	surf_lsm_v(0)%qcsws, &
	1505	surf_usm_v(0)%qcsws, &
	1506	surf_def_v(1)%qcsws, &
	1507	surf_lsm_v(1)%qcsws, &
	1508	surf_usm_v(1)%qcsws, &
	1509	surf_def_v(2)%qcsws, &
	1510	surf_lsm_v(2)%qcsws, &
	1511	surf_usm_v(2)%qcsws, &
	1512	surf_def_v(3)%qcsws, &
	1513	surf_lsm_v(3)%qcsws, &
	1514	surf_usm_v(3)%qcsws )
	1515	ELSE
	1516	!
	1517	!-- Output of averaged data
	1518	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1519	REAL( average_count_surf, KIND=wp )
	1520	surfaces%var_av(:,n_out) = 0.0_wp
	1521
	1522	ENDIF
	1523
	1524	CASE ( 'ncsws' )
	1525	!
	1526	!-- Output of instantaneous data
	1527	IF ( av == 0 ) THEN
	1528	CALL surface_output_collect( surf_def_h(0)%ncsws, &
	1529	surf_def_h(1)%ncsws, &
	1530	surf_lsm_h%ncsws, &
	1531	surf_usm_h%ncsws, &
	1532	surf_def_v(0)%ncsws, &
	1533	surf_lsm_v(0)%ncsws, &
	1534	surf_usm_v(0)%ncsws, &
	1535	surf_def_v(1)%ncsws, &
	1536	surf_lsm_v(1)%ncsws, &
	1537	surf_usm_v(1)%ncsws, &
	1538	surf_def_v(2)%ncsws, &
	1539	surf_lsm_v(2)%ncsws, &
	1540	surf_usm_v(2)%ncsws, &
	1541	surf_def_v(3)%ncsws, &
	1542	surf_lsm_v(3)%ncsws, &
	1543	surf_usm_v(3)%ncsws )
	1544	ELSE
	1545	!
	1546	!-- Output of averaged data
	1547	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1548	REAL( average_count_surf, KIND=wp )
	1549	surfaces%var_av(:,n_out) = 0.0_wp
	1550
	1551	ENDIF
	1552
	1553	CASE ( 'qrsws' )
	1554	!
	1555	!-- Output of instantaneous data
	1556	IF ( av == 0 ) THEN
	1557	CALL surface_output_collect( surf_def_h(0)%qrsws, &
	1558	surf_def_h(1)%qrsws, &
	1559	surf_lsm_h%qrsws, &
	1560	surf_usm_h%qrsws, &
	1561	surf_def_v(0)%qrsws, &
	1562	surf_lsm_v(0)%qrsws, &
	1563	surf_usm_v(0)%qrsws, &
	1564	surf_def_v(1)%qrsws, &
	1565	surf_lsm_v(1)%qrsws, &
	1566	surf_usm_v(1)%qrsws, &
	1567	surf_def_v(2)%qrsws, &
	1568	surf_lsm_v(2)%qrsws, &
	1569	surf_usm_v(2)%qrsws, &
	1570	surf_def_v(3)%qrsws, &
	1571	surf_lsm_v(3)%qrsws, &
	1572	surf_usm_v(3)%qrsws )
	1573	ELSE
	1574	!
	1575	!-- Output of averaged data
	1576	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1577	REAL( average_count_surf, KIND=wp )
	1578	surfaces%var_av(:,n_out) = 0.0_wp
	1579
	1580	ENDIF
	1581
	1582	CASE ( 'nrsws' )
	1583	!
	1584	!-- Output of instantaneous data
	1585	IF ( av == 0 ) THEN
	1586	CALL surface_output_collect( surf_def_h(0)%nrsws, &
	1587	surf_def_h(1)%nrsws, &
	1588	surf_lsm_h%nrsws, &
	1589	surf_usm_h%nrsws, &
	1590	surf_def_v(0)%nrsws, &
	1591	surf_lsm_v(0)%nrsws, &
	1592	surf_usm_v(0)%nrsws, &
	1593	surf_def_v(1)%nrsws, &
	1594	surf_lsm_v(1)%nrsws, &
	1595	surf_usm_v(1)%nrsws, &
	1596	surf_def_v(2)%nrsws, &
	1597	surf_lsm_v(2)%nrsws, &
	1598	surf_usm_v(2)%nrsws, &
	1599	surf_def_v(3)%nrsws, &
	1600	surf_lsm_v(3)%nrsws, &
	1601	surf_usm_v(3)%nrsws )
	1602	ELSE
	1603	!
	1604	!-- Output of averaged data
	1605	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1606	REAL( average_count_surf, KIND=wp )
	1607	surfaces%var_av(:,n_out) = 0.0_wp
	1608
	1609	ENDIF
	1610
	1611	CASE ( 'sasws' )
	1612	!
	1613	!-- Output of instantaneous data
	1614	IF ( av == 0 ) THEN
	1615	CALL surface_output_collect( surf_def_h(0)%sasws, &
	1616	surf_def_h(1)%sasws, &
	1617	surf_lsm_h%sasws, &
	1618	surf_usm_h%sasws, &
	1619	surf_def_v(0)%sasws, &
	1620	surf_lsm_v(0)%sasws, &
	1621	surf_usm_v(0)%sasws, &
	1622	surf_def_v(1)%sasws, &
	1623	surf_lsm_v(1)%sasws, &
	1624	surf_usm_v(1)%sasws, &
	1625	surf_def_v(2)%sasws, &
	1626	surf_lsm_v(2)%sasws, &
	1627	surf_usm_v(2)%sasws, &
	1628	surf_def_v(3)%sasws, &
	1629	surf_lsm_v(3)%sasws, &
	1630	surf_usm_v(3)%sasws )
	1631	ELSE
	1632	!
	1633	!-- Output of averaged data
	1634	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1635	REAL( average_count_surf, KIND=wp )
	1636	surfaces%var_av(:,n_out) = 0.0_wp
	1637
	1638	ENDIF
	1639
	1640	CASE ( 'q_surface' )
	1641	!
	1642	!-- Output of instantaneous data
	1643	IF ( av == 0 ) THEN
	1644	CALL surface_output_collect( surf_def_h(0)%q_surface, &
	1645	surf_def_h(1)%q_surface, &
	1646	surf_lsm_h%q_surface, &
	1647	surf_usm_h%q_surface, &
	1648	surf_def_v(0)%q_surface, &
	1649	surf_lsm_v(0)%q_surface, &
	1650	surf_usm_v(0)%q_surface, &
	1651	surf_def_v(1)%q_surface, &
	1652	surf_lsm_v(1)%q_surface, &
	1653	surf_usm_v(1)%q_surface, &
	1654	surf_def_v(2)%q_surface, &
	1655	surf_lsm_v(2)%q_surface, &
	1656	surf_usm_v(2)%q_surface, &
	1657	surf_def_v(3)%q_surface, &
	1658	surf_lsm_v(3)%q_surface, &
	1659	surf_usm_v(3)%q_surface )
	1660	ELSE
	1661	!
	1662	!-- Output of averaged data
	1663	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1664	REAL( average_count_surf, KIND=wp )
	1665	surfaces%var_av(:,n_out) = 0.0_wp
	1666
	1667	ENDIF
	1668
	1669	CASE ( 'theta_surface' )
	1670	!
	1671	!-- Output of instantaneous data
	1672	IF ( av == 0 ) THEN
	1673	CALL surface_output_collect( surf_def_h(0)%pt_surface, &
	1674	surf_def_h(1)%pt_surface, &
	1675	surf_lsm_h%pt_surface, &
	1676	surf_usm_h%pt_surface, &
	1677	surf_def_v(0)%pt_surface, &
	1678	surf_lsm_v(0)%pt_surface, &
	1679	surf_usm_v(0)%pt_surface, &
	1680	surf_def_v(1)%pt_surface, &
	1681	surf_lsm_v(1)%pt_surface, &
	1682	surf_usm_v(1)%pt_surface, &
	1683	surf_def_v(2)%pt_surface, &
	1684	surf_lsm_v(2)%pt_surface, &
	1685	surf_usm_v(2)%pt_surface, &
	1686	surf_def_v(3)%pt_surface, &
	1687	surf_lsm_v(3)%pt_surface, &
	1688	surf_usm_v(3)%pt_surface )
	1689	ELSE
	1690	!
	1691	!-- Output of averaged data
	1692	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1693	REAL( average_count_surf, KIND=wp )
	1694	surfaces%var_av(:,n_out) = 0.0_wp
	1695
	1696	ENDIF
	1697
	1698	CASE ( 'thetav_surface' )
	1699	!
	1700	!-- Output of instantaneous data
	1701	IF ( av == 0 ) THEN
	1702	CALL surface_output_collect( surf_def_h(0)%vpt_surface, &
	1703	surf_def_h(1)%vpt_surface, &
	1704	surf_lsm_h%vpt_surface, &
	1705	surf_usm_h%vpt_surface, &
	1706	surf_def_v(0)%vpt_surface, &
	1707	surf_lsm_v(0)%vpt_surface, &
	1708	surf_usm_v(0)%vpt_surface, &
	1709	surf_def_v(1)%vpt_surface, &
	1710	surf_lsm_v(1)%vpt_surface, &
	1711	surf_usm_v(1)%vpt_surface, &
	1712	surf_def_v(2)%vpt_surface, &
	1713	surf_lsm_v(2)%vpt_surface, &
	1714	surf_usm_v(2)%vpt_surface, &
	1715	surf_def_v(3)%vpt_surface, &
	1716	surf_lsm_v(3)%vpt_surface, &
	1717	surf_usm_v(3)%vpt_surface)
	1718	ELSE
	1719	!
	1720	!-- Output of averaged data
	1721	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1722	REAL( average_count_surf, KIND=wp )
	1723	surfaces%var_av(:,n_out) = 0.0_wp
	1724
	1725	ENDIF
	1726
	1727	CASE ( 'rad_net' )
	1728	!
	1729	!-- Output of instantaneous data
	1730	IF ( av == 0 ) THEN
	1731	CALL surface_output_collect( surf_def_h(0)%rad_net, &
	1732	surf_def_h(1)%rad_net, &
	1733	surf_lsm_h%rad_net, &
	1734	surf_usm_h%rad_net, &
	1735	surf_def_v(0)%rad_net, &
	1736	surf_lsm_v(0)%rad_net, &
	1737	surf_usm_v(0)%rad_net, &
	1738	surf_def_v(1)%rad_net, &
	1739	surf_lsm_v(1)%rad_net, &
	1740	surf_usm_v(1)%rad_net, &
	1741	surf_def_v(2)%rad_net, &
	1742	surf_lsm_v(2)%rad_net, &
	1743	surf_usm_v(2)%rad_net, &
	1744	surf_def_v(3)%rad_net, &
	1745	surf_lsm_v(3)%rad_net, &
	1746	surf_usm_v(3)%rad_net )
	1747	ELSE
	1748	!
	1749	!-- Output of averaged data
	1750	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1751	REAL( average_count_surf, KIND=wp )
	1752	surfaces%var_av(:,n_out) = 0.0_wp
	1753
	1754	ENDIF
	1755
	1756	CASE ( 'rad_lw_in' )
	1757	!
	1758	!-- Output of instantaneous data
	1759	IF ( av == 0 ) THEN
	1760	CALL surface_output_collect( surf_def_h(0)%rad_lw_in, &
	1761	surf_def_h(1)%rad_lw_in, &
	1762	surf_lsm_h%rad_lw_in, &
	1763	surf_usm_h%rad_lw_in, &
	1764	surf_def_v(0)%rad_lw_in, &
	1765	surf_lsm_v(0)%rad_lw_in, &
	1766	surf_usm_v(0)%rad_lw_in, &
	1767	surf_def_v(1)%rad_lw_in, &
	1768	surf_lsm_v(1)%rad_lw_in, &
	1769	surf_usm_v(1)%rad_lw_in, &
	1770	surf_def_v(2)%rad_lw_in, &
	1771	surf_lsm_v(2)%rad_lw_in, &
	1772	surf_usm_v(2)%rad_lw_in, &
	1773	surf_def_v(3)%rad_lw_in, &
	1774	surf_lsm_v(3)%rad_lw_in, &
	1775	surf_usm_v(3)%rad_lw_in )
	1776	ELSE
	1777	!
	1778	!-- Output of averaged data
	1779	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1780	REAL( average_count_surf, KIND=wp )
	1781	surfaces%var_av(:,n_out) = 0.0_wp
	1782
	1783	ENDIF
	1784
	1785	CASE ( 'rad_lw_out' )
	1786	!
	1787	!-- Output of instantaneous data
	1788	IF ( av == 0 ) THEN
	1789	CALL surface_output_collect( surf_def_h(0)%rad_lw_out, &
	1790	surf_def_h(1)%rad_lw_out, &
	1791	surf_lsm_h%rad_lw_out, &
	1792	surf_usm_h%rad_lw_out, &
	1793	surf_def_v(0)%rad_lw_out, &
	1794	surf_lsm_v(0)%rad_lw_out, &
	1795	surf_usm_v(0)%rad_lw_out, &
	1796	surf_def_v(1)%rad_lw_out, &
	1797	surf_lsm_v(1)%rad_lw_out, &
	1798	surf_usm_v(1)%rad_lw_out, &
	1799	surf_def_v(2)%rad_lw_out, &
	1800	surf_lsm_v(2)%rad_lw_out, &
	1801	surf_usm_v(2)%rad_lw_out, &
	1802	surf_def_v(3)%rad_lw_out, &
	1803	surf_lsm_v(3)%rad_lw_out, &
	1804	surf_usm_v(3)%rad_lw_out )
	1805	ELSE
	1806	!
	1807	!-- Output of averaged data
	1808	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1809	REAL( average_count_surf, KIND=wp )
	1810	surfaces%var_av(:,n_out) = 0.0_wp
	1811
	1812	ENDIF
	1813
	1814	CASE ( 'rad_sw_in' )
	1815	!
	1816	!-- Output of instantaneous data
	1817	IF ( av == 0 ) THEN
	1818	CALL surface_output_collect( surf_def_h(0)%rad_sw_in, &
	1819	surf_def_h(1)%rad_sw_in, &
	1820	surf_lsm_h%rad_sw_in, &
	1821	surf_usm_h%rad_sw_in, &
	1822	surf_def_v(0)%rad_sw_in, &
	1823	surf_lsm_v(0)%rad_sw_in, &
	1824	surf_usm_v(0)%rad_sw_in, &
	1825	surf_def_v(1)%rad_sw_in, &
	1826	surf_lsm_v(1)%rad_sw_in, &
	1827	surf_usm_v(1)%rad_sw_in, &
	1828	surf_def_v(2)%rad_sw_in, &
	1829	surf_lsm_v(2)%rad_sw_in, &
	1830	surf_usm_v(2)%rad_sw_in, &
	1831	surf_def_v(3)%rad_sw_in, &
	1832	surf_lsm_v(3)%rad_sw_in, &
	1833	surf_usm_v(3)%rad_sw_in )
	1834	ELSE
	1835	!
	1836	!-- Output of averaged data
	1837	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1838	REAL( average_count_surf, KIND=wp )
	1839	surfaces%var_av(:,n_out) = 0.0_wp
	1840
	1841	ENDIF
	1842
	1843	CASE ( 'rad_sw_out' )
	1844	!
	1845	!-- Output of instantaneous data
	1846	IF ( av == 0 ) THEN
	1847	CALL surface_output_collect( surf_def_h(0)%rad_sw_out, &
	1848	surf_def_h(1)%rad_sw_out, &
	1849	surf_lsm_h%rad_sw_out, &
	1850	surf_usm_h%rad_sw_out, &
	1851	surf_def_v(0)%rad_sw_out, &
	1852	surf_lsm_v(0)%rad_sw_out, &
	1853	surf_usm_v(0)%rad_sw_out, &
	1854	surf_def_v(1)%rad_sw_out, &
	1855	surf_lsm_v(1)%rad_sw_out, &
	1856	surf_usm_v(1)%rad_sw_out, &
	1857	surf_def_v(2)%rad_sw_out, &
	1858	surf_lsm_v(2)%rad_sw_out, &
	1859	surf_usm_v(2)%rad_sw_out, &
	1860	surf_def_v(3)%rad_sw_out, &
	1861	surf_lsm_v(3)%rad_sw_out, &
	1862	surf_usm_v(3)%rad_sw_out )
	1863	ELSE
	1864	!
	1865	!-- Output of averaged data
	1866	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1867	REAL( average_count_surf, KIND=wp )
	1868	surfaces%var_av(:,n_out) = 0.0_wp
	1869
	1870	ENDIF
	1871
	1872	CASE ( 'ghf' )
	1873	!
	1874	!-- Output of instantaneous data
	1875	IF ( av == 0 ) THEN
	1876	CALL surface_output_collect( surf_def_h(0)%ghf, &
	1877	surf_def_h(1)%ghf, &
	1878	surf_lsm_h%ghf, &
	1879	surf_usm_h%ghf, &
	1880	surf_def_v(0)%ghf, &
	1881	surf_lsm_v(0)%ghf, &
	1882	surf_usm_v(0)%ghf, &
	1883	surf_def_v(1)%ghf, &
	1884	surf_lsm_v(1)%ghf, &
	1885	surf_usm_v(1)%ghf, &
	1886	surf_def_v(2)%ghf, &
	1887	surf_lsm_v(2)%ghf, &
	1888	surf_usm_v(2)%ghf, &
	1889	surf_def_v(3)%ghf, &
	1890	surf_lsm_v(3)%ghf, &
	1891	surf_usm_v(3)%ghf )
	1892	ELSE
	1893	!
	1894	!-- Output of averaged data
	1895	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1896	REAL( average_count_surf, KIND=wp )
	1897	surfaces%var_av(:,n_out) = 0.0_wp
	1898
	1899	ENDIF
	1900
	1901	CASE ( 'r_a' )
	1902	!
	1903	!-- Output of instantaneous data
	1904	IF ( av == 0 ) THEN
	1905	CALL surface_output_collect( surf_def_h(0)%r_a, &
	1906	surf_def_h(1)%r_a, &
	1907	surf_lsm_h%r_a, &
	1908	surf_usm_h%r_a, &
	1909	surf_def_v(0)%r_a, &
	1910	surf_lsm_v(0)%r_a, &
	1911	surf_usm_v(0)%r_a, &
	1912	surf_def_v(1)%r_a, &
	1913	surf_lsm_v(1)%r_a, &
	1914	surf_usm_v(1)%r_a, &
	1915	surf_def_v(2)%r_a, &
	1916	surf_lsm_v(2)%r_a, &
	1917	surf_usm_v(2)%r_a, &
	1918	surf_def_v(3)%r_a, &
	1919	surf_lsm_v(3)%r_a, &
	1920	surf_usm_v(3)%r_a )
	1921	ELSE
	1922	!
	1923	!-- Output of averaged data
	1924	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1925	REAL( average_count_surf, KIND=wp )
	1926	surfaces%var_av(:,n_out) = 0.0_wp
	1927
	1928	ENDIF
	1929
	1930	CASE ( 'r_soil' )
	1931	!
	1932	!-- Output of instantaneous data
	1933	IF ( av == 0 ) THEN
	1934	CALL surface_output_collect( surf_def_h(0)%r_soil, &
	1935	surf_def_h(1)%r_soil, &
	1936	surf_lsm_h%r_soil, &
	1937	surf_usm_h%r_soil, &
	1938	surf_def_v(0)%r_soil, &
	1939	surf_lsm_v(0)%r_soil, &
	1940	surf_usm_v(0)%r_soil, &
	1941	surf_def_v(1)%r_soil, &
	1942	surf_lsm_v(1)%r_soil, &
	1943	surf_usm_v(1)%r_soil, &
	1944	surf_def_v(2)%r_soil, &
	1945	surf_lsm_v(2)%r_soil, &
	1946	surf_usm_v(2)%r_soil, &
	1947	surf_def_v(3)%r_soil, &
	1948	surf_lsm_v(3)%r_soil, &
	1949	surf_usm_v(3)%r_soil )
	1950	ELSE
	1951	!
	1952	!-- Output of averaged data
	1953	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1954	REAL( average_count_surf, KIND=wp )
	1955	surfaces%var_av(:,n_out) = 0.0_wp
	1956
	1957	ENDIF
	1958
	1959	CASE ( 'r_canopy' )
	1960	!
	1961	!-- Output of instantaneous data
	1962	IF ( av == 0 ) THEN
	1963	CALL surface_output_collect( surf_def_h(0)%r_canopy, &
	1964	surf_def_h(1)%r_canopy, &
	1965	surf_lsm_h%r_canopy, &
	1966	surf_usm_h%r_canopy, &
	1967	surf_def_v(0)%r_canopy, &
	1968	surf_lsm_v(0)%r_canopy, &
	1969	surf_usm_v(0)%r_canopy, &
	1970	surf_def_v(1)%r_canopy, &
	1971	surf_lsm_v(1)%r_canopy, &
	1972	surf_usm_v(1)%r_canopy, &
	1973	surf_def_v(2)%r_canopy, &
	1974	surf_lsm_v(2)%r_canopy, &
	1975	surf_usm_v(2)%r_canopy, &
	1976	surf_def_v(3)%r_canopy, &
	1977	surf_lsm_v(3)%r_canopy, &
	1978	surf_usm_v(3)%r_canopy )
	1979	ELSE
	1980	!
	1981	!-- Output of averaged data
	1982	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	1983	REAL( average_count_surf, KIND=wp )
	1984	surfaces%var_av(:,n_out) = 0.0_wp
	1985
	1986	ENDIF
	1987
	1988	CASE ( 'r_s' )
	1989	!
	1990	!-- Output of instantaneous data
	1991	IF ( av == 0 ) THEN
	1992	CALL surface_output_collect( surf_def_h(0)%r_s, &
	1993	surf_def_h(1)%r_s, &
	1994	surf_lsm_h%r_s, &
	1995	surf_usm_h%r_s, &
	1996	surf_def_v(0)%r_s, &
	1997	surf_lsm_v(0)%r_s, &
	1998	surf_usm_v(0)%r_s, &
	1999	surf_def_v(1)%r_s, &
	2000	surf_lsm_v(1)%r_s, &
	2001	surf_usm_v(1)%r_s, &
	2002	surf_def_v(2)%r_s, &
	2003	surf_lsm_v(2)%r_s, &
	2004	surf_usm_v(2)%r_s, &
	2005	surf_def_v(3)%r_s, &
	2006	surf_lsm_v(3)%r_s, &
	2007	surf_usm_v(3)%r_s )
	2008	ELSE
	2009	!
	2010	!-- Output of averaged data
	2011	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2012	REAL( average_count_surf, KIND=wp )
	2013	surfaces%var_av(:,n_out) = 0.0_wp
	2014
	2015	ENDIF
	2016
[3572]	2017	CASE ( 'rad_sw_dir' )
[3421]	2018	!
[3572]	2019	!-- Output of instantaneous data
	2020	IF ( av == 0 ) THEN
	2021	CALL surface_output_collect( surf_def_h(0)%rad_sw_dir, &
	2022	surf_def_h(1)%rad_sw_dir, &
	2023	surf_lsm_h%rad_sw_dir, &
	2024	surf_usm_h%rad_sw_dir, &
	2025	surf_def_v(0)%rad_sw_dir, &
	2026	surf_lsm_v(0)%rad_sw_dir, &
	2027	surf_usm_v(0)%rad_sw_dir, &
	2028	surf_def_v(1)%rad_sw_dir, &
	2029	surf_lsm_v(1)%rad_sw_dir, &
	2030	surf_usm_v(1)%rad_sw_dir, &
	2031	surf_def_v(2)%rad_sw_dir, &
	2032	surf_lsm_v(2)%rad_sw_dir, &
	2033	surf_usm_v(2)%rad_sw_dir, &
	2034	surf_def_v(3)%rad_sw_dir, &
	2035	surf_lsm_v(3)%rad_sw_dir, &
	2036	surf_usm_v(3)%rad_sw_dir )
	2037	ELSE
	2038	!
	2039	!-- Output of averaged data
	2040	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2041	REAL( average_count_surf, KIND=wp )
	2042	surfaces%var_av(:,n_out) = 0.0_wp
	2043
	2044	ENDIF
	2045
	2046	CASE ( 'rad_sw_dif' )
	2047	!
	2048	!-- Output of instantaneous data
	2049	IF ( av == 0 ) THEN
	2050	CALL surface_output_collect( surf_def_h(0)%rad_sw_dif, &
	2051	surf_def_h(1)%rad_sw_dif, &
	2052	surf_lsm_h%rad_sw_dif, &
	2053	surf_usm_h%rad_sw_dif, &
	2054	surf_def_v(0)%rad_sw_dif, &
	2055	surf_lsm_v(0)%rad_sw_dif, &
	2056	surf_usm_v(0)%rad_sw_dif, &
	2057	surf_def_v(1)%rad_sw_dif, &
	2058	surf_lsm_v(1)%rad_sw_dif, &
	2059	surf_usm_v(1)%rad_sw_dif, &
	2060	surf_def_v(2)%rad_sw_dif, &
	2061	surf_lsm_v(2)%rad_sw_dif, &
	2062	surf_usm_v(2)%rad_sw_dif, &
	2063	surf_def_v(3)%rad_sw_dif, &
	2064	surf_lsm_v(3)%rad_sw_dif, &
	2065	surf_usm_v(3)%rad_sw_dif )
	2066	ELSE
	2067	!
	2068	!-- Output of averaged data
	2069	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2070	REAL( average_count_surf, KIND=wp )
	2071	surfaces%var_av(:,n_out) = 0.0_wp
	2072
	2073	ENDIF
	2074
	2075	CASE ( 'rad_sw_ref' )
	2076	!
	2077	!-- Output of instantaneous data
	2078	IF ( av == 0 ) THEN
	2079	CALL surface_output_collect( surf_def_h(0)%rad_sw_ref, &
	2080	surf_def_h(1)%rad_sw_ref, &
	2081	surf_lsm_h%rad_sw_ref, &
	2082	surf_usm_h%rad_sw_ref, &
	2083	surf_def_v(0)%rad_sw_ref, &
	2084	surf_lsm_v(0)%rad_sw_ref, &
	2085	surf_usm_v(0)%rad_sw_ref, &
	2086	surf_def_v(1)%rad_sw_ref, &
	2087	surf_lsm_v(1)%rad_sw_ref, &
	2088	surf_usm_v(1)%rad_sw_ref, &
	2089	surf_def_v(2)%rad_sw_ref, &
	2090	surf_lsm_v(2)%rad_sw_ref, &
	2091	surf_usm_v(2)%rad_sw_ref, &
	2092	surf_def_v(3)%rad_sw_ref, &
	2093	surf_lsm_v(3)%rad_sw_ref, &
	2094	surf_usm_v(3)%rad_sw_ref )
	2095	ELSE
	2096	!
	2097	!-- Output of averaged data
	2098	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2099	REAL( average_count_surf, KIND=wp )
	2100	surfaces%var_av(:,n_out) = 0.0_wp
	2101
	2102	ENDIF
	2103
	2104	CASE ( 'rad_sw_res' )
	2105	!
	2106	!-- Output of instantaneous data
	2107	IF ( av == 0 ) THEN
	2108	CALL surface_output_collect( surf_def_h(0)%rad_sw_res, &
	2109	surf_def_h(1)%rad_sw_res, &
	2110	surf_lsm_h%rad_sw_res, &
	2111	surf_usm_h%rad_sw_res, &
	2112	surf_def_v(0)%rad_sw_res, &
	2113	surf_lsm_v(0)%rad_sw_res, &
	2114	surf_usm_v(0)%rad_sw_res, &
	2115	surf_def_v(1)%rad_sw_res, &
	2116	surf_lsm_v(1)%rad_sw_res, &
	2117	surf_usm_v(1)%rad_sw_res, &
	2118	surf_def_v(2)%rad_sw_res, &
	2119	surf_lsm_v(2)%rad_sw_res, &
	2120	surf_usm_v(2)%rad_sw_res, &
	2121	surf_def_v(3)%rad_sw_res, &
	2122	surf_lsm_v(3)%rad_sw_res, &
	2123	surf_usm_v(3)%rad_sw_res )
	2124	ELSE
	2125	!
	2126	!-- Output of averaged data
	2127	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2128	REAL( average_count_surf, KIND=wp )
	2129	surfaces%var_av(:,n_out) = 0.0_wp
	2130
	2131	ENDIF
	2132
	2133	CASE ( 'rad_lw_dif' )
	2134	!
	2135	!-- Output of instantaneous data
	2136	IF ( av == 0 ) THEN
	2137	CALL surface_output_collect( surf_def_h(0)%rad_lw_dif, &
	2138	surf_def_h(1)%rad_lw_dif, &
	2139	surf_lsm_h%rad_lw_dif, &
	2140	surf_usm_h%rad_lw_dif, &
	2141	surf_def_v(0)%rad_lw_dif, &
	2142	surf_lsm_v(0)%rad_lw_dif, &
	2143	surf_usm_v(0)%rad_lw_dif, &
	2144	surf_def_v(1)%rad_lw_dif, &
	2145	surf_lsm_v(1)%rad_lw_dif, &
	2146	surf_usm_v(1)%rad_lw_dif, &
	2147	surf_def_v(2)%rad_lw_dif, &
	2148	surf_lsm_v(2)%rad_lw_dif, &
	2149	surf_usm_v(2)%rad_lw_dif, &
	2150	surf_def_v(3)%rad_lw_dif, &
	2151	surf_lsm_v(3)%rad_lw_dif, &
	2152	surf_usm_v(3)%rad_lw_dif )
	2153	ELSE
	2154	!
	2155	!-- Output of averaged data
	2156	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2157	REAL( average_count_surf, KIND=wp )
	2158	surfaces%var_av(:,n_out) = 0.0_wp
	2159
	2160	ENDIF
	2161
	2162	CASE ( 'rad_lw_ref' )
	2163	!
	2164	!-- Output of instantaneous data
	2165	IF ( av == 0 ) THEN
	2166	CALL surface_output_collect( surf_def_h(0)%rad_lw_ref, &
	2167	surf_def_h(1)%rad_lw_ref, &
	2168	surf_lsm_h%rad_lw_ref, &
	2169	surf_usm_h%rad_lw_ref, &
	2170	surf_def_v(0)%rad_lw_ref, &
	2171	surf_lsm_v(0)%rad_lw_ref, &
	2172	surf_usm_v(0)%rad_lw_ref, &
	2173	surf_def_v(1)%rad_lw_ref, &
	2174	surf_lsm_v(1)%rad_lw_ref, &
	2175	surf_usm_v(1)%rad_lw_ref, &
	2176	surf_def_v(2)%rad_lw_ref, &
	2177	surf_lsm_v(2)%rad_lw_ref, &
	2178	surf_usm_v(2)%rad_lw_ref, &
	2179	surf_def_v(3)%rad_lw_ref, &
	2180	surf_lsm_v(3)%rad_lw_ref, &
	2181	surf_usm_v(3)%rad_lw_ref )
	2182	ELSE
	2183	!
	2184	!-- Output of averaged data
	2185	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2186	REAL( average_count_surf, KIND=wp )
	2187	surfaces%var_av(:,n_out) = 0.0_wp
	2188
	2189	ENDIF
	2190
	2191	CASE ( 'rad_lw_res' )
	2192	!
	2193	!-- Output of instantaneous data
	2194	IF ( av == 0 ) THEN
	2195	CALL surface_output_collect( surf_def_h(0)%rad_lw_res, &
	2196	surf_def_h(1)%rad_lw_res, &
	2197	surf_lsm_h%rad_lw_res, &
	2198	surf_usm_h%rad_lw_res, &
	2199	surf_def_v(0)%rad_lw_res, &
	2200	surf_lsm_v(0)%rad_lw_res, &
	2201	surf_usm_v(0)%rad_lw_res, &
	2202	surf_def_v(1)%rad_lw_res, &
	2203	surf_lsm_v(1)%rad_lw_res, &
	2204	surf_usm_v(1)%rad_lw_res, &
	2205	surf_def_v(2)%rad_lw_res, &
	2206	surf_lsm_v(2)%rad_lw_res, &
	2207	surf_usm_v(2)%rad_lw_res, &
	2208	surf_def_v(3)%rad_lw_res, &
	2209	surf_lsm_v(3)%rad_lw_res, &
	2210	surf_usm_v(3)%rad_lw_res )
	2211	ELSE
	2212	!
	2213	!-- Output of averaged data
	2214	surfaces%var_out(:) = surfaces%var_av(:,n_out) / &
	2215	REAL( average_count_surf, KIND=wp )
	2216	surfaces%var_av(:,n_out) = 0.0_wp
	2217
	2218	ENDIF
	2219
	2220	!
[3421]	2221	!-- Add further variables:
	2222	!-- 'css', 'cssws', 'qsws_liq', 'qsws_soil', 'qsws_veg'
	2223
	2224	END SELECT
	2225	!
	2226	!-- Write to binary file:
	2227	!-- - surfaces%points ( 3, 1-npoints )
	2228	!-- - surfaces%polygons ( 5, 1-ns )
	2229	!-- - surfaces%var_out ( 1-ns, time )
	2230	!-- - Dimension: 1-nsurfaces, 1-npoints - can be ordered consecutively
	2231	!-- - Distinguish between average and non-average data
	2232	DO i = 0, io_blocks-1
	2233	IF ( i == io_group ) THEN
[3494]	2234	WRITE ( 25+av ) LEN_TRIM( 'time' )
	2235	WRITE ( 25+av ) 'time'
	2236	WRITE ( 25+av ) simulated_time
	2237	WRITE ( 25+av ) LEN_TRIM( trimvar )
	2238	WRITE ( 25+av ) TRIM( trimvar )
[3421]	2239	WRITE ( 25+av ) surfaces%var_out
	2240	ENDIF
[3483]	2241	#if defined( __parallel )
[3421]	2242	CALL MPI_BARRIER( comm2d, ierr )
[3483]	2243	#endif
[3421]	2244	ENDDO
	2245
	2246	ENDDO
	2247
	2248	!
[3494]	2249	!-- If averaged output was written to NetCDF file, set the counter to zero
	2250	IF ( av == 1 ) average_count_surf = 0
[3421]	2251
	2252	END SUBROUTINE surface_output
	2253
	2254	!------------------------------------------------------------------------------!
	2255	! Description:
	2256	! ------------
	2257	!> Routine for controlling the data averaging.
	2258	!------------------------------------------------------------------------------!
	2259	SUBROUTINE surface_output_averaging
	2260
	2261	IMPLICIT NONE
	2262
	2263	CHARACTER(LEN=100) :: trimvar !< dummy variable for current output variable
	2264
	2265	INTEGER(iwp) :: n_out !< counter variables for surface output
	2266
	2267	n_out = 0
	2268	DO WHILE ( dosurf(1,n_out+1)(1:1) /= ' ' )
	2269
	2270	n_out = n_out + 1
	2271	trimvar = TRIM( dosurf(1,n_out) )
	2272
	2273	SELECT CASE ( trimvar )
	2274
	2275	CASE ( 'us' )
	2276	CALL surface_output_sum_up( surf_def_h(0)%us, &
	2277	surf_def_h(1)%us, &
	2278	surf_lsm_h%us, &
	2279	surf_usm_h%us, &
	2280	surf_def_v(0)%us, &
	2281	surf_lsm_v(0)%us, &
	2282	surf_usm_v(0)%us, &
	2283	surf_def_v(1)%us, &
	2284	surf_lsm_v(1)%us, &
	2285	surf_usm_v(1)%us, &
	2286	surf_def_v(2)%us, &
	2287	surf_lsm_v(2)%us, &
	2288	surf_usm_v(2)%us, &
	2289	surf_def_v(3)%us, &
	2290	surf_lsm_v(3)%us, &
	2291	surf_usm_v(3)%us, n_out )
	2292
	2293	CASE ( 'ts' )
	2294	CALL surface_output_sum_up( surf_def_h(0)%ts, &
	2295	surf_def_h(1)%ts, &
	2296	surf_lsm_h%ts, &
	2297	surf_usm_h%ts, &
	2298	surf_def_v(0)%ts, &
	2299	surf_lsm_v(0)%ts, &
	2300	surf_usm_v(0)%ts, &
	2301	surf_def_v(1)%ts, &
	2302	surf_lsm_v(1)%ts, &
	2303	surf_usm_v(1)%ts, &
	2304	surf_def_v(2)%ts, &
	2305	surf_lsm_v(2)%ts, &
	2306	surf_usm_v(2)%ts, &
	2307	surf_def_v(3)%ts, &
	2308	surf_lsm_v(3)%ts, &
	2309	surf_usm_v(3)%ts, n_out )
	2310
	2311	CASE ( 'qs' )
	2312	CALL surface_output_sum_up( surf_def_h(0)%qs, &
	2313	surf_def_h(1)%qs, &
	2314	surf_lsm_h%qs, &
	2315	surf_usm_h%qs, &
	2316	surf_def_v(0)%qs, &
	2317	surf_lsm_v(0)%qs, &
	2318	surf_usm_v(0)%qs, &
	2319	surf_def_v(1)%qs, &
	2320	surf_lsm_v(1)%qs, &
	2321	surf_usm_v(1)%qs, &
	2322	surf_def_v(2)%qs, &
	2323	surf_lsm_v(2)%qs, &
	2324	surf_usm_v(2)%qs, &
	2325	surf_def_v(3)%qs, &
	2326	surf_lsm_v(3)%qs, &
	2327	surf_usm_v(3)%qs, n_out )
	2328
	2329	CASE ( 'ss' )
	2330	CALL surface_output_sum_up( surf_def_h(0)%ss, &
	2331	surf_def_h(1)%ss, &
	2332	surf_lsm_h%ss, &
	2333	surf_usm_h%ss, &
	2334	surf_def_v(0)%ss, &
	2335	surf_lsm_v(0)%ss, &
	2336	surf_usm_v(0)%ss, &
	2337	surf_def_v(1)%ss, &
	2338	surf_lsm_v(1)%ss, &
	2339	surf_usm_v(1)%ss, &
	2340	surf_def_v(2)%ss, &
	2341	surf_lsm_v(2)%ss, &
	2342	surf_usm_v(2)%ss, &
	2343	surf_def_v(3)%ss, &
	2344	surf_lsm_v(3)%ss, &
	2345	surf_usm_v(3)%ss, n_out )
	2346
	2347	CASE ( 'qcs' )
	2348	CALL surface_output_sum_up( surf_def_h(0)%qcs, &
	2349	surf_def_h(1)%qcs, &
	2350	surf_lsm_h%qcs, &
	2351	surf_usm_h%qcs, &
	2352	surf_def_v(0)%qcs, &
	2353	surf_lsm_v(0)%qcs, &
	2354	surf_usm_v(0)%qcs, &
	2355	surf_def_v(1)%qcs, &
	2356	surf_lsm_v(1)%qcs, &
	2357	surf_usm_v(1)%qcs, &
	2358	surf_def_v(2)%qcs, &
	2359	surf_lsm_v(2)%qcs, &
	2360	surf_usm_v(2)%qcs, &
	2361	surf_def_v(3)%qcs, &
	2362	surf_lsm_v(3)%qcs, &
	2363	surf_usm_v(3)%qcs, n_out )
	2364
	2365	CASE ( 'ncs' )
	2366	CALL surface_output_sum_up( surf_def_h(0)%ncs, &
	2367	surf_def_h(1)%ncs, &
	2368	surf_lsm_h%ncs, &
	2369	surf_usm_h%ncs, &
	2370	surf_def_v(0)%ncs, &
	2371	surf_lsm_v(0)%ncs, &
	2372	surf_usm_v(0)%ncs, &
	2373	surf_def_v(1)%ncs, &
	2374	surf_lsm_v(1)%ncs, &
	2375	surf_usm_v(1)%ncs, &
	2376	surf_def_v(2)%ncs, &
	2377	surf_lsm_v(2)%ncs, &
	2378	surf_usm_v(2)%ncs, &
	2379	surf_def_v(3)%ncs, &
	2380	surf_lsm_v(3)%ncs, &
	2381	surf_usm_v(3)%ncs, n_out )
	2382
	2383	CASE ( 'qrs' )
	2384	CALL surface_output_sum_up( surf_def_h(0)%qrs, &
	2385	surf_def_h(1)%qrs, &
	2386	surf_lsm_h%qrs, &
	2387	surf_usm_h%qrs, &
	2388	surf_def_v(0)%qrs, &
	2389	surf_lsm_v(0)%qrs, &
	2390	surf_usm_v(0)%qrs, &
	2391	surf_def_v(1)%qrs, &
	2392	surf_lsm_v(1)%qrs, &
	2393	surf_usm_v(1)%qrs, &
	2394	surf_def_v(2)%qrs, &
	2395	surf_lsm_v(2)%qrs, &
	2396	surf_usm_v(2)%qrs, &
	2397	surf_def_v(3)%qrs, &
	2398	surf_lsm_v(3)%qrs, &
	2399	surf_usm_v(3)%qrs, n_out )
	2400
	2401	CASE ( 'nrs' )
	2402	CALL surface_output_sum_up( surf_def_h(0)%nrs, &
	2403	surf_def_h(1)%nrs, &
	2404	surf_lsm_h%nrs, &
	2405	surf_usm_h%nrs, &
	2406	surf_def_v(0)%nrs, &
	2407	surf_lsm_v(0)%nrs, &
	2408	surf_usm_v(0)%nrs, &
	2409	surf_def_v(1)%nrs, &
	2410	surf_lsm_v(1)%nrs, &
	2411	surf_usm_v(1)%nrs, &
	2412	surf_def_v(2)%nrs, &
	2413	surf_lsm_v(2)%nrs, &
	2414	surf_usm_v(2)%nrs, &
	2415	surf_def_v(3)%nrs, &
	2416	surf_lsm_v(3)%nrs, &
	2417	surf_usm_v(3)%nrs, n_out )
	2418
	2419	CASE ( 'ol' )
	2420	CALL surface_output_sum_up( surf_def_h(0)%ol, &
	2421	surf_def_h(1)%ol, &
	2422	surf_lsm_h%ol, &
	2423	surf_usm_h%ol, &
	2424	surf_def_v(0)%ol, &
	2425	surf_lsm_v(0)%ol, &
	2426	surf_usm_v(0)%ol, &
	2427	surf_def_v(1)%ol, &
	2428	surf_lsm_v(1)%ol, &
	2429	surf_usm_v(1)%ol, &
	2430	surf_def_v(2)%ol, &
	2431	surf_lsm_v(2)%ol, &
	2432	surf_usm_v(2)%ol, &
	2433	surf_def_v(3)%ol, &
	2434	surf_lsm_v(3)%ol, &
	2435	surf_usm_v(3)%ol, n_out )
	2436
	2437	CASE ( 'z0' )
	2438	CALL surface_output_sum_up( surf_def_h(0)%z0, &
	2439	surf_def_h(1)%z0, &
	2440	surf_lsm_h%z0, &
	2441	surf_usm_h%z0, &
	2442	surf_def_v(0)%z0, &
	2443	surf_lsm_v(0)%z0, &
	2444	surf_usm_v(0)%z0, &
	2445	surf_def_v(1)%z0, &
	2446	surf_lsm_v(1)%z0, &
	2447	surf_usm_v(1)%z0, &
	2448	surf_def_v(2)%z0, &
	2449	surf_lsm_v(2)%z0, &
	2450	surf_usm_v(2)%z0, &
	2451	surf_def_v(3)%z0, &
	2452	surf_lsm_v(3)%z0, &
	2453	surf_usm_v(3)%z0, n_out )
	2454
	2455	CASE ( 'z0h' )
	2456	CALL surface_output_sum_up( surf_def_h(0)%z0h, &
	2457	surf_def_h(1)%z0h, &
	2458	surf_lsm_h%z0h, &
	2459	surf_usm_h%z0h, &
	2460	surf_def_v(0)%z0h, &
	2461	surf_lsm_v(0)%z0h, &
	2462	surf_usm_v(0)%z0h, &
	2463	surf_def_v(1)%z0h, &
	2464	surf_lsm_v(1)%z0h, &
	2465	surf_usm_v(1)%z0h, &
	2466	surf_def_v(2)%z0h, &
	2467	surf_lsm_v(2)%z0h, &
	2468	surf_usm_v(2)%z0h, &
	2469	surf_def_v(3)%z0h, &
	2470	surf_lsm_v(3)%z0h, &
	2471	surf_usm_v(3)%z0h, n_out )
	2472
	2473	CASE ( 'z0q' )
	2474	CALL surface_output_sum_up( surf_def_h(0)%z0q, &
	2475	surf_def_h(1)%z0q, &
	2476	surf_lsm_h%z0q, &
	2477	surf_usm_h%z0q, &
	2478	surf_def_v(0)%z0q, &
	2479	surf_lsm_v(0)%z0q, &
	2480	surf_usm_v(0)%z0q, &
	2481	surf_def_v(1)%z0q, &
	2482	surf_lsm_v(1)%z0q, &
	2483	surf_usm_v(1)%z0q, &
	2484	surf_def_v(2)%z0q, &
	2485	surf_lsm_v(2)%z0q, &
	2486	surf_usm_v(2)%z0q, &
	2487	surf_def_v(3)%z0q, &
	2488	surf_lsm_v(3)%z0q, &
	2489	surf_usm_v(3)%z0q, n_out )
	2490
	2491	CASE ( 'theta1' )
	2492	CALL surface_output_sum_up( surf_def_h(0)%pt1, &
	2493	surf_def_h(1)%pt1, &
	2494	surf_lsm_h%pt1, &
	2495	surf_usm_h%pt1, &
	2496	surf_def_v(0)%pt1, &
	2497	surf_lsm_v(0)%pt1, &
	2498	surf_usm_v(0)%pt1, &
	2499	surf_def_v(1)%pt1, &
	2500	surf_lsm_v(1)%pt1, &
	2501	surf_usm_v(1)%pt1, &
	2502	surf_def_v(2)%pt1, &
	2503	surf_lsm_v(2)%pt1, &
	2504	surf_usm_v(2)%pt1, &
	2505	surf_def_v(3)%pt1, &
	2506	surf_lsm_v(3)%pt1, &
	2507	surf_usm_v(3)%pt1, n_out )
	2508
	2509	CASE ( 'qv1' )
	2510	CALL surface_output_sum_up( surf_def_h(0)%qv1, &
	2511	surf_def_h(1)%qv1, &
	2512	surf_lsm_h%qv1, &
	2513	surf_usm_h%qv1, &
	2514	surf_def_v(0)%qv1, &
	2515	surf_lsm_v(0)%qv1, &
	2516	surf_usm_v(0)%qv1, &
	2517	surf_def_v(1)%qv1, &
	2518	surf_lsm_v(1)%qv1, &
	2519	surf_usm_v(1)%qv1, &
	2520	surf_def_v(2)%qv1, &
	2521	surf_lsm_v(2)%qv1, &
	2522	surf_usm_v(2)%qv1, &
	2523	surf_def_v(3)%qv1, &
	2524	surf_lsm_v(3)%qv1, &
	2525	surf_usm_v(3)%qv1, n_out )
	2526
	2527	CASE ( 'thetav1' )
	2528	CALL surface_output_sum_up( surf_def_h(0)%vpt1, &
	2529	surf_def_h(1)%vpt1, &
	2530	surf_lsm_h%vpt1, &
	2531	surf_usm_h%vpt1, &
	2532	surf_def_v(0)%vpt1, &
	2533	surf_lsm_v(0)%vpt1, &
	2534	surf_usm_v(0)%vpt1, &
	2535	surf_def_v(1)%vpt1, &
	2536	surf_lsm_v(1)%vpt1, &
	2537	surf_usm_v(1)%vpt1, &
	2538	surf_def_v(2)%vpt1, &
	2539	surf_lsm_v(2)%vpt1, &
	2540	surf_usm_v(2)%vpt1, &
	2541	surf_def_v(3)%vpt1, &
	2542	surf_lsm_v(3)%vpt1, &
	2543	surf_usm_v(3)%vpt1, n_out )
	2544
	2545	CASE ( 'usws' )
	2546	CALL surface_output_sum_up( surf_def_h(0)%usws, &
	2547	surf_def_h(1)%usws, &
	2548	surf_lsm_h%usws, &
	2549	surf_usm_h%usws, &
	2550	surf_def_v(0)%usws, &
	2551	surf_lsm_v(0)%usws, &
	2552	surf_usm_v(0)%usws, &
	2553	surf_def_v(1)%usws, &
	2554	surf_lsm_v(1)%usws, &
	2555	surf_usm_v(1)%usws, &
	2556	surf_def_v(2)%usws, &
	2557	surf_lsm_v(2)%usws, &
	2558	surf_usm_v(2)%usws, &
	2559	surf_def_v(3)%usws, &
	2560	surf_lsm_v(3)%usws, &
	2561	surf_usm_v(3)%usws, n_out )
	2562
	2563	CASE ( 'vsws' )
	2564	CALL surface_output_sum_up( surf_def_h(0)%vsws, &
	2565	surf_def_h(1)%vsws, &
	2566	surf_lsm_h%vsws, &
	2567	surf_usm_h%vsws, &
	2568	surf_def_v(0)%vsws, &
	2569	surf_lsm_v(0)%vsws, &
	2570	surf_usm_v(0)%vsws, &
	2571	surf_def_v(1)%vsws, &
	2572	surf_lsm_v(1)%vsws, &
	2573	surf_usm_v(1)%vsws, &
	2574	surf_def_v(2)%vsws, &
	2575	surf_lsm_v(2)%vsws, &
	2576	surf_usm_v(2)%vsws, &
	2577	surf_def_v(3)%vsws, &
	2578	surf_lsm_v(3)%vsws, &
	2579	surf_usm_v(3)%vsws, n_out )
	2580
	2581	CASE ( 'shf' )
	2582	CALL surface_output_sum_up( surf_def_h(0)%shf, &
	2583	surf_def_h(1)%shf, &
	2584	surf_lsm_h%shf, &
	2585	surf_usm_h%shf, &
	2586	surf_def_v(0)%shf, &
	2587	surf_lsm_v(0)%shf, &
	2588	surf_usm_v(0)%shf, &
	2589	surf_def_v(1)%shf, &
	2590	surf_lsm_v(1)%shf, &
	2591	surf_usm_v(1)%shf, &
	2592	surf_def_v(2)%shf, &
	2593	surf_lsm_v(2)%shf, &
	2594	surf_usm_v(2)%shf, &
	2595	surf_def_v(3)%shf, &
	2596	surf_lsm_v(3)%shf, &
	2597	surf_usm_v(3)%shf, n_out )
	2598
	2599	CASE ( 'qsws' )
	2600	CALL surface_output_sum_up( surf_def_h(0)%qsws, &
	2601	surf_def_h(1)%qsws, &
	2602	surf_lsm_h%qsws, &
	2603	surf_usm_h%qsws, &
	2604	surf_def_v(0)%qsws, &
	2605	surf_lsm_v(0)%qsws, &
	2606	surf_usm_v(0)%qsws, &
	2607	surf_def_v(1)%qsws, &
	2608	surf_lsm_v(1)%qsws, &
	2609	surf_usm_v(1)%qsws, &
	2610	surf_def_v(2)%qsws, &
	2611	surf_lsm_v(2)%qsws, &
	2612	surf_usm_v(2)%qsws, &
	2613	surf_def_v(3)%qsws, &
	2614	surf_lsm_v(3)%qsws, &
	2615	surf_usm_v(3)%qsws, n_out )
	2616
	2617	CASE ( 'ssws' )
	2618	CALL surface_output_sum_up( surf_def_h(0)%ssws, &
	2619	surf_def_h(1)%ssws, &
	2620	surf_lsm_h%ssws, &
	2621	surf_usm_h%ssws, &
	2622	surf_def_v(0)%ssws, &
	2623	surf_lsm_v(0)%ssws, &
	2624	surf_usm_v(0)%ssws, &
	2625	surf_def_v(1)%ssws, &
	2626	surf_lsm_v(1)%ssws, &
	2627	surf_usm_v(1)%ssws, &
	2628	surf_def_v(2)%ssws, &
	2629	surf_lsm_v(2)%ssws, &
	2630	surf_usm_v(2)%ssws, &
	2631	surf_def_v(3)%ssws, &
	2632	surf_lsm_v(3)%ssws, &
	2633	surf_usm_v(3)%ssws, n_out )
	2634
	2635	CASE ( 'qcsws' )
	2636	CALL surface_output_sum_up( surf_def_h(0)%qcsws, &
	2637	surf_def_h(1)%qcsws, &
	2638	surf_lsm_h%qcsws, &
	2639	surf_usm_h%qcsws, &
	2640	surf_def_v(0)%qcsws, &
	2641	surf_lsm_v(0)%qcsws, &
	2642	surf_usm_v(0)%qcsws, &
	2643	surf_def_v(1)%qcsws, &
	2644	surf_lsm_v(1)%qcsws, &
	2645	surf_usm_v(1)%qcsws, &
	2646	surf_def_v(2)%qcsws, &
	2647	surf_lsm_v(2)%qcsws, &
	2648	surf_usm_v(2)%qcsws, &
	2649	surf_def_v(3)%qcsws, &
	2650	surf_lsm_v(3)%qcsws, &
	2651	surf_usm_v(3)%qcsws, n_out )
	2652
	2653	CASE ( 'ncsws' )
	2654	CALL surface_output_sum_up( surf_def_h(0)%ncsws, &
	2655	surf_def_h(1)%ncsws, &
	2656	surf_lsm_h%ncsws, &
	2657	surf_usm_h%ncsws, &
	2658	surf_def_v(0)%ncsws, &
	2659	surf_lsm_v(0)%ncsws, &
	2660	surf_usm_v(0)%ncsws, &
	2661	surf_def_v(1)%ncsws, &
	2662	surf_lsm_v(1)%ncsws, &
	2663	surf_usm_v(1)%ncsws, &
	2664	surf_def_v(2)%ncsws, &
	2665	surf_lsm_v(2)%ncsws, &
	2666	surf_usm_v(2)%ncsws, &
	2667	surf_def_v(3)%ncsws, &
	2668	surf_lsm_v(3)%ncsws, &
	2669	surf_usm_v(3)%ncsws, n_out )
	2670
	2671	CASE ( 'qrsws' )
	2672	CALL surface_output_sum_up( surf_def_h(0)%qrsws, &
	2673	surf_def_h(1)%qrsws, &
	2674	surf_lsm_h%qrsws, &
	2675	surf_usm_h%qrsws, &
	2676	surf_def_v(0)%qrsws, &
	2677	surf_lsm_v(0)%qrsws, &
	2678	surf_usm_v(0)%qrsws, &
	2679	surf_def_v(1)%qrsws, &
	2680	surf_lsm_v(1)%qrsws, &
	2681	surf_usm_v(1)%qrsws, &
	2682	surf_def_v(2)%qrsws, &
	2683	surf_lsm_v(2)%qrsws, &
	2684	surf_usm_v(2)%qrsws, &
	2685	surf_def_v(3)%qrsws, &
	2686	surf_lsm_v(3)%qrsws, &
	2687	surf_usm_v(3)%qrsws, n_out )
	2688
	2689	CASE ( 'nrsws' )
	2690	CALL surface_output_sum_up( surf_def_h(0)%nrsws, &
	2691	surf_def_h(1)%nrsws, &
	2692	surf_lsm_h%nrsws, &
	2693	surf_usm_h%nrsws, &
	2694	surf_def_v(0)%nrsws, &
	2695	surf_lsm_v(0)%nrsws, &
	2696	surf_usm_v(0)%nrsws, &
	2697	surf_def_v(1)%nrsws, &
	2698	surf_lsm_v(1)%nrsws, &
	2699	surf_usm_v(1)%nrsws, &
	2700	surf_def_v(2)%nrsws, &
	2701	surf_lsm_v(2)%nrsws, &
	2702	surf_usm_v(2)%nrsws, &
	2703	surf_def_v(3)%nrsws, &
	2704	surf_lsm_v(3)%nrsws, &
	2705	surf_usm_v(3)%nrsws, n_out )
	2706
	2707	CASE ( 'sasws' )
	2708	CALL surface_output_sum_up( surf_def_h(0)%sasws, &
	2709	surf_def_h(1)%sasws, &
	2710	surf_lsm_h%sasws, &
	2711	surf_usm_h%sasws, &
	2712	surf_def_v(0)%sasws, &
	2713	surf_lsm_v(0)%sasws, &
	2714	surf_usm_v(0)%sasws, &
	2715	surf_def_v(1)%sasws, &
	2716	surf_lsm_v(1)%sasws, &
	2717	surf_usm_v(1)%sasws, &
	2718	surf_def_v(2)%sasws, &
	2719	surf_lsm_v(2)%sasws, &
	2720	surf_usm_v(2)%sasws, &
	2721	surf_def_v(3)%sasws, &
	2722	surf_lsm_v(3)%sasws, &
	2723	surf_usm_v(3)%sasws, n_out )
	2724
	2725	CASE ( 'q_surface' )
	2726	CALL surface_output_sum_up( surf_def_h(0)%q_surface, &
	2727	surf_def_h(1)%q_surface, &
	2728	surf_lsm_h%q_surface, &
	2729	surf_usm_h%q_surface, &
	2730	surf_def_v(0)%q_surface, &
	2731	surf_lsm_v(0)%q_surface, &
	2732	surf_usm_v(0)%q_surface, &
	2733	surf_def_v(1)%q_surface, &
	2734	surf_lsm_v(1)%q_surface, &
	2735	surf_usm_v(1)%q_surface, &
	2736	surf_def_v(2)%q_surface, &
	2737	surf_lsm_v(2)%q_surface, &
	2738	surf_usm_v(2)%q_surface, &
	2739	surf_def_v(3)%q_surface, &
	2740	surf_lsm_v(3)%q_surface, &
	2741	surf_usm_v(3)%q_surface, n_out )
	2742
	2743	CASE ( 'theta_surface' )
	2744	CALL surface_output_sum_up( surf_def_h(0)%pt_surface, &
	2745	surf_def_h(1)%pt_surface, &
	2746	surf_lsm_h%pt_surface, &
	2747	surf_usm_h%pt_surface, &
	2748	surf_def_v(0)%pt_surface, &
	2749	surf_lsm_v(0)%pt_surface, &
	2750	surf_usm_v(0)%pt_surface, &
	2751	surf_def_v(1)%pt_surface, &
	2752	surf_lsm_v(1)%pt_surface, &
	2753	surf_usm_v(1)%pt_surface, &
	2754	surf_def_v(2)%pt_surface, &
	2755	surf_lsm_v(2)%pt_surface, &
	2756	surf_usm_v(2)%pt_surface, &
	2757	surf_def_v(3)%pt_surface, &
	2758	surf_lsm_v(3)%pt_surface, &
	2759	surf_usm_v(3)%pt_surface, n_out )
	2760
	2761	CASE ( 'thetav_surface' )
	2762	CALL surface_output_sum_up( surf_def_h(0)%vpt_surface, &
	2763	surf_def_h(1)%vpt_surface, &
	2764	surf_lsm_h%vpt_surface, &
	2765	surf_usm_h%vpt_surface, &
	2766	surf_def_v(0)%vpt_surface, &
	2767	surf_lsm_v(0)%vpt_surface, &
	2768	surf_usm_v(0)%vpt_surface, &
	2769	surf_def_v(1)%vpt_surface, &
	2770	surf_lsm_v(1)%vpt_surface, &
	2771	surf_usm_v(1)%vpt_surface, &
	2772	surf_def_v(2)%vpt_surface, &
	2773	surf_lsm_v(2)%vpt_surface, &
	2774	surf_usm_v(2)%vpt_surface, &
	2775	surf_def_v(3)%vpt_surface, &
	2776	surf_lsm_v(3)%vpt_surface, &
	2777	surf_usm_v(3)%vpt_surface, n_out )
	2778
	2779	CASE ( 'rad_net' )
	2780	CALL surface_output_sum_up( surf_def_h(0)%rad_net, &
	2781	surf_def_h(1)%rad_net, &
	2782	surf_lsm_h%rad_net, &
	2783	surf_usm_h%rad_net, &
	2784	surf_def_v(0)%rad_net, &
	2785	surf_lsm_v(0)%rad_net, &
	2786	surf_usm_v(0)%rad_net, &
	2787	surf_def_v(1)%rad_net, &
	2788	surf_lsm_v(1)%rad_net, &
	2789	surf_usm_v(1)%rad_net, &
	2790	surf_def_v(2)%rad_net, &
	2791	surf_lsm_v(2)%rad_net, &
	2792	surf_usm_v(2)%rad_net, &
	2793	surf_def_v(3)%rad_net, &
	2794	surf_lsm_v(3)%rad_net, &
	2795	surf_usm_v(3)%rad_net, n_out )
	2796
	2797	CASE ( 'rad_lw_in' )
	2798	CALL surface_output_sum_up( surf_def_h(0)%rad_lw_in, &
	2799	surf_def_h(1)%rad_lw_in, &
	2800	surf_lsm_h%rad_lw_in, &
	2801	surf_usm_h%rad_lw_in, &
	2802	surf_def_v(0)%rad_lw_in, &
	2803	surf_lsm_v(0)%rad_lw_in, &
	2804	surf_usm_v(0)%rad_lw_in, &
	2805	surf_def_v(1)%rad_lw_in, &
	2806	surf_lsm_v(1)%rad_lw_in, &
	2807	surf_usm_v(1)%rad_lw_in, &
	2808	surf_def_v(2)%rad_lw_in, &
	2809	surf_lsm_v(2)%rad_lw_in, &
	2810	surf_usm_v(2)%rad_lw_in, &
	2811	surf_def_v(3)%rad_lw_in, &
	2812	surf_lsm_v(3)%rad_lw_in, &
	2813	surf_usm_v(3)%rad_lw_in, n_out )
	2814
	2815	CASE ( 'rad_lw_out' )
	2816	CALL surface_output_sum_up( surf_def_h(0)%rad_lw_out, &
	2817	surf_def_h(1)%rad_lw_out, &
	2818	surf_lsm_h%rad_lw_out, &
	2819	surf_usm_h%rad_lw_out, &
	2820	surf_def_v(0)%rad_lw_out, &
	2821	surf_lsm_v(0)%rad_lw_out, &
	2822	surf_usm_v(0)%rad_lw_out, &
	2823	surf_def_v(1)%rad_lw_out, &
	2824	surf_lsm_v(1)%rad_lw_out, &
	2825	surf_usm_v(1)%rad_lw_out, &
	2826	surf_def_v(2)%rad_lw_out, &
	2827	surf_lsm_v(2)%rad_lw_out, &
	2828	surf_usm_v(2)%rad_lw_out, &
	2829	surf_def_v(3)%rad_lw_out, &
	2830	surf_lsm_v(3)%rad_lw_out, &
	2831	surf_usm_v(3)%rad_lw_out, n_out )
	2832
	2833	CASE ( 'rad_sw_in' )
	2834	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_in, &
	2835	surf_def_h(1)%rad_sw_in, &
	2836	surf_lsm_h%rad_sw_in, &
	2837	surf_usm_h%rad_sw_in, &
	2838	surf_def_v(0)%rad_sw_in, &
	2839	surf_lsm_v(0)%rad_sw_in, &
	2840	surf_usm_v(0)%rad_sw_in, &
	2841	surf_def_v(1)%rad_sw_in, &
	2842	surf_lsm_v(1)%rad_sw_in, &
	2843	surf_usm_v(1)%rad_sw_in, &
	2844	surf_def_v(2)%rad_sw_in, &
	2845	surf_lsm_v(2)%rad_sw_in, &
	2846	surf_usm_v(2)%rad_sw_in, &
	2847	surf_def_v(3)%rad_sw_in, &
	2848	surf_lsm_v(3)%rad_sw_in, &
	2849	surf_usm_v(3)%rad_sw_in, n_out )
	2850
	2851	CASE ( 'rad_sw_out' )
	2852	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_out, &
	2853	surf_def_h(1)%rad_sw_out, &
	2854	surf_lsm_h%rad_sw_out, &
	2855	surf_usm_h%rad_sw_out, &
	2856	surf_def_v(0)%rad_sw_out, &
	2857	surf_lsm_v(0)%rad_sw_out, &
	2858	surf_usm_v(0)%rad_sw_out, &
	2859	surf_def_v(1)%rad_sw_out, &
	2860	surf_lsm_v(1)%rad_sw_out, &
	2861	surf_usm_v(1)%rad_sw_out, &
	2862	surf_def_v(2)%rad_sw_out, &
	2863	surf_lsm_v(2)%rad_sw_out, &
	2864	surf_usm_v(2)%rad_sw_out, &
	2865	surf_def_v(3)%rad_sw_out, &
	2866	surf_lsm_v(3)%rad_sw_out, &
	2867	surf_usm_v(3)%rad_sw_out, n_out )
	2868
	2869	CASE ( 'ghf' )
	2870	CALL surface_output_sum_up( surf_def_h(0)%ghf, &
	2871	surf_def_h(1)%ghf, &
	2872	surf_lsm_h%ghf, &
	2873	surf_usm_h%ghf, &
	2874	surf_def_v(0)%ghf, &
	2875	surf_lsm_v(0)%ghf, &
	2876	surf_usm_v(0)%ghf, &
	2877	surf_def_v(1)%ghf, &
	2878	surf_lsm_v(1)%ghf, &
	2879	surf_usm_v(1)%ghf, &
	2880	surf_def_v(2)%ghf, &
	2881	surf_lsm_v(2)%ghf, &
	2882	surf_usm_v(2)%ghf, &
	2883	surf_def_v(3)%ghf, &
	2884	surf_lsm_v(3)%ghf, &
	2885	surf_usm_v(3)%ghf, n_out )
	2886
	2887	CASE ( 'r_a' )
	2888	CALL surface_output_sum_up( surf_def_h(0)%r_a, &
	2889	surf_def_h(1)%r_a, &
	2890	surf_lsm_h%r_a, &
	2891	surf_usm_h%r_a, &
	2892	surf_def_v(0)%r_a, &
	2893	surf_lsm_v(0)%r_a, &
	2894	surf_usm_v(0)%r_a, &
	2895	surf_def_v(1)%r_a, &
	2896	surf_lsm_v(1)%r_a, &
	2897	surf_usm_v(1)%r_a, &
	2898	surf_def_v(2)%r_a, &
	2899	surf_lsm_v(2)%r_a, &
	2900	surf_usm_v(2)%r_a, &
	2901	surf_def_v(3)%r_a, &
	2902	surf_lsm_v(3)%r_a, &
	2903	surf_usm_v(3)%r_a, n_out )
	2904
	2905	CASE ( 'r_soil' )
	2906	CALL surface_output_sum_up( surf_def_h(0)%r_soil, &
	2907	surf_def_h(1)%r_soil, &
	2908	surf_lsm_h%r_soil, &
	2909	surf_usm_h%r_soil, &
	2910	surf_def_v(0)%r_soil, &
	2911	surf_lsm_v(0)%r_soil, &
	2912	surf_usm_v(0)%r_soil, &
	2913	surf_def_v(1)%r_soil, &
	2914	surf_lsm_v(1)%r_soil, &
	2915	surf_usm_v(1)%r_soil, &
	2916	surf_def_v(2)%r_soil, &
	2917	surf_lsm_v(2)%r_soil, &
	2918	surf_usm_v(2)%r_soil, &
	2919	surf_def_v(3)%r_soil, &
	2920	surf_lsm_v(3)%r_soil, &
	2921	surf_usm_v(3)%r_soil, n_out )
	2922
	2923	CASE ( 'r_canopy' )
	2924	CALL surface_output_sum_up( surf_def_h(0)%r_canopy, &
	2925	surf_def_h(1)%r_canopy, &
	2926	surf_lsm_h%r_canopy, &
	2927	surf_usm_h%r_canopy, &
	2928	surf_def_v(0)%r_canopy, &
	2929	surf_lsm_v(0)%r_canopy, &
	2930	surf_usm_v(0)%r_canopy, &
	2931	surf_def_v(1)%r_canopy, &
	2932	surf_lsm_v(1)%r_canopy, &
	2933	surf_usm_v(1)%r_canopy, &
	2934	surf_def_v(2)%r_canopy, &
	2935	surf_lsm_v(2)%r_canopy, &
	2936	surf_usm_v(2)%r_canopy, &
	2937	surf_def_v(3)%r_canopy, &
	2938	surf_lsm_v(3)%r_canopy, &
	2939	surf_usm_v(3)%r_canopy, n_out )
	2940
	2941	CASE ( 'r_s' )
	2942	CALL surface_output_sum_up( surf_def_h(0)%r_s, &
	2943	surf_def_h(1)%r_s, &
	2944	surf_lsm_h%r_s, &
	2945	surf_usm_h%r_s, &
	2946	surf_def_v(0)%r_s, &
	2947	surf_lsm_v(0)%r_s, &
	2948	surf_usm_v(0)%r_s, &
	2949	surf_def_v(1)%r_s, &
	2950	surf_lsm_v(1)%r_s, &
	2951	surf_usm_v(1)%r_s, &
	2952	surf_def_v(2)%r_s, &
	2953	surf_lsm_v(2)%r_s, &
	2954	surf_usm_v(2)%r_s, &
	2955	surf_def_v(3)%r_s, &
	2956	surf_lsm_v(3)%r_s, &
	2957	surf_usm_v(3)%r_s, n_out )
	2958
[3572]	2959
	2960	CASE ( 'rad_sw_dir' )
	2961	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_dir, &
	2962	surf_def_h(1)%rad_sw_dir, &
	2963	surf_lsm_h%rad_sw_dir, &
	2964	surf_usm_h%rad_sw_dir, &
	2965	surf_def_v(0)%rad_sw_dir, &
	2966	surf_lsm_v(0)%rad_sw_dir, &
	2967	surf_usm_v(0)%rad_sw_dir, &
	2968	surf_def_v(1)%rad_sw_dir, &
	2969	surf_lsm_v(1)%rad_sw_dir, &
	2970	surf_usm_v(1)%rad_sw_dir, &
	2971	surf_def_v(2)%rad_sw_dir, &
	2972	surf_lsm_v(2)%rad_sw_dir, &
	2973	surf_usm_v(2)%rad_sw_dir, &
	2974	surf_def_v(3)%rad_sw_dir, &
	2975	surf_lsm_v(3)%rad_sw_dir, &
	2976	surf_usm_v(3)%rad_sw_dir, n_out )
	2977	CASE ( 'rad_sw_dif' )
	2978	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_dif, &
	2979	surf_def_h(1)%rad_sw_dif, &
	2980	surf_lsm_h%rad_sw_dif, &
	2981	surf_usm_h%rad_sw_dif, &
	2982	surf_def_v(0)%rad_sw_dif, &
	2983	surf_lsm_v(0)%rad_sw_dif, &
	2984	surf_usm_v(0)%rad_sw_dif, &
	2985	surf_def_v(1)%rad_sw_dif, &
	2986	surf_lsm_v(1)%rad_sw_dif, &
	2987	surf_usm_v(1)%rad_sw_dif, &
	2988	surf_def_v(2)%rad_sw_dif, &
	2989	surf_lsm_v(2)%rad_sw_dif, &
	2990	surf_usm_v(2)%rad_sw_dif, &
	2991	surf_def_v(3)%rad_sw_dif, &
	2992	surf_lsm_v(3)%rad_sw_dif, &
	2993	surf_usm_v(3)%rad_sw_dif, n_out )
	2994
	2995	CASE ( 'rad_sw_ref' )
	2996	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_ref, &
	2997	surf_def_h(1)%rad_sw_ref, &
	2998	surf_lsm_h%rad_sw_ref, &
	2999	surf_usm_h%rad_sw_ref, &
	3000	surf_def_v(0)%rad_sw_ref, &
	3001	surf_lsm_v(0)%rad_sw_ref, &
	3002	surf_usm_v(0)%rad_sw_ref, &
	3003	surf_def_v(1)%rad_sw_ref, &
	3004	surf_lsm_v(1)%rad_sw_ref, &
	3005	surf_usm_v(1)%rad_sw_ref, &
	3006	surf_def_v(2)%rad_sw_ref, &
	3007	surf_lsm_v(2)%rad_sw_ref, &
	3008	surf_usm_v(2)%rad_sw_ref, &
	3009	surf_def_v(3)%rad_sw_ref, &
	3010	surf_lsm_v(3)%rad_sw_ref, &
	3011	surf_usm_v(3)%rad_sw_ref, n_out )
	3012
	3013	CASE ( 'rad_sw_res' )
	3014	CALL surface_output_sum_up( surf_def_h(0)%rad_sw_res, &
	3015	surf_def_h(1)%rad_sw_res, &
	3016	surf_lsm_h%rad_sw_res, &
	3017	surf_usm_h%rad_sw_res, &
	3018	surf_def_v(0)%rad_sw_res, &
	3019	surf_lsm_v(0)%rad_sw_res, &
	3020	surf_usm_v(0)%rad_sw_res, &
	3021	surf_def_v(1)%rad_sw_res, &
	3022	surf_lsm_v(1)%rad_sw_res, &
	3023	surf_usm_v(1)%rad_sw_res, &
	3024	surf_def_v(2)%rad_sw_res, &
	3025	surf_lsm_v(2)%rad_sw_res, &
	3026	surf_usm_v(2)%rad_sw_res, &
	3027	surf_def_v(3)%rad_sw_res, &
	3028	surf_lsm_v(3)%rad_sw_res, &
	3029	surf_usm_v(3)%rad_sw_res, n_out )
	3030
	3031	CASE ( 'rad_lw_dif' )
	3032	CALL surface_output_sum_up( surf_def_h(0)%rad_lw_dif, &
	3033	surf_def_h(1)%rad_lw_dif, &
	3034	surf_lsm_h%rad_lw_dif, &
	3035	surf_usm_h%rad_lw_dif, &
	3036	surf_def_v(0)%rad_lw_dif, &
	3037	surf_lsm_v(0)%rad_lw_dif, &
	3038	surf_usm_v(0)%rad_lw_dif, &
	3039	surf_def_v(1)%rad_lw_dif, &
	3040	surf_lsm_v(1)%rad_lw_dif, &
	3041	surf_usm_v(1)%rad_lw_dif, &
	3042	surf_def_v(2)%rad_lw_dif, &
	3043	surf_lsm_v(2)%rad_lw_dif, &
	3044	surf_usm_v(2)%rad_lw_dif, &
	3045	surf_def_v(3)%rad_lw_dif, &
	3046	surf_lsm_v(3)%rad_lw_dif, &
	3047	surf_usm_v(3)%rad_lw_dif, n_out )
	3048
	3049	CASE ( 'rad_lw_ref' )
	3050	CALL surface_output_sum_up( surf_def_h(0)%rad_lw_ref, &
	3051	surf_def_h(1)%rad_lw_ref, &
	3052	surf_lsm_h%rad_lw_ref, &
	3053	surf_usm_h%rad_lw_ref, &
	3054	surf_def_v(0)%rad_lw_ref, &
	3055	surf_lsm_v(0)%rad_lw_ref, &
	3056	surf_usm_v(0)%rad_lw_ref, &
	3057	surf_def_v(1)%rad_lw_ref, &
	3058	surf_lsm_v(1)%rad_lw_ref, &
	3059	surf_usm_v(1)%rad_lw_ref, &
	3060	surf_def_v(2)%rad_lw_ref, &
	3061	surf_lsm_v(2)%rad_lw_ref, &
	3062	surf_usm_v(2)%rad_lw_ref, &
	3063	surf_def_v(3)%rad_lw_ref, &
	3064	surf_lsm_v(3)%rad_lw_ref, &
	3065	surf_usm_v(3)%rad_lw_ref, n_out )
	3066
	3067	CASE ( 'rad_lw_res' )
	3068	CALL surface_output_sum_up( surf_def_h(0)%rad_lw_res, &
	3069	surf_def_h(1)%rad_lw_res, &
	3070	surf_lsm_h%rad_lw_res, &
	3071	surf_usm_h%rad_lw_res, &
	3072	surf_def_v(0)%rad_lw_res, &
	3073	surf_lsm_v(0)%rad_lw_res, &
	3074	surf_usm_v(0)%rad_lw_res, &
	3075	surf_def_v(1)%rad_lw_res, &
	3076	surf_lsm_v(1)%rad_lw_res, &
	3077	surf_usm_v(1)%rad_lw_res, &
	3078	surf_def_v(2)%rad_lw_res, &
	3079	surf_lsm_v(2)%rad_lw_res, &
	3080	surf_usm_v(2)%rad_lw_res, &
	3081	surf_def_v(3)%rad_lw_res, &
	3082	surf_lsm_v(3)%rad_lw_res, &
	3083	surf_usm_v(3)%rad_lw_res, n_out )
	3084
[3421]	3085	END SELECT
	3086	ENDDO
	3087
	3088
	3089	END SUBROUTINE surface_output_averaging
	3090
	3091	!------------------------------------------------------------------------------!
	3092	! Description:
	3093	! ------------
	3094	!> Sum-up the surface data for average output variables.
	3095	!------------------------------------------------------------------------------!
[3494]	3096	SUBROUTINE surface_output_sum_up( var_def_h0, var_def_h1, &
	3097	var_lsm_h, var_usm_h, &
	3098	var_def_v0, var_lsm_v0, var_usm_v0, &
	3099	var_def_v1, var_lsm_v1, var_usm_v1, &
	3100	var_def_v2, var_lsm_v2, var_usm_v2, &
	3101	var_def_v3, var_lsm_v3, var_usm_v3, n_out )
[3421]	3102
	3103	IMPLICIT NONE
	3104
	3105	INTEGER(iwp) :: m !< running index for surface elements
	3106	INTEGER(iwp) :: n_out !< index for output variable
	3107	INTEGER(iwp) :: n_surf !< running index for surface elements
	3108
	3109	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_h0 !< output variable at upward-facing default-type surfaces
	3110	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_h1 !< output variable at downward-facing default-type surfaces
	3111	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_h !< output variable at upward-facing natural-type surfaces
	3112	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_h !< output variable at upward-facing urban-type surfaces
	3113	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v0 !< output variable at northward-facing default-type surfaces
	3114	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v1 !< output variable at southward-facing default-type surfaces
	3115	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v2 !< output variable at eastward-facing default-type surfaces
	3116	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v3 !< output variable at westward-facing default-type surfaces
	3117	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v0 !< output variable at northward-facing natural-type surfaces
	3118	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v1 !< output variable at southward-facing natural-type surfaces
	3119	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v2 !< output variable at eastward-facing natural-type surfaces
	3120	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v3 !< output variable at westward-facing natural-type surfaces
	3121	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v0 !< output variable at northward-facing urban-type surfaces
	3122	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v1 !< output variable at southward-facing urban-type surfaces
	3123	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v2 !< output variable at eastward-facing urban-type surfaces
	3124	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v3 !< output variable at westward-facing urban-type surfaces
	3125
	3126	!
	3127	!-- Set counter variable to zero before the variable is written to
	3128	!-- the output array.
	3129	n_surf = 0
	3130
	3131	!
	3132	!-- Write the horizontal surfaces.
	3133	!-- Before each the variable is written to the output data structure, first
	3134	!-- check if the variable for the respective surface type is defined.
	3135	!-- If a variable is not defined, skip the block and increment the counter
	3136	!-- variable by the number of surface elements of this type. Usually this
	3137	!-- is zere, however, there might be the situation that e.g. urban surfaces
	3138	!-- are defined but the respective variable is not allocated for this surface
	3139	!-- type. To write the data on the exact position, increment the counter.
	3140	IF ( ALLOCATED( var_def_h0 ) ) THEN
	3141	DO m = 1, surf_def_h(0)%ns
	3142	n_surf = n_surf + 1
	3143	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3144	+ var_def_h0(m)
	3145	ENDDO
	3146	ELSE
	3147	n_surf = n_surf + surf_def_h(0)%ns
	3148	ENDIF
	3149	IF ( ALLOCATED( var_def_h1 ) ) THEN
	3150	DO m = 1, surf_def_h(1)%ns
	3151	n_surf = n_surf + 1
	3152	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3153	+ var_def_h1(m)
	3154	ENDDO
	3155	ELSE
	3156	n_surf = n_surf + surf_def_h(1)%ns
	3157	ENDIF
	3158	IF ( ALLOCATED( var_lsm_h ) ) THEN
	3159	DO m = 1, surf_lsm_h%ns
	3160	n_surf = n_surf + 1
	3161	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3162	+ var_lsm_h(m)
	3163	ENDDO
	3164	ELSE
	3165	n_surf = n_surf + surf_lsm_h%ns
	3166	ENDIF
	3167	IF ( ALLOCATED( var_usm_h ) ) THEN
	3168	DO m = 1, surf_usm_h%ns
	3169	n_surf = n_surf + 1
	3170	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3171	+ var_usm_h(m)
	3172	ENDDO
	3173	ELSE
	3174	n_surf = n_surf + surf_usm_h%ns
	3175	ENDIF
	3176	!
	3177	!-- Write northward-facing
	3178	IF ( ALLOCATED( var_def_v0 ) ) THEN
	3179	DO m = 1, surf_def_v(0)%ns
	3180	n_surf = n_surf + 1
	3181	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3182	+ var_def_v0(m)
	3183	ENDDO
	3184	ELSE
	3185	n_surf = n_surf + surf_def_v(0)%ns
	3186	ENDIF
	3187	IF ( ALLOCATED( var_lsm_v0 ) ) THEN
	3188	DO m = 1, surf_lsm_v(0)%ns
	3189	n_surf = n_surf + 1
	3190	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3191	+ var_lsm_v0(m)
	3192	ENDDO
	3193	ELSE
	3194	n_surf = n_surf + surf_lsm_v(0)%ns
	3195	ENDIF
	3196	IF ( ALLOCATED( var_usm_v0 ) ) THEN
	3197	DO m = 1, surf_usm_v(0)%ns
	3198	n_surf = n_surf + 1
	3199	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3200	+ var_usm_v0(m)
	3201	ENDDO
	3202	ELSE
	3203	n_surf = n_surf + surf_usm_v(0)%ns
	3204	ENDIF
	3205	!
	3206	!-- Write southward-facing
	3207	IF ( ALLOCATED( var_def_v1 ) ) THEN
	3208	DO m = 1, surf_def_v(1)%ns
	3209	n_surf = n_surf + 1
	3210	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3211	+ var_def_v1(m)
	3212	ENDDO
	3213	ELSE
	3214	n_surf = n_surf + surf_def_v(1)%ns
	3215	ENDIF
	3216	IF ( ALLOCATED( var_lsm_v1 ) ) THEN
	3217	DO m = 1, surf_lsm_v(1)%ns
	3218	n_surf = n_surf + 1
	3219	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3220	+ var_lsm_v1(m)
	3221	ENDDO
	3222	ELSE
	3223	n_surf = n_surf + surf_lsm_v(1)%ns
	3224	ENDIF
	3225	IF ( ALLOCATED( var_usm_v1 ) ) THEN
	3226	DO m = 1, surf_usm_v(1)%ns
	3227	n_surf = n_surf + 1
	3228	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3229	+ var_usm_v1(m)
	3230	ENDDO
	3231	ELSE
	3232	n_surf = n_surf + surf_usm_v(1)%ns
	3233	ENDIF
	3234	!
	3235	!-- Write eastward-facing
	3236	IF ( ALLOCATED( var_def_v2 ) ) THEN
	3237	DO m = 1, surf_def_v(2)%ns
	3238	n_surf = n_surf + 1
	3239	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3240	+ var_def_v2(m)
	3241	ENDDO
	3242	ELSE
	3243	n_surf = n_surf + surf_def_v(2)%ns
	3244	ENDIF
	3245	IF ( ALLOCATED( var_lsm_v2 ) ) THEN
	3246	DO m = 1, surf_lsm_v(2)%ns
	3247	n_surf = n_surf + 1
	3248	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3249	+ var_lsm_v2(m)
	3250	ENDDO
	3251	ELSE
	3252	n_surf = n_surf + surf_lsm_v(2)%ns
	3253	ENDIF
	3254	IF ( ALLOCATED( var_usm_v2 ) ) THEN
	3255	DO m = 1, surf_usm_v(2)%ns
	3256	n_surf = n_surf + 1
	3257	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3258	+ var_usm_v2(m)
	3259	ENDDO
	3260	ELSE
	3261	n_surf = n_surf + surf_usm_v(2)%ns
	3262	ENDIF
	3263	!
	3264	!-- Write westward-facing
	3265	IF ( ALLOCATED( var_def_v3 ) ) THEN
	3266	DO m = 1, surf_def_v(3)%ns
	3267	n_surf = n_surf + 1
	3268	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3269	+ var_def_v3(m)
	3270	ENDDO
	3271	ELSE
	3272	n_surf = n_surf + surf_def_v(3)%ns
	3273	ENDIF
	3274	IF ( ALLOCATED( var_lsm_v3 ) ) THEN
	3275	DO m = 1, surf_lsm_v(3)%ns
	3276	n_surf = n_surf + 1
	3277	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3278	+ var_lsm_v3(m)
	3279	ENDDO
	3280	ELSE
	3281	n_surf = n_surf + surf_lsm_v(3)%ns
	3282	ENDIF
	3283	IF ( ALLOCATED( var_usm_v3 ) ) THEN
	3284	DO m = 1, surf_usm_v(3)%ns
	3285	n_surf = n_surf + 1
	3286	surfaces%var_av(n_surf,n_out) = surfaces%var_av(n_surf,n_out) &
	3287	+ var_usm_v3(m)
	3288	ENDDO
	3289	ELSE
	3290	n_surf = n_surf + surf_usm_v(3)%ns
	3291	ENDIF
	3292
	3293	END SUBROUTINE surface_output_sum_up
	3294
	3295	!------------------------------------------------------------------------------!
	3296	! Description:
	3297	! ------------
	3298	!> Collect the surface data from different types and different orientation.
	3299	!------------------------------------------------------------------------------!
	3300	SUBROUTINE surface_output_collect( var_def_h0, var_def_h1, &
[3494]	3301	var_lsm_h, var_usm_h, &
	3302	var_def_v0, var_lsm_v0, var_usm_v0, &
	3303	var_def_v1, var_lsm_v1, var_usm_v1, &
	3304	var_def_v2, var_lsm_v2, var_usm_v2, &
	3305	var_def_v3, var_lsm_v3, var_usm_v3 )
[3421]	3306
	3307	IMPLICIT NONE
	3308
	3309	INTEGER(iwp) :: m !< running index for surface elements
	3310	INTEGER(iwp) :: n_surf !< running index for surface elements
	3311
	3312	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_h0 !< output variable at upward-facing default-type surfaces
	3313	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_h1 !< output variable at downward-facing default-type surfaces
	3314	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_h !< output variable at upward-facing natural-type surfaces
	3315	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_h !< output variable at upward-facing urban-type surfaces
	3316	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v0 !< output variable at northward-facing default-type surfaces
	3317	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v1 !< output variable at southward-facing default-type surfaces
	3318	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v2 !< output variable at eastward-facing default-type surfaces
	3319	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_def_v3 !< output variable at westward-facing default-type surfaces
	3320	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v0 !< output variable at northward-facing natural-type surfaces
	3321	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v1 !< output variable at southward-facing natural-type surfaces
	3322	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v2 !< output variable at eastward-facing natural-type surfaces
	3323	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_lsm_v3 !< output variable at westward-facing natural-type surfaces
	3324	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v0 !< output variable at northward-facing urban-type surfaces
	3325	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v1 !< output variable at southward-facing urban-type surfaces
	3326	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v2 !< output variable at eastward-facing urban-type surfaces
	3327	REAL(wp), DIMENSION(:), ALLOCATABLE, INTENT(IN) :: var_usm_v3 !< output variable at westward-facing urban-type surfaces
	3328
	3329	!
	3330	!-- Set counter variable to zero before the variable is written to
	3331	!-- the output array.
	3332	n_surf = 0
	3333
	3334	!
	3335	!-- Write the horizontal surfaces.
	3336	!-- Before each the variable is written to the output data structure, first
	3337	!-- check if the variable for the respective surface type is defined.
	3338	!-- If a variable is not defined, skip the block and increment the counter
	3339	!-- variable by the number of surface elements of this type. Usually this
	3340	!-- is zere, however, there might be the situation that e.g. urban surfaces
	3341	!-- are defined but the respective variable is not allocated for this surface
	3342	!-- type. To write the data on the exact position, increment the counter.
	3343	IF ( ALLOCATED( var_def_h0 ) ) THEN
	3344	DO m = 1, surf_def_h(0)%ns
	3345	n_surf = n_surf + 1
	3346	surfaces%var_out(n_surf) = var_def_h0(m)
	3347	ENDDO
	3348	ELSE
	3349	n_surf = n_surf + surf_def_h(0)%ns
	3350	ENDIF
	3351	IF ( ALLOCATED( var_def_h1 ) ) THEN
	3352	DO m = 1, surf_def_h(1)%ns
	3353	n_surf = n_surf + 1
	3354	surfaces%var_out(n_surf) = var_def_h1(m)
	3355	ENDDO
	3356	ELSE
	3357	n_surf = n_surf + surf_def_h(1)%ns
	3358	ENDIF
	3359	IF ( ALLOCATED( var_lsm_h ) ) THEN
	3360	DO m = 1, surf_lsm_h%ns
	3361	n_surf = n_surf + 1
	3362	surfaces%var_out(n_surf) = var_lsm_h(m)
	3363	ENDDO
	3364	ELSE
	3365	n_surf = n_surf + surf_lsm_h%ns
	3366	ENDIF
	3367	IF ( ALLOCATED( var_usm_h ) ) THEN
	3368	DO m = 1, surf_usm_h%ns
	3369	n_surf = n_surf + 1
	3370	surfaces%var_out(n_surf) = var_usm_h(m)
	3371	ENDDO
	3372	ELSE
	3373	n_surf = n_surf + surf_usm_h%ns
	3374	ENDIF
	3375	!
	3376	!-- Write northward-facing
	3377	IF ( ALLOCATED( var_def_v0 ) ) THEN
	3378	DO m = 1, surf_def_v(0)%ns
	3379	n_surf = n_surf + 1
	3380	surfaces%var_out(n_surf) = var_def_v0(m)
	3381	ENDDO
	3382	ELSE
	3383	n_surf = n_surf + surf_def_v(0)%ns
	3384	ENDIF
	3385	IF ( ALLOCATED( var_lsm_v0 ) ) THEN
	3386	DO m = 1, surf_lsm_v(0)%ns
	3387	n_surf = n_surf + 1
	3388	surfaces%var_out(n_surf) = var_lsm_v0(m)
	3389	ENDDO
	3390	ELSE
	3391	n_surf = n_surf + surf_lsm_v(0)%ns
	3392	ENDIF
	3393	IF ( ALLOCATED( var_usm_v0 ) ) THEN
	3394	DO m = 1, surf_usm_v(0)%ns
	3395	n_surf = n_surf + 1
	3396	surfaces%var_out(n_surf) = var_usm_v0(m)
	3397	ENDDO
	3398	ELSE
	3399	n_surf = n_surf + surf_usm_v(0)%ns
	3400	ENDIF
	3401	!
	3402	!-- Write southward-facing
	3403	IF ( ALLOCATED( var_def_v1 ) ) THEN
	3404	DO m = 1, surf_def_v(1)%ns
	3405	n_surf = n_surf + 1
	3406	surfaces%var_out(n_surf) = var_def_v1(m)
	3407	ENDDO
	3408	ELSE
	3409	n_surf = n_surf + surf_def_v(1)%ns
	3410	ENDIF
	3411	IF ( ALLOCATED( var_lsm_v1 ) ) THEN
	3412	DO m = 1, surf_lsm_v(1)%ns
	3413	n_surf = n_surf + 1
	3414	surfaces%var_out(n_surf) = var_lsm_v1(m)
	3415	ENDDO
	3416	ELSE
	3417	n_surf = n_surf + surf_lsm_v(1)%ns
	3418	ENDIF
	3419	IF ( ALLOCATED( var_usm_v1 ) ) THEN
	3420	DO m = 1, surf_usm_v(1)%ns
	3421	n_surf = n_surf + 1
	3422	surfaces%var_out(n_surf) = var_usm_v1(m)
	3423	ENDDO
	3424	ELSE
	3425	n_surf = n_surf + surf_usm_v(1)%ns
	3426	ENDIF
	3427	!
	3428	!-- Write eastward-facing
	3429	IF ( ALLOCATED( var_def_v2 ) ) THEN
	3430	DO m = 1, surf_def_v(2)%ns
	3431	n_surf = n_surf + 1
	3432	surfaces%var_out(n_surf) = var_def_v2(m)
	3433	ENDDO
	3434	ELSE
	3435	n_surf = n_surf + surf_def_v(2)%ns
	3436	ENDIF
	3437	IF ( ALLOCATED( var_lsm_v2 ) ) THEN
	3438	DO m = 1, surf_lsm_v(2)%ns
	3439	n_surf = n_surf + 1
	3440	surfaces%var_out(n_surf) = var_lsm_v2(m)
	3441	ENDDO
	3442	ELSE
	3443	n_surf = n_surf + surf_lsm_v(2)%ns
	3444	ENDIF
	3445	IF ( ALLOCATED( var_usm_v2 ) ) THEN
	3446	DO m = 1, surf_usm_v(2)%ns
	3447	n_surf = n_surf + 1
	3448	surfaces%var_out(n_surf) = var_usm_v2(m)
	3449	ENDDO
	3450	ELSE
	3451	n_surf = n_surf + surf_usm_v(2)%ns
	3452	ENDIF
	3453	!
	3454	!-- Write westward-facing
	3455	IF ( ALLOCATED( var_def_v3 ) ) THEN
	3456	DO m = 1, surf_def_v(3)%ns
	3457	n_surf = n_surf + 1
	3458	surfaces%var_out(n_surf) = var_def_v3(m)
	3459	ENDDO
	3460	ELSE
	3461	n_surf = n_surf + surf_def_v(3)%ns
	3462	ENDIF
	3463	IF ( ALLOCATED( var_lsm_v3 ) ) THEN
	3464	DO m = 1, surf_lsm_v(3)%ns
	3465	n_surf = n_surf + 1
	3466	surfaces%var_out(n_surf) = var_lsm_v3(m)
	3467	ENDDO
	3468	ELSE
	3469	n_surf = n_surf + surf_lsm_v(3)%ns
	3470	ENDIF
	3471	IF ( ALLOCATED( var_usm_v3 ) ) THEN
	3472	DO m = 1, surf_usm_v(3)%ns
	3473	n_surf = n_surf + 1
	3474	surfaces%var_out(n_surf) = var_usm_v3(m)
	3475	ENDDO
	3476	ELSE
	3477	n_surf = n_surf + surf_usm_v(3)%ns
	3478	ENDIF
	3479
	3480	END SUBROUTINE surface_output_collect
	3481
	3482
	3483
	3484	! SUBROUTINE data_output_surf_driver (time)
	3485	! IMPLICIT NONE
	3486	!
	3487	! REAL(kind=wp),INTENT(IN) :: time
	3488	! CHARACTER(LEN=64) :: filename
	3489	! INTEGER,SAVE :: file_number=1
	3490	!
	3491	! if(output_mode == "legacy_VTK") then
	3492	! CALL prepare_vtk_data
	3493	!
	3494	! write(filename,'(a,i5.5,a)') 'surface_output_',file_number,'.vtk'
	3495	! if(myid == surf_io_pe) write(9,*) 'Legacy VTK surface output to file: ',trim(filename)
	3496	! file_number = file_number+1
	3497	!
	3498	! CALL write_ascii_vtk (filename, time)
	3499	! END IF
	3500	!
	3501	! RETURN
	3502	! END SUBROUTINE data_output_surf_driver
	3503
	3504
	3505	!------------------------------------------------------------------------------!
	3506	! Description:
	3507	! ------------
	3508	!> Parin for output of surface parameters
	3509	!------------------------------------------------------------------------------!
	3510	SUBROUTINE surface_output_parin
	3511
	3512	IMPLICIT NONE
	3513
	3514	CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file
	3515
	3516
	3517	NAMELIST /surface_output_parameters/ &
	3518	averaging_interval_surf, data_output_surf, &
	3519	dt_dosurf, dt_dosurf_av, &
	3520	skip_time_dosurf, skip_time_dosurf_av
	3521
	3522	line = ' '
	3523
	3524	!
	3525	!-- Try to find the namelist
	3526	REWIND ( 11 )
	3527	line = ' '
	3528	DO WHILE ( INDEX( line, '&surface_output_parameters' ) == 0 )
	3529	READ ( 11, '(A)', END=14 ) line
	3530	ENDDO
	3531	BACKSPACE ( 11 )
	3532
	3533	!
	3534	!-- Read namelist
	3535	READ ( 11, surface_output_parameters, ERR = 10 )
	3536	!
	3537	!-- Set flag that indicates that surface data output is switched on
	3538	surface_data_output = .TRUE.
	3539	GOTO 14
	3540
	3541	10 BACKSPACE( 11 )
	3542	READ( 11 , '(A)') line
	3543	CALL parin_fail_message( 'surface_output_parameters', line )
	3544
	3545	14 CONTINUE
	3546
	3547
	3548	END SUBROUTINE surface_output_parin
	3549
	3550
	3551	!------------------------------------------------------------------------------!
	3552	! Description:
	3553	! ------------
	3554	!> Check the input parameters for consistency. Further pre-process the given
	3555	!> output variables, i.e. separate them into average and non-average output
	3556	!> variables and map them onto internal output array.
	3557	!------------------------------------------------------------------------------!
	3558	SUBROUTINE surface_output_check_parameters
	3559
	3560	USE control_parameters, &
[3614]	3561	ONLY: averaging_interval, dt_data_output, message_string
[3421]	3562
	3563	IMPLICIT NONE
	3564
	3565	CHARACTER(LEN=100) :: trimvar !< dummy for single output variable
	3566
	3567	INTEGER(iwp) :: av !< id indicating average or non-average data output
	3568	INTEGER(iwp) :: ilen !< string length
	3569	INTEGER(iwp) :: n_out !< running index for number of output variables
	3570
	3571	!
	3572	!-- Check the average interval
	3573	IF ( averaging_interval_surf == 9999999.9_wp ) THEN
	3574	averaging_interval_surf = averaging_interval
	3575	ENDIF
	3576	!
	3577	!-- Set the default data-output interval dt_data_output if necessary
	3578	IF ( dt_dosurf == 9999999.9_wp ) dt_dosurf = dt_data_output
	3579	IF ( dt_dosurf_av == 9999999.9_wp ) dt_dosurf_av = dt_data_output
	3580
	3581	IF ( averaging_interval_surf > dt_dosurf_av ) THEN
	3582	WRITE( message_string, * ) 'averaging_interval_surf = ', &
	3583	averaging_interval_surf, ' must be <= dt_dosurf_av = ', &
	3584	dt_dosurf_av
	3585	CALL message( 'surface_output_check_parameters', &
	3586	'PA0087', 1, 2, 0, 6, 0 )
	3587	ENDIF
	3588
	3589	!
	3590	!-- Count number of output variables and separate output strings for
	3591	!-- average and non-average output variables.
	3592	n_out = 0
	3593	DO WHILE ( data_output_surf(n_out+1)(1:1) /= ' ' )
	3594
	3595	n_out = n_out + 1
	3596	ilen = LEN_TRIM( data_output_surf(n_out) )
	3597	trimvar = TRIM( data_output_surf(n_out) )
	3598
	3599	!
	3600	!-- Check for data averaging
	3601	av = 0
	3602	IF ( ilen > 3 ) THEN
	3603	IF ( data_output_surf(n_out)(ilen-2:ilen) == '_av' ) THEN
	3604	trimvar = data_output_surf(n_out)(1:ilen-3)
	3605	av = 1
	3606	ENDIF
	3607	ENDIF
	3608
	3609	dosurf_no(av) = dosurf_no(av) + 1
	3610	dosurf(av,dosurf_no(av)) = TRIM( trimvar )
	3611
	3612	!
	3613	!-- Check if all output variables are known
	3614	SELECT CASE ( TRIM( trimvar ) )
	3615
	3616	CASE ( 'css', 'cssws', 'qsws_liq', 'qsws_soil', 'qsws_veg' )
	3617	message_string = TRIM( trimvar ) // &
	3618	' is not yet implemented in the surface output'
	3619	CALL message( 'surface_output_check_parameters', &
	3620	'PA0087', 1, 2, 0, 6, 0 )
	3621
	3622	CASE ( 'us', 'ts', 'qs', 'ss', 'qcs', 'ncs', 'qrs', 'nrs', 'ol' )
	3623
	3624	CASE ( 'z0', 'z0h', 'z0q' )
	3625
	3626	CASE ( 'theta1', 'qv1', 'thetav1' )
	3627
	3628	CASE ( 'usws', 'vsws' )
	3629
	3630	CASE ( 'qcsws', 'ncsws', 'qrsws', 'nrsws', 'sasws' )
	3631
	3632	CASE ( 'shf', 'qsws', 'ssws' )
	3633
	3634	CASE ( 'q_surface', 'theta_surface', 'thetav_surface' )
	3635
	3636	CASE ( 'rad_net' )
	3637
[3572]	3638	CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_lw_dif', 'rad_lw_ref', &
	3639	'rad_lw_res' )
[3421]	3640
[3572]	3641	CASE ( 'rad_sw_in', 'rad_sw_out', 'rad_sw_dif', 'rad_sw_ref', &
	3642	'rad_sw_res', 'rad_sw_dir' )
	3643
[3421]	3644	CASE ( 'ghf' )
	3645
	3646	CASE ( 'r_a', 'r_canopy', 'r_soil', 'r_s' )
	3647
	3648	CASE DEFAULT
	3649	message_string = TRIM( trimvar ) // &
	3650	' is not part of the surface output'
	3651	CALL message( 'surface_output_check_parameters', &
	3652	'PA0087', 1, 2, 0, 6, 0 )
	3653	END SELECT
	3654
	3655	ENDDO
	3656
	3657	END SUBROUTINE surface_output_check_parameters
[3494]	3658
	3659
	3660	!------------------------------------------------------------------------------!
	3661	! Description:
	3662	! ------------
	3663	!> Last action.
	3664	!------------------------------------------------------------------------------!
	3665	SUBROUTINE surface_output_last_action( av )
	3666
	3667	USE control_parameters, &
	3668	ONLY: io_blocks, io_group
[3421]	3669
[3494]	3670	USE pegrid, &
	3671	ONLY: comm2d, ierr
	3672
	3673	IMPLICIT NONE
	3674
	3675	INTEGER(iwp) :: av !< id indicating average or non-average data output
	3676	INTEGER(iwp) :: i !< loop index
	3677
	3678	!
	3679	!-- Return, if nothing to output
	3680	IF ( dosurf_no(av) == 0 ) RETURN
	3681	!
	3682	!-- Open files
	3683	CALL check_open( 25+av )
	3684	!
	3685	!-- Write time coordinate
	3686	DO i = 0, io_blocks-1
	3687	IF ( i == io_group ) THEN
	3688	WRITE ( 25+av ) LEN_TRIM( 'END' )
	3689	WRITE ( 25+av ) 'END'
	3690	ENDIF
	3691	#if defined( __parallel )
	3692	CALL MPI_BARRIER( comm2d, ierr )
	3693	#endif
	3694	ENDDO
	3695
	3696	END SUBROUTINE surface_output_last_action
	3697
[3421]	3698	!--Private Subroutines (Only called from tinside his module)
	3699
	3700
	3701
	3702	! SUBROUTINE write_ascii_vtk (filename,time)
	3703	! IMPLICIT NONE
	3704	! CHARACTER(LEN=*),INTENT(IN) :: filename
	3705	! REAL(kind=wp),INTENT(IN) :: time
	3706	!
	3707	! INTEGER :: file_id=251
	3708	! INTEGER :: n
	3709	!
	3710	! IF(myid == surf_io_pe) then
	3711	! open(file_id,file=TRIM(filename)) ! Please use check_open in the final version
	3712	!
	3713	! write(file_id,'(a)') "# vtk DataFile Version 3.0"
	3714	! write(file_id,'(a,f8.2,a)') "legacy vtk File generated by PALM, simulation time = ",time," sec"
	3715	! write(file_id,'(a)') "ASCII"
	3716	!
	3717	! write(file_id,'(a)') "DATASET POLYDATA"
	3718	! write(file_id,'(a,i5,a)') "POINTS ",npo," float"
	3719	! do n=1,npo
	3720	! write (file_id,'(8f10.1)') points(1:3,n)
	3721	! end do
	3722	!
	3723	! write(file_id,'(a,8i10)') "POLYGONS ",npg, 5*npg
	3724	! do n=1,npg
	3725	! write(file_id,'(8i10)') polygons(:,n)
	3726	! end do
	3727	!
	3728	! write(file_id,'(a,i10)') "CELL_DATA ",npg
	3729	! write(file_id,'(a,i10)') "SCALARS cell_scalars float 1 "
	3730	! write(file_id,'(a,i10)') "LOOKUP_TABLE default "
	3731	!
	3732	! do n=1,npg
	3733	! write(file_id,'(f10.2)') cell_scalars(n)
	3734	! end do
	3735	!
	3736	! close (file_id)
	3737	!
	3738	! DEALLOCATE (points,polygons,cell_scalars)
	3739	! ENDIF
	3740	!
	3741	! RETURN
	3742	! END SUBROUTINE write_ascii_vtk
	3743
	3744	END MODULE surface_output_mod

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |