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

Last change on this file since 220 was 220, checked in by raasch, 13 years ago

some small bugfixes in user_module, user_read_restart_data, read_3d_binary, flow_statistics and mrun

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