Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |