Home

Context Navigation

source: palm/trunk/SOURCE/modules.f90 @ 415

Last change on this file since 415 was 412, checked in by raasch, 15 years ago
data asignments for dvrp arrays switched back to old settings due to runtime problems; type problem in mpi_waitall (poisfft_hybrid) fixed for mpi2 libraries
Property svn:keywords set to `Id`
File size: 63.5 KB

Rev	Line
[1]	1	MODULE advection
	2
	3
	4	!------------------------------------------------------------------------------!
[343]	5	! Current revisions:
[1]	6	! -----------------
[405]	7	! -var_ts: replaced by dots_max,
	8	! initial data assignments to some dvrp arrays changed due to error messages
	9	! from gfortran compiler
[411]	10	! +large_scale_subsidence, ws_vertical_gradient, ws_vertical_gradient_level,
	11	! ws_vertical_gradient_level_ind, w_subs
[392]	12	!
[410]	13	! Branch revisions:
	14	! -----------------
	15	! Masked data output: + dt_domask, mask_01~20_x\|y\|z, mask_01~20_x\|y\|z_loop,
	16	! mask_scale\|_x\|y\|z, masks, netcdf_format_mask[_av], skip_time_domask
	17	!
[392]	18	! Former revisions:
	19	! -----------------
	20	! $Id: modules.f90 412 2009-12-14 07:14:46Z raasch $
	21	!
	22	! 388 2009-09-23 09:40:33Z raasch
[388]	23	! +prho, prho_1
[366]	24	! +bc_lr_cyc, bc_ns_cyc
[345]	25	! +output_for_t0
[343]	26	! translation error of actual -> current revisions fixed
	27	! +q* in dots_label, dots_unit. increased dots_num respectively
[288]	28	! typographical error in dots_unit fixed
[264]	29	! +clip_dvrp_*, cluster_size, color_interval, dvrpsize_interval, dvrp_overlap,
[284]	30	! dvrp_total_overlap, groundplate_color, local_dvrserver_running, n*_dvrp,
	31	! interval_*_dvrp_prt, isosurface_color, particle_color, particle_dvrpsize,
	32	! topography color in dvrp_variables,
[336]	33	! vertical_particle_advection is a 1d-array,
[240]	34	! +canyon_height, canyon_width_x, canyon_width_y, canyon_wall_left,
[291]	35	! canyon_wall_south, conserve_volume_flow_mode, coupling_start_time,
	36	! dp_external, dp_level_b, dp_level_ind_b, dp_smooth, dp_smooth_factor, dpdxy,
[305]	37	! run_coupled, time_since_reference_point, u_bulk, v_bulk in control_parameters,
	38	! default value of grid_matching changed to strict
[367]	39	! +shf_av, qsws_av
[198]	40	!
[226]	41	! 217 2008-12-09 18:00:48Z letzel
	42	! +topography_grid_convention
	43	! some dvrp-variables changed to single precision, variables for dvrp-mode
	44	! pathlines added, +target_id, abort_mode, message_string
	45	!
[198]	46	! 197 2008-09-16 15:29:03Z raasch
[189]	47	! allow 100 spectra levels instead of 10 for consistency with
[198]	48	! define_netcdf_header, +canopy_heat_flux, cthf, lai,
[153]	49	! +leaf_surface_concentration, scalar_exchange_coefficient, sec, sls
[163]	50	! +hor_index_bounds, hor_index_bounds_previous_run, id_inflow, id_recycling,
[151]	51	! inflow_damping_*, mean_inflow_profiles, numprocs_previous_run, nx_on_file,
	52	! ny_on_file, offset_ocean_*, recycling_plane, recycling_width, turbulent_inflow
[146]	53	! -myid_char_14
[139]	54	!
	55	! 138 2007-11-28 10:03:58Z letzel
[138]	56	! +drag_coefficient, pch_index, lad_surface, lad_vertical_gradient,
	57	! lad_vertical_gradient_level, plant_canopy, lad, lad_s, lad_u, lad_v,
	58	! lad_w, cdc, lad_vertical_gradient_level_ind, canopy_mode
[132]	59	! +dt_sort_particles, ngp_2dh_s_inner, time_sort_particles, flags,
	60	! wall_flags_1..10, wall_humidityflux(0:4), wall_qflux(0:4),
	61	! wall_salinityflux(0:4), wall_scalarflux(0:4)
[110]	62	!
	63	! 108 2007-08-24 15:10:38Z letzel
[103]	64	! +comm_inter, constant_top_momentumflux, coupling_char, coupling_mode,
[108]	65	! coupling_mode_remote, c_u, c_v, c_w, dt_coupling, e_init, humidity_remote,
	66	! ngp_xy, nxlu, nysv, port_name, qswst_remote, terminate_coupled,
	67	! terminate_coupled_remote, time_coupling, top_momentumflux_u\|v, type_xy,
	68	! uswst, vswst
[77]	69	!
[98]	70	! 97 2007-06-21 08:23:15Z raasch
	71	! +atmos_ocean_sign, ocean, r, + salinity variables
	72	! defaults of .._vertical_gradient_levels changed from -1.0 to -9999999.9
	73	! hydro_press renamed hyp, use_pt_reference renamed use_reference
	74	!
[90]	75	! 89 2007-05-25 12:08:31Z raasch
	76	! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index,
	77	! dopr_initial_index and dopr_unit enlarged,
	78	! var_hom and var_sum renamed pr_palm
	79	!
[83]	80	! 82 2007-04-16 15:40:52Z raasch
	81	! +return_addres, return_username
	82	! Cpp-directive lcmuk renamed lc
	83	!
[77]	84	! 75 2007-03-22 09:54:05Z raasch
[72]	85	! +arrays precipitation_amount, precipitation_rate, precipitation_rate_av,
[73]	86	! rif_wall, z0_av, +arrays u_m_l, u_m_r, etc. for radiation boundary conditions,
[63]	87	! +loop_optimization, netcdf_64bit_3d, zu_s_inner, zw_w_inner, id_var_zusi_*,
	88	! id_var_zwwi_*, ts_value, u_nzb_p1_for_vfc, v_nzb_p1_for_vfc, pt_reference,
[75]	89	! use_pt_reference, precipitation_amount_interval, revision
	90	! +age_m in particle_type, moisture renamed humidity,
	91	! -data_output_ts, dots_n, uvmean_outflow, uxrp, vynp,
	92	! arrays dots_label and dots_unit now dimensioned with dots_max,
	93	! setting of palm version moved to main program
[46]	94	!
	95	! 37 2007-03-01 08:33:54Z raasch
[19]	96	! +constant_top_heatflux, top_heatflux, use_top_fluxes, +arrays for top fluxes,
	97	! +nzt_diff, default of bc_pt_t renamed "initial_gradient"
	98	! Bugfix: p is not a pointer
[1]	99	!
[3]	100	! RCS Log replace by Id keyword, revision history cleaned up
	101	!
[1]	102	! Revision 1.95 2007/02/11 13:18:30 raasch
	103	! version 3.1b (last under RCS control)
	104	!
	105	! Revision 1.1 1997/07/24 11:21:26 raasch
	106	! Initial revision
	107	!
	108	!
	109	! Description:
	110	! ------------
	111	! Definition of variables for special advection schemes
	112	!------------------------------------------------------------------------------!
	113
	114	REAL :: spl_gamma_x, spl_gamma_y
	115
	116	REAL, DIMENSION(:), ALLOCATABLE :: aex, bex, dex, eex, spl_z_x, spl_z_y
	117	REAL, DIMENSION(:,:), ALLOCATABLE :: spl_tri_x, spl_tri_y, spl_tri_zu, &
	118	spl_tri_zw
	119
	120	SAVE
	121
	122	END MODULE advection
	123
	124
	125
	126
	127	MODULE array_kind
	128
	129	!------------------------------------------------------------------------------!
	130	! Description:
	131	! ------------
	132	! Definition of type parameters (used for the definition of single or double
	133	! precision variables)
	134	!------------------------------------------------------------------------------!
	135
	136	INTEGER, PARAMETER :: dpk = SELECTED_REAL_KIND( 12 ), &
	137	spk = SELECTED_REAL_KIND( 6 )
	138
	139	SAVE
	140
	141	END MODULE array_kind
	142
	143
	144
	145
	146	MODULE arrays_3d
	147
	148	!------------------------------------------------------------------------------!
	149	! Description:
	150	! ------------
	151	! Definition of all arrays defined on the computational grid
	152	!------------------------------------------------------------------------------!
	153
	154	USE array_kind
	155
	156	REAL, DIMENSION(:), ALLOCATABLE :: &
[151]	157	ddzu, dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, km_damp_x, &
	158	km_damp_y, lad, l_grid, pt_init, q_init, rdf, sa_init, ug, u_init, &
[411]	159	u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
[1]	160
	161	REAL, DIMENSION(:,:), ALLOCATABLE :: &
[151]	162	c_u, c_v, c_w, dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg, &
	163	mean_inflow_profiles, pt_slope_ref, qs, qswst_remote, ts, us, z0
[1]	164
	165	REAL, DIMENSION(:,:), ALLOCATABLE, TARGET :: &
[95]	166	qsws_1, qsws_2, qswst_1, qswst_2, rif_1, rif_2, saswsb_1, saswst_1, &
[102]	167	shf_1, shf_2, tswst_1, tswst_2, usws_1, usws_2, uswst_1, uswst_2, &
	168	vsws_1, vsws_2, vswst_1, vswst_2
[1]	169
	170	REAL, DIMENSION(:,:), POINTER :: &
[95]	171	qsws, qsws_m, qswst, qswst_m, rif, rif_m, saswsb, saswst, shf, &
[102]	172	shf_m, tswst, tswst_m, usws, uswst, usws_m, uswst_m, vsws, vswst, &
	173	vsws_m, vswst_m
[1]	174
	175	REAL, DIMENSION(:,:,:), ALLOCATABLE :: &
[153]	176	canopy_heat_flux, cdc, d, diss, lad_s, lad_u, lad_v, lad_w, lai, &
	177	l_wall, sec, sls, tend, u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n, &
	178	v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s
[1]	179
	180	REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
	181	ql_v, ql_vp
	182
	183	REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
