Ignore:
Timestamp:
Jan 12, 2018 4:03:39 PM (4 years ago)
Author:
suehring
Message:

In case of natural- and urban-type surfaces output surfaces fluxes in W/m2

File:
1 edited

Legend:

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

    r2742 r2743  
    2525! -----------------
    2626! $Id$
     27! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
     28!
     29! 2742 2018-01-12 14:59:47Z suehring
    2730! Enable output of surface temperature
    2831!
     
    194197
    195198    USE arrays_3d,                                                             &
    196         ONLY:  dzw, e, nc, nr, p, pt, precipitation_amount, precipitation_rate,&
     199        ONLY:  dzw, e, heatflux_output_conversion, nc, nr, p, pt,              &
     200               precipitation_amount, precipitation_rate,                       &
    197201               prr, q, qc, ql, ql_c, ql_v, ql_vp, qr, rho_ocean, s, sa,        &
    198                tend, u, v, vpt, w, zu, zw
     202               tend, u, v, vpt, w, zu, zw, waterflux_output_conversion
    199203       
    200204    USE averaging
    201205       
    202206    USE cloud_parameters,                                                      &
    203         ONLY:  hyrho, l_d_cp, pt_d_t
     207        ONLY:  cp, hyrho, l_d_cp, l_v, pt_d_t
    204208               
    205209    USE control_parameters,                                                    &
     
    809813             CASE ( 'qsws*_xy' )        ! 2d-array
    810814                IF ( av == 0 ) THEN
     815!
     816!--                In case of default surfaces, clean-up flux by density.
     817!--                In case of land- and urban-surfaces, convert fluxes into
     818!--                dynamic units
    811819                   DO  m = 1, surf_def_h(0)%ns
    812820                      i = surf_def_h(0)%i(m)
    813821                      j = surf_def_h(0)%j(m)
    814                       local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m)
     822                      k = surf_def_h(0)%k(m)
     823                      local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m) *            &
     824                                            waterflux_output_conversion(k)
    815825                   ENDDO
    816826                   DO  m = 1, surf_lsm_h%ns
    817827                      i = surf_lsm_h%i(m)
    818828                      j = surf_lsm_h%j(m)
    819                       local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m)
     829                      k = surf_lsm_h%k(m)
     830                      local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) * l_v
    820831                   ENDDO
    821832                   DO  m = 1, surf_usm_h%ns
    822833                      i = surf_usm_h%i(m)
    823834                      j = surf_usm_h%j(m)
    824                       local_pf(i,j,nzb+1) = surf_usm_h%qsws(m)
     835                      k = surf_usm_h%k(m)
     836                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v
    825837                   ENDDO
    826838                ELSE
     
    901913             CASE ( 'shf*_xy' )        ! 2d-array
    902914                IF ( av == 0 ) THEN
     915!
     916!--                In case of default surfaces, clean-up flux by density.
     917!--                In case of land- and urban-surfaces, convert fluxes into
     918!--                dynamic units.
    903919                   DO  m = 1, surf_def_h(0)%ns
    904920                      i = surf_def_h(0)%i(m)
    905921                      j = surf_def_h(0)%j(m)
    906                       local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m)
     922                      k = surf_def_h(0)%k(m)
     923                      local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m) *             &
     924                                            heatflux_output_conversion(k)
    907925                   ENDDO
    908926                   DO  m = 1, surf_lsm_h%ns
    909927                      i = surf_lsm_h%i(m)
    910928                      j = surf_lsm_h%j(m)
    911                       local_pf(i,j,nzb+1) = surf_lsm_h%shf(m)
     929                      k = surf_lsm_h%k(m)
     930                      local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) * cp
    912931                   ENDDO
    913932                   DO  m = 1, surf_usm_h%ns
    914933                      i = surf_usm_h%i(m)
    915934                      j = surf_usm_h%j(m)
    916                       local_pf(i,j,nzb+1) = surf_usm_h%shf(m)
     935                      k = surf_usm_h%k(m)
     936                      local_pf(i,j,nzb+1) = surf_usm_h%shf(m) * cp
    917937                   ENDDO
    918938                ELSE
Note: See TracChangeset for help on using the changeset viewer.