Changeset 4845 for palm/trunk/SOURCE/read_restart_data_mod.f90
- Timestamp:
- Jan 18, 2021 11:15:37 AM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/read_restart_data_mod.f90 (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/read_restart_data_mod.f90
r4828 r4845 25 25 ! ----------------- 26 26 ! $Id$ 27 ! arrays removed that were required for radiation boundary conditions 28 ! 29 ! 4828 2021-01-05 11:21:41Z Giersch 27 30 ! include time_indoor into restart mechanism 28 31 ! … … 149 152 ONLY: inflow_damping_factor, mean_inflow_profiles, pt_init, & 150 153 q_init, ref_state, sa_init, s_init, u_init, ug, v_init, vg, & 151 e, kh, km, p, pt, q, ql, s, u, u_m_l, u_m_n, u_m_r, u_m_s, & 152 v, v_m_l, v_m_n, v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s 154 e, kh, km, p, pt, q, ql, s, u, v, vpt, w 153 155 154 156 USE averaging … … 2058 2060 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 2059 2061 2060 CASE ( 'u_m_l' )2061 IF ( k == 1 ) THEN2062 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2063 1:2) )2064 READ ( 13 ) tmp_3d_non_standard2065 ENDIF2066 IF ( bc_radiation_l ) THEN2067 u_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2068 ENDIF2069 2070 CASE ( 'u_m_n' )2071 IF ( k == 1 ) THEN2072 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &2073 nxl_on_file-nbgp:nxr_on_file+nbgp) )2074 READ ( 13 ) tmp_3d_non_standard2075 ENDIF2076 IF ( bc_radiation_n ) THEN2077 u_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2078 ENDIF2079 2080 CASE ( 'u_m_r' )2081 IF ( k == 1 ) THEN2082 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2083 nx-1:nx) )2084 READ ( 13 ) tmp_3d_non_standard2085 ENDIF2086 IF ( bc_radiation_r ) THEN2087 u_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2088 ENDIF2089 2090 CASE ( 'u_m_s' )2091 IF ( k == 1 ) THEN2092 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1, &2093 nxl_on_file-nbgp:nxr_on_file+nbgp) )2094 READ ( 13 ) tmp_3d_non_standard2095 ENDIF2096 IF ( bc_radiation_s ) THEN2097 u_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2098 ENDIF2099 2100 2062 CASE ( 'us_av' ) 2101 2063 IF ( .NOT. ALLOCATED( us_av ) ) THEN … … 2119 2081 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 2120 2082 2121 CASE ( 'v_m_l' )2122 IF ( k == 1 ) THEN2123 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2124 0:1) )2125 READ ( 13 ) tmp_3d_non_standard2126 ENDIF2127 IF ( bc_radiation_l ) THEN2128 v_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2129 ENDIF2130 2131 CASE ( 'v_m_n' )2132 IF ( k == 1 ) THEN2133 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &2134 nxl_on_file-nbgp:nxr_on_file+nbgp) )2135 READ ( 13 ) tmp_3d_non_standard2136 ENDIF2137 IF ( bc_radiation_n ) THEN2138 v_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2139 ENDIF2140 2141 CASE ( 'v_m_r' )2142 IF ( k == 1 ) THEN2143 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2144 nx-1:nx) )2145 READ ( 13 ) tmp_3d_non_standard2146 ENDIF2147 IF ( bc_radiation_r ) THEN2148 v_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2149 ENDIF2150 2151 CASE ( 'v_m_s' )2152 IF ( k == 1 ) THEN2153 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,1:2, &2154 nxl_on_file-nbgp:nxr_on_file+nbgp) )2155 READ ( 13 ) tmp_3d_non_standard2156 ENDIF2157 IF ( bc_radiation_s ) THEN2158 v_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2159 ENDIF2160 2161 2083 CASE ( 'vpt' ) 2162 2084 IF ( k == 1 ) READ ( 13 ) tmp_3d … … 2184 2106 w_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 2185 2107 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 2186 2187 CASE ( 'w_m_l' )2188 IF ( k == 1 ) THEN2189 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2190 0:1) )2191 READ ( 13 ) tmp_3d_non_standard2192 ENDIF2193 IF ( bc_radiation_l ) THEN2194 w_m_l(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2195 ENDIF2196 2197 CASE ( 'w_m_n' )2198 IF ( k == 1 ) THEN2199 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,ny-1:ny, &2200 nxl_on_file-nbgp:nxr_on_file+nbgp) )2201 READ ( 13 ) tmp_3d_non_standard2202 ENDIF2203 IF ( bc_radiation_n ) THEN2204 w_m_n(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2205 ENDIF2206 2207 CASE ( 'w_m_r' )2208 IF ( k == 1 ) THEN2209 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp, &2210 nx-1:nx) )2211 READ ( 13 ) tmp_3d_non_standard2212 ENDIF2213 IF ( bc_radiation_r ) THEN2214 w_m_r(:,nysc-nbgp:nync+nbgp,:) = tmp_3d_non_standard(:,nysf-nbgp:nynf+nbgp,:)2215 ENDIF2216 2217 CASE ( 'w_m_s' )2218 IF ( k == 1 ) THEN2219 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,0:1, &2220 nxl_on_file-nbgp:nxr_on_file+nbgp) )2221 READ ( 13 ) tmp_3d_non_standard2222 ENDIF2223 IF ( bc_radiation_s ) THEN2224 w_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp)2225 ENDIF2226 2108 2227 2109 CASE ( 'z0_av' ) … … 2488 2370 ENDIF 2489 2371 2490 CALL rd_mpi_io_check_array( 'u_m_l' , found = array_found )2491 IF ( array_found ) THEN2492 IF ( .NOT. ALLOCATED( u_m_l ) ) ALLOCATE( u_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2493 CALL rrd_mpi_io( 'u_m_l', u_m_l )2494 ENDIF2495 2496 CALL rd_mpi_io_check_array( 'u_m_n' , found = array_found )2497 IF ( array_found ) THEN2498 IF ( .NOT. ALLOCATED( u_m_n ) ) ALLOCATE( u_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2499 CALL rrd_mpi_io( 'u_m_n', u_m_n )2500 ENDIF2501 2502 CALL rd_mpi_io_check_array( 'u_m_r' , found = array_found )2503 IF ( array_found ) THEN2504 IF ( .NOT. ALLOCATED( u_m_r ) ) ALLOCATE( u_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2505 CALL rrd_mpi_io( 'u_m_r', u_m_r )2506 ENDIF2507 2508 CALL rd_mpi_io_check_array( 'u_m_s' , found = array_found )2509 IF ( array_found ) THEN2510 IF ( .NOT. ALLOCATED( u_m_s ) ) ALLOCATE( u_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2511 CALL rrd_mpi_io( 'u_m_s', u_m_s )2512 ENDIF2513 2514 2372 CALL rd_mpi_io_check_array( 'us_av' , found = array_found ) 2515 2373 IF ( array_found ) THEN … … 2526 2384 ENDIF 2527 2385 2528 CALL rd_mpi_io_check_array( 'v_m_l' , found = array_found )2529 IF ( array_found ) THEN2530 IF ( .NOT. ALLOCATED( v_m_l ) ) ALLOCATE( v_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2531 CALL rrd_mpi_io( 'v_m_l', v_m_l )2532 ENDIF2533 2534 CALL rd_mpi_io_check_array( 'v_m_n' , found = array_found )2535 IF ( array_found ) THEN2536 IF ( .NOT. ALLOCATED( v_m_n ) ) ALLOCATE( v_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2537 CALL rrd_mpi_io( 'v_m_n', v_m_n )2538 ENDIF2539 2540 CALL rd_mpi_io_check_array( 'v_m_r' , found = array_found )2541 IF ( array_found ) THEN2542 IF ( .NOT. ALLOCATED( v_m_r ) ) ALLOCATE( v_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2543 CALL rrd_mpi_io( 'v_m_r', v_m_r )2544 ENDIF2545 2546 CALL rd_mpi_io_check_array( 'v_m_s' , found = array_found )2547 IF ( array_found ) THEN2548 IF ( .NOT. ALLOCATED( v_m_s ) ) ALLOCATE( v_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2549 CALL rrd_mpi_io( 'v_m_s', v_m_s )2550 ENDIF2551 2552 2386 CALL rd_mpi_io_check_array( 'vpt' , found = array_found ) 2553 2387 IF ( array_found ) THEN … … 2567 2401 IF ( .NOT. ALLOCATED( w_av ) ) ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 2568 2402 CALL rrd_mpi_io( 'w_av', w_av ) 2569 ENDIF2570 2571 CALL rd_mpi_io_check_array( 'w_m_l' , found = array_found )2572 IF ( array_found ) THEN2573 IF ( .NOT. ALLOCATED( w_m_l ) ) ALLOCATE( w_m_l(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2574 CALL rrd_mpi_io( 'w_m_l', w_m_l )2575 ENDIF2576 2577 CALL rd_mpi_io_check_array( 'w_m_n' , found = array_found )2578 IF ( array_found ) THEN2579 IF ( .NOT. ALLOCATED( w_m_n ) ) ALLOCATE( w_m_n(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2580 CALL rrd_mpi_io( 'w_m_n', w_m_n )2581 ENDIF2582 2583 CALL rd_mpi_io_check_array( 'w_m_r' , found = array_found )2584 IF ( array_found ) THEN2585 IF ( .NOT. ALLOCATED( w_m_r ) ) ALLOCATE( w_m_r(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2586 CALL rrd_mpi_io( 'w_m_r', w_m_r )2587 ENDIF2588 2589 CALL rd_mpi_io_check_array( 'w_m_s' , found = array_found )2590 IF ( array_found ) THEN2591 IF ( .NOT. ALLOCATED( w_m_s ) ) ALLOCATE( w_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2592 CALL rrd_mpi_io( 'w_m_s', w_m_s )2593 2403 ENDIF 2594 2404
Note: See TracChangeset
for help on using the changeset viewer.