[388]	184	e_1, e_2, e_3, kh_1, kh_2, km_1, km_2, p, prho_1, pt_1, pt_2, pt_3, &
	185	q_1, q_2, q_3, ql_1, ql_2, rho_1, sa_1, sa_2, sa_3, u_1, u_2, u_3, &
	186	v_1, v_2, v_3, vpt_1, vpt_2, w_1, w_2, w_3
[1]	187
	188	REAL, DIMENSION(:,:,:), POINTER :: &
[388]	189	e, e_m, e_p, kh, kh_m, km, km_m, prho, pt, pt_m, pt_p, q, q_m, q_p, &
	190	ql, ql_c, rho, sa, sa_p, te_m, tpt_m, tq_m, tsa_m, tu_m, tv_m, tw_m, &
	191	u, u_m, u_p, v, v_m, v_p, vpt, vpt_m, w, w_m, w_p
[51]	192
[75]	193	REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: rif_wall
[1]	194
	195	SAVE
	196
	197	END MODULE arrays_3d
	198
	199
	200
	201
	202	MODULE averaging
	203
	204	!------------------------------------------------------------------------------!
	205	! Description:
	206	! ------------
	207	! Definition of variables needed for time-averaging of 2d/3d data
	208	!------------------------------------------------------------------------------!
	209
[354]	210	REAL, DIMENSION(:,:), ALLOCATABLE :: lwp_av, precipitation_rate_av, &
	211	qsws_av, shf_av,ts_av, us_av, z0_av
[1]	212
	213	REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
	214	e_av, p_av, pc_av, pr_av, pt_av, q_av, ql_av, ql_c_av, ql_v_av, &
[96]	215	ql_vp_av, qv_av, rho_av, s_av, sa_av, u_av, v_av, vpt_av, w_av
[1]	216
	217	END MODULE averaging
	218
	219
	220
	221
	222	MODULE cloud_parameters
	223
	224	!------------------------------------------------------------------------------!
	225	! Description:
	226	! ------------
	227	! Definition of variables and constants for cloud physics
	228	!------------------------------------------------------------------------------!
	229
	230	REAL :: b_cond, cp = 1005.0, diff_coeff_l = 0.23E-4, &
	231	effective_coll_efficiency, l_d_cp, l_d_r, l_d_rv, l_v = 2.5E+06, &
	232	mass_of_solute, molecular_weight_of_solute, &
	233	prec_time_const = 0.001, ql_crit = 0.0005, rho_l = 1.0E3, &
[95]	234	r_d = 287.0, r_v = 461.51, thermal_conductivity_l = 2.43E-2
[1]	235
[95]	236	REAL, DIMENSION(:), ALLOCATABLE :: pt_d_t, t_d_pt
[1]	237
[72]	238	REAL, DIMENSION(:,:), ALLOCATABLE :: precipitation_amount, &
	239	precipitation_rate
	240
[1]	241	SAVE
	242
	243	END MODULE cloud_parameters
	244
	245
	246
	247
	248	MODULE constants
	249
	250	!------------------------------------------------------------------------------!
	251	! Description:
	252	! ------------
	253	! Definition of general constants
	254	!------------------------------------------------------------------------------!
	255
	256	REAL :: pi = 3.141592654
	257
	258	SAVE
	259
	260	END MODULE constants
	261
	262
	263
	264
	265	MODULE control_parameters
	266
	267	!------------------------------------------------------------------------------!
	268	! Description:
	269	! ------------
	270	! Definition of parameters for program control
	271	!------------------------------------------------------------------------------!
	272
	273	TYPE plot_precision
	274	CHARACTER (LEN=6) :: variable
	275	INTEGER :: precision
	276	END TYPE plot_precision
	277
	278	TYPE(plot_precision), DIMENSION(100) :: plot_3d_precision = &
	279	(/ plot_precision( 'u', 2 ), plot_precision( 'v', 2 ), &
	280	plot_precision( 'w', 2 ), plot_precision( 'p', 5 ), &
	281	plot_precision( 'pt', 2 ), &
	282	( plot_precision( ' ', 1 ), i9 = 1,95 ) /)
	283
	284	TYPE file_status
	285	LOGICAL :: opened, opened_before
	286	END TYPE file_status
	287
	288	TYPE(file_status), DIMENSION(200) :: openfile = file_status(.FALSE.,.FALSE.)
	289
	290
[410]	291	INTEGER, PARAMETER :: mask_xyz_dimension = 100, max_masks = 20
	292
[213]	293	CHARACTER (LEN=1) :: cycle_mg = 'w', timestep_reason = ' '
	294	CHARACTER (LEN=2) :: coupling_char = ''
	295	CHARACTER (LEN=5) :: write_binary = 'false'
[305]	296	CHARACTER (LEN=6) :: grid_matching = 'strict'
[213]	297	CHARACTER (LEN=8) :: run_date, run_time
	298	CHARACTER (LEN=9) :: simulated_time_chr
[217]	299	CHARACTER (LEN=11) :: topography_grid_convention = ' '
[213]	300	CHARACTER (LEN=12) :: version = ' ', revision = ' '
[241]	301	CHARACTER (LEN=16) :: conserve_volume_flow_mode = 'default', &
	302	loop_optimization = 'default', &
[213]	303	momentum_advec = 'pw-scheme', &
	304	psolver = 'poisfft', &
	305	scalar_advec = 'pw-scheme'
	306	CHARACTER (LEN=20) :: bc_e_b = 'neumann', bc_lr = 'cyclic', &
	307	bc_ns = 'cyclic', bc_p_b = 'neumann', &
	308	bc_p_t = 'dirichlet', bc_pt_b = 'dirichlet', &
	309	bc_pt_t = 'initial_gradient', &
	310	bc_q_b = 'dirichlet', bc_q_t = 'neumann', &
	311	bc_s_b = 'dirichlet', bc_s_t = 'neumann', &
	312	bc_sa_t = 'neumann', &
	313	bc_uv_b = 'dirichlet', bc_uv_t = 'dirichlet', &
	314	canopy_mode = 'block', &
	315	coupling_mode = 'uncoupled', &
	316	coupling_mode_remote = 'uncoupled', &
	317	dissipation_1d = 'as_in_3d_model', &
	318	fft_method = 'system-specific', &
	319	mixing_length_1d = 'as_in_3d_model', &
	320	random_generator = 'numerical-recipes', &
	321	return_addres, return_username, &
	322	timestep_scheme = 'runge-kutta-3'
	323	CHARACTER (LEN=40) :: avs_data_file, topography = 'flat'
	324	CHARACTER (LEN=64) :: host = ' '
	325	CHARACTER (LEN=80) :: log_message, run_identifier
	326	CHARACTER (LEN=100) :: initializing_actions = ' ', run_description_header
	327	CHARACTER (LEN=1000) :: message_string = ' '
[1]	328
	329	CHARACTER (LEN=7), DIMENSION(100) :: do3d_comp_prec = ' '
	330	CHARACTER (LEN=10), DIMENSION(10) :: data_output_format = ' '
[354]	331	CHARACTER (LEN=11), DIMENSION(100) :: data_output = ' ', &
[1]	332	data_output_user = ' ', doav = ' '
[410]	333	CHARACTER (LEN=10), DIMENSION(100) :: &
	334	data_output_mask_01 = ' ', data_output_mask_02 = ' ', &
	335	data_output_mask_03 = ' ', data_output_mask_04 = ' ', &
	336	data_output_mask_05 = ' ', data_output_mask_06 = ' ', &
	337	data_output_mask_07 = ' ', data_output_mask_08 = ' ', &
	338	data_output_mask_09 = ' ', data_output_mask_10 = ' ', &
	339	data_output_mask_11 = ' ', data_output_mask_12 = ' ', &
	340	data_output_mask_13 = ' ', data_output_mask_14 = ' ', &
	341	data_output_mask_15 = ' ', data_output_mask_16 = ' ', &
	342	data_output_mask_17 = ' ', data_output_mask_18 = ' ', &
	343	data_output_mask_19 = ' ', data_output_mask_20 = ' ', &
	344	data_output_mask_user_01 = ' ', data_output_mask_user_02 = ' ', &
	345	data_output_mask_user_03 = ' ', data_output_mask_user_04 = ' ', &
	346	data_output_mask_user_05 = ' ', data_output_mask_user_06 = ' ', &
	347	data_output_mask_user_07 = ' ', data_output_mask_user_08 = ' ', &
	348	data_output_mask_user_09 = ' ', data_output_mask_user_10 = ' ', &
	349	data_output_mask_user_11 = ' ', data_output_mask_user_12 = ' ', &
	350	data_output_mask_user_13 = ' ', data_output_mask_user_14 = ' ', &
	351	data_output_mask_user_15 = ' ', data_output_mask_user_16 = ' ', &
	352	data_output_mask_user_17 = ' ', data_output_mask_user_18 = ' ', &
	353	data_output_mask_user_19 = ' ', data_output_mask_user_20 = ' '
	354
[87]	355	CHARACTER (LEN=10), DIMENSION(300) :: data_output_pr = ' '
	356	CHARACTER (LEN=10), DIMENSION(200) :: data_output_pr_user = ' '
[410]	357	CHARACTER (LEN=20), DIMENSION(50) :: netcdf_precision = ' '
[1]	358
[410]	359	CHARACTER (LEN=10), DIMENSION(max_masks,0:1,100) :: domask = ' '
[1]	360	CHARACTER (LEN=10), DIMENSION(0:1,100) :: do2d = ' ', do3d = ' '
	361
[216]	362	INTEGER :: abort_mode = 1, average_count_pr = 0, average_count_sp = 0, &
[1]	363	average_count_3d = 0, current_timestep_number = 0, &
	364	dist_range = 0, disturbance_level_ind_b, &
	365	disturbance_level_ind_t, doav_n = 0, dopr_n = 0, &
	366	dopr_time_count = 0, dopts_time_count = 0, &
[48]	367	dosp_time_count = 0, dots_time_count = 0, &
[1]	368	do2d_xy_n = 0, do2d_xz_n = 0, do2d_yz_n = 0, do3d_avs_n = 0, &
[240]	369	dp_level_ind_b = 0, &
[1]	370	dvrp_filecount = 0, dz_stretch_level_index, gamma_mg, &
	371	grid_level, ibc_e_b, ibc_p_b, ibc_p_t, ibc_pt_b, ibc_pt_t, &
