Changeset 4518 for palm/trunk/SOURCE/read_restart_data_mod.f90
- Timestamp:
- May 4, 2020 3:44:28 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/read_restart_data_mod.f90
r4517 r4518 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move input of diagnostic output quantities to doq_rrd_local 28 ! 29 ! 4517 2020-05-03 14:29:30Z raasch 27 30 ! qsurf and ssurf added 28 31 ! … … 130 133 USE cpulog, & 131 134 ONLY: cpu_log, log_point_s 132 133 USE diagnostic_output_quantities_mod, &134 ONLY: pt_2m_av, &135 ti_av, &136 u_center_av, &137 uu_av, &138 uv_10m_av, &139 v_center_av, &140 vv_av, &141 wspeed_av, &142 ww_av143 135 144 136 USE grid_variables, & … … 1773 1765 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1774 1766 1775 CASE ( 'pt_2m_av' )1776 IF ( .NOT. ALLOCATED( pt_2m_av ) ) THEN1777 ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )1778 ENDIF1779 IF ( k == 1 ) READ ( 13 ) tmp_2d1780 pt_2m_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1781 tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1782 1783 1767 CASE ( 'q' ) 1784 1768 IF ( k == 1 ) READ ( 13 ) tmp_3d … … 1895 1879 tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1896 1880 1897 CASE ( 'ti_av' )1898 IF ( .NOT. ALLOCATED( ti_av ) ) THEN1899 ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )1900 ENDIF1901 IF ( k == 1 ) THEN1902 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &1903 nxl_on_file:nxr_on_file) )1904 READ ( 13 ) tmp_3d_non_standard1905 ENDIF1906 ti_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)1907 1908 1881 CASE ( 'ts_av' ) 1909 1882 IF ( .NOT. ALLOCATED( ts_av ) ) THEN … … 1934 1907 u_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 1935 1908 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1936 1937 CASE ( 'u_center_av' )1938 IF ( .NOT. ALLOCATED( u_center_av ) ) THEN1939 ALLOCATE( u_center_av(nzb:nzt+1,nys:nyn,nxl:nxr) )1940 ENDIF1941 IF ( k == 1 ) THEN1942 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &1943 nxl_on_file:nxr_on_file) )1944 READ ( 13 ) tmp_3d_non_standard1945 ENDIF1946 u_center_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)1947 1948 CASE ( 'uu_av' )1949 IF ( .NOT. ALLOCATED( uu_av ) ) THEN1950 ALLOCATE( uu_av(nzb:nzt+1,nys:nyn,nxl:nxr) )1951 ENDIF1952 IF ( k == 1 ) THEN1953 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &1954 nxl_on_file:nxr_on_file) )1955 READ ( 13 ) tmp_3d_non_standard1956 ENDIF1957 uu_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)1958 1959 CASE ( 'uv_10m_av' )1960 IF ( .NOT. ALLOCATED( uv_10m_av ) ) THEN1961 ALLOCATE( uv_10m_av(nysg:nyng,nxlg:nxrg) )1962 ENDIF1963 IF ( k == 1 ) READ ( 13 ) tmp_2d1964 uv_10m_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1965 tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1966 1909 1967 1910 CASE ( 'u_m_l' ) … … 2026 1969 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 2027 1970 2028 CASE ( 'v_center_av' )2029 IF ( .NOT. ALLOCATED( v_center_av ) ) THEN2030 ALLOCATE( v_center_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2031 ENDIF2032 IF ( k == 1 ) THEN2033 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &2034 nxl_on_file:nxr_on_file) )2035 READ ( 13 ) tmp_3d_non_standard2036 ENDIF2037 v_center_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)2038 2039 CASE ( 'vv_av' )2040 IF ( .NOT. ALLOCATED( vv_av ) ) THEN2041 ALLOCATE( vv_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2042 ENDIF2043 IF ( k == 1 ) THEN2044 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &2045 nxl_on_file:nxr_on_file) )2046 READ ( 13 ) tmp_3d_non_standard2047 ENDIF2048 vv_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)2049 2050 1971 CASE ( 'v_m_l' ) 2051 1972 IF ( k == 1 ) THEN … … 2114 2035 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 2115 2036 2116 CASE ( 'ww_av' )2117 IF ( .NOT. ALLOCATED( ww_av ) ) THEN2118 ALLOCATE( ww_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2119 ENDIF2120 IF ( k == 1 ) THEN2121 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &2122 nxl_on_file:nxr_on_file) )2123 READ ( 13 ) tmp_3d_non_standard2124 ENDIF2125 ww_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)2126 2127 2037 CASE ( 'w_m_l' ) 2128 2038 IF ( k == 1 ) THEN … … 2164 2074 w_m_s(:,:,nxlc-nbgp:nxrc+nbgp) = tmp_3d_non_standard(:,:,nxlf-nbgp:nxrf+nbgp) 2165 2075 ENDIF 2166 2167 CASE ( 'wspeed_av' )2168 IF ( .NOT. ALLOCATED( wspeed_av ) ) THEN2169 ALLOCATE( wspeed_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2170 ENDIF2171 IF ( k == 1 ) THEN2172 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, &2173 nxl_on_file:nxr_on_file) )2174 READ ( 13 ) tmp_3d_non_standard2175 ENDIF2176 wspeed_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf)2177 2076 2178 2077 CASE ( 'z0_av' ) … … 2326 2225 ENDIF 2327 2226 2328 CALL rd_mpi_io_check_array( 'pt_2m_av' , found = array_found )2329 IF ( array_found ) THEN2330 IF ( .NOT. ALLOCATED( pt_2m_av ) ) ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )2331 CALL rrd_mpi_io( 'pt_2m_av', pt_2m_av )2332 ENDIF2333 2334 2227 CALL rd_mpi_io_check_array( 'q' , found = array_found ) 2335 2228 IF ( array_found ) THEN … … 2425 2318 ENDIF 2426 2319 2427 CALL rd_mpi_io_check_array( 'ti_av' , found = array_found )2428 IF ( array_found ) THEN2429 IF ( .NOT. ALLOCATED( ti_av ) ) ALLOCATE( ti_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )2430 CALL rrd_mpi_io( 'ti_av', ti_av )2431 ENDIF2432 2433 2320 CALL rd_mpi_io_check_array( 'ts_av' , found = array_found ) 2434 2321 IF ( array_found ) THEN … … 2451 2338 ENDIF 2452 2339 2453 CALL rd_mpi_io_check_array( 'u_center_av' , found = array_found )2454 IF ( array_found ) THEN2455 IF ( .NOT. ALLOCATED( u_center_av ) ) ALLOCATE( u_center_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2456 CALL rrd_mpi_io( 'u_center_av', u_center_av )2457 ENDIF2458 2459 CALL rd_mpi_io_check_array( 'uu_av' , found = array_found )2460 IF ( array_found ) THEN2461 IF ( .NOT. ALLOCATED( uu_av ) ) ALLOCATE( uu_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2462 CALL rrd_mpi_io( 'uu_av', uu_av )2463 ENDIF2464 2465 CALL rd_mpi_io_check_array( 'uv_10m_av' , found = array_found )2466 IF ( array_found ) THEN2467 IF ( .NOT. ALLOCATED( uv_10m_av ) ) ALLOCATE( uv_10m_av(nysg:nyng,nxlg:nxrg) )2468 CALL rrd_mpi_io( 'uv_10m_av', uv_10m_av )2469 ENDIF2470 2471 2340 CALL rd_mpi_io_check_array( 'u_m_l' , found = array_found ) 2472 2341 IF ( array_found ) THEN … … 2507 2376 ENDIF 2508 2377 2509 CALL rd_mpi_io_check_array( 'v_center_av' , found = array_found )2510 IF ( array_found ) THEN2511 IF ( .NOT. ALLOCATED( v_center_av ) ) ALLOCATE( v_center_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2512 CALL rrd_mpi_io( 'v_center_av', v_center_av )2513 ENDIF2514 2515 CALL rd_mpi_io_check_array( 'vv_av' , found = array_found )2516 IF ( array_found ) THEN2517 IF ( .NOT. ALLOCATED( vv_av ) ) ALLOCATE( vv_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2518 CALL rrd_mpi_io( 'vv_av', vv_av )2519 ENDIF2520 2521 2378 CALL rd_mpi_io_check_array( 'v_m_l' , found = array_found ) 2522 2379 IF ( array_found ) THEN … … 2562 2419 ENDIF 2563 2420 2564 CALL rd_mpi_io_check_array( 'ww_av' , found = array_found )2565 IF ( array_found ) THEN2566 IF ( .NOT. ALLOCATED( ww_av ) ) ALLOCATE( w_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2567 CALL rrd_mpi_io( 'ww_av', ww_av )2568 ENDIF2569 2570 2421 CALL rd_mpi_io_check_array( 'w_m_l' , found = array_found ) 2571 2422 IF ( array_found ) THEN … … 2590 2441 IF ( .NOT. ALLOCATED( w_m_s ) ) ALLOCATE( w_m_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 2591 2442 CALL rrd_mpi_io( 'w_m_s', w_m_s ) 2592 ENDIF2593 2594 CALL rd_mpi_io_check_array( 'wspeed_av' , found = array_found )2595 IF ( array_found ) THEN2596 IF ( .NOT. ALLOCATED( wspeed_av ) ) ALLOCATE( wspeed_av(nzb:nzt+1,nys:nyn,nxl:nxr) )2597 CALL rrd_mpi_io( 'wspeed_av', wspeed_av )2598 2443 ENDIF 2599 2444
Note: See TracChangeset
for help on using the changeset viewer.