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

Added support for RRTMG radiation code

File:
1 edited

Legend:

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

    r1556 r1585  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added support for RRTMG
    2323!
    2424! Former revisions:
     
    154154
    155155    USE radiation_model_mod,                                                   &
    156         ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av
     156        ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av, rad_sw_out,       &
     157               rad_sw_out_av, rad_lw_in, rad_lw_in_av, rad_lw_out,             &
     158               rad_lw_out_av
    157159
    158160    IMPLICIT NONE
     
    913915                level_z(nzb+1) = zu(nzb+1)
    914916
    915              CASE ( 'rad_sw_in*_xy' )        ! 2d-array
    916                 IF ( av == 0 ) THEN
    917                    DO  i = nxlg, nxrg
    918                       DO  j = nysg, nyng
    919                          local_pf(i,j,nzb+1) =  rad_sw_in(j,i)
    920                       ENDDO
    921                    ENDDO
    922                 ELSE
    923                    DO  i = nxlg, nxrg
    924                       DO  j = nysg, nyng
    925                          local_pf(i,j,nzb+1) =  rad_sw_in_av(j,i)
    926                       ENDDO
    927                    ENDDO
    928                 ENDIF
    929                 resorted = .TRUE.
    930                 two_d = .TRUE.
    931                 level_z(nzb+1) = zu(nzb+1)
     917
     918             CASE ( 'rad_lw_in_xy', 'rad_lw_in_xz', 'rad_lw_in_yz' )
     919                IF ( av == 0 )  THEN
     920                   to_be_resorted => rad_lw_in
     921                ELSE
     922                   to_be_resorted => rad_lw_in_av
     923                ENDIF
     924
     925             CASE ( 'rad_lw_out_xy', 'rad_lw_out_xz', 'rad_lw_out_yz' )
     926                IF ( av == 0 )  THEN
     927                   to_be_resorted => rad_lw_out
     928                ELSE
     929                   to_be_resorted => rad_lw_out_av
     930                ENDIF
     931
     932             CASE ( 'rad_sw_in_xy', 'rad_sw_in_xz', 'rad_sw_in_yz' )
     933                IF ( av == 0 )  THEN
     934                   to_be_resorted => rad_sw_in
     935                ELSE
     936                   to_be_resorted => rad_sw_in_av
     937                ENDIF
     938
     939             CASE ( 'rad_sw_out_xy', 'rad_sw_out_xz', 'rad_sw_out_yz' )
     940                IF ( av == 0 )  THEN
     941                   to_be_resorted => rad_sw_out
     942                ELSE
     943                   to_be_resorted => rad_sw_out_av
     944                ENDIF
    932945
    933946             CASE ( 'rho_xy', 'rho_xz', 'rho_yz' )
Note: See TracChangeset for help on using the changeset viewer.