Home

Context Navigation

source: palm/trunk/SOURCE/read_var_list.f90 @ 1016

Last change on this file since 1016 was 1015, checked in by raasch, 12 years ago
Starting with changes required for GPU optimization. OpenACC statements for using NVIDIA GPUs added. Adjustment of mixing length to the Prandtl mixing length at first grid point above ground removed. mask array is set zero for ghost boundaries
Property svn:keywords set to `Id`
File size: 32.4 KB

Rev	Line
[1]	1	SUBROUTINE read_var_list
	2
	3	!------------------------------------------------------------------------------!
[257]	4	! Current revisions:
[622]	5	! ------------------
[1015]	6	! -adjust_mixing_length
[601]	7	!
	8	! Former revisions:
	9	! -----------------
	10	! $Id: read_var_list.f90 1015 2012-09-27 09:23:24Z raasch $
	11	!
[1004]	12	! 1003 2012-09-14 14:35:53Z raasch
	13	! -grid_matching
	14	!
[1002]	15	! 1001 2012-09-13 14:08:46Z raasch
	16	! -cut_spline_overshoot, dt_fixed, last_dt_change, long_filter_factor,
	17	! overshoot_limit_, ups_limit_
	18	!
[979]	19	! 978 2012-08-09 08:28:32Z fricke
	20	! -km_damp_max, outflow_damping_width
	21	! +pt_damping_factor, pt_damping_width
	22	! +z0h_factor
	23	!
[941]	24	! 940 2012-07-09 14:31:00Z raasch
	25	! +neutral
	26	!
[928]	27	! 927 2012-06-06 19:15:04Z raasch
	28	! +masking_method
	29	!
[850]	30	! 849 2012-03-15 10:35:09Z raasch
	31	! first_call_advec_particles renamed first_call_lpm
	32	!
[826]	33	! 824 2012-02-17 09:09:57Z raasch
	34	! +curvature_solution_effects
	35	!
[627]	36	! 626 2010-12-10 13:04:12Z suehring
	37	! idum replaced by cdum in read_parts_of_var_list
	38	!
[623]	39	! 622 2010-12-10 08:08:13Z raasch
	40	! +collective_wait
	41	!
[601]	42	! 600 2010-11-24 16:10:51Z raasch
[600]	43	! +call_psolver_at_all_substeps, cfl_factor, cycle_mg, mg_cycles,
	44	! mg_switch_to_pe0_level, ngsrb, nsor, omega_sor, psolver,
	45	! rayleigh_damping_factor, rayleigh_damping_height, residual_limit
	46	! in routine skip_var_list (end of this file), variable ldum is replaced
	47	! by cdum(LEN=1), because otherwise read errors (too few data on file)
	48	! appear due to one of the additional parameters (cycle_mg) which are now
	49	! stored on the restart file
[392]	50	!
[592]	51	! 591 2010-10-28 06:35:52Z helmke
	52	! remove print command
	53	!
[588]	54	! 587 2010-10-27 08:36:51Z helmke
	55	! +time_domask
	56	!
[581]	57	! 580 2010-10-05 13:59:11Z heinze
	58	! Renaming of ws_vertical_gradient to subs_vertical_gradient,
	59	! ws_vertical_gradient_level to subs_vertical_gradient_level and
	60	! ws_vertical_gradient_level_ind to subs_vertical_gradient_level_i
	61	!
[482]	62	! 411 2009-12-11 14:15:58Z heinze
	63	! +large_scale_subsidence, ws_vertical_gradient, ws_vertical_gradient_level,
	64	! ws_vertical_gradient_level_ind
	65	!
[392]	66	! 345 2009-07-01 14:37:56Z heinze
[345]	67	! +output_for_t0
[336]	68	! dt_fixed is read into a dummy variable.
[257]	69	! Output of messages replaced by message handling routine.
[240]	70	! +canyon_height, canyon_width_x, canyon_width_y, canyon_wall_left,
[291]	71	! canyon_wall_south, conserve_volume_flow_mode, coupling_start_time,
	72	! dp_external, dp_level_b, dp_smooth, dpdxy, run_coupled,
	73	! time_since_reference_point, topography_grid_convention, u_bulk, v_bulk
[198]	74	!
[226]	75	! 216 2008-11-25 07:12:43Z raasch
	76	! limitations for nx_on_file, ny_on_file removed (read_parts_of_var_list)
	77	!
[198]	78	! 173 2008-05-23 20:39:38Z raasch
[153]	79	! +cthf, leaf_surface_concentration, scalar_exchange_coefficient
[151]	80	! +numprocs_previous_run, hor_index_bounds_previous_run, inflow_damping_factor,
	81	! inflow_damping_height, inflow_damping_width, mean_inflow_profiles,
	82	! recycling_width, turbulent_inflow,
[145]	83	! -cross_ts_*, npex, npey,
	84	! hom_sum, volume_flow_area, volume_flow_initial moved from
[146]	85	! read_3d_binary to here,
[147]	86	! routines read_parts_of_var_list and skip_var_list added at the end
[138]	87	!
[139]	88	! 138 2007-11-28 10:03:58Z letzel
	89	! +canopy_mode, drag_coefficient, lad, lad_surface, lad_vertical_gradient,
	90	! lad_vertical_gradient_level, lad_vertical_gradient_level_ind, pch_index,
	91	! plant_canopy, time_sort_particles
	92	!
[110]	93	! 102 2007-07-27 09:09:17Z raasch
	94	! +time_coupling, top_momentumflux_u\|v
	95	!
