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/average_3d_data.f90

    r1556 r1585  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Adapted for RRTMG
    2323!
    2424! Former revisions:
     
    8888
    8989    USE radiation_model_mod,                                                   &
    90         ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av
     90        ONLY:  rad_net, rad_net_av, rad_lw_in, rad_lw_in_av, rad_lw_out,       &
     91               rad_lw_out_av, rad_sw_in, rad_sw_in_av, rad_sw_out,             &
     92               rad_sw_out_av
    9193
    9294    IMPLICIT NONE
     
    349351             ENDDO
    350352
    351          CASE ( 'rad_sw_in*' )
    352              DO  i = nxlg, nxrg
    353                 DO  j = nysg, nyng
    354                    rad_sw_in_av(j,i) = rad_sw_in_av(j,i) / REAL( average_count_3d, KIND=wp )
    355                 ENDDO
    356              ENDDO
    357 
    358353         CASE ( 'rad_net*' )
    359354             DO  i = nxlg, nxrg
    360355                DO  j = nysg, nyng
    361356                   rad_net_av(j,i) = rad_net_av(j,i) / REAL( average_count_3d, KIND=wp )
     357                ENDDO
     358             ENDDO
     359
     360          CASE ( 'rad_lw_in' )
     361             DO  i = nxlg, nxrg
     362                DO  j = nysg, nyng
     363                   DO  k = nzb, nzt+1
     364                      rad_lw_in_av(k,j,i) = rad_lw_in_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     365                   ENDDO
     366                ENDDO
     367             ENDDO
     368
     369          CASE ( 'rad_lw_out' )
     370             DO  i = nxlg, nxrg
     371                DO  j = nysg, nyng
     372                   DO  k = nzb, nzt+1
     373                      rad_lw_out_av(k,j,i) = rad_lw_out_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     374                   ENDDO
     375                ENDDO
     376             ENDDO
     377
     378          CASE ( 'rad_sw_in' )
     379             DO  i = nxlg, nxrg
     380                DO  j = nysg, nyng
     381                   DO  k = nzb, nzt+1
     382                      rad_sw_in_av(k,j,i) = rad_sw_in_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     383                   ENDDO
     384                ENDDO
     385             ENDDO
     386
     387          CASE ( 'rad_sw_out' )
     388             DO  i = nxlg, nxrg
     389                DO  j = nysg, nyng
     390                   DO  k = nzb, nzt+1
     391                      rad_sw_out_av(k,j,i) = rad_sw_out_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     392                   ENDDO
    362393                ENDDO
    363394             ENDDO
Note: See TracChangeset for help on using the changeset viewer.