Home

Context Navigation

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

Last change on this file since 2284 was 2269, checked in by suehring, 8 years ago
Enable restarts with USM with different number of PEs; some bugfixes in new surface structure in USM; formatting adjustments and descriptions in surface_mod
Property svn:keywords set to `Id`
File size: 45.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |