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

Last change on this file since 207 was 198, checked in by raasch, 16 years ago

file headers updated for the next release 3.5

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