Home

Context Navigation

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

Last change on this file since 2849 was 2817, checked in by knoop, 7 years ago
Preliminary gust module interface implemented
Property svn:keywords set to `Id`
File size: 51.3 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |