Changeset 2977 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Apr 17, 2018 10:27:57 AM (6 years ago)

Author:

kanani

Message:

Fixes for radiative transfer model

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Implement changes from branch radiation (r2948-2971) with minor modifications,
+! plus some formatting.
+! (moh.hefny)
+! Added flag to check the existence of vertical urban/land surfaces, required
+! to activate RTM
+! 
+! 2970 2018-04-13 15:09:23Z suehring
 ! Remove un-necessary initialization of surface elements in old large-scale 
 ! forcing mode
@@ -450 +457 @@
     INTEGER(iwp) ::  ns_v_on_file(0:3)                       !< total number of vertical surfaces with the same facing, required for writing restart data 
 
+    LOGICAL ::  vertical_surfaces_exist = .FALSE.   !< flag indicating that there are vertical urban/land surfaces 
+                                                    !< in the domain (required to activiate RTM)
+
 
     SAVE
@@ -488 +498 @@
     PUBLIC bc_h, ind_pav_green, ind_veg_wall, ind_wat_win, ns_h_on_file,       &
            ns_v_on_file, surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v,       &
-           surf_usm_h, surf_usm_v, surf_type
+           surf_usm_h, surf_usm_v, surf_type, vertical_surfaces_exist
 !
 !-- Public subroutines and functions
@@ -615 +625 @@
     SUBROUTINE init_surface_arrays
 
+
+       USE pegrid
+
+
        IMPLICIT NONE
 
@@ -628 +642 @@
        INTEGER(iwp), DIMENSION(0:3) ::  num_lsm_v !< number of vertically-aligned natural surfaces 
        INTEGER(iwp), DIMENSION(0:3) ::  num_usm_v !< number of vertically-aligned urban surfaces 
+
+       INTEGER(iwp)              ::  num_surf_v_l !< number of vertically-aligned local urban/land surfaces
+       INTEGER(iwp)              ::  num_surf_v   !< number of vertically-aligned total urban/land surfaces
 
        LOGICAL ::  building                       !< flag indicating building grid point
@@ -920 +937 @@
                                                nys, nyn, nxl, nxr )
        ENDDO
-
+!
+!--    Set the flag for the existence of vertical urban/land surfaces
+       num_surf_v_l = 0
+       DO  l = 0, 3
+          num_surf_v_l = num_surf_v_l + surf_usm_v(l)%ns + surf_lsm_v(l)%ns
+       ENDDO
+
+#if defined( __parallel )
+       CALL MPI_ALLREDUCE( num_surf_v_l, num_surf_v, 1, MPI_INTEGER, MPI_SUM, comm2d, ierr)
+#else
+       num_surf_v = num_surf_v_l
+#endif
+       IF ( num_surf_v > 0 ) vertical_surfaces_exist = .TRUE.
+        
 
     END SUBROUTINE init_surface_arrays