Changeset 4039 for palm/trunk/SOURCE/read_restart_data_mod.f90
- Timestamp:
- Jun 18, 2019 10:32:41 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/read_restart_data_mod.f90
r4017 r4039 25 25 ! ----------------- 26 26 ! $Id$ 27 ! input of uu_av, vv_av, ww_av added 28 ! 29 ! 4017 2019-06-06 12:16:46Z schwenkel 27 30 ! bugfix for r3998, allocation of 3d temporary arrays of various dimensions revised 28 31 ! … … 111 114 !> 112 115 !> @todo: Revise max_pr_cs (profiles for chemistry) 116 !> @todo: Modularize reading of restart data for diagnostic quantities, which 117 !> is not possible with the current module-interface structure 113 118 !------------------------------------------------------------------------------! 114 119 MODULE read_restart_data_mod … … 135 140 136 141 USE diagnostic_output_quantities_mod, & 137 ONLY: ti_av 142 ONLY: ti_av, uu_av, vv_av, ww_av 138 143 139 144 USE grid_variables, & … … 1553 1558 u_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 1554 1559 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1560 1561 CASE ( 'uu_av' ) 1562 IF ( .NOT. ALLOCATED( uu_av ) ) THEN 1563 ALLOCATE( uu_av(nzb:nzt+1,nys:nyn,nxl:nxr) ) 1564 ENDIF 1565 IF ( k == 1 ) THEN 1566 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, & 1567 nxl_on_file:nxr_on_file) ) 1568 READ ( 13 ) tmp_3d_non_standard 1569 ENDIF 1570 uu_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf) 1555 1571 1556 1572 CASE ( 'u_m_l' ) … … 1614 1630 v_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 1615 1631 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1632 1633 CASE ( 'vv_av' ) 1634 IF ( .NOT. ALLOCATED( vv_av ) ) THEN 1635 ALLOCATE( vv_av(nzb:nzt+1,nys:nyn,nxl:nxr) ) 1636 ENDIF 1637 IF ( k == 1 ) THEN 1638 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, & 1639 nxl_on_file:nxr_on_file) ) 1640 READ ( 13 ) tmp_3d_non_standard 1641 ENDIF 1642 vv_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf) 1616 1643 1617 1644 CASE ( 'v_m_l' ) … … 1680 1707 w_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 1681 1708 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1709 1710 CASE ( 'ww_av' ) 1711 IF ( .NOT. ALLOCATED( ww_av ) ) THEN 1712 ALLOCATE( ww_av(nzb:nzt+1,nys:nyn,nxl:nxr) ) 1713 ENDIF 1714 IF ( k == 1 ) THEN 1715 ALLOCATE( tmp_3d_non_standard(nzb:nzt+1,nys_on_file:nyn_on_file, & 1716 nxl_on_file:nxr_on_file) ) 1717 READ ( 13 ) tmp_3d_non_standard 1718 ENDIF 1719 ww_av(:,nysc:nync,nxlc:nxrc) = tmp_3d_non_standard(:,nysf:nynf,nxlf:nxrf) 1682 1720 1683 1721 CASE ( 'w_m_l' )
Note: See TracChangeset
for help on using the changeset viewer.