Changeset 1738 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Dec 18, 2015 1:56:05 PM (8 years ago)

Author:

raasch

Message:

bugfixes for calculations in statistical regions which do not contain grid points in the lowest vertical levels, mean surface level height considered in the calculation of the characteristic vertical velocity

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (10 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! calculate mean surface level height for each statistic region
 !
 ! Former revisions:
@@ -273 +273 @@
     
     USE statistics,                                                            &
-        ONLY:  hom, hom_sum, pr_palm, rmask, spectrum_x, spectrum_y,           &
-               statistic_regions, sums, sums_divnew_l, sums_divold_l, sums_l,  &
-               sums_l_l, sums_up_fraction_l, sums_wsts_bc_l, ts_value,         &
-               var_d, weight_pres, weight_substep 
+        ONLY:  hom, hom_sum, mean_surface_level_height, pr_palm, rmask,        &
+               spectrum_x, spectrum_y, statistic_regions, sums, sums_divnew_l, &
+               sums_divold_l, sums_l, sums_l_l, sums_up_fraction_l,            &
+               sums_wsts_bc_l, ts_value, var_d, weight_pres, weight_substep
  
     USE surface_layer_fluxes_mod,                                              &
@@ -299 +299 @@
     REAL(wp), DIMENSION(1:2) ::  volume_flow_initial_l  !<
 
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  mean_surface_level_height_l    !<
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_l    !<
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_tmp  !<
@@ -306 +307 @@
 !
 !-- Allocate arrays
-    ALLOCATE( ngp_2dh(0:statistic_regions), ngp_2dh_l(0:statistic_regions), &
+    ALLOCATE( mean_surface_level_height(0:statistic_regions),               &
+              mean_surface_level_height_l(0:statistic_regions),             &
+              ngp_2dh(0:statistic_regions), ngp_2dh_l(0:statistic_regions), &
               ngp_3d(0:statistic_regions),                                  &
               ngp_3d_inner(0:statistic_regions),                            &
@@ -1820 +1823 @@
 !
 !-- Compute total sum of active mask grid points
+!-- and the mean surface level height for each statistic region
 !-- ngp_2dh: number of grid points of a horizontal cross section through the
 !--          total domain
@@ -1834 +1838 @@
     ngp_sums          = ( nz + 2 ) * ( pr_palm + max_pr_user )
 
+    mean_surface_level_height   = 0.0_wp
+    mean_surface_level_height_l = 0.0_wp
+
     DO  sr = 0, statistic_regions
        DO  i = nxl, nxr
@@ -1841 +1848 @@
 !--             All xy-grid points
                 ngp_2dh_l(sr) = ngp_2dh_l(sr) + 1
+                mean_surface_level_height_l(sr) = mean_surface_level_height_l(sr) + &
+                                                  zw(nzb_s_inner(j,i))
 !
 !--             xy-grid points above topography
@@ -1873 +1882 @@
                         MPI_SUM, comm2d, ierr )
     ngp_3d_inner = INT( ngp_3d_inner_tmp, KIND = SELECTED_INT_KIND( 18 ) )
+    IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+    CALL MPI_ALLREDUCE( mean_surface_level_height_l(0),                       &
+                        mean_surface_level_height(0), sr, MPI_REAL,           &
+                        MPI_SUM, comm2d, ierr )
+    mean_surface_level_height = mean_surface_level_height / REAL( ngp_2dh )
 #else
     ngp_2dh         = ngp_2dh_l
@@ -1878 +1892 @@
     ngp_2dh_s_inner = ngp_2dh_s_inner_l
     ngp_3d_inner    = INT( ngp_3d_inner_l, KIND = SELECTED_INT_KIND( 18 ) )
+    mean_surface_level_height = mean_surface_level_height_l / REAL( ngp_2dh_l )
 #endif
 
@@ -1892 +1907 @@
     ngp_2dh_s_inner = MAX( 1, ngp_2dh_s_inner(:,:) ) 
 
-    DEALLOCATE( ngp_2dh_l, ngp_2dh_outer_l, ngp_3d_inner_l, ngp_3d_inner_tmp )
+    DEALLOCATE( mean_surface_level_height_l, ngp_2dh_l, ngp_2dh_outer_l,      &
+                ngp_3d_inner_l, ngp_3d_inner_tmp )
 
     CALL location_message( 'leaving init_3d_model', .TRUE. )