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

Last change on this file since 449 was 449, checked in by raasch, 12 years ago

branch revision comments from Marcus (rev 410) replaced by normal revision comments

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