Changeset 4502 for palm/trunk/SOURCE/surface_data_output_mod.f90
- Timestamp:
- Apr 17, 2020 4:14:16 PM (5 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/surface_data_output_mod.f90 (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_data_output_mod.f90
r4500 r4502 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Implementation of ice microphysics 28 ! 29 ! 4500 2020-04-17 10:12:45Z suehring 27 30 ! - Correct output of ground/wall heat flux at USM surfaces 28 31 ! - Add conversion factor to heat and momentum-flux output … … 1695 1698 ENDIF 1696 1699 1700 CASE ( 'qis' ) 1701 ! 1702 !-- Output of instantaneous data 1703 IF ( av == 0 ) THEN 1704 CALL surface_data_output_collect( surf_def_h(0)%qis, & 1705 surf_def_h(1)%qis, & 1706 surf_lsm_h%qis, & 1707 surf_usm_h%qis, & 1708 surf_def_v(0)%qis, & 1709 surf_lsm_v(0)%qis, & 1710 surf_usm_v(0)%qis, & 1711 surf_def_v(1)%qis, & 1712 surf_lsm_v(1)%qis, & 1713 surf_usm_v(1)%qis, & 1714 surf_def_v(2)%qis, & 1715 surf_lsm_v(2)%qis, & 1716 surf_usm_v(2)%qis, & 1717 surf_def_v(3)%qis, & 1718 surf_lsm_v(3)%qis, & 1719 surf_usm_v(3)%qis ) 1720 ELSE 1721 ! 1722 !-- Output of averaged data 1723 surfaces%var_out(:) = surfaces%var_av(:,n_out) / & 1724 REAL( average_count_surf, KIND=wp ) 1725 surfaces%var_av(:,n_out) = 0.0_wp 1726 1727 ENDIF 1728 1729 CASE ( 'nis' ) 1730 ! 1731 !-- Output of instantaneous data 1732 IF ( av == 0 ) THEN 1733 CALL surface_data_output_collect( surf_def_h(0)%nis, & 1734 surf_def_h(1)%nis, & 1735 surf_lsm_h%nis, & 1736 surf_usm_h%nis, & 1737 surf_def_v(0)%nis, & 1738 surf_lsm_v(0)%nis, & 1739 surf_usm_v(0)%nis, & 1740 surf_def_v(1)%nis, & 1741 surf_lsm_v(1)%nis, & 1742 surf_usm_v(1)%nis, & 1743 surf_def_v(2)%nis, & 1744 surf_lsm_v(2)%nis, & 1745 surf_usm_v(2)%nis, & 1746 surf_def_v(3)%nis, & 1747 surf_lsm_v(3)%nis, & 1748 surf_usm_v(3)%nis ) 1749 ELSE 1750 ! 1751 !-- Output of averaged data 1752 surfaces%var_out(:) = surfaces%var_av(:,n_out) / & 1753 REAL( average_count_surf, KIND=wp ) 1754 surfaces%var_av(:,n_out) = 0.0_wp 1755 1756 ENDIF 1757 1697 1758 CASE ( 'qrs' ) 1698 1759 ! … … 2153 2214 surf_lsm_v(3)%ncsws, & 2154 2215 surf_usm_v(3)%ncsws ) 2216 ELSE 2217 ! 2218 !-- Output of averaged data 2219 surfaces%var_out(:) = surfaces%var_av(:,n_out) / & 2220 REAL( average_count_surf, KIND=wp ) 2221 surfaces%var_av(:,n_out) = 0.0_wp 2222 2223 ENDIF 2224 2225 2226 CASE ( 'qisws' ) 2227 ! 2228 !-- Output of instantaneous data 2229 IF ( av == 0 ) THEN 2230 CALL surface_data_output_collect( surf_def_h(0)%qisws, & 2231 surf_def_h(1)%qisws, & 2232 surf_lsm_h%qisws, & 2233 surf_usm_h%qisws, & 2234 surf_def_v(0)%qisws, & 2235 surf_lsm_v(0)%qisws, & 2236 surf_usm_v(0)%qisws, & 2237 surf_def_v(1)%qisws, & 2238 surf_lsm_v(1)%qisws, & 2239 surf_usm_v(1)%qisws, & 2240 surf_def_v(2)%qisws, & 2241 surf_lsm_v(2)%qisws, & 2242 surf_usm_v(2)%qisws, & 2243 surf_def_v(3)%qisws, & 2244 surf_lsm_v(3)%qisws, & 2245 surf_usm_v(3)%qisws ) 2246 ELSE 2247 ! 2248 !-- Output of averaged data 2249 surfaces%var_out(:) = surfaces%var_av(:,n_out) / & 2250 REAL( average_count_surf, KIND=wp ) 2251 surfaces%var_av(:,n_out) = 0.0_wp 2252 2253 ENDIF 2254 2255 CASE ( 'nisws' ) 2256 ! 2257 !-- Output of instantaneous data 2258 IF ( av == 0 ) THEN 2259 CALL surface_data_output_collect( surf_def_h(0)%nisws, & 2260 surf_def_h(1)%nisws, & 2261 surf_lsm_h%nisws, & 2262 surf_usm_h%nisws, & 2263 surf_def_v(0)%nisws, & 2264 surf_lsm_v(0)%nisws, & 2265 surf_usm_v(0)%nisws, & 2266 surf_def_v(1)%nisws, & 2267 surf_lsm_v(1)%nisws, & 2268 surf_usm_v(1)%nisws, & 2269 surf_def_v(2)%nisws, & 2270 surf_lsm_v(2)%nisws, & 2271 surf_usm_v(2)%nisws, & 2272 surf_def_v(3)%nisws, & 2273 surf_lsm_v(3)%nisws, & 2274 surf_usm_v(3)%nisws ) 2155 2275 ELSE 2156 2276 ! … … 3120 3240 surf_usm_v(3)%ncs, n_out ) 3121 3241 3242 CASE ( 'qis' ) 3243 CALL surface_data_output_sum_up( surf_def_h(0)%qis, & 3244 surf_def_h(1)%qis, & 3245 surf_lsm_h%qis, & 3246 surf_usm_h%qis, & 3247 surf_def_v(0)%qis, & 3248 surf_lsm_v(0)%qis, & 3249 surf_usm_v(0)%qis, & 3250 surf_def_v(1)%qis, & 3251 surf_lsm_v(1)%qis, & 3252 surf_usm_v(1)%qis, & 3253 surf_def_v(2)%qis, & 3254 surf_lsm_v(2)%qis, & 3255 surf_usm_v(2)%qis, & 3256 surf_def_v(3)%qis, & 3257 surf_lsm_v(3)%qis, & 3258 surf_usm_v(3)%qrs, n_out ) 3259 3260 CASE ( 'nis' ) 3261 CALL surface_data_output_sum_up( surf_def_h(0)%nis, & 3262 surf_def_h(1)%nis, & 3263 surf_lsm_h%nis, & 3264 surf_usm_h%nis, & 3265 surf_def_v(0)%nis, & 3266 surf_lsm_v(0)%nis, & 3267 surf_usm_v(0)%nis, & 3268 surf_def_v(1)%nis, & 3269 surf_lsm_v(1)%nis, & 3270 surf_usm_v(1)%nis, & 3271 surf_def_v(2)%nis, & 3272 surf_lsm_v(2)%nis, & 3273 surf_usm_v(2)%nis, & 3274 surf_def_v(3)%nis, & 3275 surf_lsm_v(3)%nis, & 3276 surf_usm_v(3)%nis, n_out ) 3277 3122 3278 CASE ( 'qrs' ) 3123 3279 CALL surface_data_output_sum_up( surf_def_h(0)%qrs, & … … 3411 3567 surf_lsm_v(3)%ncsws, & 3412 3568 surf_usm_v(3)%ncsws, n_out ) 3569 3570 CASE ( 'qisws' ) 3571 CALL surface_data_output_sum_up( surf_def_h(0)%qisws, & 3572 surf_def_h(1)%qisws, & 3573 surf_lsm_h%qisws, & 3574 surf_usm_h%qisws, & 3575 surf_def_v(0)%qisws, & 3576 surf_lsm_v(0)%qisws, & 3577 surf_usm_v(0)%qisws, & 3578 surf_def_v(1)%qisws, & 3579 surf_lsm_v(1)%qisws, & 3580 surf_usm_v(1)%qisws, & 3581 surf_def_v(2)%qisws, & 3582 surf_lsm_v(2)%qisws, & 3583 surf_usm_v(2)%qisws, & 3584 surf_def_v(3)%qisws, & 3585 surf_lsm_v(3)%qisws, & 3586 surf_usm_v(3)%qisws, n_out ) 3587 3588 CASE ( 'nisws' ) 3589 CALL surface_data_output_sum_up( surf_def_h(0)%nisws, & 3590 surf_def_h(1)%nisws, & 3591 surf_lsm_h%nisws, & 3592 surf_usm_h%nisws, & 3593 surf_def_v(0)%nisws, & 3594 surf_lsm_v(0)%nisws, & 3595 surf_usm_v(0)%nisws, & 3596 surf_def_v(1)%nisws, & 3597 surf_lsm_v(1)%nisws, & 3598 surf_usm_v(1)%nisws, & 3599 surf_def_v(2)%nisws, & 3600 surf_lsm_v(2)%nisws, & 3601 surf_usm_v(2)%nisws, & 3602 surf_def_v(3)%nisws, & 3603 surf_lsm_v(3)%nisws, & 3604 surf_usm_v(3)%nisws, n_out ) 3413 3605 3414 3606 CASE ( 'qrsws' ) … … 4533 4725 unit = 'm/s' 4534 4726 4535 CASE ( 'ss', 'qcs', 'ncs', 'q rs', 'nrs' )4727 CASE ( 'ss', 'qcs', 'ncs', 'qis', 'nis', 'qrs', 'nrs' ) 4536 4728 unit = '1' 4537 4729 … … 4545 4737 unit = 'm2/s2' 4546 4738 4547 CASE ( 'qcsws', 'ncsws', 'q rsws', 'nrsws', 'sasws' )4739 CASE ( 'qcsws', 'ncsws', 'qisws', 'nisws', 'qrsws', 'nrsws', 'sasws' ) 4548 4740 4549 4741 CASE ( 'shf' )
Note: See TracChangeset
for help on using the changeset viewer.
