Ignore:
Timestamp:
Apr 30, 2015 7:05:52 AM (6 years ago)
Author:
maronga
Message:

Added support for RRTMG radiation code

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r1556 r1585  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Adapted for RRTMG
    2323!
    2424! Former revisions:
     
    112112
    113113    USE radiation_model_mod,                                                   &
    114         ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av
     114        ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av, rad_sw_out,       &
     115               rad_sw_out_av, rad_lw_in, rad_lw_in_av, rad_lw_out,             &
     116               rad_lw_out_av
    115117
    116118    IMPLICIT NONE
     
    326328                rad_net_av = 0.0_wp
    327329
    328              CASE ( 'rad_sw_in*' )
     330             CASE ( 'rad_lw_in' )
     331                IF ( .NOT. ALLOCATED( rad_lw_in_av ) )  THEN
     332                   ALLOCATE( rad_lw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     333                ENDIF
     334                rad_lw_in_av = 0.0_wp
     335
     336             CASE ( 'rad_lw_out' )
     337                IF ( .NOT. ALLOCATED( rad_lw_out_av ) )  THEN
     338                   ALLOCATE( rad_lw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     339                ENDIF
     340                rad_lw_out_av = 0.0_wp
     341
     342             CASE ( 'rad_sw_in' )
    329343                IF ( .NOT. ALLOCATED( rad_sw_in_av ) )  THEN
    330                    ALLOCATE( rad_sw_in_av(nysg:nyng,nxlg:nxrg) )
     344                   ALLOCATE( rad_sw_in_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    331345                ENDIF
    332346                rad_sw_in_av = 0.0_wp
     347
     348             CASE ( 'rad_sw_out' )
     349                IF ( .NOT. ALLOCATED( rad_sw_out_av ) )  THEN
     350                   ALLOCATE( rad_sw_out_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     351                ENDIF
     352                rad_sw_out_av = 0.0_wp
    333353
    334354             CASE ( 'rho' )
     
    732752             ENDDO
    733753
    734           CASE ( 'rad_sw_in*' )
    735              DO  i = nxlg, nxrg
    736                 DO  j = nysg, nyng
    737                    rad_sw_in_av(j,i) = rad_sw_in_av(j,i) + rad_sw_in(j,i)
     754          CASE ( 'rad_lw_in' )
     755             DO  i = nxlg, nxrg
     756                DO  j = nysg, nyng
     757                   DO  k = nzb, nzt+1
     758                      rad_lw_in_av(k,j,i) = rad_lw_in_av(k,j,i) + rad_lw_in(k,j,i)
     759                   ENDDO
     760                ENDDO
     761             ENDDO
     762
     763          CASE ( 'rad_lw_out' )
     764             DO  i = nxlg, nxrg
     765                DO  j = nysg, nyng
     766                   DO  k = nzb, nzt+1
     767                      rad_lw_out_av(k,j,i) = rad_lw_out_av(k,j,i) + rad_lw_out(k,j,i)
     768                   ENDDO
     769                ENDDO
     770             ENDDO
     771
     772
     773          CASE ( 'rad_sw_in' )
     774             DO  i = nxlg, nxrg
     775                DO  j = nysg, nyng
     776                   DO  k = nzb, nzt+1
     777                      rad_sw_in_av(k,j,i) = rad_sw_in_av(k,j,i) + rad_sw_in(k,j,i)
     778                   ENDDO
     779                ENDDO
     780             ENDDO
     781
     782          CASE ( 'rad_sw_out' )
     783             DO  i = nxlg, nxrg
     784                DO  j = nysg, nyng
     785                   DO  k = nzb, nzt+1
     786                      rad_sw_out_av(k,j,i) = rad_sw_out_av(k,j,i) + rad_sw_out(k,j,i)
     787                   ENDDO
    738788                ENDDO
    739789             ENDDO
Note: See TracChangeset for help on using the changeset viewer.