[98]	96	! 95 2007-06-02 16:48:38Z raasch
	97	! +bc_sa_t, ocean, sa_init, sa_surface, sa_vertical_gradient,
	98	! sa_vertical_gradient_level, bottom/top_salinity_flux
	99	!
[90]	100	! 87 2007-05-22 15:46:47Z raasch
	101	! +max_pr_user (version 3.1), var_hom renamed pr_palm
	102	!
[77]	103	! 75 2007-03-22 09:54:05Z raasch
	104	! +loop_optimization, pt_reference, moisture renamed humidity
	105	!
[39]	106	! 20 2007-02-26 00:12:32Z raasch
	107	! +top_heatflux, use_top_fluxes
	108	!
[3]	109	! RCS Log replace by Id keyword, revision history cleaned up
	110	!
[1]	111	! Revision 1.34 2006/08/22 14:14:27 raasch
	112	! +dz_max
	113	!
	114	! Revision 1.1 1998/03/18 20:18:48 raasch
	115	! Initial revision
	116	!
	117	!
	118	! Description:
	119	! ------------
[143]	120	! Reading values of global control variables from restart-file (binary format)
[1]	121	!------------------------------------------------------------------------------!
	122
	123	USE arrays_3d
	124	USE averaging
[824]	125	USE cloud_parameters
[87]	126	USE control_parameters
[1]	127	USE grid_variables
	128	USE indices
	129	USE model_1d
[345]	130	USE netcdf_control
[116]	131	USE particle_attributes
[1]	132	USE pegrid
	133	USE profil_parameter
	134	USE statistics
	135
	136	IMPLICIT NONE
	137
	138	CHARACTER (LEN=10) :: binary_version, version_on_file
	139	CHARACTER (LEN=30) :: variable_chr
[87]	140
[336]	141	LOGICAL :: ldum
[143]	142
[336]	143
[1]	144	CALL check_open( 13 )
	145
	146	!
	147	!-- Make version number check first
	148	READ ( 13 ) version_on_file
[600]	149	binary_version = '3.5'
[1]	150	IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN
[257]	151	WRITE( message_string, * ) 'version mismatch concerning control ', &
	152	'variables', &
	153	'&version on file = "', &
	154	TRIM( version_on_file ), '"', &
	155	'&version on program = "', &
	156	TRIM( binary_version ), '"'
	157	CALL message( 'read_var_list', 'PA0296', 1, 2, 0, 6, 0 )
[1]	158	ENDIF
	159
	160	!
[145]	161	!-- Read number of PEs and horizontal index bounds of all PEs used in previous
	162	!-- run
	163	READ ( 13 ) variable_chr
	164	IF ( TRIM( variable_chr ) /= 'numprocs' ) THEN
[274]	165	WRITE( message_string, * ) 'numprocs not found in data from prior ', &
	166	'run on PE ', myid
[257]	167	CALL message( 'read_var_list', 'PA0297', 1, 2, 0, 6, 0 )
[145]	168	ENDIF
	169	READ ( 13 ) numprocs_previous_run
	170
[146]	171	IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) ) THEN
	172	ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
	173	ENDIF
[145]	174
	175	READ ( 13 ) variable_chr
	176	IF ( TRIM( variable_chr ) /= 'hor_index_bounds' ) THEN
[257]	177	WRITE( message_string, * ) 'hor_index_bounds not found in data from ', &
	178	'prior run on PE ', myid
	179	CALL message( 'read_var_list', 'PA0298', 1, 2, 0, 6, 0 )
[145]	180	ENDIF
	181	READ ( 13 ) hor_index_bounds_previous_run
	182
	183	!
[1]	184	!-- Read vertical number of gridpoints and number of different areas used
	185	!-- for computing statistics. Allocate arrays depending on these values,
	186	!-- which are needed for the following read instructions.
	187	READ ( 13 ) variable_chr
	188	IF ( TRIM( variable_chr ) /= 'nz' ) THEN
[257]	189	WRITE( message_string, * ) 'nz not found in data from prior run on PE ',&
	190	myid
	191	CALL message( 'read_var_list', 'PA0299', 1, 2, 0, 6, 0 )
[1]	192	ENDIF
	193	READ ( 13 ) nz
[87]	194
[1]	195	READ ( 13 ) variable_chr
[87]	196	IF ( TRIM( variable_chr ) /= 'max_pr_user' ) THEN
[257]	197	WRITE( message_string, * ) 'max_pr_user not found in data from ', &
[87]	198	'prior run on PE ', myid
[257]	199	CALL message( 'read_var_list', 'PA0300', 1, 2, 0, 6, 0 )
[87]	200	ENDIF
[159]	201	READ ( 13 ) max_pr_user ! This value is checked against the number of
	202	! user profiles given for the current run
	203	! in routine user_parin (it has to match)
[87]	204
	205	READ ( 13 ) variable_chr
[1]	206	IF ( TRIM( variable_chr ) /= 'statistic_regions' ) THEN
[257]	207	WRITE( message_string, * ) 'statistic_regions not found in data from ', &
[1]	208	'prior run on PE ', myid
[257]	209	CALL message( 'read_var_list', 'PA0301', 1, 2, 0, 6, 0 )
[1]	210	ENDIF
	211	READ ( 13 ) statistic_regions
[146]	212	IF ( .NOT. ALLOCATED( ug ) ) THEN
[147]	213	ALLOCATE( lad(0:nz+1), ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1), &
	214	v_init(0:nz+1), pt_init(0:nz+1), q_init(0:nz+1), &
	215	sa_init(0:nz+1), &
[146]	216	hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions), &
	217	hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
	218	ENDIF
[1]	219
	220	!
	221	!-- Now read all control parameters:
	222	!-- Caution: When the following read instructions have been changed, the
	223	!-- ------- version number stored in the variable binary_version has to be
	224	!-- increased. The same changes must also be done in write_var_list.
	225	READ ( 13 ) variable_chr
	226	DO WHILE ( TRIM( variable_chr ) /= '* end *' )
	227
	228	SELECT CASE ( TRIM( variable_chr ) )
	229
	230	CASE ( 'advected_distance_x' )
	231	READ ( 13 ) advected_distance_x
	232	CASE ( 'advected_distance_y' )
	233	READ ( 13 ) advected_distance_y
	234	CASE ( 'alpha_surface' )
	235	READ ( 13 ) alpha_surface
	236	CASE ( 'average_count_pr' )
	237	READ ( 13 ) average_count_pr
	238	CASE ( 'average_count_sp' )
	239	READ ( 13 ) average_count_sp
	240	CASE ( 'average_count_3d' )
	241	READ ( 13 ) average_count_3d
	242	CASE ( 'bc_e_b' )
	243	READ ( 13 ) bc_e_b
	244	CASE ( 'bc_lr' )
	245	READ ( 13 ) bc_lr
	246	CASE ( 'bc_ns' )
	247	READ ( 13 ) bc_ns
	248	CASE ( 'bc_p_b' )
	249	READ ( 13 ) bc_p_b
	250	CASE ( 'bc_p_t' )
	251	READ ( 13 ) bc_p_t
	252	CASE ( 'bc_pt_b' )
	253	READ ( 13 ) bc_pt_b
	254	CASE ( 'bc_pt_t' )
	255	READ ( 13 ) bc_pt_t
	256	CASE ( 'bc_pt_t_val' )
	257	READ ( 13 ) bc_pt_t_val
	258	CASE ( 'bc_q_b' )
	259	READ ( 13 ) bc_q_b
	260	CASE ( 'bc_q_t' )
	261	READ ( 13 ) bc_q_t
	262	CASE ( 'bc_q_t_val' )
	263	READ ( 13 ) bc_q_t_val
	264	CASE ( 'bc_s_b' )
	265	READ ( 13 ) bc_s_b
	266	CASE ( 'bc_s_t' )
	267	READ ( 13 ) bc_s_t
[95]	268	CASE ( 'bc_sa_t' )
	269	READ ( 13 ) bc_sa_t
[1]	270	CASE ( 'bc_uv_b' )
	271	READ ( 13 ) bc_uv_b
	272	CASE ( 'bc_uv_t' )
	273	READ ( 13 ) bc_uv_t
[95]	274	CASE ( 'bottom_salinityflux' )
	275	READ ( 13 ) bottom_salinityflux
[1]	276	CASE ( 'building_height' )
	277	READ ( 13 ) building_height
	278	CASE ( 'building_length_x' )
	279	READ ( 13 ) building_length_x
	280	CASE ( 'building_length_y' )
	281	READ ( 13 ) building_length_y
	282	CASE ( 'building_wall_left' )
	283	READ ( 13 ) building_wall_left
	284	CASE ( 'building_wall_south' )
	285	READ ( 13 ) building_wall_south
[600]	286	CASE ( 'call_psolver_at_all_substeps' )
	287	READ ( 13 ) call_psolver_at_all_substeps
[138]	288	CASE ( 'canopy_mode' )
	289	READ ( 13 ) canopy_mode
[240]	290	CASE ( 'canyon_height' )
	291	READ ( 13 ) canyon_height
	292	CASE ( 'canyon_width_x' )
	293	READ ( 13 ) canyon_width_x
	294	CASE ( 'canyon_width_y' )
	295	READ ( 13 ) canyon_width_y
	296	CASE ( 'canyon_wall_left' )
	297	READ ( 13 ) canyon_wall_left
	298	CASE ( 'canyon_wall_south' )
	299	READ ( 13 ) canyon_wall_south
[600]	300	CASE ( 'cfl_factor' )
	301	READ ( 13 ) cfl_factor
[1]	302	CASE ( 'cloud_droplets' )
	303	READ ( 13 ) cloud_droplets
	304	CASE ( 'cloud_physics' )
	305	READ ( 13 ) cloud_physics
[622]	306	CASE ( 'collective_wait' )
	307	READ ( 13 ) collective_wait
[1]	308	CASE ( 'conserve_volume_flow' )
	309	READ ( 13 ) conserve_volume_flow
[241]	310	CASE ( 'conserve_volume_flow_mode' )
	311	READ ( 13 ) conserve_volume_flow_mode
[291]	312	CASE ( 'coupling_start_time' )
	313	READ ( 13 ) coupling_start_time
[153]	314	CASE ( 'cthf' )
	315	READ ( 13 ) cthf
[1]	316	CASE ( 'current_timestep_number' )
	317	READ ( 13 ) current_timestep_number
[824]	318	CASE ( 'curvature_solution_effects' )
	319	READ ( 13 ) curvature_solution_effects
[600]	320	CASE ( 'cycle_mg' )
	321	READ ( 13 ) cycle_mg
[1]	322	CASE ( 'damp_level_1d' )
	323	READ ( 13 ) damp_level_1d
	324	CASE ( 'dissipation_1d' )
	325	READ ( 13 ) dissipation_1d
[240]	326	CASE ( 'dp_external' )
	327	READ ( 13 ) dp_external
	328	CASE ( 'dp_level_b' )
	329	READ ( 13 ) dp_level_b
	330	CASE ( 'dp_smooth' )
	331	READ ( 13 ) dp_smooth
	332	CASE ( 'dpdxy' )
	333	READ ( 13 ) dpdxy
[138]	334	CASE ( 'drag_coefficient' )
	335	READ ( 13 ) drag_coefficient
[1]	336	CASE ( 'dt_pr_1d' )
	337	READ ( 13 ) dt_pr_1d
	338	CASE ( 'dt_run_control_1d' )
	339	READ ( 13 ) dt_run_control_1d
	340	CASE ( 'dt_3d' )
	341	READ ( 13 ) dt_3d
	342	CASE ( 'dvrp_filecount' )
	343	READ ( 13 ) dvrp_filecount
	344	CASE ( 'dx' )
	345	READ ( 13 ) dx
	346	CASE ( 'dy' )
	347	READ ( 13 ) dy
	348	CASE ( 'dz' )
	349	READ ( 13 ) dz
	350	CASE ( 'dz_max' )
	351	READ ( 13 ) dz_max
	352	CASE ( 'dz_stretch_factor' )
	353	READ ( 13 ) dz_stretch_factor
	354	CASE ( 'dz_stretch_level' )
	355	READ ( 13 ) dz_stretch_level
	356	CASE ( 'e_min' )
	357	READ ( 13 ) e_min
	358	CASE ( 'end_time_1d' )
	359	READ ( 13 ) end_time_1d
	360	CASE ( 'fft_method' )
	361	READ ( 13 ) fft_method
[849]	362	CASE ( 'first_call_lpm' )
	363	READ ( 13 ) first_call_lpm
[1]	364	CASE ( 'galilei_transformation' )
	365	READ ( 13 ) galilei_transformation
	366	CASE ( 'hom' )
	367	READ ( 13 ) hom
[143]	368	CASE ( 'hom_sum' )
	369	READ ( 13 ) hom_sum
[151]	370	CASE ( 'humidity' )
	371	READ ( 13 ) humidity
	372	CASE ( 'inflow_damping_factor' )
	373	IF ( .NOT. ALLOCATED( inflow_damping_factor ) ) THEN
	374	ALLOCATE( inflow_damping_factor(0:nz+1) )
	375	ENDIF
	376	READ ( 13 ) inflow_damping_factor
	377	CASE ( 'inflow_damping_height' )
	378	READ ( 13 ) inflow_damping_height
	379	CASE ( 'inflow_damping_width' )
	380	READ ( 13 ) inflow_damping_width
[1]	381	CASE ( 'inflow_disturbance_begin' )
	382	READ ( 13 ) inflow_disturbance_begin
	383	CASE ( 'inflow_disturbance_end' )
	384	READ ( 13 ) inflow_disturbance_end
	385	CASE ( 'km_constant' )
	386	READ ( 13 ) km_constant
[138]	387	CASE ( 'lad' )
	388	READ ( 13 ) lad
	389	CASE ( 'lad_surface' )
	390	READ ( 13 ) lad_surface
	391	CASE ( 'lad_vertical_gradient' )
	392	READ ( 13 ) lad_vertical_gradient
	393	CASE ( 'lad_vertical_gradient_level' )
	394	READ ( 13 ) lad_vertical_gradient_level
	395	CASE ( 'lad_vertical_gradient_level_in' )
	396	READ ( 13 ) lad_vertical_gradient_level_ind
[411]	397	CASE ( 'large_scale_subsidence' )
	398	READ ( 13 ) large_scale_subsidence
[153]	399	CASE ( 'leaf_surface_concentration' )
	400	READ ( 13 ) leaf_surface_concentration
[63]	401	CASE ( 'loop_optimization' )
	402	READ ( 13 ) loop_optimization
[927]	403	CASE ( 'masking_method' )
	404	READ ( 13 ) masking_method
[151]	405	CASE ( 'mean_inflow_profiles' )
	406	IF ( .NOT. ALLOCATED( mean_inflow_profiles ) ) THEN
	407	ALLOCATE( mean_inflow_profiles(0:nz+1,5) )
	408	ENDIF
	409	READ ( 13 ) mean_inflow_profiles
[600]	410	CASE ( 'mg_cycles' )
	411	READ ( 13 ) mg_cycles
	412	CASE ( 'mg_switch_to_pe0_level' )
	413	READ ( 13 ) mg_switch_to_pe0_level
[1]	414	CASE ( 'mixing_length_1d' )
	415	READ ( 13 ) mixing_length_1d
	416	CASE ( 'momentum_advec' )
	417	READ ( 13 ) momentum_advec
	418	CASE ( 'netcdf_precision' )
	419	READ ( 13 ) netcdf_precision
[940]	420	CASE ( 'neutral' )
	421	READ ( 13 ) neutral
[600]	422	CASE ( 'ngsrb' )
	423	READ ( 13 ) ngsrb
	424	CASE ( 'nsor' )
	425	READ ( 13 ) nsor
[1]	426	CASE ( 'nsor_ini' )
	427	READ ( 13 ) nsor_ini
	428	CASE ( 'nx' )
	429	READ ( 13 ) nx
[147]	430	nx_on_file = nx
[1]	431	CASE ( 'ny' )
	432	READ ( 13 ) ny
[147]	433	ny_on_file = ny
[94]	434	CASE ( 'ocean' )
	435	READ ( 13 ) ocean
[1]	436	CASE ( 'old_dt' )
	437	READ ( 13 ) old_dt
	438	CASE ( 'omega' )
	439	READ ( 13 ) omega
[600]	440	CASE ( 'omega_sor' )
	441	READ ( 13 ) omega_sor
[345]	442	CASE ( 'output_for_t0' )
	443	READ (13) output_for_t0
[1]	444	CASE ( 'passive_scalar' )
	445	READ ( 13 ) passive_scalar
[138]	446	CASE ( 'pch_index' )
	447	READ ( 13 ) pch_index
[1]	448	CASE ( 'phi' )
	449	READ ( 13 ) phi
[138]	450	CASE ( 'plant_canopy' )
	451	READ ( 13 ) plant_canopy
[1]	452	CASE ( 'prandtl_layer' )
	453	READ ( 13 ) prandtl_layer
[600]	454	CASE ( 'prandtl_number' )
	455	READ ( 13 ) prandtl_number
[1]	456	CASE ( 'precipitation' )
	457	READ ( 13 ) precipitation
[600]	458	CASE ( 'psolver' )
	459	READ ( 13 ) psolver
[978]	460	CASE ( 'pt_damping_factor' )
	461	READ ( 13 ) pt_damping_factor
	462	CASE ( 'pt_damping_width' )
	463	READ ( 13 ) pt_damping_width
[1]	464	CASE ( 'pt_init' )
	465	READ ( 13 ) pt_init
[57]	466	CASE ( 'pt_reference' )
	467	READ ( 13 ) pt_reference
[1]	468	CASE ( 'pt_surface' )
	469	READ ( 13 ) pt_surface
	470	CASE ( 'pt_surface_initial_change' )
	471	READ ( 13 ) pt_surface_initial_change
	472	CASE ( 'pt_vertical_gradient' )
	473	READ ( 13 ) pt_vertical_gradient
	474	CASE ( 'pt_vertical_gradient_level' )
	475	READ ( 13 ) pt_vertical_gradient_level
	476	CASE ( 'pt_vertical_gradient_level_ind' )
	477	READ ( 13 ) pt_vertical_gradient_level_ind
	478	CASE ( 'q_init' )
	479	READ ( 13 ) q_init
	480	CASE ( 'q_surface' )
	481	READ ( 13 ) q_surface
	482	CASE ( 'q_surface_initial_change' )
	483	READ ( 13 ) q_surface_initial_change
	484	CASE ( 'q_vertical_gradient' )
	485	READ ( 13 ) q_vertical_gradient
	486	CASE ( 'q_vertical_gradient_level' )
	487	READ ( 13 ) q_vertical_gradient_level
	488	CASE ( 'q_vertical_gradient_level_ind' )
	489	READ ( 13 ) q_vertical_gradient_level_ind
	490	CASE ( 'radiation' )
	491	READ ( 13 ) radiation
	492	CASE ( 'random_generator' )
	493	READ ( 13 ) random_generator
	494	CASE ( 'random_heatflux' )
	495	READ ( 13 ) random_heatflux
[600]	496	CASE ( 'rayleigh_damping_factor' )
	497	READ ( 13 ) rayleigh_damping_factor
	498	CASE ( 'rayleigh_damping_height' )
	499	READ ( 13 ) rayleigh_damping_height
[151]	500	CASE ( 'recycling_width' )
	501	READ ( 13 ) recycling_width
[600]	502	CASE ( 'residual_limit' )
	503	READ ( 13 ) residual_limit
[1]	504	CASE ( 'rif_max' )
	505	READ ( 13 ) rif_max
	506	CASE ( 'rif_min' )
	507	READ ( 13 ) rif_min
	508	CASE ( 'roughness_length' )
	509	READ ( 13 ) roughness_length
	510	CASE ( 'runnr' )
	511	READ ( 13 ) runnr
[291]	512	CASE ( 'run_coupled' )
	513	READ ( 13 ) run_coupled
[94]	514	CASE ( 'sa_init' )
	515	READ ( 13 ) sa_init
	516	CASE ( 'sa_surface' )
	517	READ ( 13 ) sa_surface
	518	CASE ( 'sa_vertical_gradient' )
	519	READ ( 13 ) sa_vertical_gradient
	520	CASE ( 'sa_vertical_gradient_level' )
	521	READ ( 13 ) sa_vertical_gradient_level
[1]	522	CASE ( 'scalar_advec' )
	523	READ ( 13 ) scalar_advec
[153]	524	CASE ( 'scalar_exchange_coefficient' )
	525	READ ( 13 ) scalar_exchange_coefficient
[1]	526	CASE ( 'simulated_time' )
	527	READ ( 13 ) simulated_time
	528	CASE ( 'surface_heatflux' )
	529	READ ( 13 ) surface_heatflux
	530	CASE ( 'surface_pressure' )
	531	READ ( 13 ) surface_pressure
	532	CASE ( 'surface_scalarflux' )
	533	READ ( 13 ) surface_scalarflux
	534	CASE ( 'surface_waterflux' )
	535	READ ( 13 ) surface_waterflux
	536	CASE ( 's_surface' )
	537	READ ( 13 ) s_surface
	538	CASE ( 's_surface_initial_change' )
	539	READ ( 13 ) s_surface_initial_change
	540	CASE ( 's_vertical_gradient' )
	541	READ ( 13 ) s_vertical_gradient
	542	CASE ( 's_vertical_gradient_level' )
	543	READ ( 13 ) s_vertical_gradient_level
[102]	544	CASE ( 'time_coupling' )
	545	READ ( 13 ) time_coupling
[1]	546	CASE ( 'time_disturb' )
	547	READ ( 13 ) time_disturb
[587]	548	CASE ( 'time_dopr' )
[1]	549	READ ( 13 ) time_dopr
[587]	550	CASE ( 'time_domask' )
	551	READ ( 13 ) time_domask
[1]	552	CASE ( 'time_dopr_av' )
	553	READ ( 13 ) time_dopr_av
	554	CASE ( 'time_dopr_listing' )
	555	READ ( 13 ) time_dopr_listing
	556	CASE ( 'time_dopts' )
	557	READ ( 13 ) time_dopts
	558	CASE ( 'time_dosp' )
	559	READ ( 13 ) time_dosp
	560	CASE ( 'time_dots' )
	561	READ ( 13 ) time_dots
	562	CASE ( 'time_do2d_xy' )
	563	READ ( 13 ) time_do2d_xy
	564	CASE ( 'time_do2d_xz' )
	565	READ ( 13 ) time_do2d_xz
	566	CASE ( 'time_do2d_yz' )
	567	READ ( 13 ) time_do2d_yz
	568	CASE ( 'time_do3d' )
	569	READ ( 13 ) time_do3d
	570	CASE ( 'time_do_av' )
	571	READ ( 13 ) time_do_av
	572	CASE ( 'time_do_sla' )
	573	READ ( 13 ) time_do_sla
	574	CASE ( 'time_dvrp' )
	575	READ ( 13 ) time_dvrp
	576	CASE ( 'time_restart' )
	577	READ ( 13 ) time_restart
	578	CASE ( 'time_run_control' )
	579	READ ( 13 ) time_run_control
[291]	580	CASE ( 'time_since_reference_point' )
	581	READ ( 13 ) time_since_reference_point
[116]	582	CASE ( 'time_sort_particles' )
	583	READ ( 13 ) time_sort_particles
[1]	584	CASE ( 'timestep_scheme' )
	585	READ ( 13 ) timestep_scheme
	586	CASE ( 'topography' )
	587	READ ( 13 ) topography
[256]	588	CASE ( 'topography_grid_convention' )
	589	READ ( 13 ) topography_grid_convention
[19]	590	CASE ( 'top_heatflux' )
	591	READ ( 13 ) top_heatflux
[102]	592	CASE ( 'top_momentumflux_u' )
	593	READ ( 13 ) top_momentumflux_u
	594	CASE ( 'top_momentumflux_v' )
	595	READ ( 13 ) top_momentumflux_v
[94]	596	CASE ( 'top_salinityflux' )
	597	READ ( 13 ) top_salinityflux
[1]	598	CASE ( 'tsc' )
	599	READ ( 13 ) tsc
[151]	600	CASE ( 'turbulent_inflow' )
	601	READ ( 13 ) turbulent_inflow
[241]	602	CASE ( 'u_bulk' )
	603	READ ( 13 ) u_bulk
[1]	604	CASE ( 'u_init' )
	605	READ ( 13 ) u_init
	606	CASE ( 'u_max' )
	607	READ ( 13 ) u_max
	608	CASE ( 'u_max_ijk' )
	609	READ ( 13 ) u_max_ijk
	610	CASE ( 'ug' )
	611	READ ( 13 ) ug
	612	CASE ( 'ug_surface' )
	613	READ ( 13 ) ug_surface
	614	CASE ( 'ug_vertical_gradient' )
	615	READ ( 13 ) ug_vertical_gradient
	616	CASE ( 'ug_vertical_gradient_level' )
	617	READ ( 13 ) ug_vertical_gradient_level
	618	CASE ( 'ug_vertical_gradient_level_ind' )
	619	READ ( 13 ) ug_vertical_gradient_level_ind
	620	CASE ( 'use_surface_fluxes' )
	621	READ ( 13 ) use_surface_fluxes
[20]	622	CASE ( 'use_top_fluxes' )
	623	READ ( 13 ) use_top_fluxes
[1]	624	CASE ( 'use_ug_for_galilei_tr' )
	625	READ ( 13 ) use_ug_for_galilei_tr
	626	CASE ( 'use_upstream_for_tke' )
	627	READ ( 13 ) use_upstream_for_tke
[241]	628	CASE ( 'v_bulk' )
	629	READ ( 13 ) v_bulk
[1]	630	CASE ( 'v_init' )
	631	READ ( 13 ) v_init
	632	CASE ( 'v_max' )
	633	READ ( 13 ) v_max
	634	CASE ( 'v_max_ijk' )
	635	READ ( 13 ) v_max_ijk
	636	CASE ( 'vg' )
	637	READ ( 13 ) vg
	638	CASE ( 'vg_surface' )
	639	READ ( 13 ) vg_surface
	640	CASE ( 'vg_vertical_gradient' )
	641	READ ( 13 ) vg_vertical_gradient
	642	CASE ( 'vg_vertical_gradient_level' )
	643	READ ( 13 ) vg_vertical_gradient_level
	644	CASE ( 'vg_vertical_gradient_level_ind' )
	645	READ ( 13 ) vg_vertical_gradient_level_ind
[143]	646	CASE ( 'volume_flow_area' )
	647	READ ( 13 ) volume_flow_area
	648	CASE ( 'volume_flow_initial' )
	649	READ ( 13 ) volume_flow_initial
[1]	650	CASE ( 'wall_adjustment' )
	651	READ ( 13 ) wall_adjustment
[580]	652	CASE ( 'subs_vertical_gradient' )
	653	READ ( 13 ) subs_vertical_gradient
	654	CASE ( 'subs_vertical_gradient_level' )
	655	READ ( 13 ) subs_vertical_gradient_level
	656	CASE ( 'subs_vertical_gradient_level_i' )
	657	READ ( 13 ) subs_vertical_gradient_level_i
[1]	658	CASE ( 'w_max' )
	659	READ ( 13 ) w_max
	660	CASE ( 'w_max_ijk' )
	661	READ ( 13 ) w_max_ijk
[978]	662	CASE ( 'z0h_factor' )
	663	READ ( 13 ) z0h_factor
[1]	664
	665	CASE DEFAULT
[257]	666	WRITE( message_string, * ) 'unknown variable named "', &
	667	TRIM( variable_chr ), '" found in', &
	668	' data from prior run on PE ', myid
	669	CALL message( 'read_var_list', 'PA0302', 1, 2, 0, 6, 0 )
[1]	670	END SELECT
	671	!
	672	!-- Read next string
	673	READ ( 13 ) variable_chr
	674
	675	ENDDO
	676
	677
[146]	678	END SUBROUTINE read_var_list
[143]	679
[146]	680
	681
[147]	682	SUBROUTINE read_parts_of_var_list
	683
	684	!------------------------------------------------------------------------------!
	685	! Description:
	686	! ------------
	687	! Skipping the global control variables from restart-file (binary format)
	688	! except some informations needed when reading restart data from a previous
[151]	689	! run which used a smaller total domain or/and a different domain decomposition.
[147]	690	!------------------------------------------------------------------------------!
	691
[151]	692	USE arrays_3d
[147]	693	USE control_parameters
	694	USE indices
	695	USE pegrid
	696	USE statistics
	697
	698	IMPLICIT NONE
	699
	700	CHARACTER (LEN=10) :: version_on_file
	701	CHARACTER (LEN=30) :: variable_chr
[626]	702	CHARACTER (LEN=1) :: cdum
[147]	703
	704	INTEGER :: idum, max_pr_user_on_file, nz_on_file, &
[170]	705	statistic_regions_on_file, tmp_mpru, tmp_sr
[147]	706
[151]	707	REAL, DIMENSION(:,:,:), ALLOCATABLE :: hom_sum_on_file
[147]	708	REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: hom_on_file
	709
	710
	711	CALL check_open( 13 )
	712
	713	WRITE (9,*) 'rpovl: after check open 13'
	714	CALL local_flush( 9 )
	715	READ ( 13 ) version_on_file
	716
	717	!
	718	!-- Read number of PEs and horizontal index bounds of all PEs used in previous
	719	!-- run
	720	READ ( 13 ) variable_chr
	721	IF ( TRIM( variable_chr ) /= 'numprocs' ) THEN
[274]	722	WRITE( message_string, * ) 'numprocs not found in data from prior ', &
	723	'run on PE ', myid
[257]	724	CALL message( 'read_parts_of_var_list', 'PA0297', 1, 2, 0, 6, 0 )
[147]	725	ENDIF
	726	READ ( 13 ) numprocs_previous_run
	727
	728	IF ( .NOT. ALLOCATED( hor_index_bounds_previous_run ) ) THEN
	729	ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
	730	ENDIF
	731
	732	READ ( 13 ) variable_chr
	733	IF ( TRIM( variable_chr ) /= 'hor_index_bounds' ) THEN
[257]	734	WRITE( message_string, * ) 'hor_index_bounds not found in data from ', &
	735	'prior run on PE ', myid
	736	CALL message( 'read_parts_of_var_list', 'PA0298', 1, 2, 0, 6, 0 )
[147]	737	ENDIF
	738	READ ( 13 ) hor_index_bounds_previous_run
	739
	740	!
	741	!-- Read vertical number of gridpoints and number of different areas used
	742	!-- for computing statistics. Allocate arrays depending on these values,
	743	!-- which are needed for the following read instructions.
	744	READ ( 13 ) variable_chr
	745	IF ( TRIM( variable_chr ) /= 'nz' ) THEN
[257]	746	message_string = 'nz not found in restart data file'
	747	CALL message( 'read_parts_of_var_list', 'PA0303', 1, 2, 0, 6, 0 )
[147]	748	ENDIF
	749	READ ( 13 ) nz_on_file
	750	IF ( nz_on_file /= nz ) THEN
[257]	751	WRITE( message_string, * ) 'mismatch concerning number of ', &
[274]	752	'gridpoints along z', &
	753	'&nz on file = "', nz_on_file, '"', &
	754	'&nz from run = "', nz, '"'
[257]	755	CALL message( 'read_parts_of_var_list', 'PA0304', 1, 2, 0, 6, 0 )
[147]	756	ENDIF
	757
	758	READ ( 13 ) variable_chr
	759	IF ( TRIM( variable_chr ) /= 'max_pr_user' ) THEN
[257]	760	message_string = 'max_pr_user not found in restart data file'
	761	CALL message( 'read_parts_of_var_list', 'PA0305', 1, 2, 0, 6, 0 )
[147]	762	ENDIF
	763	READ ( 13 ) max_pr_user_on_file
[170]	764	IF ( max_pr_user_on_file /= max_pr_user ) THEN
[274]	765	WRITE( message_string, * ) 'number of user profiles on res', &
	766	'tart data file differs from the current ', &
	767	'run&max_pr_user on file = "', &
	768	max_pr_user_on_file, '"', &
	769	'&max_pr_user from run = "', &
[257]	770	max_pr_user, '"'
	771	CALL message( 'read_parts_of_var_list', 'PA0306', 0, 0, 0, 6, 0 )
[170]	772	tmp_mpru = MIN( max_pr_user_on_file, max_pr_user )
	773	ELSE
	774	tmp_mpru = max_pr_user
[147]	775	ENDIF
	776
	777	READ ( 13 ) variable_chr
	778	IF ( TRIM( variable_chr ) /= 'statistic_regions' ) THEN
[257]	779	message_string = 'statistic_regions not found in restart data file'
	780	CALL message( 'read_parts_of_var_list', 'PA0307', 1, 2, 0, 6, 0 )
[147]	781	ENDIF
	782	READ ( 13 ) statistic_regions_on_file
[173]	783	IF ( statistic_regions_on_file /= statistic_regions ) THEN
[274]	784	WRITE( message_string, * ) 'statistic regions on restart data file ',&
	785	'differ from the current run', &
	786	'&statistic regions on file = "', &
	787	statistic_regions_on_file, '"', &
	788	'&statistic regions from run = "', &
	789	statistic_regions, '"', &
[257]	790	'&statistic data may be lost!'
	791	CALL message( 'read_parts_of_var_list', 'PA0308', 0, 1, 0, 6, 0 )
[170]	792	tmp_sr = MIN( statistic_regions_on_file, statistic_regions )
	793	ELSE
	794	tmp_sr = statistic_regions
[147]	795	ENDIF
	796
	797
	798	!
	799	!-- Now read and check some control parameters and skip the rest
	800	WRITE (9,*) 'wpovl: begin reading variables'
	801	CALL local_flush( 9 )
	802	READ ( 13 ) variable_chr
	803
	804	DO WHILE ( TRIM( variable_chr ) /= '* end *' )
	805
	806	SELECT CASE ( TRIM( variable_chr ) )
	807
[151]	808	CASE ( 'average_count_pr' )
	809	READ ( 13 ) average_count_pr
[257]	810	IF ( average_count_pr /= 0 ) THEN
[274]	811	WRITE( message_string, * ) 'inflow profiles not temporally ', &
	812	'averaged. &Averaging will be done now using ', &
	813	average_count_pr, ' samples.'
	814	CALL message( 'read_parts_of_var_list', 'PA0309', &
	815	0, 1, 0, 6, 0 )
[151]	816	ENDIF
	817
[147]	818	CASE ( 'hom' )
	819	ALLOCATE( hom_on_file(0:nz+1,2,pr_palm+max_pr_user_on_file, &
	820	0:statistic_regions_on_file) )
	821	READ ( 13 ) hom_on_file
[170]	822	hom(:,:,1:pr_palm+tmp_mpru,0:tmp_sr) = &
	823	hom_on_file(:,:,1:pr_palm+tmp_mpru,0:tmp_sr)
[147]	824	DEALLOCATE( hom_on_file )
	825
[151]	826	CASE ( 'hom_sum' )
	827	ALLOCATE( hom_sum_on_file(0:nz+1,pr_palm+max_pr_user_on_file, &
	828	0:statistic_regions_on_file) )
	829	READ ( 13 ) hom_sum_on_file
[170]	830	hom_sum(:,1:pr_palm+tmp_mpru,0:tmp_sr) = &
	831	hom_sum_on_file(:,1:pr_palm+tmp_mpru,0:tmp_sr)
[151]	832	DEALLOCATE( hom_sum_on_file )
	833
[147]	834	CASE ( 'nx' )
	835	READ ( 13 ) nx_on_file
	836
	837	CASE ( 'ny' )
	838	READ ( 13 ) ny_on_file
	839
	840
	841	CASE DEFAULT
	842
[626]	843	READ ( 13 ) cdum
[147]	844
	845	END SELECT
	846
	847	READ ( 13 ) variable_chr
	848
	849	ENDDO
	850
[151]	851	!
	852	!-- Calculate the temporal average of vertical profiles, if neccessary
	853	IF ( average_count_pr /= 0 ) THEN
	854	hom_sum = hom_sum / REAL( average_count_pr )
	855	ENDIF
[147]	856
[151]	857
[147]	858	END SUBROUTINE read_parts_of_var_list
	859
	860
	861
[146]	862	SUBROUTINE skip_var_list
	863
	864	!------------------------------------------------------------------------------!
	865	! Description:
	866	! ------------
	867	! Skipping the global control variables from restart-file (binary format)
	868	!------------------------------------------------------------------------------!
	869
[147]	870	IMPLICIT NONE
	871
[146]	872	CHARACTER (LEN=10) :: version_on_file
	873	CHARACTER (LEN=30) :: variable_chr
	874
[600]	875	CHARACTER (LEN=1) :: cdum
[146]	876
	877
[147]	878	WRITE (9,*) 'skipvl #1'
	879	CALL local_flush( 9 )
[146]	880	READ ( 13 ) version_on_file
	881
[147]	882	WRITE (9,*) 'skipvl before variable_chr'
	883	CALL local_flush( 9 )
[146]	884	READ ( 13 ) variable_chr
[147]	885	WRITE (9,*) 'skipvl after variable_chr'
	886	CALL local_flush( 9 )
[146]	887
	888	DO WHILE ( TRIM( variable_chr ) /= '* end *' )
	889
[147]	890	WRITE (9,*) 'skipvl chr = ', variable_chr
	891	CALL local_flush( 9 )
[600]	892	READ ( 13 ) cdum
[146]	893	READ ( 13 ) variable_chr
	894
	895	ENDDO
[147]	896	WRITE (9,*) 'skipvl last'
	897	CALL local_flush( 9 )
[146]	898
	899
	900	END SUBROUTINE skip_var_list

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |