Home

Context Navigation

source: palm/trunk/SOURCE/read_3d_binary.f90 @ 236

Last change on this file since 236 was 226, checked in by raasch, 16 years ago
preparations for the next release
Property svn:keywords set to `Id`
File size: 40.8 KB

Rev	Line
[1]	1	SUBROUTINE read_3d_binary
	2
	3	!------------------------------------------------------------------------------!
	4	! Actual revisions:
	5	! -----------------
[226]	6	!
	7	!
	8	! Former revisions:
	9	! -----------------
	10	! $Id: read_3d_binary.f90 226 2009-02-02 07:39:34Z letzel $
	11	!
	12	! 220 2008-12-18 07:00:36Z raasch
[216]	13	! reading mechanism completely revised (subdomain/total domain size can vary
	14	! arbitrarily between current and previous run)
	15	! Bugfix: reading of spectrum_x\|y from restart files ignored if total numbers
	16	! of grid points do not match
[198]	17	!
	18	! 150 2008-02-29 08:19:58Z raasch
[146]	19	! Files from which restart data are to be read are determined and subsequently
[147]	20	! opened. The total domain on the restart file is allowed to be smaller than
	21	! the current total domain. In this case it will be periodically mapped on the
	22	! current domain (needed for recycling method).
[146]	23	! +call of user_read_restart_data, -dopr_time_count,
[145]	24	! hom_sum, volume_flow_area, volume_flow_initial moved to read_var_list,
[143]	25	! reading of old profil parameters (cross_..., dopr_crossindex, profile_***)
[146]	26	! removed, initialization of spectrum_x\|y removed
[1]	27	!
[110]	28	! 102 2007-07-27 09:09:17Z raasch
	29	! +uswst, uswst_m, vswst, vswst_m
	30	!
[98]	31	! 96 2007-06-04 08:07:41Z raasch
	32	! +rho_av, sa, sa_av, saswsb, saswst
	33	!
[77]	34	! 73 2007-03-20 08:33:14Z raasch
	35	! +precipitation_amount, precipitation_rate_av, rif_wall, u_m_l, u_m_r, etc.,
	36	! z0_av
	37	!
[39]	38	! 19 2007-02-23 04:53:48Z raasch
[77]	39	! +qswst, qswst_m, tswst, tswst_m
[39]	40	!
[3]	41	! RCS Log replace by Id keyword, revision history cleaned up
	42	!
[1]	43	! Revision 1.4 2006/08/04 15:02:32 raasch
	44	! +iran, iran_part
	45	!
	46	! Revision 1.1 2004/04/30 12:47:27 raasch
	47	! Initial revision
	48	!
	49	!
	50	! Description:
	51	! ------------
	52	! Binary input of variables and arrays from restart file
	53	!------------------------------------------------------------------------------!
	54
	55	USE arrays_3d
	56	USE averaging
[72]	57	USE cloud_parameters
[1]	58	USE control_parameters
	59	USE cpulog
	60	USE indices
	61	USE interfaces
	62	USE particle_attributes
	63	USE pegrid
	64	USE profil_parameter
	65	USE random_function_mod
	66	USE statistics
	67
	68	IMPLICIT NONE
	69
[146]	70	CHARACTER (LEN=5) :: myid_char_save
[1]	71	CHARACTER (LEN=10) :: binary_version, version_on_file
	72	CHARACTER (LEN=20) :: field_chr
	73
[216]	74	INTEGER :: files_to_be_opened, i, j, k, myid_on_file, &
[147]	75	numprocs_on_file, nxlc, nxlf, nxlpr, nxl_on_file, nxrc, nxrf, &
	76	nxrpr, nxr_on_file, nync, nynf, nynpr, nyn_on_file, nysc, &
	77	nysf, nyspr, nys_on_file, nzb_on_file, nzt_on_file, offset_x, &
[216]	78	offset_y, shift_x, shift_y
[1]	79
[216]	80	INTEGER, DIMENSION(numprocs_previous_run) :: file_list, overlap_count
[146]	81
[216]	82	INTEGER, DIMENSION(numprocs_previous_run,1000) :: nxlfa, nxrfa, nynfa, &
	83	nysfa, offset_xa, &
	84	offset_ya
	85	REAL :: rdummy
	86
[146]	87	REAL, DIMENSION(:,:), ALLOCATABLE :: tmp_2d
[216]	88	REAL, DIMENSION(:,:,:), ALLOCATABLE :: tmp_3d, tmp_3dwul, tmp_3dwun, &
	89	tmp_3dwur, tmp_3dwus, tmp_3dwvl, &
	90	tmp_3dwvn, tmp_3dwvr, tmp_3dwvs, &
	91	tmp_3dwwl, tmp_3dwwn, tmp_3dwwr, &
	92	tmp_3dwws
[146]	93	REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: tmp_4d
	94
	95
[1]	96	!
[146]	97	!-- Read data from previous model run.
[1]	98	CALL cpu_log( log_point_s(14), 'read_3d_binary', 'start' )
	99
	100	!
[146]	101	!-- Check which of the restart files contain data needed for the subdomain
	102	!-- of this PE
	103	files_to_be_opened = 0
[143]	104
[146]	105	DO i = 1, numprocs_previous_run
[216]	106	!
	107	!-- Store array bounds of the previous run ("pr") in temporary scalars
[147]	108	nxlpr = hor_index_bounds_previous_run(1,i-1)
	109	nxrpr = hor_index_bounds_previous_run(2,i-1)
	110	nyspr = hor_index_bounds_previous_run(3,i-1)
	111	nynpr = hor_index_bounds_previous_run(4,i-1)
	112
[143]	113	!
[147]	114	!-- Determine the offsets. They may be non-zero in case that the total domain
	115	!-- on file is smaller than the current total domain.
	116	offset_x = ( nxl / ( nx_on_file + 1 ) ) * ( nx_on_file + 1 )
	117	offset_y = ( nys / ( ny_on_file + 1 ) ) * ( ny_on_file + 1 )
	118
	119	!
[216]	120	!-- Start with this offset and then check, if the subdomain on file
	121	!-- matches another time(s) in the current subdomain by shifting it
	122	!-- for nx_on_file+1, ny_on_file+1 respectively
[146]	123
[216]	124	shift_y = 0
	125	j = 0 ! counter for the number of files to be opened
[146]	126
[216]	127	DO WHILE ( nyspr+shift_y <= nyn-offset_y .AND. &
	128	nynpr+shift_y >= nys-offset_y )
[147]	129
[216]	130	shift_x = 0
[146]	131
[216]	132	DO WHILE ( nxlpr+shift_x <= nxr-offset_x .AND. &
	133	nxrpr+shift_x >= nxl-offset_x )
[147]	134
[216]	135	j = j +1
[147]	136
[216]	137	IF ( j > 1000 ) THEN
[147]	138	!
[216]	139	!-- Array bound exceeded
	140	PRINT*, '+++ read_3d_binary: data from subdomain of previous', &
	141	' run mapped more than 1000 times'
	142	#if defined( __parallel )
	143	CALL MPI_ABORT( comm2d, 9999, ierr )
	144	#else
	145	STOP
	146	#endif
	147	ENDIF
[147]	148
[216]	149	IF ( j == 1 ) THEN
	150	files_to_be_opened = files_to_be_opened + 1
	151	file_list(files_to_be_opened) = i-1
	152	ENDIF
[147]	153
[216]	154	offset_xa(files_to_be_opened,j) = offset_x + shift_x
	155	offset_ya(files_to_be_opened,j) = offset_y + shift_y
[147]	156
	157	!
	158	!-- Index bounds of overlapping data
[216]	159	nxlfa(files_to_be_opened,j) = MAX( nxl-offset_x-shift_x, nxlpr )
	160	nxrfa(files_to_be_opened,j) = MIN( nxr-offset_x-shift_x, nxrpr )
	161	nysfa(files_to_be_opened,j) = MAX( nys-offset_y, nyspr )
	162	nynfa(files_to_be_opened,j) = MIN( nyn-offset_y, nynpr )
[147]	163
[216]	164	shift_x = shift_x + ( nx_on_file + 1 )
[147]	165
[216]	166	ENDDO
[147]	167
[216]	168	shift_y = shift_y + ( ny_on_file + 1 )
[147]	169
[216]	170	ENDDO
[147]	171
[216]	172	IF ( j > 0 ) overlap_count(files_to_be_opened) = j
[147]	173
	174	!
[216]	175	!-- Test output, to be removed later
	176	IF ( j > 0 ) THEN
	177	WRITE (9,) '** reading from file: ', i, j, ' times'
	178	WRITE (9,*) ' nxl = ', nxl, ' nxr = ', nxr, ' nys = ', nys, ' nyn = ', nyn
	179	WRITE (9,*) ' '
	180	DO k = 1, j
	181	WRITE (9,*) 'k = ', k
	182	WRITE (9,'(6(A,I2))') 'nxlfa = ', nxlfa(files_to_be_opened,k), &
	183	' nxrfa = ', nxrfa(files_to_be_opened,k), &
	184	' offset_xa = ', offset_xa(files_to_be_opened,k), &
	185	' nysfa = ', nysfa(files_to_be_opened,k), &
	186	' nynfa = ', nynfa(files_to_be_opened,k), &
	187	' offset_ya = ', offset_ya(files_to_be_opened,k)
	188	ENDDO
[147]	189	CALL local_flush( 9 )
	190	ENDIF
	191
[146]	192	ENDDO
	193
[147]	194	!
	195	!-- Save the id-string of the current process, since myid_char may now be used
	196	!-- to open files created by PEs with other id.
[146]	197	myid_char_save = myid_char
	198
[147]	199	!
	200	!-- Test output (remove later)
[146]	201	DO i = 1, numprocs_previous_run
	202	WRITE (9,*) 'i=',i-1, ' ibs= ',hor_index_bounds_previous_run(1:4,i-1)
	203	ENDDO
	204	CALL local_flush( 9 )
	205
	206	IF ( files_to_be_opened /= 1 .OR. numprocs /= numprocs_previous_run ) &
	207	THEN
	208	PRINT, '** number of PEs or virtual PE-grid changed in restart run'
	209	PRINT*, ' PE', myid, ' will read from files ', &
	210	file_list(1:files_to_be_opened)
[1]	211	ENDIF
	212
	213	!
[146]	214	!-- Read data from all restart files determined above
	215	DO i = 1, files_to_be_opened
[1]	216
[146]	217	j = file_list(i)
	218	!
	219	!-- Set the filename (underscore followed by four digit processor id)
	220	WRITE (myid_char,'(''_'',I4.4)') j
	221	WRITE (9,*) 'myid=',myid,' opening file "',myid_char,'"'
	222	CALL local_flush( 9 )
[1]	223
[146]	224	!
	225	!-- Open the restart file. If this file has been created by PE0 (_0000),
	226	!-- the global variables at the beginning of the file have to be skipped
	227	!-- first.
	228	CALL check_open( 13 )
[147]	229	WRITE (9,*) 'before skipping'
	230	CALL local_flush( 9 )
[146]	231	IF ( j == 0 ) CALL skip_var_list
[147]	232	WRITE (9,*) 'skipping done'
	233	CALL local_flush( 9 )
[1]	234
[146]	235	!
	236	!-- First compare the version numbers
	237	READ ( 13 ) version_on_file
	238	binary_version = '3.1'
	239	IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN
	240	IF ( myid == 0 ) THEN
	241	PRINT*, '+++ init_3d_model: version mismatch concerning data ', &
	242	'from prior run'
	243	PRINT*, ' version on file = "', TRIM( version_on_file ),&
	244	'"'
	245	PRINT*, ' version in program = "', TRIM( binary_version ), &
	246	'"'
	247	ENDIF
	248	CALL local_stop
	249	ENDIF
	250
	251	!
	252	!-- Read number of processors, processor-id, and array ranges.
	253	!-- Compare the array ranges with those stored in the index bound array.
	254	READ ( 13 ) numprocs_on_file, myid_on_file, nxl_on_file, nxr_on_file, &
	255	nys_on_file, nyn_on_file, nzb_on_file, nzt_on_file
	256
	257	IF ( nxl_on_file /= hor_index_bounds_previous_run(1,j) ) THEN
	258	PRINT*, '+++ read_3d_binary: problem with index bound nxl on ', &
	259	' restart file "', myid_char, '"'
	260	PRINT*, ' nxl = ', nxl_on_file, ' but it should be'
	261	PRINT*, ' = ', hor_index_bounds_previous_run(1,j)
	262	PRINT*, ' from the index bound information array'
[1]	263	#if defined( __parallel )
[146]	264	CALL MPI_ABORT( comm2d, 9999, ierr )
[1]	265	#else
[146]	266	CALL local_stop
[1]	267	#endif
[146]	268	ENDIF
[1]	269
[146]	270	IF ( nxr_on_file /= hor_index_bounds_previous_run(2,j) ) THEN
	271	PRINT*, '+++ read_3d_binary: problem with index bound nxr on ', &
	272	' restart file "', myid_char, '"'
	273	PRINT*, ' nxr = ', nxr_on_file, ' but it should be'
	274	PRINT*, ' = ', hor_index_bounds_previous_run(2,j)
	275	PRINT*, ' from the index bound information array'
[1]	276	#if defined( __parallel )
[146]	277	CALL MPI_ABORT( comm2d, 9999, ierr )
[1]	278	#else
[146]	279	CALL local_stop
[1]	280	#endif
[146]	281	ENDIF
[1]	282
[146]	283	IF ( nys_on_file /= hor_index_bounds_previous_run(3,j) ) THEN
	284	PRINT*, '+++ read_3d_binary: problem with index bound nys on ', &
	285	' restart file "', myid_char, '"'
	286	PRINT*, ' nys = ', nys_on_file, ' but it should be'
	287	PRINT*, ' = ', hor_index_bounds_previous_run(3,j)
	288	PRINT*, ' from the index bound information array'
[1]	289	#if defined( __parallel )
[146]	290	CALL MPI_ABORT( comm2d, 9999, ierr )
[1]	291	#else
[146]	292	CALL local_stop
[1]	293	#endif
[146]	294	ENDIF
[1]	295
[146]	296	IF ( nyn_on_file /= hor_index_bounds_previous_run(4,j) ) THEN
	297	PRINT*, '+++ read_3d_binary: problem with index bound nyn on ', &
	298	' restart file "', myid_char, '"'
	299	PRINT*, ' nyn = ', nyn_on_file, ' but it should be'
	300	PRINT*, ' = ', hor_index_bounds_previous_run(4,j)
	301	PRINT*, ' from the index bound information array'
[1]	302	#if defined( __parallel )
[146]	303	CALL MPI_ABORT( comm2d, 9999, ierr )
[1]	304	#else
[146]	305	CALL local_stop
[1]	306	#endif
[146]	307	ENDIF
[1]	308
[146]	309	IF ( nzb_on_file /= nzb ) THEN
	310	PRINT*, '+++ read_3d_binary: mismatch between actual data and data '
	311	PRINT*, ' from prior run on PE ', myid
	312	PRINT*, ' nzb on file = ', nzb_on_file
	313	PRINT*, ' nzb = ', nzb
	314	CALL local_stop
	315	ENDIF
[1]	316
[146]	317	IF ( nzt_on_file /= nzt ) THEN
	318	PRINT*, '+++ read_3d_binary: mismatch between actual data and data '
	319	PRINT*, ' from prior run on PE ', myid
	320	PRINT*, ' nzt on file = ', nzt_on_file
	321	PRINT*, ' nzt = ', nzt
	322	CALL local_stop
	323	ENDIF
[1]	324
	325	!
[146]	326	!-- Allocate temporary arrays sized as the arrays on the restart file
	327	ALLOCATE( tmp_2d(nys_on_file-1:nyn_on_file+1, &
	328	nxl_on_file-1:nxr_on_file+1), &
	329	tmp_3d(nzb:nzt+1,nys_on_file-1:nyn_on_file+1, &
	330	nxl_on_file-1:nxr_on_file+1) )
[1]	331
	332	!
[146]	333	!-- Read arrays
	334	!-- ATTENTION: If the following read commands have been altered, the
	335	!-- ---------- version number of the variable binary_version must be altered,
	336	!-- too. Furthermore, the output list of arrays in write_3d_binary
	337	!-- must also be altered accordingly.
	338	READ ( 13 ) field_chr
	339	DO WHILE ( TRIM( field_chr ) /= '* end *' )
[1]	340
[216]	341	!
	342	!-- Map data on file as often as needed (data are read only for k=1)
	343	DO k = 1, overlap_count(i)
[1]	344
[216]	345	!
	346	!-- Get the index range of the subdomain on file which overlap with the
	347	!-- current subdomain
	348	nxlf = nxlfa(i,k)
	349	nxlc = nxlfa(i,k) + offset_xa(i,k)
	350	nxrf = nxrfa(i,k)
	351	nxrc = nxrfa(i,k) + offset_xa(i,k)
	352	nysf = nysfa(i,k)
	353	nysc = nysfa(i,k) + offset_ya(i,k)
	354	nynf = nynfa(i,k)
	355	nync = nynfa(i,k) + offset_ya(i,k)
	356
	357
	358	WRITE (9,*) 'var = ', field_chr
	359	CALL local_flush( 9 )
	360
	361	SELECT CASE ( TRIM( field_chr ) )
	362
	363	CASE ( 'e' )
	364	IF ( k == 1 ) READ ( 13 ) tmp_3d
	365	e(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	366	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[1]	367
[216]	368	CASE ( 'e_av' )
	369	IF ( .NOT. ALLOCATED( e_av ) ) THEN
	370	ALLOCATE( e_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	371	ENDIF
	372	IF ( k == 1 ) READ ( 13 ) tmp_3d
	373	e_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	374	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	375
[216]	376	CASE ( 'e_m' )
	377	IF ( k == 1 ) READ ( 13 ) tmp_3d
	378	e_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	379	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	380
[216]	381	CASE ( 'iran' ) ! matching random numbers is still unresolved
	382	! issue
	383	IF ( k == 1 ) READ ( 13 ) iran, iran_part
[146]	384
[216]	385	CASE ( 'kh' )
	386	IF ( k == 1 ) READ ( 13 ) tmp_3d
	387	kh(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	388	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	389
[216]	390	CASE ( 'kh_m' )
	391	IF ( k == 1 ) READ ( 13 ) tmp_3d
	392	kh_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	393	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	394
[216]	395	CASE ( 'km' )
	396	IF ( k == 1 ) READ ( 13 ) tmp_3d
	397	km(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	398	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	399
[216]	400	CASE ( 'km_m' )
	401	IF ( k == 1 ) READ ( 13 ) tmp_3d
	402	km_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	403	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	404
[216]	405	CASE ( 'lwp_av' )
	406	IF ( .NOT. ALLOCATED( lwp_av ) ) THEN
	407	ALLOCATE( lwp_av(nys-1:nyn+1,nxl-1:nxr+1) )
	408	ENDIF
	409	IF ( k == 1 ) READ ( 13 ) tmp_2d
	410	lwp_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	411	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	412
[216]	413	CASE ( 'p' )
	414	IF ( k == 1 ) READ ( 13 ) tmp_3d
	415	p(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	416	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	417
[216]	418	CASE ( 'p_av' )
	419	IF ( .NOT. ALLOCATED( p_av ) ) THEN
	420	ALLOCATE( p_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	421	ENDIF
	422	IF ( k == 1 ) READ ( 13 ) tmp_3d
	423	p_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	424	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	425
[216]	426	CASE ( 'pc_av' )
	427	IF ( .NOT. ALLOCATED( pc_av ) ) THEN
	428	ALLOCATE( pc_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	429	ENDIF
	430	IF ( k == 1 ) READ ( 13 ) tmp_3d
	431	pc_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	432	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	433
[216]	434	CASE ( 'pr_av' )
	435	IF ( .NOT. ALLOCATED( pr_av ) ) THEN
	436	ALLOCATE( pr_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	437	ENDIF
	438	IF ( k == 1 ) READ ( 13 ) tmp_3d
	439	pr_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	440	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	441
[216]	442	CASE ( 'precipitation_amount' )
	443	IF ( k == 1 ) READ ( 13 ) tmp_2d
	444	precipitation_amount(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	445	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	446
[216]	447	CASE ( 'precipitation_rate_a' )
	448	IF ( .NOT. ALLOCATED( precipitation_rate_av ) ) THEN
	449	ALLOCATE( precipitation_rate_av(nys-1:nyn+1,nxl-1:nxr+1) )
	450	ENDIF
	451	IF ( k == 1 ) READ ( 13 ) tmp_2d
	452	precipitation_rate_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	453	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	454
[216]	455	CASE ( 'pt' )
	456	IF ( k == 1 ) READ ( 13 ) tmp_3d
	457	pt(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	458	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	459
[216]	460	CASE ( 'pt_av' )
	461	IF ( .NOT. ALLOCATED( pt_av ) ) THEN
	462	ALLOCATE( pt_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	463	ENDIF
	464	IF ( k == 1 ) READ ( 13 ) tmp_3d
	465	pt_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	466	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	467
[216]	468	CASE ( 'pt_m' )
	469	IF ( k == 1 ) READ ( 13 ) tmp_3d
	470	pt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	471	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	472
[216]	473	CASE ( 'q' )
	474	IF ( k == 1 ) READ ( 13 ) tmp_3d
	475	q(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	476	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	477
[216]	478	CASE ( 'q_av' )
	479	IF ( .NOT. ALLOCATED( q_av ) ) THEN
	480	ALLOCATE( q_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	481	ENDIF
	482	IF ( k == 1 ) READ ( 13 ) tmp_3d
	483	q_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	484	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	485
[216]	486	CASE ( 'q_m' )
	487	IF ( k == 1 ) READ ( 13 ) tmp_3d
	488	q_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	489	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	490
[216]	491	CASE ( 'ql' )
	492	IF ( k == 1 ) READ ( 13 ) tmp_3d
	493	ql(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	494	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	495
[216]	496	CASE ( 'ql_av' )
	497	IF ( .NOT. ALLOCATED( ql_av ) ) THEN
	498	ALLOCATE( ql_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	499	ENDIF
	500	IF ( k == 1 ) READ ( 13 ) tmp_3d
	501	ql_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	502	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	503
[216]	504	CASE ( 'ql_c_av' )
	505	IF ( .NOT. ALLOCATED( ql_c_av ) ) THEN
	506	ALLOCATE( ql_c_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	507	ENDIF
	508	IF ( k == 1 ) READ ( 13 ) tmp_3d
	509	ql_c_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	510	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	511
[216]	512	CASE ( 'ql_v_av' )
	513	IF ( .NOT. ALLOCATED( ql_v_av ) ) THEN
	514	ALLOCATE( ql_v_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	515	ENDIF
	516	IF ( k == 1 ) READ ( 13 ) tmp_3d
	517	ql_v_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	518	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	519
[216]	520	CASE ( 'ql_vp_av' )
	521	IF ( .NOT. ALLOCATED( ql_vp_av ) ) THEN
	522	ALLOCATE( ql_vp_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	523	ENDIF
	524	IF ( k == 1 ) READ ( 13 ) tmp_3d
	525	ql_vp_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	526	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	527
[216]	528	CASE ( 'qs' )
	529	IF ( k == 1 ) READ ( 13 ) tmp_2d
	530	qs(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	531	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	532
[216]	533	CASE ( 'qsws' )
	534	IF ( k == 1 ) READ ( 13 ) tmp_2d
	535	qsws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	536	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	537
[216]	538	CASE ( 'qsws_m' )
	539	IF ( k == 1 ) READ ( 13 ) tmp_2d
	540	qsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	541	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	542
[216]	543	CASE ( 'qswst' )
	544	IF ( k == 1 ) READ ( 13 ) tmp_2d
	545	qswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	546	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	547
[216]	548	CASE ( 'qswst_m' )
	549	IF ( k == 1 ) READ ( 13 ) tmp_2d
	550	qswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	551	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	552
[216]	553	CASE ( 'qv_av' )
	554	IF ( .NOT. ALLOCATED( qv_av ) ) THEN
	555	ALLOCATE( qv_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	556	ENDIF
	557	IF ( k == 1 ) READ ( 13 ) tmp_3d
	558	qv_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	559	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	560
[216]	561	CASE ( 'random_iv' ) ! still unresolved issue
	562	IF ( k == 1 ) READ ( 13 ) random_iv
	563	IF ( k == 1 ) READ ( 13 ) random_iy
[146]	564
[216]	565	CASE ( 'rho_av' )
	566	IF ( .NOT. ALLOCATED( rho_av ) ) THEN
	567	ALLOCATE( rho_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	568	ENDIF
	569	IF ( k == 1 ) READ ( 13 ) tmp_3d
	570	rho_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	571	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	572
[216]	573	CASE ( 'rif' )
	574	IF ( k == 1 ) READ ( 13 ) tmp_2d
	575	rif(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	576	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	577
[216]	578	CASE ( 'rif_m' )
	579	IF ( k == 1 ) READ ( 13 ) tmp_2d
	580	rif_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	581	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	582
[216]	583	CASE ( 'rif_wall' )
	584	IF ( k == 1 ) THEN
	585	ALLOCATE( tmp_4d(nzb:nzt+1,nys_on_file-1:nyn_on_file+1, &
	586	nxl_on_file-1:nxr_on_file+1,1:4) )
	587	READ ( 13 ) tmp_4d
	588	ENDIF
	589	rif_wall(:,nysc-1:nync+1,nxlc-1:nxrc+1,:) = &
[147]	590	tmp_4d(:,nysf-1:nynf+1,nxlf-1:nxrf+1,:)
[146]	591
[216]	592	CASE ( 's_av' )
	593	IF ( .NOT. ALLOCATED( s_av ) ) THEN
	594	ALLOCATE( s_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	595	ENDIF
	596	IF ( k == 1 ) READ ( 13 ) tmp_3d
	597	s_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	598	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	599
[216]	600	CASE ( 'sa' )
	601	IF ( k == 1 ) READ ( 13 ) tmp_3d
	602	sa(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	603	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	604
[216]	605	CASE ( 'sa_av' )
	606	IF ( .NOT. ALLOCATED( sa_av ) ) THEN
	607	ALLOCATE( sa_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	608	ENDIF
	609	IF ( k == 1 ) READ ( 13 ) tmp_3d
	610	sa_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	611	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	612
[216]	613	CASE ( 'saswsb' )
	614	IF ( k == 1 ) READ ( 13 ) tmp_2d
	615	saswsb(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	616	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	617
[216]	618	CASE ( 'saswst' )
	619	IF ( k == 1 ) READ ( 13 ) tmp_2d
	620	saswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	621	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	622
[216]	623	CASE ( 'shf' )
	624	IF ( k == 1 ) READ ( 13 ) tmp_2d
	625	shf(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	626	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	627
[216]	628	CASE ( 'shf_m' )
	629	IF ( k == 1 ) READ ( 13 ) tmp_2d
	630	shf_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	631	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	632
[216]	633	CASE ( 'spectrum_x' )
	634	IF ( k == 1 ) THEN
	635	IF ( nx_on_file /= nx ) THEN
	636	IF ( myid == 0 ) THEN
	637	PRINT*, '+++ WARNING: read_3d_binary: spectrum_x', &
	638	' on restart file ignored because'
	639	PRINT*, ' total numbers of grid points (nx) ', &
	640	'do not match'
	641	ENDIF
	642	READ ( 13 ) rdummy
	643	ELSE
	644	READ ( 13 ) spectrum_x
	645	ENDIF
	646	ENDIF
[146]	647
[216]	648	CASE ( 'spectrum_y' )
	649	IF ( k == 1 ) THEN
	650	IF ( ny_on_file /= ny ) THEN
	651	IF ( myid == 0 ) THEN
	652	PRINT*, '+++ WARNING: read_3d_binary: spectrum_y', &
	653	' on restart file ignored because'
	654	PRINT*, ' total numbers of grid points (ny) ', &
	655	'do not match'
	656	ENDIF
	657	READ ( 13 ) rdummy
	658	ELSE
	659	READ ( 13 ) spectrum_y
	660	ENDIF
	661	ENDIF
[146]	662
[216]	663	CASE ( 'ts' )
	664	IF ( k == 1 ) READ ( 13 ) tmp_2d
	665	ts(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	666	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	667
[216]	668	CASE ( 'ts_av' )
	669	IF ( .NOT. ALLOCATED( ts_av ) ) THEN
	670	ALLOCATE( ts_av(nys-1:nyn+1,nxl-1:nxr+1) )
	671	ENDIF
	672	IF ( k == 1 ) READ ( 13 ) tmp_2d
	673	ts_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	674	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	675
[216]	676	CASE ( 'tswst' )
	677	IF ( k == 1 ) READ ( 13 ) tmp_2d
	678	tswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	679	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	680
[216]	681	CASE ( 'tswst_m' )
	682	IF ( k == 1 ) READ ( 13 ) tmp_2d
	683	tswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	684	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	685
[216]	686	CASE ( 'u' )
	687	IF ( k == 1 ) READ ( 13 ) tmp_3d
	688	u(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	689	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	690
[216]	691	CASE ( 'u_av' )
	692	IF ( .NOT. ALLOCATED( u_av ) ) THEN
	693	ALLOCATE( u_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	694	ENDIF
	695	IF ( k == 1 ) READ ( 13 ) tmp_3d
	696	u_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	697	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	698
[216]	699	CASE ( 'u_m' )
	700	IF ( k == 1 ) READ ( 13 ) tmp_3d
	701	u_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	702	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	703
[216]	704	CASE ( 'u_m_l' )
	705	IF ( k == 1 ) THEN
	706	ALLOCATE( tmp_3dwul(nzb:nzt+1, &
	707	nys_on_file-1:nyn_on_file+1,1:2) )
	708	READ ( 13 ) tmp_3dwul
	709	ENDIF
	710	IF ( outflow_l ) THEN
	711	u_m_l(:,nysc-1:nync+1,:) = tmp_3dwul(:,nysf-1:nynf+1,:)
	712	ENDIF
[146]	713
[216]	714	CASE ( 'u_m_n' )
	715	IF ( k == 1 ) THEN
	716	ALLOCATE( tmp_3dwun(nzb:nzt+1,ny-1:ny, &
	717	nxl_on_file-1:nxr_on_file+1) )
	718	READ ( 13 ) tmp_3dwun
	719	ENDIF
	720	IF ( outflow_n ) THEN
	721	u_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dwun(:,:,nxlf-1:nxrf+1)
	722	ENDIF
[146]	723
[216]	724	CASE ( 'u_m_r' )
	725	IF ( k == 1 ) THEN
	726	ALLOCATE( tmp_3dwur(nzb:nzt+1,&
	727	nys_on_file-1:nyn_on_file+1,nx-1:nx) )
	728	READ ( 13 ) tmp_3dwur
	729	ENDIF
	730	IF ( outflow_r ) THEN
	731	u_m_r(:,nysc-1:nync+1,:) = tmp_3dwur(:,nysf-1:nynf+1,:)
	732	ENDIF
[146]	733
[216]	734	CASE ( 'u_m_s' )
	735	IF ( k == 1 ) THEN
	736	ALLOCATE( tmp_3dwus(nzb:nzt+1,0:1, &
	737	nxl_on_file-1:nxr_on_file+1) )
	738	READ ( 13 ) tmp_3dwus
	739	ENDIF
	740	IF ( outflow_s ) THEN
	741	u_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dwus(:,:,nxlf-1:nxrf+1)
	742	ENDIF
[146]	743
[216]	744	CASE ( 'us' )
	745	IF ( k == 1 ) READ ( 13 ) tmp_2d
	746	us(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	747	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	748
[216]	749	CASE ( 'usws' )
	750	IF ( k == 1 ) READ ( 13 ) tmp_2d
	751	usws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	752	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	753
[216]	754	CASE ( 'uswst' )
	755	IF ( k == 1 ) READ ( 13 ) tmp_2d
	756	uswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	757	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	758
[216]	759	CASE ( 'usws_m' )
	760	IF ( k == 1 ) READ ( 13 ) tmp_2d
	761	usws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	762	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	763
[216]	764	CASE ( 'uswst_m' )
	765	IF ( k == 1 ) READ ( 13 ) tmp_2d
	766	uswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	767	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	768
[216]	769	CASE ( 'us_av' )
	770	IF ( .NOT. ALLOCATED( us_av ) ) THEN
	771	ALLOCATE( us_av(nys-1:nyn+1,nxl-1:nxr+1) )
	772	ENDIF
	773	IF ( k == 1 ) READ ( 13 ) tmp_2d
	774	us_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	775	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	776
[216]	777	CASE ( 'v' )
	778	IF ( k == 1 ) READ ( 13 ) tmp_3d
	779	v(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	780	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	781
[216]	782	CASE ( 'v_av' )
	783	IF ( .NOT. ALLOCATED( v_av ) ) THEN
	784	ALLOCATE( v_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	785	ENDIF
	786	IF ( k == 1 ) READ ( 13 ) tmp_3d
	787	v_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	788	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	789
[216]	790	CASE ( 'v_m' )
	791	IF ( k == 1 ) READ ( 13 ) tmp_3d
	792	v_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	793	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	794
[216]	795	CASE ( 'v_m_l' )
	796	IF ( k == 1 ) THEN
	797	ALLOCATE( tmp_3dwvl(nzb:nzt+1,&
	798	nys_on_file-1:nyn_on_file+1,0:1) )
	799	READ ( 13 ) tmp_3dwvl
	800	ENDIF
	801	IF ( outflow_l ) THEN
	802	v_m_l(:,nysc-1:nync+1,:) = tmp_3dwvl(:,nysf-1:nynf+1,:)
	803	ENDIF
[146]	804
[216]	805	CASE ( 'v_m_n' )
	806	IF ( k == 1 ) THEN
	807	ALLOCATE( tmp_3dwvn(nzb:nzt+1,ny-1:ny, &
	808	nxl_on_file-1:nxr_on_file+1) )
	809	READ ( 13 ) tmp_3dwvn
	810	ENDIF
	811	IF ( outflow_n ) THEN
	812	v_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dwvn(:,:,nxlf-1:nxrf+1)
	813	ENDIF
[146]	814
[216]	815	CASE ( 'v_m_r' )
	816	IF ( k == 1 ) THEN
	817	ALLOCATE( tmp_3dwvr(nzb:nzt+1,&
	818	nys_on_file-1:nyn_on_file+1,nx-1:nx) )
	819	READ ( 13 ) tmp_3dwvr
	820	ENDIF
	821	IF ( outflow_r ) THEN
	822	v_m_r(:,nysc-1:nync+1,:) = tmp_3dwvr(:,nysf-1:nynf+1,:)
	823	ENDIF
[146]	824
[216]	825	CASE ( 'v_m_s' )
	826	IF ( k == 1 ) THEN
	827	ALLOCATE( tmp_3dwvs(nzb:nzt+1,1:2, &
	828	nxl_on_file-1:nxr_on_file+1) )
	829	READ ( 13 ) tmp_3dwvs
	830	ENDIF
	831	IF ( outflow_s ) THEN
	832	v_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dwvs(:,:,nxlf-1:nxrf+1)
	833	ENDIF
[146]	834
[216]	835	CASE ( 'vpt' )
	836	IF ( k == 1 ) READ ( 13 ) tmp_3d
	837	vpt(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	838	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	839
[216]	840	CASE ( 'vpt_av' )
	841	IF ( .NOT. ALLOCATED( vpt_av ) ) THEN
	842	ALLOCATE( vpt_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	843	ENDIF
	844	IF ( k == 1 ) READ ( 13 ) tmp_3d
	845	vpt_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	846	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	847
[216]	848	CASE ( 'vpt_m' )
	849	IF ( k == 1 ) READ ( 13 ) tmp_3d
	850	vpt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	851	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	852
[216]	853	CASE ( 'vsws' )
	854	IF ( k == 1 ) READ ( 13 ) tmp_2d
	855	vsws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	856	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	857
[216]	858	CASE ( 'vswst' )
	859	IF ( k == 1 ) READ ( 13 ) tmp_2d
	860	vswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	861	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	862
[216]	863	CASE ( 'vsws_m' )
	864	IF ( k == 1 ) READ ( 13 ) tmp_2d
	865	vsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	866	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	867
[216]	868	CASE ( 'vswst_m' )
	869	IF ( k == 1 ) READ ( 13 ) tmp_2d
	870	vswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	871	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	872
[216]	873	CASE ( 'w' )
	874	IF ( k == 1 ) READ ( 13 ) tmp_3d
	875	w(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	876	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	877
[216]	878	CASE ( 'w_av' )
	879	IF ( .NOT. ALLOCATED( w_av ) ) THEN
	880	ALLOCATE( w_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	881	ENDIF
	882	IF ( k == 1 ) READ ( 13 ) tmp_3d
	883	w_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	884	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	885
[216]	886	CASE ( 'w_m' )
	887	IF ( k == 1 ) READ ( 13 ) tmp_3d
	888	w_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	889	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	890
[216]	891	CASE ( 'w_m_l' )
	892	IF ( k == 1 ) THEN
	893	ALLOCATE( tmp_3dwwl(nzb:nzt+1,&
	894	nys_on_file-1:nyn_on_file+1,0:1) )
	895	READ ( 13 ) tmp_3dwwl
	896	ENDIF
	897	IF ( outflow_l ) THEN
	898	w_m_l(:,nysc-1:nync+1,:) = tmp_3dwwl(:,nysf-1:nynf+1,:)
	899	ENDIF
[146]	900
[216]	901	CASE ( 'w_m_n' )
	902	IF ( k == 1 ) THEN
	903	ALLOCATE( tmp_3dwwn(nzb:nzt+1,ny-1:ny, &
	904	nxl_on_file-1:nxr_on_file+1) )
	905	READ ( 13 ) tmp_3dwwn
	906	ENDIF
	907	IF ( outflow_n ) THEN
	908	w_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dwwn(:,:,nxlf-1:nxrf+1)
	909	ENDIF
[146]	910
[216]	911	CASE ( 'w_m_r' )
	912	IF ( k == 1 ) THEN
	913	ALLOCATE( tmp_3dwwr(nzb:nzt+1,&
	914	nys_on_file-1:nyn_on_file+1,nx-1:nx) )
	915	READ ( 13 ) tmp_3dwwr
	916	ENDIF
	917	IF ( outflow_r ) THEN
	918	w_m_r(:,nysc-1:nync+1,:) = tmp_3dwwr(:,nysf-1:nynf+1,:)
	919	ENDIF
[146]	920
[216]	921	CASE ( 'w_m_s' )
	922	IF ( k == 1 ) THEN
	923	ALLOCATE( tmp_3dwws(nzb:nzt+1,0:1, &
	924	nxl_on_file-1:nxr_on_file+1) )
	925	READ ( 13 ) tmp_3dwws
	926	ENDIF
	927	IF ( outflow_s ) THEN
	928	w_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dwws(:,:,nxlf-1:nxrf+1)
	929	ENDIF
	930	DEALLOCATE( tmp_3dwws )
[146]	931
[216]	932	CASE ( 'z0' )
	933	IF ( k == 1 ) READ ( 13 ) tmp_2d
	934	z0(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	935	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	936
[216]	937	CASE ( 'z0_av' )
	938	IF ( .NOT. ALLOCATED( z0_av ) ) THEN
	939	ALLOCATE( z0_av(nys-1:nyn+1,nxl-1:nxr+1) )
	940	ENDIF
	941	IF ( k == 1 ) READ ( 13 ) tmp_2d
	942	z0_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
[147]	943	tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
[146]	944
[216]	945	CASE DEFAULT
	946	PRINT*, '+++ read_3d_binary: unknown field named "', &
	947	TRIM( field_chr ), '" found in'
	948	PRINT*, ' data from prior run on PE ',myid
	949	CALL local_stop
[146]	950
[216]	951	END SELECT
	952
	953	ENDDO ! overlap loop
	954
[1]	955	!
[216]	956	!-- Deallocate arrays needed for specific variables only
	957	IF ( ALLOCATED( tmp_3dwul ) ) DEALLOCATE( tmp_3dwul )
	958	IF ( ALLOCATED( tmp_3dwun ) ) DEALLOCATE( tmp_3dwun )
	959	IF ( ALLOCATED( tmp_3dwur ) ) DEALLOCATE( tmp_3dwur )
	960	IF ( ALLOCATED( tmp_3dwus ) ) DEALLOCATE( tmp_3dwus )
	961	IF ( ALLOCATED( tmp_3dwvl ) ) DEALLOCATE( tmp_3dwvl )
	962	IF ( ALLOCATED( tmp_3dwvn ) ) DEALLOCATE( tmp_3dwvn )
	963	IF ( ALLOCATED( tmp_3dwvr ) ) DEALLOCATE( tmp_3dwvr )
	964	IF ( ALLOCATED( tmp_3dwvs ) ) DEALLOCATE( tmp_3dwvs )
	965	IF ( ALLOCATED( tmp_3dwwl ) ) DEALLOCATE( tmp_3dwwl )
	966	IF ( ALLOCATED( tmp_3dwwn ) ) DEALLOCATE( tmp_3dwwn )
	967	IF ( ALLOCATED( tmp_3dwwr ) ) DEALLOCATE( tmp_3dwwr )
	968	IF ( ALLOCATED( tmp_3dwws ) ) DEALLOCATE( tmp_3dwws )
	969	IF ( ALLOCATED( tmp_4d ) ) DEALLOCATE( tmp_4d )
	970
	971	!
[146]	972	!-- Read next character string
	973	READ ( 13 ) field_chr
[1]	974
[146]	975	ENDDO ! loop over variables
[1]	976
	977	!
[146]	978	!-- Read user-defined restart data
[220]	979	CALL user_read_restart_data( i, nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
	980	nynfa, nyn_on_file, nysfa, nys_on_file, &
	981	offset_xa, offset_ya, overlap_count(i), &
[216]	982	tmp_2d, tmp_3d )
[145]	983
	984	!
[146]	985	!-- Close the restart file
	986	CALL close_file( 13 )
	987
	988	DEALLOCATE( tmp_2d, tmp_3d )
	989
	990	ENDDO ! loop over restart files
	991
	992
	993	!
	994	!-- Restore the original filename for the restart file to be written
	995	myid_char = myid_char_save
	996
	997
	998	!
[1]	999	!-- End of time measuring for reading binary data
	1000	CALL cpu_log( log_point_s(14), 'read_3d_binary', 'stop' )
	1001
	1002	END SUBROUTINE read_3d_binary

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |