Changeset 3176

Timestamp:

Jul 26, 2018 5:12:48 PM (6 years ago)

Author:

suehring

Message:

Major bugfix in calculation of ol and ts at building roofs; bugfix in restart data for surface elements; in 2D data output, mask latent heat flux at urban-type surfaces

Location:

palm/trunk/SOURCE

Files:

• data_output_2d.f90 (modified) (3 diffs)
• sum_up_3d_data.f90 (modified) (2 diffs)
• surface_mod.f90 (modified) (2 diffs)
• time_integration.f90 (modified) (1 diff)
• urban_surface_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove output of latent heat flux at urban-surfaces and set fill values 
+! instead 
+! 
+! 3052 2018-05-31 06:11:20Z raasch
 ! Do not open FORTRAN binary files in case of parallel netCDF I/O
 ! 
@@ -933 +937 @@
              CASE ( 'qsws*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN
+                   local_pf(:,:,nzb+1) = REAL( fill_value, KIND = wp )
 !
 !--                In case of default surfaces, clean-up flux by density.
-!--                In case of land- and urban-surfaces, convert fluxes into
+!--                In case of land-surfaces, convert fluxes into
 !--                dynamic units
                    DO  m = 1, surf_def_h(0)%ns
@@ -949 +954 @@
                       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)
-                      k = surf_usm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v
                    ENDDO
                 ELSE

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove output of latent heat flux at urban-surfaces and set fill values 
+! instead 
+! 
+! 3173 2018-07-26 12:55:23Z suehring
 ! Bugfix for last commit
 ! 
@@ -894 +898 @@
                          qsws_av(j,i) = qsws_av(j,i) +                         &
                                          surf_lsm_h%qsws(m) * l_v
-                      ELSEIF ( match_usm )  THEN
-                         m = surf_usm_h%end_index(j,i)
-                         qsws_av(j,i) = qsws_av(j,i) +                         &
-                                         surf_usm_h%qsws(m) * l_v
                       ENDIF
                    ENDDO

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix in restarts for surface elements.
+! 
+! 3152 2018-07-19 13:26:52Z suehring
 ! q_surface added
 ! 
@@ -2690 +2693 @@
                          surf_h(l)%pt_surface(mm(l)) = surf_lsm_h%pt_surface(m)
                       IF ( ALLOCATED( surf_def_h(l)%q_surface ) )              &
-                         surf_h(l)%q_surface(mm(l)) = surf_def_h(l)%q_surface(m)
+                         surf_h(l)%q_surface(mm(l)) = surf_lsm_h%q_surface(m)
                       IF ( ALLOCATED( surf_def_h(l)%vpt_surface ) )            &
-                         surf_h(l)%vpt_surface(mm(l)) = surf_def_h(l)%vpt_surface(m)
+                         surf_h(l)%vpt_surface(mm(l)) = surf_lsm_h%vpt_surface(m)
                       IF ( ALLOCATED( surf_lsm_h%usws ) )                      &
                          surf_h(0)%usws(mm(0))    = surf_lsm_h%usws(m)

palm/trunk/SOURCE/time_integration.f90

@@ -1015 +1015 @@
              IF (urban_surface) THEN
                 CALL cpu_log( log_point(74), 'urban_surface', 'start' )
+                
                 CALL usm_surface_energy_balance
                 IF ( usm_material_model )  THEN

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix, update virtual potential surface temparture, else heat fluxes on 
+! roofs might become unphysical
+! 
+! 3152 2018-07-19 13:26:52Z suehring
 ! Initialize q_surface, which might be used in surface_layer_fluxes
 ! 
@@ -7438 +7442 @@
                            surf_usm_h%tt_surface_green_m(m)
 !
-!--        Store surface temperature        
+!--        Store surface temperature on pt_surface. Further, in case humidity is used
+!--        store also vpt_surface, which is, due to the lack of moisture on roofs simply
+!--        assumed to be the surface temperature.
            surf_usm_h%pt_surface(m) = ( surf_usm_h%frac(ind_veg_wall,m) * t_surf_h_p(m)   &
                                + surf_usm_h%frac(ind_wat_win,m) * t_surf_window_h_p(m)   &
                                + surf_usm_h%frac(ind_pav_green,m) * t_surf_green_h_p(m) )  &
                                / exn(k)
+                               
+           IF ( humidity )  surf_usm_h%vpt_surface(m) =                        &
+                                                   surf_usm_h%pt_surface(m)
                        
 !--        calculate true tendency
@@ -7637 +7646 @@
                                 surf_usm_v(l)%tt_surface_green_m(m)
 !
-!--           Store surface temperature        
+!--           Store surface temperature. Further, in case humidity is used
+!--           store also vpt_surface, which is, due to the lack of moisture on roofs simply
+!--           assumed to be the surface temperature.     
               surf_usm_v(l)%pt_surface(m) =  ( surf_usm_v(l)%frac(ind_veg_wall,m) * t_surf_v_p(l)%t(m)  &
                                       + surf_usm_v(l)%frac(ind_wat_win,m) * t_surf_window_v_p(l)%t(m)  &
                                       + surf_usm_v(l)%frac(ind_pav_green,m) * t_surf_green_v_p(l)%t(m) ) &
                                       / exn(k) 
+                                      
+              IF ( humidity )  surf_usm_v(l)%vpt_surface(m) =                  &
+                                                    surf_usm_v(l)%pt_surface(m)
             
 !--           calculate true tendency