Changeset 2743 for palm/trunk/SOURCE/data_output_2d.f90

Timestamp:

Jan 12, 2018 4:03:39 PM (6 years ago)

Author:

suehring

Message:

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

File:

• palm/trunk/SOURCE/data_output_2d.f90 (modified) (4 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
+! 
+! 2742 2018-01-12 14:59:47Z suehring
 ! Enable output of surface temperature
 ! 
@@ -194 +197 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dzw, e, nc, nr, p, pt, precipitation_amount, precipitation_rate,&
+        ONLY:  dzw, e, heatflux_output_conversion, nc, nr, p, pt,              &
+               precipitation_amount, precipitation_rate,                       &
                prr, q, qc, ql, ql_c, ql_v, ql_vp, qr, rho_ocean, s, sa,        &
-               tend, u, v, vpt, w, zu, zw
+               tend, u, v, vpt, w, zu, zw, waterflux_output_conversion
         
     USE averaging
         
     USE cloud_parameters,                                                      &
-        ONLY:  hyrho, l_d_cp, pt_d_t
+        ONLY:  cp, hyrho, l_d_cp, l_v, pt_d_t
                
     USE control_parameters,                                                    &
@@ -809 +813 @@
              CASE ( 'qsws*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN
+!
+!--                In case of default surfaces, clean-up flux by density.
+!--                In case of land- and urban-surfaces, convert fluxes into
+!--                dynamic units
                    DO  m = 1, surf_def_h(0)%ns
                       i = surf_def_h(0)%i(m)
                       j = surf_def_h(0)%j(m)
-                      local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m)
+                      k = surf_def_h(0)%k(m)
+                      local_pf(i,j,nzb+1) = surf_def_h(0)%qsws(m) *            &
+                                            waterflux_output_conversion(k)
                    ENDDO
                    DO  m = 1, surf_lsm_h%ns
                       i = surf_lsm_h%i(m)
                       j = surf_lsm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m)
+                      k = surf_lsm_h%k(m)
+                      local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) * l_v
                    ENDDO
                    DO  m = 1, surf_usm_h%ns
                       i = surf_usm_h%i(m)
                       j = surf_usm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m)
+                      k = surf_usm_h%k(m)
+                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v
                    ENDDO
                 ELSE
@@ -901 +913 @@
              CASE ( 'shf*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN
+!
+!--                In case of default surfaces, clean-up flux by density.
+!--                In case of land- and urban-surfaces, convert fluxes into
+!--                dynamic units.
                    DO  m = 1, surf_def_h(0)%ns
                       i = surf_def_h(0)%i(m)
                       j = surf_def_h(0)%j(m)
-                      local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m)
+                      k = surf_def_h(0)%k(m)
+                      local_pf(i,j,nzb+1) = surf_def_h(0)%shf(m) *             &
+                                            heatflux_output_conversion(k)
                    ENDDO
                    DO  m = 1, surf_lsm_h%ns
                       i = surf_lsm_h%i(m)
                       j = surf_lsm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_lsm_h%shf(m)
+                      k = surf_lsm_h%k(m)
+                      local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) * cp
                    ENDDO
                    DO  m = 1, surf_usm_h%ns
                       i = surf_usm_h%i(m)
                       j = surf_usm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%shf(m)
+                      k = surf_usm_h%k(m)
+                      local_pf(i,j,nzb+1) = surf_usm_h%shf(m) * cp
                    ENDDO
                 ELSE