Changeset 1976 for palm/trunk/SOURCE/sum_up_3d_data.f90
- Timestamp:
- Jul 27, 2016 1:28:04 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/sum_up_3d_data.f90
r1973 r1976 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! Radiation actions are now done directly in the respective module 22 22 ! 23 23 ! Former revisions: … … 138 138 139 139 USE radiation_model_mod, & 140 ONLY: rad_net, rad_net_av, rad_sw_in, rad_sw_in_av, rad_sw_out, & 141 rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr, & 142 rad_sw_hr_av, rad_lw_in, rad_lw_in_av, rad_lw_out, & 143 rad_lw_out_av, rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr, & 144 rad_lw_hr_av 140 ONLY: radiation, radiation_3d_data_averaging 145 141 146 142 … … 291 287 qv_av = 0.0_wp 292 288 293 CASE ( 'rad_net*' )294 IF ( .NOT. ALLOCATED( rad_net_av ) ) THEN295 ALLOCATE( rad_net_av(nysg:nyng,nxlg:nxrg) )296 ENDIF297 rad_net_av = 0.0_wp298 299 CASE ( 'rad_lw_in' )300 IF ( .NOT. ALLOCATED( rad_lw_in_av ) ) THEN301 ALLOCATE( rad_lw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )302 ENDIF303 rad_lw_in_av = 0.0_wp304 305 CASE ( 'rad_lw_out' )306 IF ( .NOT. ALLOCATED( rad_lw_out_av ) ) THEN307 ALLOCATE( rad_lw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )308 ENDIF309 rad_lw_out_av = 0.0_wp310 311 CASE ( 'rad_lw_cs_hr' )312 IF ( .NOT. ALLOCATED( rad_lw_cs_hr_av ) ) THEN313 ALLOCATE( rad_lw_cs_hr_av(nzb+1:nzt+1,nysg:nyng,nxlg:nxrg) )314 ENDIF315 rad_lw_cs_hr_av = 0.0_wp316 317 CASE ( 'rad_lw_hr' )318 IF ( .NOT. ALLOCATED( rad_lw_hr_av ) ) THEN319 ALLOCATE( rad_lw_hr_av(nzb+1:nzt+1,nysg:nyng,nxlg:nxrg) )320 ENDIF321 rad_lw_hr_av = 0.0_wp322 323 CASE ( 'rad_sw_in' )324 IF ( .NOT. ALLOCATED( rad_sw_in_av ) ) THEN325 ALLOCATE( rad_sw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )326 ENDIF327 rad_sw_in_av = 0.0_wp328 329 CASE ( 'rad_sw_out' )330 IF ( .NOT. ALLOCATED( rad_sw_out_av ) ) THEN331 ALLOCATE( rad_sw_out_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )332 ENDIF333 rad_sw_out_av = 0.0_wp334 335 CASE ( 'rad_sw_cs_hr' )336 IF ( .NOT. ALLOCATED( rad_sw_cs_hr_av ) ) THEN337 ALLOCATE( rad_sw_cs_hr_av(nzb+1:nzt+1,nysg:nyng,nxlg:nxrg) )338 ENDIF339 rad_sw_cs_hr_av = 0.0_wp340 341 CASE ( 'rad_sw_hr' )342 IF ( .NOT. ALLOCATED( rad_sw_hr_av ) ) THEN343 ALLOCATE( rad_sw_hr_av(nzb+1:nzt+1,nysg:nyng,nxlg:nxrg) )344 ENDIF345 rad_sw_hr_av = 0.0_wp346 347 289 CASE ( 'rho' ) 348 290 IF ( .NOT. ALLOCATED( rho_av ) ) THEN … … 429 371 IF ( land_surface ) THEN 430 372 CALL lsm_3d_data_averaging( 'allocate', doav(ii) ) 373 ENDIF 374 375 ! 376 !-- Radiation quantity 377 IF ( radiation ) THEN 378 CALL radiation_3d_data_averaging( 'allocate', doav(ii) ) 431 379 ENDIF 432 380 … … 655 603 ENDDO 656 604 657 CASE ( 'rad_net*' )658 DO i = nxlg, nxrg659 DO j = nysg, nyng660 rad_net_av(j,i) = rad_net_av(j,i) + rad_net(j,i)661 ENDDO662 ENDDO663 664 CASE ( 'rad_lw_in' )665 DO i = nxlg, nxrg666 DO j = nysg, nyng667 DO k = nzb, nzt+1668 rad_lw_in_av(k,j,i) = rad_lw_in_av(k,j,i) + rad_lw_in(k,j,i)669 ENDDO670 ENDDO671 ENDDO672 673 CASE ( 'rad_lw_out' )674 DO i = nxlg, nxrg675 DO j = nysg, nyng676 DO k = nzb, nzt+1677 rad_lw_out_av(k,j,i) = rad_lw_out_av(k,j,i) + rad_lw_out(k,j,i)678 ENDDO679 ENDDO680 ENDDO681 682 CASE ( 'rad_lw_cs_hr' )683 DO i = nxlg, nxrg684 DO j = nysg, nyng685 DO k = nzb, nzt+1686 rad_lw_cs_hr_av(k,j,i) = rad_lw_cs_hr_av(k,j,i) + rad_lw_cs_hr(k,j,i)687 ENDDO688 ENDDO689 ENDDO690 691 CASE ( 'rad_lw_hr' )692 DO i = nxlg, nxrg693 DO j = nysg, nyng694 DO k = nzb, nzt+1695 rad_lw_hr_av(k,j,i) = rad_lw_hr_av(k,j,i) + rad_lw_hr(k,j,i)696 ENDDO697 ENDDO698 ENDDO699 700 CASE ( 'rad_sw_in' )701 DO i = nxlg, nxrg702 DO j = nysg, nyng703 DO k = nzb, nzt+1704 rad_sw_in_av(k,j,i) = rad_sw_in_av(k,j,i) + rad_sw_in(k,j,i)705 ENDDO706 ENDDO707 ENDDO708 709 CASE ( 'rad_sw_out' )710 DO i = nxlg, nxrg711 DO j = nysg, nyng712 DO k = nzb, nzt+1713 rad_sw_out_av(k,j,i) = rad_sw_out_av(k,j,i) + rad_sw_out(k,j,i)714 ENDDO715 ENDDO716 ENDDO717 718 CASE ( 'rad_sw_cs_hr' )719 DO i = nxlg, nxrg720 DO j = nysg, nyng721 DO k = nzb, nzt+1722 rad_sw_cs_hr_av(k,j,i) = rad_sw_cs_hr_av(k,j,i) + rad_sw_cs_hr(k,j,i)723 ENDDO724 ENDDO725 ENDDO726 727 CASE ( 'rad_sw_hr' )728 DO i = nxlg, nxrg729 DO j = nysg, nyng730 DO k = nzb, nzt+1731 rad_sw_hr_av(k,j,i) = rad_sw_hr_av(k,j,i) + rad_sw_hr(k,j,i)732 ENDDO733 ENDDO734 ENDDO735 736 605 CASE ( 'rho' ) 737 606 DO i = nxlg, nxrg … … 854 723 855 724 ! 725 !-- Radiation quantity 726 IF ( radiation ) THEN 727 CALL radiation_3d_data_averaging( 'sum', doav(ii) ) 728 ENDIF 729 730 ! 856 731 !-- User-defined quantity 857 732 CALL user_3d_data_averaging( 'sum', doav(ii) )
Note: See TracChangeset
for help on using the changeset viewer.