Changeset 3994 for palm/trunk/SOURCE/time_integration.f90
- Timestamp:
- May 22, 2019 6:08:09 PM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/time_integration.f90
r3988 r3994 25 25 ! ----------------- 26 26 ! $Id$ 27 ! output of turbulence intensity added 28 ! 29 ! 3988 2019-05-22 11:32:37Z kanani 27 30 ! Implement steerable output interval for virtual measurements 28 31 ! … … 561 564 USE date_and_time_mod, & 562 565 ONLY: calc_date_and_time, hour_call_emis, hour_of_year 566 567 USE diagnostic_output_quantities_mod, & 568 ONLY: diagnostic_output_quantities_calculate, & 569 timestep_number_at_prev_calc 563 570 564 571 USE flight_mod, & … … 1655 1662 time_do_sla = time_do_sla + dt_3d 1656 1663 IF ( time_do_sla >= dt_averaging_input ) THEN 1664 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1665 CALL diagnostic_output_quantities_calculate 1666 1657 1667 CALL sum_up_3d_data 1658 1668 average_count_3d = average_count_3d + 1 … … 1733 1743 !-- 2d-data output (cross-sections) 1734 1744 IF ( time_do2d_xy >= dt_do2d_xy .AND. time_since_reference_point >= skip_time_do2d_xy ) THEN 1745 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1746 CALL diagnostic_output_quantities_calculate 1747 1735 1748 CALL data_output_2d( 'xy', 0 ) 1736 1749 time_do2d_xy = MOD( time_do2d_xy, MAX( dt_do2d_xy, dt_3d ) ) 1737 1750 ENDIF 1738 1751 IF ( time_do2d_xz >= dt_do2d_xz .AND. time_since_reference_point >= skip_time_do2d_xz ) THEN 1752 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1753 1754 CALL diagnostic_output_quantities_calculate 1739 1755 CALL data_output_2d( 'xz', 0 ) 1740 1756 time_do2d_xz = MOD( time_do2d_xz, MAX( dt_do2d_xz, dt_3d ) ) 1741 1757 ENDIF 1742 1758 IF ( time_do2d_yz >= dt_do2d_yz .AND. time_since_reference_point >= skip_time_do2d_yz ) THEN 1759 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1760 CALL diagnostic_output_quantities_calculate 1761 1743 1762 CALL data_output_2d( 'yz', 0 ) 1744 1763 time_do2d_yz = MOD( time_do2d_yz, MAX( dt_do2d_yz, dt_3d ) ) … … 1748 1767 !-- 3d-data output (volume data) 1749 1768 IF ( time_do3d >= dt_do3d .AND. time_since_reference_point >= skip_time_do3d ) THEN 1769 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1770 CALL diagnostic_output_quantities_calculate 1771 1750 1772 CALL data_output_3d( 0 ) 1751 1773 time_do3d = MOD( time_do3d, MAX( dt_do3d, dt_3d ) ) … … 1757 1779 IF ( time_domask(mid) >= dt_domask(mid) & 1758 1780 .AND. time_since_reference_point >= skip_time_domask(mid) ) THEN 1781 IF ( current_timestep_number > timestep_number_at_prev_calc ) & 1782 CALL diagnostic_output_quantities_calculate 1783 1759 1784 CALL data_output_mask( 0 ) 1760 1785 time_domask(mid) = MOD( time_domask(mid), MAX( dt_domask(mid), dt_3d ) )
Note: See TracChangeset
for help on using the changeset viewer.