Changeset 3274 for palm/trunk/SOURCE/sum_up_3d_data.f90
- Timestamp:
- Sep 24, 2018 3:42:55 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/sum_up_3d_data.f90
r3241 r3274 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Modularization of all bulk cloud physics code components 28 ! 29 ! 3241 2018-09-12 15:02:00Z raasch 27 30 ! unused variables removed 28 31 ! … … 204 207 ONLY: dzw, e, heatflux_output_conversion, nc, nr, p, prr, pt, & 205 208 q, qc, ql, ql_c, ql_v, qr, rho_ocean, s, sa, u, v, vpt, w, & 206 waterflux_output_conversion 209 waterflux_output_conversion, d_exner 207 210 208 211 USE averaging, & 209 ONLY: e_av, ghf_av, lpt_av, lwp_av, nc_av, & 210 nr_av, & 211 ol_av, p_av, pc_av, pr_av, prr_av, pt_av, q_av, qc_av, ql_av, & 212 ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, r_a_av, & 212 ONLY: e_av, ghf_av, lpt_av, lwp_av, & 213 ol_av, p_av, pc_av, pr_av, pt_av, q_av, ql_av, & 214 ql_c_av, ql_v_av, ql_vp_av, qsws_av, qv_av, r_a_av, & 213 215 rho_ocean_av, s_av, sa_av, shf_av, ssws_av, ts_av, tsurf_av, & 214 216 u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av 215 217 218 USE basic_constants_and_equations_mod, & 219 ONLY: c_p, lv_d_cp, l_v 220 221 USE bulk_cloud_model_mod, & 222 ONLY: bulk_cloud_model, bcm_3d_data_averaging 223 216 224 USE chemistry_model_mod, & 217 225 ONLY: chem_3d_data_averaging 218 226 219 USE cloud_parameters, &220 ONLY: cp, l_d_cp, l_v, pt_d_t221 222 227 USE control_parameters, & 223 ONLY: air_chemistry, average_count_3d, cloud_physics, doav, doav_n,&228 ONLY: air_chemistry, average_count_3d, doav, doav_n, & 224 229 land_surface, rho_surface, urban_surface, uv_exposure, & 225 230 varnamelength … … 323 328 lwp_av = 0.0_wp 324 329 325 CASE ( 'nc' )326 IF ( .NOT. ALLOCATED( nc_av ) ) THEN327 ALLOCATE( nc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )328 ENDIF329 nc_av = 0.0_wp330 331 CASE ( 'nr' )332 IF ( .NOT. ALLOCATED( nr_av ) ) THEN333 ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )334 ENDIF335 nr_av = 0.0_wp336 337 330 CASE ( 'ol*' ) 338 331 IF ( .NOT. ALLOCATED( ol_av ) ) THEN … … 359 352 pr_av = 0.0_wp 360 353 361 CASE ( 'prr' )362 IF ( .NOT. ALLOCATED( prr_av ) ) THEN363 ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )364 ENDIF365 prr_av = 0.0_wp366 367 354 CASE ( 'pt' ) 368 355 IF ( .NOT. ALLOCATED( pt_av ) ) THEN … … 377 364 q_av = 0.0_wp 378 365 379 CASE ( 'qc' )380 IF ( .NOT. ALLOCATED( qc_av ) ) THEN381 ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )382 ENDIF383 qc_av = 0.0_wp384 385 366 CASE ( 'ql' ) 386 367 IF ( .NOT. ALLOCATED( ql_av ) ) THEN … … 406 387 ENDIF 407 388 ql_vp_av = 0.0_wp 408 409 CASE ( 'qr' )410 IF ( .NOT. ALLOCATED( qr_av ) ) THEN411 ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )412 ENDIF413 qr_av = 0.0_wp414 389 415 390 CASE ( 'qsws*' ) … … 532 507 !-- Turbulence closure module 533 508 CALL tcm_3d_data_averaging( 'allocate', doav(ii) ) 509 510 ! 511 !-- Microphysics module quantities 512 IF ( bulk_cloud_model ) THEN 513 CALL bcm_3d_data_averaging( 'allocate', doav(ii) ) 514 ENDIF 534 515 535 516 ! … … 654 635 ENDIF 655 636 656 CASE ( 'nc' )657 IF ( ALLOCATED( nc_av ) ) THEN658 DO i = nxlg, nxrg659 DO j = nysg, nyng660 DO k = nzb, nzt+1661 nc_av(k,j,i) = nc_av(k,j,i) + nc(k,j,i)662 ENDDO663 ENDDO664 ENDDO665 ENDIF666 667 CASE ( 'nr' )668 IF ( ALLOCATED( nr_av ) ) THEN669 DO i = nxlg, nxrg670 DO j = nysg, nyng671 DO k = nzb, nzt+1672 nr_av(k,j,i) = nr_av(k,j,i) + nr(k,j,i)673 ENDDO674 ENDDO675 ENDDO676 ENDIF677 678 637 CASE ( 'ol*' ) 679 638 IF ( ALLOCATED( ol_av ) ) THEN … … 758 717 ENDIF 759 718 760 CASE ( 'prr' )761 IF ( ALLOCATED( prr_av ) ) THEN762 DO i = nxlg, nxrg763 DO j = nysg, nyng764 DO k = nzb, nzt+1765 prr_av(k,j,i) = prr_av(k,j,i) + prr(k,j,i)766 ENDDO767 ENDDO768 ENDDO769 ENDIF770 771 719 CASE ( 'pt' ) 772 720 IF ( ALLOCATED( pt_av ) ) THEN 773 IF ( .NOT. cloud_physics) THEN721 IF ( .NOT. bulk_cloud_model ) THEN 774 722 DO i = nxlg, nxrg 775 723 DO j = nysg, nyng … … 783 731 DO j = nysg, nyng 784 732 DO k = nzb, nzt+1 785 pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + l _d_cp * &786 pt_d_t(k) * ql(k,j,i)733 pt_av(k,j,i) = pt_av(k,j,i) + pt(k,j,i) + lv_d_cp * & 734 d_exner(k) * ql(k,j,i) 787 735 ENDDO 788 736 ENDDO … … 797 745 DO k = nzb, nzt+1 798 746 q_av(k,j,i) = q_av(k,j,i) + q(k,j,i) 799 ENDDO800 ENDDO801 ENDDO802 ENDIF803 804 CASE ( 'qc' )805 IF ( ALLOCATED( qc_av ) ) THEN806 DO i = nxlg, nxrg807 DO j = nysg, nyng808 DO k = nzb, nzt+1809 qc_av(k,j,i) = qc_av(k,j,i) + qc(k,j,i)810 747 ENDDO 811 748 ENDDO … … 862 799 ENDIF 863 800 ENDDO 864 ENDDO865 ENDDO866 ENDDO867 ENDIF868 869 CASE ( 'qr' )870 IF ( ALLOCATED( qr_av ) ) THEN871 DO i = nxlg, nxrg872 DO j = nysg, nyng873 DO k = nzb, nzt+1874 qr_av(k,j,i) = qr_av(k,j,i) + qr(k,j,i)875 801 ENDDO 876 802 ENDDO … … 1001 927 m = surf_lsm_h%end_index(j,i) 1002 928 shf_av(j,i) = shf_av(j,i) + & 1003 surf_lsm_h%shf(m) * c p929 surf_lsm_h%shf(m) * c_p 1004 930 ELSEIF ( match_usm ) THEN 1005 931 m = surf_usm_h%end_index(j,i) 1006 932 shf_av(j,i) = shf_av(j,i) + & 1007 surf_usm_h%shf(m) * c p933 surf_usm_h%shf(m) * c_p 1008 934 ENDIF 1009 935 ENDDO … … 1265 1191 1266 1192 ! 1193 !-- Microphysics module quantities 1194 IF ( bulk_cloud_model ) THEN 1195 CALL bcm_3d_data_averaging( 'sum', doav(ii) ) 1196 ENDIF 1197 1198 ! 1267 1199 !-- Land surface quantity 1268 1200 IF ( land_surface ) THEN
Note: See TracChangeset
for help on using the changeset viewer.