Home

Context Navigation

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

Last change on this file since 1274 was 1116, checked in by hoffmann, 12 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 45.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |