Ignore:
Timestamp:
Oct 26, 2015 4:17:44 PM (6 years ago)
Author:
maronga
Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:
1 edited

Legend:

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

    r1683 r1691  
    1414! PALM. If not, see <http://www.gnu.org/licenses/>.
    1515!
    16 ! Copyright 1997-2014 Leibniz Universitaet Hannover
     16! Copyright 1997-2015 Leibniz Universitaet Hannover
    1717!--------------------------------------------------------------------------------!
    1818!
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Added output of Obukhov length and radiative heating rates for RRTMG.
    2222!
    2323! Former revisions:
     
    9696    USE radiation_model_mod,                                                   &
    9797        ONLY:  rad_net, rad_net_av, rad_lw_in, rad_lw_in_av, rad_lw_out,       &
    98                rad_lw_out_av, rad_sw_in, rad_sw_in_av, rad_sw_out,             &
    99                rad_sw_out_av
     98               rad_lw_out_av, rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr,        &
     99               rad_lw_hr_av, rad_sw_in, rad_sw_in_av, rad_sw_out,              &
     100               rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr,        &
     101               rad_sw_hr_av
     102
    100103
    101104    IMPLICIT NONE
     
    213216             ENDDO
    214217
     218         CASE ( 'ol*' )
     219             DO  i = nxlg, nxrg
     220                DO  j = nysg, nyng
     221                   ol_av(j,i) = ol_av(j,i) / REAL( average_count_3d, KIND=wp )
     222                ENDDO
     223             ENDDO
     224
    215225          CASE ( 'p' )
    216226             DO  i = nxlg, nxrg
     
    383393             ENDDO
    384394
     395          CASE ( 'rad_lw_cs_hr' )
     396             DO  i = nxlg, nxrg
     397                DO  j = nysg, nyng
     398                   DO  k = nzb, nzt+1
     399                      rad_lw_cs_hr_av(k,j,i) = rad_lw_cs_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     400                   ENDDO
     401                ENDDO
     402             ENDDO
     403
     404          CASE ( 'rad_lw_hr' )
     405             DO  i = nxlg, nxrg
     406                DO  j = nysg, nyng
     407                   DO  k = nzb, nzt+1
     408                      rad_lw_hr_av(k,j,i) = rad_lw_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     409                   ENDDO
     410                ENDDO
     411             ENDDO
     412
    385413          CASE ( 'rad_sw_in' )
    386414             DO  i = nxlg, nxrg
     
    397425                   DO  k = nzb, nzt+1
    398426                      rad_sw_out_av(k,j,i) = rad_sw_out_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     427                   ENDDO
     428                ENDDO
     429             ENDDO
     430
     431          CASE ( 'rad_sw_cs_hr' )
     432             DO  i = nxlg, nxrg
     433                DO  j = nysg, nyng
     434                   DO  k = nzb, nzt+1
     435                      rad_sw_cs_hr_av(k,j,i) = rad_sw_cs_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     436                   ENDDO
     437                ENDDO
     438             ENDDO
     439
     440          CASE ( 'rad_sw_hr' )
     441             DO  i = nxlg, nxrg
     442                DO  j = nysg, nyng
     443                   DO  k = nzb, nzt+1
     444                      rad_sw_hr_av(k,j,i) = rad_sw_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    399445                   ENDDO
    400446                ENDDO
Note: See TracChangeset for help on using the changeset viewer.