Changeset 2743
- Timestamp:
- Jan 12, 2018 4:03:39 PM (7 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
check_parameters.f90 (modified) (3 diffs)
-
data_output_2d.f90 (modified) (4 diffs)
-
sum_up_3d_data.f90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r2742 r2743 25 25 ! ----------------- 26 26 ! $Id$ 27 ! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2. 28 ! 29 ! 2742 2018-01-12 14:59:47Z suehring 27 30 ! Enable output of surface temperature 28 31 ! … … 1047 1050 1048 1051 ! 1049 !-- When the land surface model is used, the flux output must be dynamic. 1050 IF ( land_surface ) THEN 1052 !-- When the land- or urban-surface model is used, the flux output must be 1053 !-- dynamic. 1054 IF ( land_surface .OR. urban_surface ) THEN 1051 1055 flux_output_mode = 'dynamic' 1052 1056 ENDIF … … 3351 3355 IF ( TRIM( var ) == 'z0*' ) unit = 'm' 3352 3356 IF ( TRIM( var ) == 'z0h*' ) unit = 'm' 3357 ! 3358 !-- Output of surface latent and sensible heat flux will be in W/m2 3359 !-- in case of natural- and urban-type surfaces, even if 3360 !-- flux_output_mode is set to kinematic units. 3361 IF ( land_surface .OR. urban_surface ) THEN 3362 IF ( TRIM( var ) == 'shf*' ) unit = 'W/m2' 3363 IF ( TRIM( var ) == 'qsws*' ) unit = 'W/m2' 3364 ENDIF 3353 3365 3354 3366 CASE DEFAULT -
palm/trunk/SOURCE/data_output_2d.f90
r2742 r2743 25 25 ! ----------------- 26 26 ! $Id$ 27 ! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2. 28 ! 29 ! 2742 2018-01-12 14:59:47Z suehring 27 30 ! Enable output of surface temperature 28 31 ! … … 194 197 195 198 USE arrays_3d, & 196 ONLY: dzw, e, nc, nr, p, pt, precipitation_amount, precipitation_rate,& 199 ONLY: dzw, e, heatflux_output_conversion, nc, nr, p, pt, & 200 precipitation_amount, precipitation_rate, & 197 201 prr, q, qc, ql, ql_c, ql_v, ql_vp, qr, rho_ocean, s, sa, & 198 tend, u, v, vpt, w, zu, zw 202 tend, u, v, vpt, w, zu, zw, waterflux_output_conversion 199 203 200 204 USE averaging 201 205 202 206 USE cloud_parameters, & 203 ONLY: hyrho, l_d_cp, pt_d_t207 ONLY: cp, hyrho, l_d_cp, l_v, pt_d_t 204 208 205 209 USE control_parameters, & … … 809 813 CASE ( 'qsws*_xy' ) ! 2d-array 810 814 IF ( av == 0 ) THEN 815 ! 816 !-- In case of default surfaces, clean-up flux by density. 817 !-- In case of land- and urban-surfaces, convert fluxes into 818 !-- dynamic units 811 819 DO m = 1, surf_def_h(0)%ns 812 820 i = surf_def_h(0)%i(m) 813 821 j = surf_def_h(0)%j(m) 814 local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m) 822 k = surf_def_h(0)%k(m) 823 local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m) * & 824 waterflux_output_conversion(k) 815 825 ENDDO 816 826 DO m = 1, surf_lsm_h%ns 817 827 i = surf_lsm_h%i(m) 818 828 j = surf_lsm_h%j(m) 819 local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) 829 k = surf_lsm_h%k(m) 830 local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) * l_v 820 831 ENDDO 821 832 DO m = 1, surf_usm_h%ns 822 833 i = surf_usm_h%i(m) 823 834 j = surf_usm_h%j(m) 824 local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) 835 k = surf_usm_h%k(m) 836 local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v 825 837 ENDDO 826 838 ELSE … … 901 913 CASE ( 'shf*_xy' ) ! 2d-array 902 914 IF ( av == 0 ) THEN 915 ! 916 !-- In case of default surfaces, clean-up flux by density. 917 !-- In case of land- and urban-surfaces, convert fluxes into 918 !-- dynamic units. 903 919 DO m = 1, surf_def_h(0)%ns 904 920 i = surf_def_h(0)%i(m) 905 921 j = surf_def_h(0)%j(m) 906 local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m) 922 k = surf_def_h(0)%k(m) 923 local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m) * & 924 heatflux_output_conversion(k) 907 925 ENDDO 908 926 DO m = 1, surf_lsm_h%ns 909 927 i = surf_lsm_h%i(m) 910 928 j = surf_lsm_h%j(m) 911 local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) 929 k = surf_lsm_h%k(m) 930 local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) * cp 912 931 ENDDO 913 932 DO m = 1, surf_usm_h%ns 914 933 i = surf_usm_h%i(m) 915 934 j = surf_usm_h%j(m) 916 local_pf(i,j,nzb+1) = surf_usm_h%shf(m) 935 k = surf_usm_h%k(m) 936 local_pf(i,j,nzb+1) = surf_usm_h%shf(m) * cp 917 937 ENDDO 918 938 ELSE -
palm/trunk/SOURCE/sum_up_3d_data.f90
r2742 r2743 25 25 ! ----------------- 26 26 ! $Id$ 27 ! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2. 28 ! 29 ! 2742 2018-01-12 14:59:47Z suehring 27 30 ! Enable output of surface temperature 28 31 ! … … 157 160 158 161 USE arrays_3d, & 159 ONLY: dzw, e, nc, nr, p, pt, precipitation_rate, q, qc, ql, ql_c, & 160 ql_v, qr, rho_ocean, s, sa, u, v, vpt, w 162 ONLY: dzw, e, heatflux_output_conversion, nc, nr, p, pt, & 163 precipitation_rate, q, qc, ql, ql_c, & 164 ql_v, qr, rho_ocean, s, sa, u, v, vpt, w, & 165 waterflux_output_conversion 161 166 162 167 USE averaging, & … … 173 178 174 179 USE cloud_parameters, & 175 ONLY: l_d_cp, pt_d_t180 ONLY: cp, l_d_cp, l_v, pt_d_t 176 181 177 182 USE control_parameters, & … … 740 745 741 746 CASE ( 'qsws*' ) 747 ! 748 !-- In case of default surfaces, clean-up flux by density. 749 !-- In case of land- and urban-surfaces, convert fluxes into 750 !-- dynamic units. 742 751 DO m = 1, surf_def_h(0)%ns 743 752 i = surf_def_h(0)%i(m) 744 753 j = surf_def_h(0)%j(m) 745 qsws_av(j,i) = qsws_av(j,i) + surf_def_h(0)%qsws(m) 754 qsws_av(j,i) = qsws_av(j,i) + surf_def_h(0)%qsws(m) * & 755 waterflux_output_conversion(k) 746 756 ENDDO 747 757 DO m = 1, surf_lsm_h%ns 748 758 i = surf_lsm_h%i(m) 749 759 j = surf_lsm_h%j(m) 750 qsws_av(j,i) = qsws_av(j,i) + surf_lsm_h%qsws(m) 760 qsws_av(j,i) = qsws_av(j,i) + surf_lsm_h%qsws(m) * l_v 751 761 ENDDO 752 762 DO m = 1, surf_usm_h%ns 753 763 i = surf_usm_h%i(m) 754 764 j = surf_usm_h%j(m) 755 qsws_av(j,i) = qsws_av(j,i) + surf_usm_h%qsws(m) 765 qsws_av(j,i) = qsws_av(j,i) + surf_usm_h%qsws(m) * l_v 756 766 ENDDO 757 767 … … 813 823 814 824 CASE ( 'shf*' ) 825 ! 826 !-- In case of default surfaces, clean-up flux by density. 827 !-- In case of land- and urban-surfaces, convert fluxes into 828 !-- dynamic units. 815 829 DO m = 1, surf_def_h(0)%ns 816 830 i = surf_def_h(0)%i(m) 817 831 j = surf_def_h(0)%j(m) 818 shf_av(j,i) = shf_av(j,i) + surf_def_h(0)%shf(m) 832 shf_av(j,i) = shf_av(j,i) + surf_def_h(0)%shf(m) * & 833 heatflux_output_conversion(k) 819 834 ENDDO 820 835 DO m = 1, surf_lsm_h%ns 821 836 i = surf_lsm_h%i(m) 822 837 j = surf_lsm_h%j(m) 823 shf_av(j,i) = shf_av(j,i) + surf_lsm_h%shf(m) 838 shf_av(j,i) = shf_av(j,i) + surf_lsm_h%shf(m) * cp 824 839 ENDDO 825 840 DO m = 1, surf_usm_h%ns 826 841 i = surf_usm_h%i(m) 827 842 j = surf_usm_h%j(m) 828 shf_av(j,i) = shf_av(j,i) + surf_usm_h%shf(m) 843 shf_av(j,i) = shf_av(j,i) + surf_usm_h%shf(m) * cp 829 844 ENDDO 830 845
Note: See TracChangeset
for help on using the changeset viewer.
