Home

Context Navigation

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

Last change on this file since 289 was 288, checked in by heinze, 16 years ago
NetCDF unit attribute in timeseries output added, typographical error in modules fixed
Property svn:keywords set to `Id`
File size: 54.4 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |