Home

Context Navigation

source: palm/trunk/SOURCE/parin.f90 @ 1524

Last change on this file since 1524 was 1497, checked in by maronga, 10 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 24.6 KB

Line
1	SUBROUTINE parin
2
3	!--------------------------------------------------------------------------------!
4	! This file is part of PALM.
5	!
6	! PALM is free software: you can redistribute it and/or modify it under the terms
7	! of the GNU General Public License as published by the Free Software Foundation,
8	! either version 3 of the License, or (at your option) any later version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2014 Leibniz Universitaet Hannover
18	!--------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! -----------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: parin.f90 1497 2014-12-02 17:28:07Z keck $
27	!
28	! 1496 2014-12-02 17:25:50Z maronga
29	! Renamed: "radiation -> "cloud_top_radiation"
30	!
31	! 1484 2014-10-21 10:53:05Z kanani
32	! Changes due to new module structure of the plant canopy model:
33	! canopy-model related parameters moved to new package canopy_par in
34	! subroutine package_parin
35	!
36	! 1429 2014-07-15 12:53:45Z knoop
37	! +ensemble_member_nr to prepare the random_generator for ensemble runs
38	!
39	! 1402 2014-05-09 14:25:13Z raasch
40	! location messages modified, batch_job included in envpar-NAMELIST
41	!
42	! 1384 2014-05-02 14:31:06Z raasch
43	! location messages added
44	!
45	! 1365 2014-04-22 15:03:56Z boeske
46	! Usage of large scale forcing enabled:
47	! +use_subsidence_tendencies
48	!
49	! 1361 2014-04-16 15:17:48Z hoffmann
50	! +call_microphysics_at_all_substeps
51	!
52	! 1359 2014-04-11 17:15:14Z hoffmann
53	! REAL constants provided with KIND-attribute
54	!
55	! 1353 2014-04-08 15:21:23Z heinze
56	! REAL constants provided with KIND-attribute
57	!
58	! 1327 2014-03-21 11:00:16Z raasch
59	! -data_output_format, do3d_compress, do3d_comp_prec
60	!
61	! 1320 2014-03-20 08:40:49Z raasch
62	! ONLY-attribute added to USE-statements,
63	! kind-parameters added to all INTEGER and REAL declaration statements,
64	! kinds are defined in new module kinds,
65	! old module precision_kind is removed,
66	! revision history before 2012 removed,
67	! comment fields (!:) to be used for variable explanations added to
68	! all variable declaration statements
69	!
70	! 1318 2014-03-17 13:35:16Z raasch
71	! +cpu_log_barrierwait in d3par
72	!
73	! 1301 2014-03-06 13:29:46Z heinze
74	! +large_scale_subsidence
75	!
76	! 1241 2013-10-30 11:36:58Z heinze
77	! +nudging
78	! +large_scale_forcing
79	!
80	! 1216 2013-08-26 09:31:42Z raasch
81	! +transpose_compute_overlap in inipar
82	!
83	! 1195 2013-07-01 12:27:57Z heinze
84	! Bugfix: allocate ref_state
85	!
86	! 1179 2013-06-14 05:57:58Z raasch
87	! +reference_state in inipar
88	!
89	! 1159 2013-05-21 11:58:22Z fricke
90	! +use_cmax
91	!
92	! 1128 2013-04-12 06:19:32Z raasch
93	! +background_communication in inipar
94	!
95	! 1115 2013-03-26 18:16:16Z hoffmann
96	! unused variables removed
97	!
98	! 1092 2013-02-02 11:24:22Z raasch
99	! unused variables removed
100	!
101	! 1065 2012-11-22 17:42:36Z hoffmann
102	! +nc, c_sedimentation, limiter_sedimentation, turbulence
103	! -mu_constant, mu_constant_value
104	!
105	! 1053 2012-11-13 17:11:03Z hoffmann
106	! necessary expansions according to the two new prognostic equations (nr, qr)
107	! of the two-moment cloud physics scheme and steering parameters:
108	! +_init, _surface, _surface_initial_change, _vertical_gradient,
109	! +_vertical_gradient_level, surface_waterflux_,
110	! +cloud_scheme, drizzle, mu_constant, mu_constant_value, ventilation_effect
111	!
112	! 1036 2012-10-22 13:43:42Z raasch
113	! code put under GPL (PALM 3.9)
114	!
115	! 1015 2012-09-27 09:23:24Z raasch
116	! -adjust_mixing_length
117	!
118	! 1003 2012-09-14 14:35:53Z raasch
119	! -grid_matching
120	!
121	! 1001 2012-09-13 14:08:46Z raasch
122	! -cut_spline_overshoot, long_filter_factor, overshoot_limit_, ups_limit_
123	!
124	! 996 2012-09-07 10:41:47Z raasch
125	! -use_prior_plot1d_parameters
126	!
127	! 978 2012-08-09 08:28:32Z fricke
128	! -km_damp_max, outflow_damping_width
129	! +pt_damping_factor, pt_damping_width
130	! +z0h_factor
131	!
132	! 964 2012-07-26 09:14:24Z raasch
133	! -cross_normalized_x, cross_normalized_y, cross_xtext, z_max_do1d,
134	! z_max_do1d_normalized
135	!
136	! 940 2012-07-09 14:31:00Z raasch
137	! +neutral in inipar
138	!
139	! 927 2012-06-06 19:15:04Z raasch
140	! +masking_method in inipar
141	!
142	! 824 2012-02-17 09:09:57Z raasch
143	! +curvature_solution_effects in inipar
144	!
145	! 809 2012-01-30 13:32:58Z maronga
146	! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
147	!
148	! 807 2012-01-25 11:53:51Z maronga
149	! New cpp directive "__check" implemented which is used by check_namelist_files
150	!
151	! Revision 1.1 1997/07/24 11:22:50 raasch
152	! Initial revision
153	!
154	!
155	! Description:
156	! ------------
157	! This subroutine reads variables controling the run from the NAMELIST files
158	!------------------------------------------------------------------------------!
159
160	USE arrays_3d, &
161	ONLY: pt_init, q_init, ref_state, sa_init, ug, u_init, v_init, &
162	vg
163
164	USE cloud_parameters, &
165	ONLY: c_sedimentation, curvature_solution_effects, &
166	limiter_sedimentation, nc_const, ventilation_effect
167
168	USE control_parameters, &
169	ONLY: alpha_surface, averaging_interval, averaging_interval_pr, &
170	bc_e_b, bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, &
171	bc_q_b, bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t, &
172	bottom_salinityflux, building_height, building_length_x, &
173	building_length_y, building_wall_left, building_wall_south, &
174	call_microphysics_at_all_substeps, call_psolver_at_all_substeps,&
175	canyon_height, &
176	canyon_width_x, canyon_width_y, canyon_wall_left, &
177	canyon_wall_south, cfl_factor, &
178	cloud_droplets, cloud_physics, cloud_scheme, &
179	cloud_top_radiation, conserve_volume_flow, &
180	conserve_volume_flow_mode, coupling_start_time, &
181	create_disturbances, cycle_mg, data_output, data_output_masks, &
182	data_output_pr, data_output_2d_on_each_pe, &
183	disturbance_amplitude, disturbance_energy_limit, &
184	disturbance_level_b, disturbance_level_t, dissipation_1d, &
185	do2d_at_begin, do3d_at_begin, &
186	dp_external, dp_level_b, dp_smooth, dpdxy, &
187	drizzle, dt, dt_averaging_input, dt_averaging_input_pr, &
188	dt_coupling, dt_data_output, dt_data_output_av, dt_disturb, &
189	dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy, &
190	dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart, &
191	dt_run_control, dz, dz_max, dz_stretch_factor, dz_stretch_level,&
192	end_time, ensemble_member_nr, e_init, e_min, fft_method, &
193	force_print_header, &
194	galilei_transformation, host, humidity, inflow_damping_height, &
195	inflow_damping_width, inflow_disturbance_begin, &
196	inflow_disturbance_end, initializing_actions, io_blocks, &
197	io_group, km_constant, &
198	large_scale_forcing, &
199	large_scale_subsidence, &
200	loop_optimization, masking_method, mask_scale_x, mask_scale_y, &
201	mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop, &
202	mask_y_loop, mask_z_loop, maximum_cpu_time_allowed, &
203	maximum_parallel_io_streams, max_pr_user, message_string, &
204	mg_cycles, mg_switch_to_pe0_level, mixing_length_1d, &
205	momentum_advec, netcdf_data_format, netcdf_precision, neutral, &
206	ngsrb, normalizing_region, nsor, nsor_ini, nudging, ocean, &
207	omega, omega_sor, passive_scalar, phi, nz_do3d, &
208	prandtl_layer, prandtl_number, precipitation, &
209	precipitation_amount_interval, psolver, pt_damping_factor, &
210	pt_damping_width, pt_reference, pt_surface, &
211	pt_surface_initial_change, pt_vertical_gradient, &
212	pt_vertical_gradient_level, q_surface, &
213	q_surface_initial_change, q_vertical_gradient, &
214	q_vertical_gradient_level, random_generator, &
215	random_heatflux, rayleigh_damping_factor, &
216	rayleigh_damping_height, recycling_width, reference_state, &
217	residual_limit, restart_time, return_addres, return_username, &
218	revision, rif_max, rif_min, roughness_length, runnr, &
219	run_identifier, sa_surface, sa_vertical_gradient, &
220	sa_vertical_gradient_level, scalar_advec, &
221	scalar_rayleigh_damping, &
222	section_xy, section_xz, section_yz, skip_time_data_output, &
223	skip_time_data_output_av, skip_time_dopr, skip_time_do2d_xy, &
224	skip_time_do2d_xz, skip_time_do2d_yz, skip_time_do3d, &
225	skip_time_domask, subs_vertical_gradient, &
226	subs_vertical_gradient_level, surface_heatflux, &
227	surface_pressure, surface_scalarflux, surface_waterflux, &
228	synchronous_exchange,s_surface, s_surface_initial_change, &
229	s_vertical_gradient, s_vertical_gradient_level, &
230	termination_time_needed, timestep_scheme, topography, &
231	topography_grid_convention, top_heatflux, top_momentumflux_u, &
232	top_momentumflux_v, top_salinityflux, &
233	transpose_compute_overlap, turbulence, turbulent_inflow, &
234	ug_surface, ug_vertical_gradient, ug_vertical_gradient_level, &
235	use_subsidence_tendencies, use_surface_fluxes, use_cmax, &
236	use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke, &
237	uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient,&
238	vg_vertical_gradient_level, v_bulk, v_profile, &
239	wall_adjustment, wall_heatflux, wall_humidityflux, &
240	wall_scalarflux, write_binary, z0h_factor, z_max_do2d
241
242	USE cpulog, &
243	ONLY: cpu_log_barrierwait
244
245	USE dvrp_variables, &
246	ONLY: local_dvrserver_running
247
248	USE grid_variables, &
249	ONLY: dx, dy
250
251	USE indices, &
252	ONLY: nx, ny, nz
253
254	USE model_1d, &
255	ONLY: damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
256
257	USE pegrid
258
259	USE profil_parameter, &
260	ONLY: cross_profiles, cross_ts_uymax, cross_ts_uymin, &
261	profile_columns, profile_rows
262
263	USE progress_bar, &
264	ONLY : batch_job
265
266	USE statistics, &
267	ONLY: hom, hom_sum, pr_palm, region, statistic_regions
268
269	IMPLICIT NONE
270
271	INTEGER(iwp) :: i !:
272
273
274	NAMELIST /inipar/ alpha_surface, background_communication, bc_e_b, bc_lr, &
275	bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, &
276	bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t, &
277	bottom_salinityflux, building_height, building_length_x, &
278	building_length_y, building_wall_left, building_wall_south, &
279	call_psolver_at_all_substeps, call_microphysics_at_all_substeps, &
280	canyon_height, &
281	canyon_width_x, canyon_width_y, canyon_wall_left, &
282	canyon_wall_south, c_sedimentation, cfl_factor, cloud_droplets, &
283	cloud_physics, cloud_scheme, cloud_top_radiation, &
284	collective_wait, conserve_volume_flow, &
285	conserve_volume_flow_mode, coupling_start_time, &
286	curvature_solution_effects, cycle_mg, damp_level_1d, &
287	dissipation_1d, &
288	dp_external, dp_level_b, dp_smooth, dpdxy, &
289	drizzle, dt, dt_pr_1d, dt_run_control_1d, dx, dy, dz, dz_max, &
290	dz_stretch_factor, dz_stretch_level, end_time_1d, &
291	ensemble_member_nr, &
292	e_init, e_min, fft_method, galilei_transformation, humidity, &
293	inflow_damping_height, inflow_damping_width, &
294	inflow_disturbance_begin, inflow_disturbance_end, &
295	initializing_actions, km_constant, &
296	large_scale_forcing, large_scale_subsidence, &
297	limiter_sedimentation, &
298	loop_optimization, masking_method, mg_cycles, &
299	mg_switch_to_pe0_level, mixing_length_1d, momentum_advec, &
300	nc_const, netcdf_precision, neutral, ngsrb, &
301	nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor, &
302	passive_scalar, phi, prandtl_layer, &
303	prandtl_number, precipitation, psolver, pt_damping_factor, &
304	pt_damping_width, pt_reference, pt_surface, &
305	pt_surface_initial_change, pt_vertical_gradient, &
306	pt_vertical_gradient_level, q_surface, q_surface_initial_change, &
307	q_vertical_gradient, q_vertical_gradient_level, &
308	random_generator, random_heatflux, &
309	rayleigh_damping_factor, rayleigh_damping_height, recycling_width,&
310	reference_state, residual_limit, &
311	rif_max, rif_min, roughness_length, sa_surface, &
312	sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec, &
313	scalar_rayleigh_damping, &
314	statistic_regions, subs_vertical_gradient, &
315	subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
316	surface_scalarflux, surface_waterflux, &
317	s_surface, s_surface_initial_change, s_vertical_gradient, &
318	s_vertical_gradient_level, timestep_scheme, &
319	topography, topography_grid_convention, top_heatflux, &
320	top_momentumflux_u, top_momentumflux_v, top_salinityflux, &
321	transpose_compute_overlap, turbulence, turbulent_inflow, &
322	use_subsidence_tendencies, ug_surface, ug_vertical_gradient, &
323	ug_vertical_gradient_level, use_surface_fluxes, use_cmax, &
324	use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke, &
325	uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient, &
326	vg_vertical_gradient_level, v_bulk, v_profile, ventilation_effect,&
327	wall_adjustment, wall_heatflux, wall_humidityflux, &
328	wall_scalarflux, z0h_factor
329
330	NAMELIST /d3par/ averaging_interval, averaging_interval_pr, &
331	cpu_log_barrierwait, create_disturbances, &
332	cross_profiles, cross_ts_uymax, cross_ts_uymin, &
333	data_output, data_output_masks, &
334	data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, &
335	disturbance_energy_limit, disturbance_level_b, &
336	disturbance_level_t, do2d_at_begin, do3d_at_begin, &
337	dt, dt_averaging_input, dt_averaging_input_pr, &
338	dt_coupling, dt_data_output, dt_data_output_av, dt_disturb, &
339	dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy, &
340	dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart, &
341	dt_run_control,end_time, force_print_header, mask_scale_x, &
342	mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop, &
343	mask_y_loop, mask_z_loop, netcdf_data_format, normalizing_region, &
344	npex, npey, nz_do3d, precipitation_amount_interval, &
345	profile_columns, profile_rows, restart_time, section_xy, &
346	section_xz, section_yz, skip_time_data_output, &
347	skip_time_data_output_av, skip_time_dopr, skip_time_do2d_xy, &
348	skip_time_do2d_xz, skip_time_do2d_yz, skip_time_do3d, &
349	skip_time_domask, synchronous_exchange, termination_time_needed, &
350	z_max_do2d
351
352
353	NAMELIST /envpar/ batch_job, host, local_dvrserver_running, &
354	maximum_cpu_time_allowed, maximum_parallel_io_streams, &
355	revision, return_addres, return_username, &
356	run_identifier, tasks_per_node, write_binary
357
358	!
359	!-- First read values of environment variables (this NAMELIST file is
360	!-- generated by mrun)
361	CALL location_message( 'reading environment parameters from ENVPAR', .FALSE. )
362	OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', ERR=30 )
363	READ ( 90, envpar, ERR=31, END=32 )
364	CLOSE ( 90 )
365	CALL location_message( 'finished', .TRUE. )
366
367	!
368	!-- Calculate the number of groups into which parallel I/O is split.
369	!-- The default for files which are opened by all PEs (or where each
370	!-- PE opens his own independent file) is, that all PEs are doing input/output
371	!-- in parallel at the same time. This might cause performance or even more
372	!-- severe problems depending on the configuration of the underlying file
373	!-- system.
374	!-- First, set the default:
375	IF ( maximum_parallel_io_streams == -1 .OR. &
376	maximum_parallel_io_streams > numprocs ) THEN
377	maximum_parallel_io_streams = numprocs
378	ENDIF
379	!
380	!-- Now calculate the number of io_blocks and the io_group to which the
381	!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
382	!-- for all PEs belonging to the same group.
383	!-- These settings are repeated in init_pegrid for the communicator comm2d,
384	!-- which is not available here
385	io_blocks = numprocs / maximum_parallel_io_streams
386	io_group = MOD( myid+1, io_blocks )
387
388	CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )
389	!
390	!-- Data is read in parallel by groups of PEs
391	DO i = 0, io_blocks-1
392	IF ( i == io_group ) THEN
393
394	!
395	!-- Open the NAMELIST-file which is send with this job
396	CALL check_open( 11 )
397
398	!
399	!-- Read the control parameters for initialization.
400	!-- The namelist "inipar" must be provided in the NAMELIST-file.
401	READ ( 11, inipar, ERR=10, END=11 )
402
403	#if defined ( __check )
404	!
405	!-- In case of a namelist file check, &inipar from the p3d file is
406	!-- used. The p3d file here must be closed and the p3df file for reading
407	!-- 3dpar is opened.
408	IF ( check_restart == 1 ) THEN
409	CALL close_file( 11 )
410	check_restart = 2
411	CALL check_open( 11 )
412	initializing_actions = 'read_restart_data'
413	ENDIF
414	#endif
415	GOTO 12
416
417	10 message_string = 'errors in \$inipar &or no \$inipar-namelist ' // &
418	'found (CRAY-machines only)'
419	CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
420
421	11 message_string = 'no \$inipar-namelist found'
422	CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
423
424	!
425	!-- If required, read control parameters from restart file (produced by
426	!-- a prior run). All PEs are reading from file created by PE0 (see
427	!-- check_open)
428	12 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN
429	#if ! defined ( __check )
430	CALL read_var_list
431	!
432	!-- The restart file will be reopened when reading the subdomain data
433	CALL close_file( 13 )
434
435	!
436	!-- Increment the run count
437	runnr = runnr + 1
438	#endif
439	ENDIF
440
441	!
442	!-- Definition of names of areas used for computing statistics. They must
443	!-- be defined at this place, because they are allowed to be redefined by
444	!-- the user in user_parin.
445	region = 'total domain'
446
447	!
448	!-- Read runtime parameters given by the user for this run (namelist
449	!-- "d3par"). The namelist "d3par" can be omitted. In that case, default
450	!-- values are used for the parameters.
451	READ ( 11, d3par, END=20 )
452
453	!
454	!-- Read control parameters for optionally used model software packages
455	20 CALL package_parin
456
457	!
458	!-- Read user-defined variables
459	CALL user_parin
460
461	!
462	!-- Check in case of initial run, if the grid point numbers are well
463	!-- defined and allocate some arrays which are already needed in
464	!-- init_pegrid or check_parameters. During restart jobs, these arrays
465	!-- will be allocated in read_var_list. All other arrays are allocated
466	!-- in init_3d_model.
467	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
468
469	IF ( nx <= 0 ) THEN
470	WRITE( message_string, * ) 'no value or wrong value given', &
471	' for nx: nx=', nx
472	CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
473	ENDIF
474	IF ( ny <= 0 ) THEN
475	WRITE( message_string, * ) 'no value or wrong value given', &
476	' for ny: ny=', ny
477	CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
478	ENDIF
479	IF ( nz <= 0 ) THEN
480	WRITE( message_string, * ) 'no value or wrong value given', &
481	' for nz: nz=', nz
482	CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
483	ENDIF
484	!
485	!-- ATTENTION: in case of changes to the following statement please
486	!-- also check the allocate statement in routine read_var_list
487	ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), &
488	ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1), &
489	u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1), &
490	hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions), &
491	hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
492
493	hom = 0.0_wp
494
495	ENDIF
496
497	!
498	!-- NAMELIST-file is not needed anymore
499	CALL close_file( 11 )
500
501	ENDIF
502	#if defined( __parallel ) && ! ( __check )
503	CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
504	#endif
505	ENDDO
506
507	CALL location_message( 'finished', .TRUE. )
508
509	RETURN
510
511	30 message_string = 'local file ENVPAR not found' // &
512	'&some variables for steering may not be properly set'
513	CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
514	RETURN
515
516	31 message_string = 'errors in local file ENVPAR' // &
517	'&some variables for steering may not be properly set'
518	CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
519	RETURN
520
521	32 message_string = 'no envpar-NAMELIST found in local file ENVPAR' // &
522	'&some variables for steering may not be properly set'
523	CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
524
525	END SUBROUTINE parin

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |