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

Timestamp:

Aug 8, 2019 8:00:47 PM (6 years ago)

Author:

suehring

Message:

Framework to input single surface variables independent on land/urban-surface model via static input file provided. This way, input and initialization of heterogeneous roughness was already added.

File:

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

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Generic routine to initialize single surface properties added
+! 
+! 4104 2019-07-17 17:08:20Z suehring
 ! Bugfix, initialization of index space for boundary data structure accidantly
 ! run over ghost points, causing a segmentation fault.
@@ -642 +645 @@
     END INTERFACE init_bc
 
+    INTERFACE init_single_surface_properties
+       MODULE PROCEDURE init_single_surface_properties
+    END INTERFACE init_single_surface_properties
+    
     INTERFACE init_surfaces
        MODULE PROCEDURE init_surfaces
@@ -685 +692 @@
 !
 !-- Public subroutines and functions
-    PUBLIC get_topography_top_index, get_topography_top_index_ji, init_bc, init_surfaces,          &
-           init_surface_arrays, surface_rrd_local, surface_restore_elements, surface_wrd_local,    &
-           surface_last_actions
+    PUBLIC get_topography_top_index,                                                               &
+           get_topography_top_index_ji,                                                            &
+           init_bc,                                                                                &
+           init_single_surface_properties,                                                         &
+           init_surfaces,                                                                          &
+           init_surface_arrays,                                                                    &
+           surface_last_actions,                                                                   &
+           surface_rrd_local,                                                                      &
+           surface_restore_elements,                                                               &
+           surface_wrd_local
 
 #if defined( _OPENACC )
@@ -3035 +3049 @@
     END SUBROUTINE init_surfaces
 
+! Description:
+! ------------
+!> Initialize single surface properties from 2D input arrays
+!------------------------------------------------------------------------------!
+    SUBROUTINE init_single_surface_properties( var_surf, var_2d,               &
+                                               ns, fill_value,                 &
+                                               index_space_i,                  &
+                                               index_space_j                   &
+                                             )
+    
+       INTEGER(iwp) ::  m  !< running index over surface elements
+       INTEGER(iwp) ::  ns !< number of surface elements in var_surf
+
+       INTEGER(iwp), DIMENSION(1:ns) ::  index_space_i !< grid indices along x direction where surface properties should be defined
+       INTEGER(iwp), DIMENSION(1:ns) ::  index_space_j !< grid indices along y direction where surface properties should be defined 
+       
+       REAL(wp) ::  fill_value !< fill value in var_2d
+       
+       REAL(wp), DIMENSION(1:ns) ::  var_surf !< 1D surface variable that should be initialized
+       REAL(wp), DIMENSION(nys:nyn,nxl:nxr) ::  var_2d !< input variable
+
+       DO  m = 1, ns
+          IF ( var_2d(index_space_j(m),index_space_i(m)) /= fill_value )  THEN
+             var_surf(m) = var_2d(index_space_j(m),index_space_i(m))
+          ENDIF
+       ENDDO
+       
+    END SUBROUTINE init_single_surface_properties
 
 !------------------------------------------------------------------------------!