Changeset 216 for palm


Ignore:
Timestamp:
Nov 25, 2008 7:12:43 AM (16 years ago)
Author:
raasch
Message:

reading mechanism for restart files completely revised

Location:
palm/trunk/SOURCE
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r215 r216  
    2121User interface has been split into single files.
    2222
     23Reading mechanism from restart files completely revised. The subdomain/total
     24domain size can vary arbitrarily between current and previous run.
     25(read_var_list, read_3d_binary, user_read_restart_data)
     26
    2327Precompilation mechanism (mbuild) completely revised: the source tarfile remains in
    2428the source directory. One depository per block (given in the config-file) is
     
    3236by this routine step by step within the next revisions.
    3337
    34 Makefile, mbuild, mrun, check_parameters, local_stop, modules
     38Neumann boundary condition at k=nzb is explicitly set for better reading,
     39although this has been already done in boundary_conds. (advec_s_bc)
     40
     41Origin of the xy-coordinate system shifted from the center of the first
     42grid cell (indices i=0, j=0) to the south-left corner of this cell. (netcdf)
     43
     44Makefile, mbuild, mrun, advec_s_bc, check_parameters, local_stop, modules, netcdf, read_var_list, read_3d_binary, user_read_restart_data
    3545
    3646New: handle_palm_message, user_3d_data_averaging user_actions user_advec_particles user_check_data_output user_check_data_output_pr user_data_output_2d user_data_output_3d user_data_output_dvrp user_define_netcdf_grid user_dvrp_coltab user_header user_init user_init_3d_model user_init_grid user_init_particles user_init_plant_canopy user_last_actions user_module user_parin user_particle_attributes user_read_restart_data user_spectra user_statistics
    3747
    3848Removed: user_interface
     49
    3950
    4051Errors:
     
    5768have been started (data_output_2d)
    5869
     70Bugfix: reading of spectrum_x|y from restart files ignored if total numbers
     71of grid points do not match (read_3d_binary)
    5972
    60 advec_particles, combine_plot_fields, data_output_2d, header, production_e
     73advec_particles, combine_plot_fields, data_output_2d, header, production_e, read_3d_binary
  • palm/trunk/SOURCE/advec_s_bc.f90

    r98 r216  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Neumann boundary condition at k=nzb is explicitly set for better reading,
     7! although this has been already done in boundary_conds
    78!
    89! Former revisions:
     
    728729          DO  i = nxl, nxr
    729730             DO  j = nys, nyn
     731                sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
    730732                sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
    731733                sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
     
    766768       DO  i = nxl, nxr
    767769          DO  j = nys, nyn
     770             sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
    768771             sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
    769772             sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
     
    788791       DO  i = nxl, nxr
    789792          DO  j = nys, nyn
     793             sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
    790794             sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
    791795             sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
     
    823827       DO  i = nxl, nxr
    824828          DO  j = nys, nyn
     829             sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
    825830             sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
    826831             sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
  • palm/trunk/SOURCE/modules.f90

    r213 r216  
    301301    CHARACTER (LEN=10), DIMENSION(0:1,100) ::  do2d = ' ', do3d = ' '
    302302
    303     INTEGER ::  abort_mode = 0, average_count_pr = 0, average_count_sp = 0, &
     303    INTEGER ::  abort_mode = 1, average_count_pr = 0, average_count_sp = 0, &
    304304                average_count_3d = 0, current_timestep_number = 0, &
    305305                dist_range = 0, disturbance_level_ind_b, &
  • palm/trunk/SOURCE/netcdf.f90

    r198 r216  
    77! Current revisions:
    88! ------------------
    9 !
     9! Origin of the xy-coordinate system shifted from the center of the first
     10! grid cell (indices i=0, j=0) to the south-left corner of this cell.
    1011!
    1112! Former revisions:
     
    457458
    458459!
    459 !--       Write data for x and xu axis (shifted by -dx/2)
     460!--       Write data for x (shifted by +dx/2) and xu axis
    460461          ALLOCATE( netcdf_data(0:nx+1) )
    461462
    462463          DO  i = 0, nx+1
    463              netcdf_data(i) = i * dx
     464             netcdf_data(i) = ( i + 0.5 ) * dx
    464465          ENDDO
    465466
     
    469470
    470471          DO  i = 0, nx+1
    471              netcdf_data(i) = ( i - 0.5 ) * dx
     472             netcdf_data(i) = i * dx
    472473          ENDDO
    473474
     
    480481
    481482!
    482 !--       Write data for y and yv axis (shifted by -dy/2)
     483!--       Write data for y (shifted by +dy/2) and yv axis
    483484          ALLOCATE( netcdf_data(0:ny+1) )
    484485
    485486          DO  i = 0, ny+1
    486              netcdf_data(i) = i * dy
     487             netcdf_data(i) = ( i + 0.5 ) * dy
    487488          ENDDO
    488489
     
    492493
    493494          DO  i = 0, ny+1
    494              netcdf_data(i) = ( i - 0.5 ) * dy
     495             netcdf_data(i) = i * dy
    495496          ENDDO
    496497
     
    10341035
    10351036!
    1036 !--       Write data for x and xu axis (shifted by -dx/2)
     1037!--       Write data for x (shifted by +dx/2) and xu axis
    10371038          ALLOCATE( netcdf_data(0:nx+1) )
    10381039
    10391040          DO  i = 0, nx+1
    1040              netcdf_data(i) = i * dx
     1041             netcdf_data(i) = ( i + 0.5 ) * dx
    10411042          ENDDO
    10421043
     
    10461047
    10471048          DO  i = 0, nx+1
    1048              netcdf_data(i) = ( i - 0.5 ) * dx
     1049             netcdf_data(i) = i * dx
    10491050          ENDDO
    10501051
     
    10571058
    10581059!
    1059 !--       Write data for y and yv axis (shifted by -dy/2)
     1060!--       Write data for y (shifted by +dy/2) and yv axis
    10601061          ALLOCATE( netcdf_data(0:ny+1) )
    10611062
    10621063          DO  i = 0, ny+1
    1063              netcdf_data(i) = i * dy
     1064             netcdf_data(i) = ( i + 0.5 ) * dy
    10641065          ENDDO
    10651066
     
    10691070
    10701071          DO  i = 0, ny+1
    1071              netcdf_data(i) = ( i - 0.5 ) * dy
     1072             netcdf_data(i) = i * dy
    10721073          ENDDO
    10731074
     
    15311532
    15321533!
    1533 !--       Write y_xz data
     1534!--       Write y_xz data (shifted by +dy/2)
     1535          DO  i = 1, ns
     1536             IF( section(i,2) == -1 )  THEN
     1537                netcdf_data(i) = -1.0  ! section averaged along y
     1538             ELSE
     1539                netcdf_data(i) = ( section(i,2) + 0.5 ) * dy
     1540             ENDIF
     1541          ENDDO
     1542          nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_y_xz(av), netcdf_data, &
     1543                                  start = (/ 1 /), count = (/ ns /) )
     1544          IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 163 )
     1545
     1546!
     1547!--       Write yv_xz data
    15341548          DO  i = 1, ns
    15351549             IF( section(i,2) == -1 )  THEN
     
    15391553             ENDIF
    15401554          ENDDO
    1541           nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_y_xz(av), netcdf_data, &
    1542                                   start = (/ 1 /), count = (/ ns /) )
    1543           IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 163 )
    1544 
    1545 !
    1546 !--       Write yv_xz data
    1547           DO  i = 1, ns
    1548              IF( section(i,2) == -1 )  THEN
    1549                 netcdf_data(i) = -1.0  ! section averaged along y
    1550              ELSE
    1551                 netcdf_data(i) = ( section(i,2) - 0.5 ) * dy
    1552              ENDIF
    1553           ENDDO
    15541555          nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_yv_xz(av), &
    15551556                                  netcdf_data, start = (/ 1 /),    &
     
    15681569
    15691570!
    1570 !--       Write data for x and xu axis (shifted by -dx/2)
     1571!--       Write data for x (shifted by +dx/2) and xu axis
    15711572          ALLOCATE( netcdf_data(0:nx+1) )
    15721573
    15731574          DO  i = 0, nx+1
    1574              netcdf_data(i) = i * dx
     1575             netcdf_data(i) = ( i + 0.5 ) * dx
    15751576          ENDDO
    15761577
     
    15801581
    15811582          DO  i = 0, nx+1
    1582              netcdf_data(i) = ( i - 0.5 ) * dx
     1583             netcdf_data(i) = i * dx
    15831584          ENDDO
    15841585
     
    20422043
    20432044!
    2044 !--       Write x_yz data
     2045!--       Write x_yz data (shifted by +dx/2)
     2046          DO  i = 1, ns
     2047             IF( section(i,3) == -1 )  THEN
     2048                netcdf_data(i) = -1.0  ! section averaged along x
     2049             ELSE
     2050                netcdf_data(i) = ( section(i,3) + 0.5 ) * dx
     2051             ENDIF
     2052          ENDDO
     2053          nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_x_yz(av), netcdf_data, &
     2054                                  start = (/ 1 /), count = (/ ns /) )
     2055          IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 202 )
     2056
     2057!
     2058!--       Write x_yz data (xu grid)
    20452059          DO  i = 1, ns
    20462060             IF( section(i,3) == -1 )  THEN
     
    20502064             ENDIF
    20512065          ENDDO
    2052           nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_x_yz(av), netcdf_data, &
    2053                                   start = (/ 1 /), count = (/ ns /) )
    2054           IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 202 )
    2055 
    2056 !
    2057 !--       Write x_yz data (xu grid)
    2058           DO  i = 1, ns
    2059              IF( section(i,3) == -1 )  THEN
    2060                 netcdf_data(i) = -1.0  ! section averaged along x
    2061              ELSE
    2062                 netcdf_data(i) = (section(i,3)-0.5) * dx
    2063              ENDIF
    2064           ENDDO
    20652066          nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_xu_yz(av), netcdf_data, &
    20662067                                  start = (/ 1 /), count = (/ ns /) )
     
    20782079
    20792080!
    2080 !--       Write data for y and yv axis (shifted by -dy/2)
     2081!--       Write data for y (shifted by +dy/2) and yv axis
    20812082          ALLOCATE( netcdf_data(0:ny+1) )
    20822083
    20832084          DO  j = 0, ny+1
    2084              netcdf_data(j) = j * dy
     2085             netcdf_data(j) = ( j + 0.5 ) * dy
    20852086          ENDDO
    20862087
     
    20902091
    20912092          DO  j = 0, ny+1
    2092              netcdf_data(j) = ( j - 0.5 ) * dy
     2093             netcdf_data(j) = j * dy
    20932094          ENDDO
    20942095
  • palm/trunk/SOURCE/read_3d_binary.f90

    r198 r216  
    44! Actual revisions:
    55! -----------------
    6 !
     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
    710!
    811! Former revisions:
     
    6669    CHARACTER (LEN=20) ::  field_chr
    6770
    68     INTEGER ::  files_to_be_opened, i, j, myid_on_file,                       &
     71    INTEGER ::  files_to_be_opened, i, j, k, myid_on_file,                    &
    6972                numprocs_on_file, nxlc, nxlf, nxlpr, nxl_on_file, nxrc, nxrf, &
    7073                nxrpr, nxr_on_file, nync, nynf, nynpr, nyn_on_file, nysc,     &
    7174                nysf, nyspr, nys_on_file, nzb_on_file, nzt_on_file, offset_x, &
    72                 offset_y
    73 
    74     INTEGER, DIMENSION(numprocs_previous_run*4) ::  file_list, nxlfa, nxrfa, &
    75                 nynfa, nysfa, offset_xa, offset_ya
     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
    7683
    7784    REAL, DIMENSION(:,:), ALLOCATABLE     ::  tmp_2d
    78     REAL, DIMENSION(:,:,:), ALLOCATABLE   ::  tmp_3d, tmp_3dw
     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
    7990    REAL, DIMENSION(:,:,:,:), ALLOCATABLE ::  tmp_4d
    8091
     
    90101
    91102    DO  i = 1, numprocs_previous_run
    92 
     103!
     104!--    Store array bounds of the previous run ("pr") in temporary scalars
    93105       nxlpr = hor_index_bounds_previous_run(1,i-1)
    94106       nxrpr = hor_index_bounds_previous_run(2,i-1)
     
    103115
    104116!
    105 !--    Only data which overlap with the current subdomain have to be read
    106        IF ( nxlpr <= nxr-offset_x  .AND.  nxrpr >= nxl-offset_x  .AND.  &
    107             nyspr <= nyn-offset_y  .AND.  nynpr >= nys-offset_y )  THEN
    108 
    109           files_to_be_opened = files_to_be_opened + 1
    110           file_list(files_to_be_opened) = i-1
    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 )
    117 
    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)
     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
    130186          CALL local_flush( 9 )
    131           offset_xa(files_to_be_opened) = offset_x
    132           offset_ya(files_to_be_opened) = offset_y
    133 
    134        ENDIF
    135 
    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 
    262187       ENDIF
    263188
     
    403328
    404329!
    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)
    415 
    416 !
    417330!--    Read arrays
    418331!--    ATTENTION: If the following read commands have been altered, the
     
    423336       DO  WHILE ( TRIM( field_chr ) /= '*** end ***' )
    424337
    425           WRITE (9,*) 'var = ', field_chr
    426           CALL local_flush( 9 )
    427           SELECT CASE ( TRIM( field_chr ) )
    428 
    429              CASE ( 'e' )
    430                 READ ( 13 )  tmp_3d
    431                 e(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    432                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    433 
    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) = &
    440                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    441 
    442              CASE ( 'e_m' )
    443                 READ ( 13 )  tmp_3d
    444                 e_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    445                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    453                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    454 
    455              CASE ( 'kh_m' )
    456                 READ ( 13 )  tmp_3d
    457                 kh_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    458                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    459 
    460              CASE ( 'km' )
    461                 READ ( 13 )  tmp_3d
    462                 km(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    463                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    464 
    465              CASE ( 'km_m' )
    466                 READ ( 13 )  tmp_3d
    467                 km_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    468                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    476                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    477 
    478              CASE ( 'p' )
    479                 READ ( 13 )  tmp_3d
    480                 p(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    481                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    489                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    497                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    505                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    506 
    507              CASE ( 'precipitation_amount' )
    508                 READ ( 13 )  tmp_2d
    509                 precipitation_amount(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    510                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    518                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    519 
    520              CASE ( 'pt' )
    521                 READ ( 13 )  tmp_3d
    522                 pt(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    523                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    531                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    532 
    533              CASE ( 'pt_m' )
    534                 READ ( 13 )  tmp_3d
    535                 pt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    536                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    537 
    538              CASE ( 'q' )
    539                 READ ( 13 )  tmp_3d
    540                 q(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    541                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    549                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    550 
    551              CASE ( 'q_m' )
    552                 READ ( 13 )  tmp_3d
    553                 q_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    554                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    555 
    556              CASE ( 'ql' )
    557                 READ ( 13 )  tmp_3d
    558                 ql(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    559                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    567                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    575                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    583                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    591                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    592 
    593              CASE ( 'qs' )
    594                 READ ( 13 )  tmp_2d
    595                 qs(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    596                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    597 
    598              CASE ( 'qsws' )
    599                 READ ( 13 )  tmp_2d
    600                 qsws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    601                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    602 
    603              CASE ( 'qsws_m' )
    604                 READ ( 13 )  tmp_2d
    605                 qsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    606                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    607 
    608              CASE ( 'qswst' )
    609                 READ ( 13 )  tmp_2d
    610                 qswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    611                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    612 
    613              CASE ( 'qswst_m' )
    614                 READ ( 13 )  tmp_2d
    615                 qswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    616                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    624                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    636                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    637 
    638              CASE ( 'rif' )
    639                 READ ( 13 )  tmp_2d
    640                 rif(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    641                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    642 
    643              CASE ( 'rif_m' )
    644                 READ ( 13 )  tmp_2d
    645                 rif_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    646                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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,:) = &
     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,:) = &
    653587                                         tmp_4d(:,nysf-1:nynf+1,nxlf-1:nxrf+1,:)
    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) = &
    662                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    663 
    664              CASE ( 'sa' )
    665                 READ ( 13 )  tmp_3d
    666                 sa(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    667                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    675                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    676 
    677              CASE ( 'saswsb' )
    678                 READ ( 13 )  tmp_2d
    679                 saswsb(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    680                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    681 
    682              CASE ( 'saswst' )
    683                 READ ( 13 )  tmp_2d
    684                 saswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    685                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    686 
    687              CASE ( 'shf' )
    688                 READ ( 13 )  tmp_2d
    689                 shf(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    690                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    691 
    692              CASE ( 'shf_m' )
    693                 READ ( 13 )  tmp_2d
    694                 shf_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    695                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    706                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    714                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    715 
    716              CASE ( 'tswst' )
    717                 READ ( 13 )  tmp_2d
    718                 tswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    719                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    720 
    721              CASE ( 'tswst_m' )
    722                 READ ( 13 )  tmp_2d
    723                 tswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    724                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    725 
    726              CASE ( 'u' )
    727                 READ ( 13 )  tmp_3d
    728                 u(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    729                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    737                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    738 
    739              CASE ( 'u_m' )
    740                 READ ( 13 )  tmp_3d
    741                 u_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    742                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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
    748                    u_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    757                    u_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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
    766                    u_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    775                    u_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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) = &
    782                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    783 
    784              CASE ( 'usws' )
    785                 READ ( 13 )  tmp_2d
    786                 usws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    787                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    788 
    789              CASE ( 'uswst' )
    790                 READ ( 13 )  tmp_2d
    791                 uswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    792                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    793 
    794              CASE ( 'usws_m' )
    795                 READ ( 13 )  tmp_2d
    796                 usws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    797                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    798 
    799              CASE ( 'uswst_m' )
    800                 READ ( 13 )  tmp_2d
    801                 uswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    802                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    810                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    811 
    812              CASE ( 'v' )
    813                 READ ( 13 )  tmp_3d
    814                 v(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    815                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    823                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    824 
    825              CASE ( 'v_m' )
    826                 READ ( 13 )  tmp_3d
    827                 v_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    828                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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
    834                    v_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    843                    v_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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
    852                    v_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    861                    v_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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) = &
    868                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    876                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    877 
    878              CASE ( 'vpt_m' )
    879                 READ ( 13 )  tmp_3d
    880                 vpt_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    881                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    882 
    883              CASE ( 'vsws' )
    884                 READ ( 13 )  tmp_2d
    885                 vsws(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    886                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    887 
    888              CASE ( 'vswst' )
    889                 READ ( 13 )  tmp_2d
    890                 vswst(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    891                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    892 
    893              CASE ( 'vsws_m' )
    894                 READ ( 13 )  tmp_2d
    895                 vsws_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    896                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    897 
    898              CASE ( 'vswst_m' )
    899                 READ ( 13 )  tmp_2d
    900                 vswst_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
    901                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    902 
    903              CASE ( 'w' )
    904                 READ ( 13 )  tmp_3d
    905                 w(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    906                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    914                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    915 
    916              CASE ( 'w_m' )
    917                 READ ( 13 )  tmp_3d
    918                 w_m(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    919                                           tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    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
    925                    w_m_l(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    934                    w_m_n(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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
    943                    w_m_r(:,nysc-1:nync+1,:) = tmp_3dw(:,nysf-1:nynf+1,:)
    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
    952                    w_m_s(:,:,nxlc-1:nxrc+1) = tmp_3dw(:,:,nxlf-1:nxrf+1)
    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) = &
    959                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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) = &
    967                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
    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
     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
    976968!
    977969!--       Read next character string
     
    982974!
    983975!--    Read user-defined restart data
    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 )
     976       CALL user_read_restart_data( 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 )
    987980
    988981!
  • palm/trunk/SOURCE/read_var_list.f90

    r198 r216  
    44! Actual revisions:
    55! -----------------
    6 !
     6! limitations for nx_on_file, ny_on_file removed (read_parts_of_var_list)
    77!
    88! Former revisions:
     
    725725          CASE ( 'nx' )
    726726             READ ( 13 )  nx_on_file
    727              IF ( nx_on_file < ( nxr - nxl ) )  THEN
    728                 PRINT*, '+++ read_parts_of_var_list: total domain along x on', &
    729                              ' restart file is smaller than current subdomain'
    730                 PRINT*, '                            nx on file = ', nx_on_file
    731                 PRINT*, '                            nxr - nxl  = ', nxr - nxl
    732                 CALL local_stop
    733              ENDIF
    734727
    735728          CASE ( 'ny' )
    736729             READ ( 13 )  ny_on_file
    737              IF ( ny_on_file < ( nyn - nys ) )  THEN
    738                 PRINT*, '+++ read_parts_of_var_list: total domain along y on', &
    739                              ' restart file is smaller than current subdomain'
    740                 PRINT*, '                            ny on file = ', ny_on_file
    741                 PRINT*, '                            nyn - nys  = ', nyn - nys
    742                 CALL local_stop
    743              ENDIF
    744730
    745731
  • palm/trunk/SOURCE/user_read_restart_data.f90

    r211 r216  
    1  SUBROUTINE user_read_restart_data( nxlc, nxlf, nxl_on_file, nxrc, nxrf,  &
    2                                     nxr_on_file, nync, nynf, nyn_on_file, &
    3                                     nysc, nysf, nys_on_file, tmp_2d, tmp_3d )
     1 SUBROUTINE user_read_restart_data( nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
     2                                    nynfa, nyn_on_file, nysfa, nys_on_file, &
     3                                    offset_xa, offset_ya, overlap_count,    &
     4                                    tmp_2d, tmp_3d )
    45
    56!------------------------------------------------------------------------------!
    67! Actual revisions:
    78! -----------------
    8 ! Former file user_interface.f90 split into one file per subroutine
     9! reading mechanism revised (subdomain/total domain size can vary arbitrarily
     10! between current and previous run),
     11! former file user_interface.f90 split into one file per subroutine
    912!
    1013! Former revisions:
     
    3033    CHARACTER (LEN=20) :: field_char
    3134
    32     INTEGER ::  nxlc, nxlf, nxl_on_file, nxrc, nxrf, nxr_on_file, nync, nynf, &
    33                 nyn_on_file, nysc, nysf, nys_on_file
     35    INTEGER ::  k, nxlc, nxlf, nxl_on_file, nxrc, nxrf, nxr_on_file, nync, &
     36                nynf, nyn_on_file, nysc, nysf, nys_on_file, overlap_count
     37
     38    INTEGER, DIMENSION(numprocs_previous_run,1000) ::  nxlfa, nxrfa, nynfa, &
     39                                                       nysfa, offset_xa, &
     40                                                       offset_ya
    3441
    3542    REAL, DIMENSION(nys_on_file-1:nyn_on_file+1,nxl_on_file-1:nxr_on_file+1) ::&
     
    4855!       DO  WHILE ( TRIM( field_char ) /= '*** end user ***' )
    4956!
    50 !          SELECT CASE ( TRIM( field_char ) )
     57!          DO  k = 1, overlap_count
    5158!
    52 !             CASE ( 'u2_av' )
    53 !                IF ( .NOT. ALLOCATED( u2_av ) ) THEN
    54 !                   ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    55 !                ENDIF
    56 !                READ ( 13 )  tmp_3d
    57 !                u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
     59!             nxlf = nxlfa(i,k)
     60!             nxlc = nxlfa(i,k) + offset_xa(i,k)
     61!             nxrf = nxrfa(i,k)
     62!             nxrc = nxrfa(i,k) + offset_xa(i,k)
     63!             nysf = nysfa(i,k)
     64!             nysc = nysfa(i,k) + offset_ya(i,k)
     65!             nynf = nynfa(i,k)
     66!             nync = nynfa(i,k) + offset_ya(i,k)
     67!
     68!
     69!             SELECT CASE ( TRIM( field_char ) )
     70!
     71!                CASE ( 'u2_av' )
     72!                   IF ( .NOT. ALLOCATED( u2_av ) ) THEN
     73!                      ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     74!                   ENDIF
     75!                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     76!                   u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    5877!                                          tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
    5978!
    60 !             CASE DEFAULT
    61 !                PRINT*, '+++ user_init: unknown variable named "', &
    62 !                        TRIM( field_char ), '" found in'
    63 !                PRINT*, '               data from prior run on PE ', myid
    64 !                CALL local_stop
     79!                CASE DEFAULT
     80!                   PRINT*, '+++ user_init: unknown variable named "', &
     81!                           TRIM( field_char ), '" found in'
     82!                   PRINT*, '               data from prior run on PE ', myid
     83!                   CALL local_stop
    6584!
    66 !          END SELECT
     85!             END SELECT
     86!
     87!          ENDDO
    6788!
    6889!          READ ( 13 )  field_char
Note: See TracChangeset for help on using the changeset viewer.