[95]	372	ibc_q_b, ibc_q_t, ibc_sa_t, ibc_uv_b, ibc_uv_t, &
[1]	373	inflow_disturbance_begin = -1, inflow_disturbance_end = -1, &
	374	intermediate_timestep_count, intermediate_timestep_count_max, &
[410]	375	iran = -1234567, last_dt_change = 0, masks = 0, &
	376	maximum_grid_level, &
[87]	377	max_pr_user = 0, mgcycles = 0, mg_cycles = -1, &
[410]	378	mg_switch_to_pe0_level = 0, mid, ngsrb = 2, nsor = 20, &
[87]	379	nsor_ini = 100, n_sor, normalizing_region = 0, &
[89]	380	nz_do1d, nz_do3d = -9999, outflow_damping_width = -1, &
[240]	381	pch_index = 0, prt_time_count = 0, recycling_plane, runnr = 0, &
[151]	382	skip_do_avs = 0, terminate_coupled = 0, &
	383	terminate_coupled_remote = 0, timestep_count = 0
[1]	384
	385	INTEGER :: dist_nxl(0:1), dist_nxr(0:1), dist_nyn(0:1), dist_nys(0:1), &
	386	do2d_no(0:1) = 0, do2d_xy_time_count(0:1), &
	387	do2d_xz_time_count(0:1), do2d_yz_time_count(0:1), &
	388	do3d_no(0:1) = 0, do3d_time_count(0:1), &
[410]	389	domask_no(max_masks,0:1) = 0, domask_time_count(max_masks,0:1),&
[138]	390	lad_vertical_gradient_level_ind(10) = -9999, &
[410]	391	mask_size(max_masks,3) = -1, mask_size_l(max_masks,3) = -1, &
	392	mask_start_l(max_masks,3) = -1, &
	393	nc_format_mask(1:max_masks,0:1), &
	394	netcdf_format_mask(1:max_masks) = 1, &
	395	netcdf_format_mask_av(1:max_masks) = 1, &
[1]	396	pt_vertical_gradient_level_ind(10) = -9999, &
	397	q_vertical_gradient_level_ind(10) = -9999, &
[94]	398	sa_vertical_gradient_level_ind(10) = -9999, &
[1]	399	section(100,3), section_xy(100) = -9999, &
	400	section_xz(100) = -9999, section_yz(100) = -9999, &
	401	ug_vertical_gradient_level_ind(10) = -9999, &
[411]	402	vg_vertical_gradient_level_ind(10) = -9999, &
	403	ws_vertical_gradient_level_ind(10) = -9999
[1]	404
	405	INTEGER, DIMENSION(:), ALLOCATABLE :: grid_level_count
	406
[410]	407	INTEGER, DIMENSION(:,:), ALLOCATABLE :: mask_i, mask_j, mask_k
	408	INTEGER, DIMENSION(:,:), ALLOCATABLE :: &
	409	mask_i_global, mask_j_global, mask_k_global
	410
[1]	411	LOGICAL :: adjust_mixing_length = .FALSE., avs_output = .FALSE., &
[366]	412	bc_lr_cyc =.TRUE., bc_ns_cyc = .TRUE., &
[1]	413	call_psolver_at_all_substeps = .TRUE., &
	414	cloud_droplets = .FALSE., cloud_physics = .FALSE., &
	415	conserve_volume_flow = .FALSE., constant_diffusion = .FALSE., &
[19]	416	constant_heatflux = .TRUE., constant_top_heatflux = .TRUE., &
[103]	417	constant_top_momentumflux = .FALSE., &
[95]	418	constant_top_salinityflux = .TRUE., &
[19]	419	constant_waterflux = .TRUE., create_disturbances = .TRUE., &
	420	cut_spline_overshoot = .TRUE., &
[1]	421	data_output_2d_on_each_pe = .TRUE., do2d_at_begin = .FALSE., &
	422	do3d_at_begin = .FALSE., do3d_compress = .FALSE., &
[240]	423	do_sum = .FALSE., dp_external = .FALSE., dp_smooth = .FALSE., &
	424	dt_changed = .FALSE., dt_fixed = .FALSE., &
[1]	425	disturbance_created = .FALSE., &
	426	first_call_advec_particles = .TRUE., &
[410]	427	force_print_header = .FALSE., format_parallel_io = .FALSE., &
	428	galilei_transformation = .FALSE.,&
[108]	429	humidity = .FALSE., humidity_remote = .FALSE., &
	430	inflow_l = .FALSE., inflow_n = .FALSE., &
[411]	431	inflow_r = .FALSE., inflow_s = .FALSE., &
	432	iso2d_output = .FALSE., large_scale_subsidence = .FALSE., &
[75]	433	mg_switch_to_pe0 = .FALSE., &
[1]	434	netcdf_output = .FALSE., netcdf_64bit = .FALSE., &
[94]	435	netcdf_64bit_3d = .TRUE., ocean = .FALSE., &
[1]	436	outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., &
	437	outflow_s = .FALSE., passive_scalar = .FALSE., &
[138]	438	plant_canopy = .FALSE., &
[1]	439	prandtl_layer = .TRUE., precipitation = .FALSE., &
	440	profil_output = .FALSE., radiation = .FALSE., &
	441	random_heatflux = .FALSE., run_control_header = .FALSE., &
[291]	442	run_coupled = .TRUE., sloping_surface = .FALSE., &
	443	stop_dt = .FALSE., terminate_run = .FALSE., &
	444	turbulent_inflow = .FALSE., &
[151]	445	use_prior_plot1d_parameters = .FALSE., use_reference = .FALSE.,&
	446	use_surface_fluxes = .FALSE., use_top_fluxes = .FALSE., &
	447	use_ug_for_galilei_tr = .TRUE., use_upstream_for_tke = .FALSE.,&
	448	wall_adjustment = .TRUE.
[1]	449
	450	LOGICAL :: data_output_xy(0:1) = .FALSE., data_output_xz(0:1) = .FALSE., &
	451	data_output_yz(0:1) = .FALSE.
	452
	453	REAL :: advected_distance_x = 0.0, advected_distance_y = 0.0, &
	454	alpha_surface = 0.0, asselin_filter_factor = 0.1, &
[97]	455	atmos_ocean_sign = 1.0, &
[1]	456	averaging_interval = 0.0, averaging_interval_pr = 9999999.9, &
	457	averaging_interval_sp = 9999999.9, bc_pt_t_val, bc_q_t_val, &
[95]	458	bottom_salinityflux = 0.0, &
[1]	459	building_height = 50.0, building_length_x = 50.0, &
	460	building_length_y = 50.0, building_wall_left = 9999999.9, &
[240]	461	building_wall_south = 9999999.9, canyon_height = 50.0, &
	462	canyon_width_x = 9999999.9, canyon_width_y = 9999999.9, &
	463	canyon_wall_left = 9999999.9, canyon_wall_south = 9999999.9, &
[291]	464	cthf = 0.0, cfl_factor = -1.0, cos_alpha_surface, &
	465	coupling_start_time, disturbance_amplitude = 0.25, &
[97]	466	disturbance_energy_limit = 0.01, &
	467	disturbance_level_b = -9999999.9, &
	468	disturbance_level_t = -9999999.9, &
[240]	469	dp_level_b = 0.0, drag_coefficient = 0.0, &
[97]	470	dt = -1.0, dt_averaging_input = 0.0, &
[102]	471	dt_averaging_input_pr = 9999999.9, dt_coupling = 9999999.9, &
	472	dt_data_output = 9999999.9, &
[1]	473	dt_data_output_av = 9999999.9, dt_disturb = 9999999.9, &
	474	dt_dopr = 9999999.9, dt_dopr_listing = 9999999.9, &
	475	dt_dopts = 9999999.9, dt_dosp = 9999999.9, dt_dots = 9999999.9, &
	476	dt_do2d_xy = 9999999.9, dt_do2d_xz = 9999999.9, &
	477	dt_do2d_yz = 9999999.9, dt_do3d = 9999999.9, dt_dvrp = 9999999.9, &
	478	dt_max = 20.0, dt_prel = 9999999.9, dt_restart = 9999999.9, &
	479	dt_run_control = 60.0, dt_3d = -1.0, dz = -1.0, &
	480	dz_max = 9999999.9, dz_stretch_factor = 1.08, &
[108]	481	dz_stretch_level = 100000.0, e_init = 0.0, e_min = 0.0, &
	482	end_time = 0.0, &
[151]	483	f = 0.0, fs = 0.0, g = 9.81, inflow_damping_height = 9999999.9, &
	484	inflow_damping_width = 9999999.9, kappa = 0.4, km_constant = -1.0,&
[153]	485	km_damp_max = -1.0, lad_surface = 0.0, &
	486	leaf_surface_concentration = 0.0, long_filter_factor = 0.0, &
[410]	487	mask_scale_x = 1.0, mask_scale_y = 1.0, mask_scale_z = 1.0, &
	488	maximum_cpu_time_allowed = 0.0, &
	489	molecular_viscosity = 1.461E-5, &
[1]	490	old_dt = 1.0E-10, omega = 7.29212E-5, omega_sor = 1.8, &
	491	overshoot_limit_e = 0.0, overshoot_limit_pt = 0.0, &
	492	overshoot_limit_u = 0.0, overshoot_limit_v = 0.0, &
	493	overshoot_limit_w = 0.0, particle_maximum_age = 9999999.9, &
[72]	494	phi = 55.0, prandtl_number = 1.0, &
[97]	495	precipitation_amount_interval = 9999999.9, prho_reference, &
	496	pt_reference = 9999999.9, pt_slope_offset = 0.0, &
	497	pt_surface = 300.0, pt_surface_initial_change = 0.0, &
	498	q_surface = 0.0, q_surface_initial_change = 0.0, &
	499	rayleigh_damping_factor = -1.0, rayleigh_damping_height = -1.0, &
[151]	500	recycling_width = 9999999.9, residual_limit = 1.0E-4, &
	501	restart_time = 9999999.9, rho_reference, rho_surface, &
[241]	502	rif_max = 1.0, rif_min = -5.0, roughness_length = 0.1, &
	503	sa_surface = 35.0, scalar_exchange_coefficient = 0.0, &
[95]	504	simulated_time = 0.0, simulated_time_at_begin, sin_alpha_surface, &
[1]	505	skip_time_data_output = 0.0, skip_time_data_output_av = 9999999.9,&
	506	skip_time_dopr = 9999999.9, skip_time_dosp = 9999999.9, &
	507	skip_time_do2d_xy = 9999999.9, skip_time_do2d_xz = 9999999.9, &
	508	skip_time_do2d_yz = 9999999.9, skip_time_do3d = 9999999.9, &
	509	surface_heatflux = 9999999.9, surface_pressure = 1013.25, &
	510	surface_scalarflux = 0.0, surface_waterflux = 0.0, &
	511	s_surface = 0.0, s_surface_initial_change = 0.0, &
[102]	512	termination_time_needed = -1.0, time_coupling = 0.0, &
[241]	513	time_disturb = 0.0, time_dopr = 0.0, time_dopr_av = 0.0, &
	514	time_dopr_listing = 0.0, time_dopts = 0.0, time_dosp = 0.0, &
	515	time_dosp_av = 0.0, time_dots = 0.0, time_do2d_xy = 0.0, &
	516	time_do2d_xz = 0.0, time_do2d_yz = 0.0, time_do3d = 0.0, &
	517	time_do_av = 0.0, time_do_sla = 0.0, time_dvrp = 0.0, &
	518	time_prel = 0.0, time_restart = 9999999.9, time_run_control = 0.0,&
[291]	519	time_since_reference_point, top_heatflux = 9999999.9, &
	520	top_momentumflux_u = 9999999.9, &
[103]	521	top_momentumflux_v = 9999999.9, top_salinityflux = 9999999.9, &
[241]	522	ug_surface = 0.0, u_bulk = 0.0, u_gtrans = 0.0, &
[19]	523	ups_limit_e = 0.0, ups_limit_pt = 0.0, ups_limit_u = 0.0, &
	524	ups_limit_v = 0.0, ups_limit_w = 0.0, vg_surface = 0.0, &
[241]	525	v_bulk = 0.0, v_gtrans = 0.0, wall_adjustment_factor = 1.8, &
	526	z_max_do1d = -1.0, z_max_do1d_normalized = -1.0, z_max_do2d = -1.0
[1]	527
	528	REAL :: do2d_xy_last_time(0:1) = -1.0, do2d_xz_last_time(0:1) = -1.0, &
[240]	529	do2d_yz_last_time(0:1) = -1.0, dpdxy(1:2) = 0.0, &
[410]	530	dt_domask(max_masks) = 9999999.9, lad_vertical_gradient(10) = 0.0,&
[138]	531	lad_vertical_gradient_level(10) = -9999999.9, &
[410]	532	mask_scale(3), &
[138]	533	pt_vertical_gradient(10) = 0.0, &
[94]	534	pt_vertical_gradient_level(10) = -9999999.9, &
[1]	535	q_vertical_gradient(10) = 0.0, &
	536	q_vertical_gradient_level(10) = -1.0, &
	537	s_vertical_gradient(10) = 0.0, &
[95]	538	s_vertical_gradient_level(10) = -1.0, &
[94]	539	sa_vertical_gradient(10) = 0.0, &
[410]	540	sa_vertical_gradient_level(10) = -9999999.9, &
	541	skip_time_domask(max_masks) = 9999999.9, threshold(20) = 0.0, &
	542	time_domask(max_masks) = 0.0, &
[1]	543	tsc(10) = (/ 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /), &
	544	ug_vertical_gradient(10) = 0.0, &
[97]	545	ug_vertical_gradient_level(10) = -9999999.9, &
[1]	546	vg_vertical_gradient(10) = 0.0, &
[97]	547	vg_vertical_gradient_level(10) = -9999999.9, &
[1]	548	volume_flow(1:2) = 0.0, volume_flow_area(1:2) = 0.0, &
[129]	549	volume_flow_initial(1:2) = 0.0, wall_heatflux(0:4) = 0.0, &
	550	wall_humidityflux(0:4) = 0.0, wall_qflux(0:4) = 0.0, &
[411]	551	wall_salinityflux(0:4) = 0.0, wall_scalarflux(0:4) = 0.0, &
	552	ws_vertical_gradient(10) = 0.0, &
	553	ws_vertical_gradient_level(10) = -9999999.9
[1]	554
[240]	555	REAL, DIMENSION(:), ALLOCATABLE :: dp_smooth_factor
[1]	556
[410]	557	REAL, DIMENSION(mask_xyz_dimension) :: &
	558	mask_01_x = -1.0, mask_02_x = -1.0, mask_03_x = -1.0, mask_04_x = -1.0, &
	559	mask_05_x = -1.0, mask_06_x = -1.0, mask_07_x = -1.0, mask_08_x = -1.0, &
	560	mask_09_x = -1.0, mask_10_x = -1.0, mask_11_x = -1.0, mask_12_x = -1.0, &
	561	mask_13_x = -1.0, mask_14_x = -1.0, mask_15_x = -1.0, mask_16_x = -1.0, &
	562	mask_17_x = -1.0, mask_18_x = -1.0, mask_19_x = -1.0, mask_20_x = -1.0, &
	563	mask_01_y = -1.0, mask_02_y = -1.0, mask_03_y = -1.0, mask_04_y = -1.0, &
	564	mask_05_y = -1.0, mask_06_y = -1.0, mask_07_y = -1.0, mask_08_y = -1.0, &
	565	mask_09_y = -1.0, mask_10_y = -1.0, mask_11_y = -1.0, mask_12_y = -1.0, &
	566	mask_13_y = -1.0, mask_14_y = -1.0, mask_15_y = -1.0, mask_16_y = -1.0, &
	567	mask_17_y = -1.0, mask_18_y = -1.0, mask_19_y = -1.0, mask_20_y = -1.0, &
	568	mask_01_z = -1.0, mask_02_z = -1.0, mask_03_z = -1.0, mask_04_z = -1.0, &
	569	mask_05_z = -1.0, mask_06_z = -1.0, mask_07_z = -1.0, mask_08_z = -1.0, &
	570	mask_09_z = -1.0, mask_10_z = -1.0, mask_11_z = -1.0, mask_12_z = -1.0, &
	571	mask_13_z = -1.0, mask_14_z = -1.0, mask_15_z = -1.0, mask_16_z = -1.0, &
	572	mask_17_z = -1.0, mask_18_z = -1.0, mask_19_z = -1.0, mask_20_z = -1.0
	573	REAL, DIMENSION(3) :: &
	574	mask_01_x_loop = (/ -1.0, -1.0, 0.0 /), &
	575	mask_02_x_loop = (/ -1.0, -1.0, 0.0 /), &
	576	mask_03_x_loop = (/ -1.0, -1.0, 0.0 /), &
	577	mask_04_x_loop = (/ -1.0, -1.0, 0.0 /), &
	578	mask_05_x_loop = (/ -1.0, -1.0, 0.0 /), &
	579	mask_06_x_loop = (/ -1.0, -1.0, 0.0 /), &
	580	mask_07_x_loop = (/ -1.0, -1.0, 0.0 /), &
	581	mask_08_x_loop = (/ -1.0, -1.0, 0.0 /), &
	582	mask_09_x_loop = (/ -1.0, -1.0, 0.0 /), &
	583	mask_10_x_loop = (/ -1.0, -1.0, 0.0 /), &
	584	mask_11_x_loop = (/ -1.0, -1.0, 0.0 /), &
	585	mask_12_x_loop = (/ -1.0, -1.0, 0.0 /), &
	586	mask_13_x_loop = (/ -1.0, -1.0, 0.0 /), &
	587	mask_14_x_loop = (/ -1.0, -1.0, 0.0 /), &
	588	mask_15_x_loop = (/ -1.0, -1.0, 0.0 /), &
	589	mask_16_x_loop = (/ -1.0, -1.0, 0.0 /), &
	590	mask_17_x_loop = (/ -1.0, -1.0, 0.0 /), &
	591	mask_18_x_loop = (/ -1.0, -1.0, 0.0 /), &
	592	mask_19_x_loop = (/ -1.0, -1.0, 0.0 /), &
	593	mask_20_x_loop = (/ -1.0, -1.0, 0.0 /), &
	594	mask_01_y_loop = (/ -1.0, -1.0, 0.0 /), &
	595	mask_02_y_loop = (/ -1.0, -1.0, 0.0 /), &
	596	mask_03_y_loop = (/ -1.0, -1.0, 0.0 /), &
	597	mask_04_y_loop = (/ -1.0, -1.0, 0.0 /), &
	598	mask_05_y_loop = (/ -1.0, -1.0, 0.0 /), &
	599	mask_06_y_loop = (/ -1.0, -1.0, 0.0 /), &
	600	mask_07_y_loop = (/ -1.0, -1.0, 0.0 /), &
	601	mask_08_y_loop = (/ -1.0, -1.0, 0.0 /), &
	602	mask_09_y_loop = (/ -1.0, -1.0, 0.0 /), &
	603	mask_10_y_loop = (/ -1.0, -1.0, 0.0 /), &
	604	mask_11_y_loop = (/ -1.0, -1.0, 0.0 /), &
	605	mask_12_y_loop = (/ -1.0, -1.0, 0.0 /), &
	606	mask_13_y_loop = (/ -1.0, -1.0, 0.0 /), &
	607	mask_14_y_loop = (/ -1.0, -1.0, 0.0 /), &
	608	mask_15_y_loop = (/ -1.0, -1.0, 0.0 /), &
	609	mask_16_y_loop = (/ -1.0, -1.0, 0.0 /), &
	610	mask_17_y_loop = (/ -1.0, -1.0, 0.0 /), &
	611	mask_18_y_loop = (/ -1.0, -1.0, 0.0 /), &
	612	mask_19_y_loop = (/ -1.0, -1.0, 0.0 /), &
	613	mask_20_y_loop = (/ -1.0, -1.0, 0.0 /), &
	614	mask_01_z_loop = (/ -1.0, -1.0, 0.0 /), &
	615	mask_02_z_loop = (/ -1.0, -1.0, 0.0 /), &
	616	mask_03_z_loop = (/ -1.0, -1.0, 0.0 /), &
	617	mask_04_z_loop = (/ -1.0, -1.0, 0.0 /), &
	618	mask_05_z_loop = (/ -1.0, -1.0, 0.0 /), &
	619	mask_06_z_loop = (/ -1.0, -1.0, 0.0 /), &
	620	mask_07_z_loop = (/ -1.0, -1.0, 0.0 /), &
	621	mask_08_z_loop = (/ -1.0, -1.0, 0.0 /), &
	622	mask_09_z_loop = (/ -1.0, -1.0, 0.0 /), &
	623	mask_10_z_loop = (/ -1.0, -1.0, 0.0 /), &
	624	mask_11_z_loop = (/ -1.0, -1.0, 0.0 /), &
	625	mask_12_z_loop = (/ -1.0, -1.0, 0.0 /), &
	626	mask_13_z_loop = (/ -1.0, -1.0, 0.0 /), &
	627	mask_14_z_loop = (/ -1.0, -1.0, 0.0 /), &
	628	mask_15_z_loop = (/ -1.0, -1.0, 0.0 /), &
	629	mask_16_z_loop = (/ -1.0, -1.0, 0.0 /), &
	630	mask_17_z_loop = (/ -1.0, -1.0, 0.0 /), &
	631	mask_18_z_loop = (/ -1.0, -1.0, 0.0 /), &
	632	mask_19_z_loop = (/ -1.0, -1.0, 0.0 /), &
	633	mask_20_z_loop = (/ -1.0, -1.0, 0.0 /)
	634	!
	635	!-- internal mask arrays ("mask,dimension,selection")
	636	REAL, DIMENSION(:,:,:), ALLOCATABLE :: mask, mask_loop
[240]	637
[1]	638	SAVE
	639
	640	END MODULE control_parameters
	641
	642
	643
	644
	645	MODULE cpulog
	646
	647	!------------------------------------------------------------------------------!
	648	! Description:
	649	! ------------
	650	! Definition of variables for cpu-time measurements
	651	!------------------------------------------------------------------------------!
	652
	653	REAL :: initial_wallclock_time
	654
	655	TYPE logpoint
	656	REAL :: isum, ivect, mean, mtime, mtimevec, sum, vector
	657	INTEGER :: counts
	658	CHARACTER (LEN=20) :: place
	659	END TYPE logpoint
	660
	661	TYPE(logpoint), DIMENSION(100) :: log_point = logpoint( 0.0, 0.0, 0.0, &
	662	0.0, 0.0, 0.0, 0.0, 0, ' ' ), &
	663	log_point_s = logpoint( 0.0, 0.0, 0.0, &
	664	0.0, 0.0, 0.0, 0.0, 0, ' ' )
	665
	666	SAVE
	667
	668	END MODULE cpulog
	669
	670
	671
	672
	673	MODULE dvrp_variables
	674
	675	!------------------------------------------------------------------------------!
	676	! Description:
	677	! ------------
	678	! Definition of variables used with dvrp-software
	679	!------------------------------------------------------------------------------!
	680
[264]	681	CHARACTER (LEN=10) :: dvrp_output = 'rtsp', particle_color = 'none', &
	682	particle_dvrpsize = 'none'
[1]	683
	684	CHARACTER (LEN=20), DIMENSION(10) :: mode_dvrp = &
	685	(/ ( ' ', i9 = 1,10 ) /)
	686
	687	CHARACTER (LEN=80) :: dvrp_directory = 'default', &
	688	dvrp_file = 'default', &
	689	dvrp_host = 'origin.rvs.uni-hannover.de', &
	690	dvrp_password = '********', &
	691	dvrp_username = ' '
	692
[264]	693	INTEGER :: cluster_size = 1, dvrp_colortable_entries = 4, &
	694	dvrp_colortable_entries_prt = 22, islice_dvrp, &
[242]	695	nx_dvrp, nxl_dvrp, nxr_dvrp, ny_dvrp, nyn_dvrp, nys_dvrp, &
	696	nz_dvrp, pathlines_fadeintime = 5, pathlines_fadeouttime = 5, &
	697	pathlines_linecount = 1000, pathlines_maxhistory = 40, &
	698	pathlines_wavecount = 10, pathlines_wavetime = 50, &
	699	vc_gradient_normals = 0, vc_mode = 0, vc_size_x = 2, &
	700	vc_size_y = 2, vc_size_z = 2
[1]	701
	702	INTEGER, DIMENSION(10) :: slicer_position_dvrp
	703
[242]	704	LOGICAL :: cyclic_dvrp = .FALSE., dvrp_overlap, dvrp_total_overlap, &
[260]	705	local_dvrserver_running, lock_steering_update = .FALSE., &
[1]	706	use_seperate_pe_for_dvrp_output = .FALSE.
	707
[242]	708	REAL :: clip_dvrp_l = 9999999.9, clip_dvrp_n = 9999999.9, &
	709	clip_dvrp_r = 9999999.9, clip_dvrp_s = 9999999.9, &
	710	superelevation = 1.0, superelevation_x = 1.0, &
[210]	711	superelevation_y = 1.0, vc_alpha = 38.0
[1]	712
[264]	713	REAL, DIMENSION(2) :: color_interval = (/ 0.0, 1.0 /), &
	714	dvrpsize_interval = (/ 0.0, 1.0 /)
	715
[284]	716	REAL, DIMENSION(3) :: groundplate_color = (/ 0.0, 0.6, 0.0 /), &
	717	topography_color = (/ 0.8, 0.7, 0.6 /)
	718
[412]	719	#if defined( __decalpha )
[405]	720	REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp = RESHAPE( (/ &
[1]	721	-1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
	722	-1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
[405]	723	-1.0, 1.0, -1.0, 1.0 /), (/ 2, 10 /) )
[1]	724
[405]	725	REAL, DIMENSION(3,10) :: isosurface_color = RESHAPE( (/ &
[284]	726	0.9, 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, &
	727	0.1, 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, &
[405]	728	0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, &
	729	0.1, 0.4, 0.1 /), (/ 3, 10 /) )
[284]	730
[264]	731	REAL(4), DIMENSION(2,100) :: interval_values_dvrp, interval_h_dvrp = &
[405]	732	RESHAPE( (/ 270.0, 225.0, 225.0, 180.0, &
	733	70.0, 25.0, 25.0, -25.0, &
	734	( 0.0, i9 = 1, 192 ) /), &
	735	(/ 2, 100 /) ), &
[210]	736	interval_l_dvrp = 0.5, interval_s_dvrp = 1.0,&
[264]	737	interval_a_dvrp = 0.0, &
	738	interval_values_dvrp_prt, &
[405]	739	interval_h_dvrp_prt = RESHAPE( &
[264]	740	(/ 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, &
	741	25.0, -25.0, ( 0.0, i9 = 1, 192 ) /), &
[405]	742	(/ 2, 100 /) ), &
[264]	743	interval_l_dvrp_prt = 0.5, &
	744	interval_s_dvrp_prt = 1.0, &
	745	interval_a_dvrp_prt = 0.0
[412]	746	#else
	747	REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp
[1]	748
[412]	749	REAL, DIMENSION(3,10) :: isosurface_color
[284]	750
[412]	751	REAL(4), DIMENSION(2,100) :: interval_values_dvrp, &
	752	interval_values_dvrp_prt, interval_h_dvrp, &
	753	interval_h_dvrp_prt, interval_l_dvrp = 0.5, &
	754	interval_l_dvrp_prt = 0.5, interval_s_dvrp = 1.0, &
	755	interval_s_dvrp_prt = 1.0, interval_a_dvrp = 0.0, &
	756	interval_a_dvrp_prt = 0.0
[1]	757
[412]	758	DATA slicer_range_limits_dvrp / -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
	759	-1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
	760	-1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
	761	-1.0, 1.0 /
[1]	762
[412]	763	DATA isosurface_color / 0.9, 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, &
	764	0.1, 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, &
	765	0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, &
	766	0.1, 0.4, 0.1 /
[284]	767
[412]	768	DATA interval_h_dvrp / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, &
	769	25.0, -25.0, 192 * 0.0 /
[271]	770
[412]	771	DATA interval_h_dvrp_prt / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, &
	772	25.0, -25.0, 192 * 0.0 /
	773	#endif
[1]	774
[210]	775	REAL(4), DIMENSION(:), ALLOCATABLE :: xcoor_dvrp, ycoor_dvrp, zcoor_dvrp
[1]	776
	777	TYPE steering
	778	CHARACTER (LEN=20) :: name
[210]	779	REAL(4) :: min, max
[1]	780	INTEGER :: imin, imax
	781	END TYPE steering
	782
	783	TYPE(steering), DIMENSION(:), ALLOCATABLE :: steering_dvrp
	784
	785	SAVE
	786
	787	END MODULE dvrp_variables
	788
	789
	790
	791
	792	MODULE grid_variables
	793
	794	!------------------------------------------------------------------------------!
	795	! Description:
	796	! ------------
	797	! Definition of grid spacings
	798	!------------------------------------------------------------------------------!
	799
	800	REAL :: ddx, ddx2, dx = 1.0, dx2, ddy, ddy2, dy = 1.0, dy2
	801
	802	REAL, DIMENSION(:), ALLOCATABLE :: ddx2_mg, ddy2_mg
	803
[48]	804	REAL, DIMENSION(:,:), ALLOCATABLE :: fwxm, fwxp, fwym, fwyp, fxm, fxp, &
	805	fym, fyp, wall_e_x, wall_e_y, &
	806	wall_u, wall_v, wall_w_x, wall_w_y, &
	807	zu_s_inner, zw_w_inner
[1]	808
	809	SAVE
	810
	811	END MODULE grid_variables
	812
	813
	814
	815
	816	MODULE indices
	817
	818	!------------------------------------------------------------------------------!
	819	! Description:
	820	! ------------
	821	! Definition of array bounds and number of gridpoints
	822	!------------------------------------------------------------------------------!
	823
[147]	824	INTEGER :: ngp_sums, nnx, nx = 0, nxa, nxl, nxlu, nxr, nxra, nx_on_file, &
	825	nny, ny = 0, nya, nyn, nyna, nys, nysv, ny_on_file, nnz, &
	826	nz = 0, nza, nzb, nzb_diff, nzt, nzta, nzt_diff
[1]	827
	828	INTEGER, DIMENSION(:), ALLOCATABLE :: &
	829	ngp_2dh, ngp_3d, ngp_3d_inner, &
	830	nnx_pe, nny_pe, nxl_mg, nxr_mg, nyn_mg, nys_mg, nzt_mg
	831
	832	INTEGER, DIMENSION(:,:), ALLOCATABLE :: &
[132]	833	ngp_2dh_outer, ngp_2dh_s_inner, mg_loc_ind, nzb_diff_s_inner, &
	834	nzb_diff_s_outer, nzb_diff_u, nzb_diff_v, nzb_inner, nzb_outer,&
	835	nzb_s_inner, nzb_s_outer, nzb_u_inner, nzb_u_outer, &
	836	nzb_v_inner, nzb_v_outer, nzb_w_inner, nzb_w_outer, nzb_2d
[1]	837
[114]	838	INTEGER, DIMENSION(:,:,:), POINTER :: flags
	839
	840	INTEGER, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_1, &
	841	wall_flags_2, wall_flags_3, wall_flags_4, wall_flags_5, &
	842	wall_flags_6, wall_flags_7, wall_flags_8, wall_flags_9, &
	843	wall_flags_10
	844
[1]	845	SAVE
	846
	847	END MODULE indices
	848
	849
	850
	851
	852	MODULE interfaces
	853
	854	!------------------------------------------------------------------------------!
	855	! Description:
	856	! ------------
	857	! Interfaces for special subroutines which use optional parameters
	858	!------------------------------------------------------------------------------!
	859
	860	INTERFACE
	861
	862	SUBROUTINE cpu_log( log_event, place, modus, barrierwait )
	863
	864	USE cpulog
	865
	866	CHARACTER (LEN=*) :: modus, place
	867	CHARACTER (LEN=*), OPTIONAL :: barrierwait
	868	TYPE(logpoint) :: log_event
	869
	870	END SUBROUTINE cpu_log
	871
	872	END INTERFACE
	873
	874
	875
	876	INTERFACE
	877
	878	SUBROUTINE global_min_max ( i1, i2, j1, j2, k1, k2, feld, mode, wert, &
	879	wert_ijk, wert1, wert1_ijk )
	880	CHARACTER (LEN=*), INTENT(IN) :: mode
	881	INTEGER, INTENT(IN) :: i1, i2, j1, j2, k1, k2
	882	INTEGER :: wert_ijk(3)
	883	INTEGER, OPTIONAL :: wert1_ijk(3)
	884	REAL :: wert
	885	REAL, OPTIONAL :: wert1
	886	REAL, INTENT(IN) :: feld(i1:i2,j1:j2,k1:k2)
	887
	888	END SUBROUTINE global_min_max
	889
	890	END INTERFACE
	891
	892	SAVE
	893
	894	END MODULE interfaces
	895
	896
	897
	898	MODULE pointer_interfaces
	899
	900	!------------------------------------------------------------------------------!
	901	! Description:
	902	! ------------
	903	! Interfaces for subroutines with pointer arguments called in
	904	! prognostic_equations
	905	!------------------------------------------------------------------------------!
	906
	907	INTERFACE
	908
	909	SUBROUTINE advec_s_bc( sk, sk_char )
	910
	911	CHARACTER (LEN=*), INTENT(IN) :: sk_char
	912	REAL, DIMENSION(:,:,:), POINTER :: sk
	913
	914	END SUBROUTINE advec_s_bc
	915
	916	END INTERFACE
	917
	918
	919	SAVE
	920
	921	END MODULE pointer_interfaces
	922
	923
	924
	925
	926	MODULE model_1d
	927
	928	!------------------------------------------------------------------------------!
	929	! Description:
	930	! ------------
	931	! Definition of variables for the 1D-model
	932	!------------------------------------------------------------------------------!
	933
	934	INTEGER :: current_timestep_number_1d = 0, damp_level_ind_1d, &
	935	last_dt_change_1d = 0
	936
	937	LOGICAL :: run_control_header_1d = .FALSE., stop_dt_1d = .FALSE.
	938
	939	REAL :: damp_level_1d = -1.0, dt_1d = 60.0, dt_max_1d = 300.0, &
	940	dt_pr_1d = 9999999.9, dt_run_control_1d = 60.0, &
	941	end_time_1d = 864000.0, old_dt_1d = 1.0E-10, &
	942	qs1d, simulated_time_1d = 0.0, time_pr_1d = 0.0, &
	943	time_run_control_1d = 0.0, ts1d, us1d, usws1d, usws1d_m, &
	944	vsws1d, vsws1d_m, z01d
	945
	946
	947	REAL, DIMENSION(:), ALLOCATABLE :: e1d, e1d_m, e1d_p, kh1d, kh1d_m, km1d, &
	948	km1d_m, l_black, l1d, l1d_m, rif1d, &
	949	te_e, te_em, te_u, te_um, te_v, te_vm, &
	950	u1d, u1d_m, u1d_p, v1d, v1d_m, v1d_p
	951
	952	SAVE
	953
	954	END MODULE model_1d
	955
	956
	957
	958
	959	MODULE netcdf_control
	960
	961	!------------------------------------------------------------------------------!
	962	! Description:
	963	! ------------
	964	! Definition of parameters and variables for netcdf control.
	965	!------------------------------------------------------------------------------!
	966
[410]	967	USE control_parameters
[1]	968	#if defined( __netcdf )
	969	USE netcdf
	970	#endif
	971
[48]	972	INTEGER, PARAMETER :: dopr_norm_num = 7, dopts_num = 26, dots_max = 100, &
[1]	973	replace_num = 6
	974
[343]	975	INTEGER :: dots_num = 23
[48]	976
[1]	977	CHARACTER, DIMENSION( replace_num ) :: &
	978	replace_char = (/ '''', '"', '*', '/', '(', ')' /), &
	979	replace_by = (/ 'p' , 'p', 's', 'o', '_', '_' /)
	980
	981	CHARACTER (LEN=6), DIMENSION(dopr_norm_num) :: dopr_norm_names = &
	982	(/ 'wpt0 ', 'ws2 ', 'tsw2 ', 'ws3 ', 'ws2tsw', 'wstsw2', &
	983	'z_i ' /)
	984
	985	CHARACTER (LEN=6), DIMENSION(dopr_norm_num) :: dopr_norm_longnames = &
	986	(/ 'wpt0 ', 'w2 ', 'tw2 ', 'w3 ', 'w2tw', 'wt*w2', &
	987	'z_i ' /)
	988
	989	CHARACTER (LEN=7), DIMENSION(dopts_num) :: dopts_label = &
	990	(/ 'tnpt ', 'x_ ', 'y_ ', 'z_ ', 'z_abs ', 'u ', &
	991	'v ', 'w ', 'u" ', 'v" ', 'w" ', 'npt_up ', &
	992	'w_up ', 'w_down ', 'npt_max', 'npt_min', 'x2 ', 'y2 ', &
	993	'z2 ', 'u2 ', 'v2 ', 'w2 ', 'u"2 ', 'v"2 ', &
	994	'w"2 ', 'npt*2 ' /)
	995
	996	CHARACTER (LEN=7), DIMENSION(dopts_num) :: dopts_unit = &
	997	(/ 'number ', 'm ', 'm ', 'm ', 'm ', 'm/s ', &
	998	'm/s ', 'm/s ', 'm/s ', 'm/s ', 'm/s ', 'number ', &
	999	'm/s ', 'm/s ', 'number ', 'number ', 'm2 ', 'm2 ', &
	1000	'm2 ', 'm2/s2 ', 'm2/s2 ', 'm2/s2 ', 'm2/s2 ', 'm2/s2 ', &
	1001	'm2/s2 ', 'number2' /)
	1002
[48]	1003	CHARACTER (LEN=7), DIMENSION(dots_max) :: dots_label = &
[1]	1004	(/ 'E ', 'E* ', 'dt ', 'u* ', 'th* ', 'umax ', &
	1005	'vmax ', 'wmax ', 'div_new', 'div_old', 'z_i_wpt', 'z_i_pt ', &
	1006	'w* ', 'w"pt"0 ', 'w"pt" ', 'wpt ', 'pt(0) ', 'pt(zp) ', &
[343]	1007	'w"u"0 ', 'w"v"0 ', 'w"q"0 ', 'mo_L ', 'q* ', &
[394]	1008	( 'unknown', i9 = 1, dots_max-23 ) /)
[1]	1009
[48]	1010	CHARACTER (LEN=7), DIMENSION(dots_max) :: dots_unit = &
[1]	1011	(/ 'm2/s2 ', 'm2/s2 ', 's ', 'm/s ', 'K ', 'm/s ', &
	1012	'm/s ', 'm/s ', 's-1 ', 's-1 ', 'm ', 'm ', &
[288]	1013	'm/s ', 'K m/s ', 'K m/s ', 'K m/s ', 'K ', 'K ', &
[343]	1014	'm2/s2 ', 'm2/s2 ', 'kg m/s ', 'm ', 'kg/kg ', &
[394]	1015	( 'unknown', i9 = 1, dots_max-23 ) /)
[1]	1016
[87]	1017	CHARACTER (LEN=7), DIMENSION(300) :: dopr_unit = 'unknown'
[1]	1018
	1019	CHARACTER (LEN=7), DIMENSION(0:1,100) :: do2d_unit, do3d_unit
	1020
	1021	CHARACTER (LEN=16), DIMENSION(25) :: prt_var_names = &
	1022	(/ 'pt_age ', 'pt_dvrp_size ', 'pt_origin_x ', &
	1023	'pt_origin_y ', 'pt_origin_z ', 'pt_radius ', &
	1024	'pt_speed_x ', 'pt_speed_y ', 'pt_speed_z ', &
	1025	'pt_weight_factor', 'pt_x ', 'pt_y ', &
	1026	'pt_z ', 'pt_color ', 'pt_group ', &
	1027	'pt_tailpoints ', 'pt_tail_id ', 'pt_density_ratio', &
	1028	'pt_exp_arg ', 'pt_exp_term ', 'not_used ', &
	1029	'not_used ', 'not_used ', 'not_used ', &
	1030	'not_used ' /)
	1031
	1032	CHARACTER (LEN=16), DIMENSION(25) :: prt_var_units = &
	1033	(/ 'seconds ', 'meters ', 'meters ', &
	1034	'meters ', 'meters ', 'meters ', &
	1035	'm/s ', 'm/s ', 'm/s ', &
	1036	'factor ', 'meters ', 'meters ', &
	1037	'meters ', 'none ', 'none ', &
	1038	'none ', 'none ', 'ratio ', &
	1039	'none ', 'none ', 'not_used ', &
	1040	'not_used ', 'not_used ', 'not_used ', &
	1041	'not_used ' /)
	1042
	1043	INTEGER :: id_dim_prtnum, id_dim_time_pr, id_dim_time_prt, &
	1044	id_dim_time_pts, id_dim_time_sp, id_dim_time_ts, id_dim_x_sp, &
	1045	id_dim_y_sp, id_dim_zu_sp, id_dim_zw_sp, id_set_pr, &
	1046	id_set_prt, id_set_pts, id_set_sp, id_set_ts, id_var_prtnum, &
	1047	id_var_rnop_prt, id_var_time_pr, id_var_time_prt, &
	1048	id_var_time_pts, id_var_time_sp, id_var_time_ts, id_var_x_sp, &
	1049	id_var_y_sp, id_var_zu_sp, id_var_zw_sp, nc_stat
	1050
	1051	INTEGER, DIMENSION(0:1) :: id_dim_time_xy, id_dim_time_xz, &
	1052	id_dim_time_yz, id_dim_time_3d, id_dim_x_xy, id_dim_xu_xy, &
	1053	id_dim_x_xz, id_dim_xu_xz, id_dim_x_yz, id_dim_xu_yz, &
	1054	id_dim_x_3d, id_dim_xu_3d, id_dim_y_xy, id_dim_yv_xy, &
	1055	id_dim_y_xz, id_dim_yv_xz, id_dim_y_yz, id_dim_yv_yz, &
	1056	id_dim_y_3d, id_dim_yv_3d, id_dim_zu_xy, id_dim_zu1_xy, &
	1057	id_dim_zu_xz, id_dim_zu_yz, id_dim_zu_3d, id_dim_zw_xy, &
	1058	id_dim_zw_xz, id_dim_zw_yz, id_dim_zw_3d, id_set_xy, &
	1059	id_set_xz, id_set_yz, id_set_3d, id_var_ind_x_yz, &
	1060	id_var_ind_y_xz, id_var_ind_z_xy, id_var_time_xy, &
	1061	id_var_time_xz, id_var_time_yz, id_var_time_3d, id_var_x_xy, &
	1062	id_var_xu_xy, id_var_x_xz, id_var_xu_xz, id_var_x_yz, &
	1063	id_var_xu_yz, id_var_x_3d, id_var_xu_3d, id_var_y_xy, &
	1064	id_var_yv_xy, id_var_y_xz, id_var_yv_xz, id_var_y_yz, &
[48]	1065	id_var_yv_yz, id_var_y_3d, id_var_yv_3d, id_var_zusi_xy, &
	1066	id_var_zusi_3d, id_var_zu_xy, id_var_zu1_xy, id_var_zu_xz, &
	1067	id_var_zu_yz, id_var_zu_3d, id_var_zwwi_xy, id_var_zwwi_3d, &
[1]	1068	id_var_zw_xy, id_var_zw_xz, id_var_zw_yz, id_var_zw_3d
	1069
[410]	1070	INTEGER, DIMENSION(10) :: id_var_dospx, id_var_dospy
[1]	1071	INTEGER, DIMENSION(20) :: id_var_prt
[410]	1072	INTEGER, DIMENSION(50) :: nc_precision
[1]	1073	INTEGER, DIMENSION(dopr_norm_num) :: id_var_norm_dopr
	1074
	1075	INTEGER, DIMENSION(dopts_num,0:10) :: id_var_dopts
	1076	INTEGER, DIMENSION(0:1,100) :: id_var_do2d, id_var_do3d
	1077	INTEGER, DIMENSION(100,0:9) :: id_dim_z_pr, id_var_dopr, &
	1078	id_var_z_pr
[48]	1079	INTEGER, DIMENSION(dots_max,0:9) :: id_var_dots
[1]	1080
[410]	1081	!
	1082	!-- masked output
	1083	CHARACTER (LEN=7), DIMENSION(max_masks,0:1,100) :: domask_unit
	1084
[345]	1085	LOGICAL :: output_for_t0 = .FALSE.
[1]	1086
[410]	1087	INTEGER, DIMENSION(1:max_masks,0:1) :: id_dim_time_mask, id_dim_x_mask, &
	1088	id_dim_xu_mask, id_dim_y_mask, id_dim_yv_mask, id_dim_zu_mask, &
	1089	id_dim_zw_mask, &
	1090	id_set_mask, &
	1091	id_var_time_mask, id_var_x_mask, id_var_xu_mask, &
	1092	id_var_y_mask, id_var_yv_mask, id_var_zu_mask, id_var_zw_mask, &
	1093	id_var_zusi_mask, id_var_zwwi_mask
	1094
	1095	INTEGER, DIMENSION(1:max_masks,0:1,100) :: id_var_domask
	1096
[1]	1097	SAVE
	1098
	1099	END MODULE netcdf_control
	1100
	1101
	1102
	1103	MODULE particle_attributes
	1104
	1105	!------------------------------------------------------------------------------!
	1106	! Description:
	1107	! ------------
	1108	! Definition of variables used to compute particle transport
	1109	!------------------------------------------------------------------------------!
	1110
	1111	USE array_kind
	1112
	1113	CHARACTER (LEN=15) :: bc_par_lr = 'cyclic', bc_par_ns = 'cyclic', &
	1114	bc_par_b = 'reflect', bc_par_t = 'absorb'
	1115
	1116	#if defined( __parallel )
	1117	INTEGER :: mpi_particle_type
	1118	#endif
	1119	INTEGER :: ibc_par_lr, ibc_par_ns, ibc_par_b, ibc_par_t, &
	1120	iran_part = -1234567, maximum_number_of_particles = 1000, &
	1121	maximum_number_of_tailpoints = 100, &
	1122	maximum_number_of_tails = 0, &
	1123	number_of_initial_particles = 0, number_of_particles = 0, &
	1124	number_of_particle_groups = 1, number_of_tails = 0, &
[150]	1125	number_of_initial_tails = 0, offset_ocean_nzt = 0, &
	1126	offset_ocean_nzt_m1 = 0, particles_per_point = 1, &
[1]	1127	particle_file_count = 0, skip_particles_for_tail = 100, &
	1128	total_number_of_particles, total_number_of_tails = 0
	1129
	1130	INTEGER, PARAMETER :: max_number_of_particle_groups = 10
	1131
	1132	INTEGER, DIMENSION(:), ALLOCATABLE :: new_tail_id
	1133	INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: prt_count, prt_start_index
	1134
	1135	LOGICAL :: particle_advection = .FALSE., random_start_position = .FALSE., &
	1136	read_particles_from_restartfile = .TRUE., &
	1137	uniform_particles = .TRUE., use_particle_tails = .FALSE., &
	1138	use_sgs_for_particles = .FALSE., &
	1139	write_particle_statistics = .FALSE.
[336]	1140
	1141	LOGICAL, DIMENSION(max_number_of_particle_groups) :: &
	1142	vertical_particle_advection = .TRUE.
	1143
[1]	1144	LOGICAL, DIMENSION(:), ALLOCATABLE :: particle_mask, tail_mask
	1145
[116]	1146	REAL :: c_0 = 3.0, dt_min_part = 0.0002, dt_sort_particles = 0.0, &
[1]	1147	dt_write_particle_data = 9999999.9, dvrp_psize = 9999999.9, &
	1148	end_time_prel = 9999999.9, initial_weighting_factor = 1.0, &
	1149	maximum_tailpoint_age = 100000.0, &
	1150	minimum_tailpoint_distance = 0.0, &
	1151	particle_advection_start = 0.0, sgs_wfu_part = 0.3333333, &
	1152	sgs_wfv_part = 0.3333333, sgs_wfw_part = 0.3333333, &
[116]	1153	time_sort_particles = 0.0, time_write_particle_data = 0.0
[1]	1154
	1155	REAL, DIMENSION(max_number_of_particle_groups) :: &
	1156	density_ratio = 9999999.9, pdx = 9999999.9, pdy = 9999999.9, &
	1157	pdz = 9999999.9, psb = 9999999.9, psl = 9999999.9, &
	1158	psn = 9999999.9, psr = 9999999.9, pss = 9999999.9, &
	1159	pst = 9999999.9, radius = 9999999.9
	1160
	1161	REAL, DIMENSION(:,:,:), ALLOCATABLE :: particle_tail_coordinates
	1162
	1163
	1164	TYPE particle_type
	1165	SEQUENCE
[57]	1166	REAL :: age, age_m, dt_sum, dvrp_psize, e_m, origin_x, origin_y, &
	1167	origin_z, radius, speed_x, speed_x_sgs, speed_y, &
	1168	speed_y_sgs, speed_z, speed_z_sgs, weight_factor, x, y, z
[1]	1169	INTEGER :: color, group, tailpoints, tail_id
	1170	END TYPE particle_type
	1171
	1172	TYPE(particle_type), DIMENSION(:), ALLOCATABLE :: initial_particles, &
	1173	particles
	1174
	1175	TYPE particle_groups_type
	1176	SEQUENCE
	1177	REAL :: density_ratio, radius, exp_arg, exp_term
	1178	END TYPE particle_groups_type
	1179
	1180	TYPE(particle_groups_type), DIMENSION(max_number_of_particle_groups) ::&
	1181	particle_groups
	1182
	1183	SAVE
	1184
	1185	END MODULE particle_attributes
	1186
	1187
	1188
	1189
	1190
	1191	MODULE pegrid
	1192
	1193	!------------------------------------------------------------------------------!
	1194	! Description:
	1195	! ------------
	1196	! Definition of variables which define processor topology and the exchange of
	1197	! ghost point layers. This modules must be placed in all routines which contain
	1198	! MPI-calls.
	1199	!------------------------------------------------------------------------------!
	1200
	1201	#if defined( __parallel )
[82]	1202	#if defined( __lc )
[1]	1203	USE MPI
	1204	#else
	1205	INCLUDE "mpif.h"
	1206	#endif
	1207	#endif
[146]	1208	CHARACTER(LEN=5) :: myid_char = ''
[206]	1209	INTEGER :: id_inflow = 0, id_recycling = 0, myid = 0, &
	1210	target_id, npex = -1, npey = -1, numprocs = 1, &
	1211	numprocs_previous_run = -1, &
[1]	1212	tasks_per_node = -9999, threads_per_task = 1
	1213
[145]	1214	INTEGER, DIMENSION(:,:), ALLOCATABLE :: hor_index_bounds, &
	1215	hor_index_bounds_previous_run
	1216
[1]	1217	#if defined( __parallel )
[102]	1218	#if defined( __mpi2 )
	1219	CHARACTER (LEN=MPI_MAX_PORT_NAME) :: port_name
	1220	#endif
	1221
	1222	INTEGER :: comm1dx, comm1dy, comm2d, comm_inter, comm_palm, ierr, myidx, &
	1223	myidy, ndim = 2, ngp_xy, ngp_y, pleft, pnorth, pright, psouth, &
[1]	1224	sendrecvcount_xy, sendrecvcount_yz, sendrecvcount_zx, &
	1225	sendrecvcount_zyd, sendrecvcount_yxd, &
[102]	1226	type_x, type_x_int, type_xy
[1]	1227
[102]	1228	INTEGER :: ibuf(12), pcoord(2), pdims(2), status(MPI_STATUS_SIZE)
	1229
[1]	1230	INTEGER, DIMENSION(:), ALLOCATABLE :: ngp_yz, type_xz
	1231
	1232	LOGICAL :: reorder = .TRUE.
	1233	LOGICAL, DIMENSION(2) :: cyclic = (/ .TRUE. , .TRUE. /), &
	1234	remain_dims
	1235	#endif
	1236
	1237	SAVE
	1238
	1239	END MODULE pegrid
	1240
	1241
	1242
	1243
	1244	MODULE profil_parameter
	1245
	1246	!------------------------------------------------------------------------------!
	1247	! Description:
	1248	! ------------
	1249	! Definition of variables which control PROFIL-output
	1250	!------------------------------------------------------------------------------!
	1251
	1252	INTEGER, PARAMETER :: crmax = 100
	1253
	1254	CHARACTER (LEN=10), DIMENSION(100) :: dopr_label = ' '
	1255
	1256	CHARACTER (LEN=10), DIMENSION(crmax) :: cross_normalized_x = ' ', &
	1257	cross_normalized_y = ' '
	1258
	1259	CHARACTER (LEN=20), DIMENSION(20) :: cross_ts_profiles = &
	1260	(/ ' E E* ', ' dt ', &
	1261	' u* w* ', ' th* ', &
	1262	' umax vmax wmax ', ' div_old div_new ', &
	1263	' z_i_wpt z_i_pt ', ' w"pt"0 w"pt" wpt ', &
	1264	' pt(0) pt(zp) ', ' splux spluy spluz ', &
	1265	' L ', &
	1266	( ' ', i9 = 1, 9 ) /)
	1267
	1268	CHARACTER (LEN=40), DIMENSION(crmax) :: cross_xtext = &
	1269	(/ 'windspeed in ms>->1 ', &
	1270	'pot. temperature in K ', &
	1271	'heat flux in K ms>->1 ', &
	1272	'momentum flux in m>2s>2 ', &
	1273	'eddy diffusivity in m>2s>->1', &
	1274	'mixing length in m ', &
	1275	( ' ', i9 = 1, 94 ) /)
	1276
	1277	CHARACTER (LEN=100), DIMENSION(crmax) :: cross_profiles = &
	1278	(/ ' u v ', &
	1279	' pt ', &
	1280	' w"pt" wpt wptBC wpt wptBC ', &
	1281	' w"u" wu wu w"v" wv wv ', &
	1282	' km kh ', &
	1283	' l ', &
	1284	( ' ', i9 = 1, 94 ) /)
	1285
	1286	INTEGER :: profile_columns = 3, profile_rows = 2, profile_number = 0
	1287
	1288	INTEGER :: cross_linecolors(100,crmax) = 1, &
	1289	cross_linestyles(100,crmax) = 0, &
	1290	cross_profile_numbers(100,crmax) = 0, &
	1291	cross_pnc_local(crmax), cross_profile_number_count(crmax) = 0, &
	1292	cross_ts_numbers(crmax,crmax) = 0, &
	1293	cross_ts_number_count(crmax) = 0, dopr_crossindex(100) = 0, &
[87]	1294	dopr_index(300) = 0, dopr_initial_index(300) = 0, &
[1]	1295	dots_crossindex(100) = 0, dots_index(100) = 0, &
	1296	linecolors(10) = (/ 2, 3, 4, 5, 7, 8, 12, 15, 16, 23 /), &
	1297	linestyles(11) = (/ 0, 7, 3, 10, 4, 1, 9, 2, 5, 8, 6 /)
	1298
	1299
	1300	REAL :: cross_normx_factor(100,crmax) = 1.0, &
	1301	cross_normy_factor(100,crmax) = 1.0, &
	1302	cross_ts_uymax(20) = &
	1303	(/ 999.999, 999.999, 999.999, 999.999, 999.999, &
	1304	999.999, 999.999, 999.999, 999.999, 999.999, &
	1305	999.999, 999.999, 999.999, 999.999, 999.999, &
	1306	999.999, 999.999, 999.999, 999.999, 999.999 /),&
	1307	cross_ts_uymax_computed(20) = 999.999, &
	1308	cross_ts_uymin(20) = &
	1309	(/ 999.999, 999.999, 999.999, -5.000, 999.999, &
	1310	999.999, 0.000, 999.999, 999.999, 999.999, &
	1311	999.999, 999.999, 999.999, 999.999, 999.999, &
	1312	999.999, 999.999, 999.999, 999.999, 999.999 /),&
	1313	cross_ts_uymin_computed(20) = 999.999, &
	1314	cross_uxmax(crmax) = 0.0, cross_uxmax_computed(crmax) = -1.0, &
	1315	cross_uxmax_normalized(crmax) = 0.0, &
	1316	cross_uxmax_normalized_computed(crmax) = -1.0, &
	1317	cross_uxmin(crmax) = 0.0, cross_uxmin_computed(crmax) = 1.0, &
	1318	cross_uxmin_normalized(crmax) = 0.0, &
	1319	cross_uxmin_normalized_computed(crmax) = 1.0, &
	1320	cross_uymax(crmax), cross_uymin(crmax)
	1321
	1322	SAVE
	1323
	1324	END MODULE profil_parameter
	1325
	1326
	1327
	1328
	1329	MODULE spectrum
	1330
	1331	!------------------------------------------------------------------------------!
	1332	! Description:
	1333	! ------------
	1334	! Definition of quantities used for computing spectra
	1335	!------------------------------------------------------------------------------!
	1336
	1337	CHARACTER (LEN=6), DIMENSION(1:5) :: header_char = (/ 'PS(u) ', 'PS(v) ',&
	1338	'PS(w) ', 'PS(pt)', 'PS(q) ' /)
	1339	CHARACTER (LEN=2), DIMENSION(10) :: spectra_direction = 'x'
	1340	CHARACTER (LEN=10), DIMENSION(10) :: data_output_sp = ' '
	1341	CHARACTER (LEN=25), DIMENSION(1:5) :: utext_char = &
	1342	(/ '-power spectrum of u ', &
	1343	'-power spectrum of v ', &
	1344	'-power spectrum of w ', &
	1345	'-power spectrum of ^1185 ', &
	1346	'-power spectrum of q ' /)
	1347	CHARACTER (LEN=39), DIMENSION(1:5) :: ytext_char = &
	1348	(/ 'k ^2236 ^2566^2569<u(k) in m>2s>->2 ', &
	1349	'k ^2236 ^2566^2569<v(k) in m>2s>->2 ', &
	1350	'k ^2236 ^2566^2569<w(k) in m>2s>->2 ', &
	1351	'k ^2236 ^2566^2569<^1185(k) in m>2s>->2', &
	1352	'k ^2236 ^2566^2569<q(k) in m>2s>->2 ' /)
	1353
	1354	INTEGER :: klist_x = 0, klist_y = 0, n_sp_x = 0, n_sp_y = 0
	1355
[189]	1356	INTEGER :: comp_spectra_level(100) = 999999, &
	1357	lstyles(100) = (/ 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1358	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1359	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1360	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1361	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1362	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1363	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1364	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1365	0, 7, 3, 10, 1, 4, 9, 2, 6, 8, &
	1366	0, 7, 3, 10, 1, 4, 9, 2, 6, 8 /), &
	1367	plot_spectra_level(100) = 999999
[1]	1368
	1369	REAL :: time_to_start_sp = 0.0
	1370
	1371	SAVE
	1372
	1373	END MODULE spectrum
	1374
	1375
	1376
	1377
	1378	MODULE statistics
	1379
	1380	!------------------------------------------------------------------------------!
	1381	! Description:
	1382	! ------------
	1383	! Definition of statistical quantities, e.g. global sums
	1384	!------------------------------------------------------------------------------!
	1385
	1386	CHARACTER (LEN=40) :: region(0:9)
[394]	1387	INTEGER :: pr_palm = 80, statistic_regions = 0
[318]	1388	INTEGER :: u_max_ijk(3), v_max_ijk(3), w_max_ijk(3)
[1]	1389	LOGICAL :: flow_statistics_called = .FALSE.
[318]	1390	REAL :: u_max, v_max, w_max
[1]	1391	REAL, DIMENSION(:), ALLOCATABLE :: sums_divnew_l, sums_divold_l
[48]	1392	REAL, DIMENSION(:,:), ALLOCATABLE :: sums, sums_wsts_bc_l, ts_value
[1]	1393	REAL, DIMENSION(:,:,:), ALLOCATABLE :: hom_sum, rmask, spectrum_x, &
	1394	spectrum_y, sums_l, sums_l_l, &
	1395	sums_up_fraction_l
	1396	REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: hom
	1397
	1398	SAVE
	1399
	1400	END MODULE statistics
	1401
	1402
	1403
	1404
	1405	MODULE transpose_indices
	1406
	1407	!------------------------------------------------------------------------------!
	1408	! Description:
	1409	! ------------
	1410	! Definition of indices for transposed arrays
	1411	!------------------------------------------------------------------------------!
	1412
	1413	INTEGER :: nxl_y, nxl_yd, nxl_z, nxr_y, nxr_ya, nxr_yd, nxr_yda, nxr_z, &
	1414	nxr_za, nyn_x, nyn_xa, nyn_z, nyn_za, nys_x, nys_z, nzb_x, &
	1415	nzb_y, nzb_yd, nzt_x, nzt_xa, nzt_y, nzt_ya, nzt_yd, nzt_yda
	1416
	1417
	1418	SAVE
	1419
	1420	END MODULE transpose_indices

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |