Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |