Changeset 2938 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Mar 27, 2018 3:52:42 PM (6 years ago)

Author:

suehring

Message:

Nesting in RANS-LES and RANS-RANS mode enabled; synthetic turbulence generator at all lateral boundaries in nesting or non-cyclic forcing mode; revised Inifor initialization in nesting mode

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Initialization of soil moisture and temperature via Inifor-provided data also
+! in nested child domains, even if no dynamic input file is available for 
+! child domain. 1D soil profiles are received from parent model.  
+! 
+! 2936 2018-03-27 14:49:27Z suehring
 ! renamed lsm_par to land_surface_parameters. Bugfix in message calls
 ! 
@@ -2365 +2370 @@
        USE control_parameters,                                                 &
            ONLY:  message_string
+
+       USE indices,                                                            &
+           ONLY:  nx, ny
+
+       USE pmc_interface,                                                      &
+           ONLY:  nested_run
     
        IMPLICIT NONE
+
+       LOGICAL      ::  init_soil_dynamically_in_child !< flag controlling initialization of soil in child domains
 
        INTEGER(iwp) ::  i                       !< running index
@@ -2384 +2397 @@
 
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  bound, bound_root_fr  !< temporary arrays for storing index bounds
+       REAL(wp), DIMENSION(:), ALLOCATABLE ::  pr_soil_init !< temporary array used for averaging soil profiles
 
 !
@@ -3894 +3908 @@
        IF (  TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
 !
+!--       In case of nested runs, check if soil is initialized via dynamic
+!--       input file in root domain and distribute this information
+!--       to all embedded child domains. This case, soil moisture and
+!--       temperature will be initialized via root domain.  
+          init_soil_dynamically_in_child = .FALSE. 
+          IF ( nested_run )  THEN
+#if defined( __parallel )
+             CALL MPI_ALLREDUCE( init_3d%from_file_tsoil  .OR.                 &
+                                 init_3d%from_file_msoil,                      &
+                                 init_soil_dynamically_in_child,               &
+                                 1, MPI_LOGICAL, MPI_LOR, MPI_COMM_WORLD, ierr )
+#endif
+          ENDIF
+
+!
 !--       First, initialize soil temperature and moisture. 
 !--       According to the initialization for surface and soil parameters, 
@@ -3923 +3952 @@
           ENDDO
 !
-!--       Level 2, if soil moisture and/or temperature  are 
-!--       provided from file, interpolate / extrapolate the provided data
-!--       onto the respective soil layers. Please note, both, zs as well as 
-!--       init_3d%z_soil indicate a depth with positive values, so that no 
-!--       distinction between atmosphere is required concerning interpolation.
-!--       Start with soil moisture
+!--       Initialization of soil moisture and temperature from file. 
+!--       In case of nested runs, only root parent reads dynamic input file. 
+!--       This case, transfer respective soil information provide
+!--       by dynamic input file from root parent domain onto all other domains.
+          IF ( init_soil_dynamically_in_child )  THEN
+!
+!--          Child domains will be only initialized with horizontally 
+!--          averaged soil profiles in parent domain (for sake of simplicity). 
+!--          If required, average soil data on root parent domain before 
+!--          distribute onto child domains.
+             IF ( init_3d%from_file_msoil  .AND.  init_3d%lod_msoil == 2 )     &
+             THEN
+                ALLOCATE( pr_soil_init(0:init_3d%nzs-1) )
+
+                DO  k = 0, init_3d%nzs-1
+                   pr_soil_init(k) = SUM( init_3d%msoil(k,nys:nyn,nxl:nxr)  )
+                ENDDO
+!
+!--             Allocate 1D array for soil-moisture profile (will not be 
+!--             allocated in lod==2 case). 
+                ALLOCATE( init_3d%msoil_init(0:init_3d%nzs-1) )
+                init_3d%msoil_init = 0.0_wp
+#if defined( __parallel )
+                CALL MPI_ALLREDUCE( pr_soil_init(0), init_3d%msoil_init(0),    &
+                                    SIZE(pr_soil_init),                        &
+                                    MPI_REAL, MPI_SUM, comm2d, ierr )
+#endif
+                init_3d%msoil_init = init_3d%msoil_init /                      &
+                                        REAL( ( nx + 1 ) * ( ny + 1), KIND=wp )
+                DEALLOCATE( pr_soil_init )
+             ENDIF
+             IF ( init_3d%from_file_tsoil  .AND.  init_3d%lod_tsoil == 2 )  THEN
+                ALLOCATE( pr_soil_init(0:init_3d%nzs-1) )
+
+                DO  k = 0, init_3d%nzs-1
+                   pr_soil_init(k) = SUM( init_3d%tsoil(k,nys:nyn,nxl:nxr)  )
+                ENDDO
+!
+!--             Allocate 1D array for soil-temperature profile (will not be 
+!--             allocated in lod==2 case). 
+                ALLOCATE( init_3d%tsoil_init(0:init_3d%nzs-1) )
+                init_3d%tsoil_init = 0.0_wp
+#if defined( __parallel )
+                CALL MPI_ALLREDUCE( pr_soil_init(0), init_3d%tsoil_init(0),    &
+                                    SIZE(pr_soil_init),                        &
+                                    MPI_REAL, MPI_SUM, comm2d, ierr )
+#endif
+                init_3d%tsoil_init = init_3d%tsoil_init /                      &
+                                        REAL( ( nx + 1 ) * ( ny + 1), KIND=wp )
+                DEALLOCATE( pr_soil_init )
+
+             ENDIF
+
+#if defined( __parallel )
+!
+!--          Distribute soil grid information on file from root to all childs.
+!--          Only process with rank 0 sends the information. 
+             CALL MPI_BCAST( init_3d%nzs,    1,                                &
+                             MPI_INTEGER, 0, MPI_COMM_WORLD, ierr )
+
+             IF ( .NOT.  ALLOCATED( init_3d%z_soil ) )                         &
+                ALLOCATE( init_3d%z_soil(1:init_3d%nzs) )
+
+             CALL MPI_BCAST( init_3d%z_soil, SIZE(init_3d%z_soil),             &
+                             MPI_REAL, 0, MPI_COMM_WORLD, ierr )
+#endif
+!
+!--          ALLOCATE arrays on child domains and set control attributes.
+             IF ( .NOT. init_3d%from_file_msoil )  THEN
+                ALLOCATE( init_3d%msoil_init(0:init_3d%nzs-1) )
+                init_3d%lod_msoil = 1
+                init_3d%from_file_msoil = .TRUE.
+             ENDIF
+             IF ( .NOT. init_3d%from_file_tsoil )  THEN
+                ALLOCATE( init_3d%tsoil_init(0:init_3d%nzs-1) )
+                init_3d%lod_tsoil = 1
+                init_3d%from_file_tsoil = .TRUE.
+             ENDIF
+
+
+#if defined( __parallel )
+!
+!--          Distribute soil profiles from root to all childs
+             CALL MPI_BCAST( init_3d%msoil_init, SIZE(init_3d%msoil_init),     &
+                             MPI_REAL, 0, MPI_COMM_WORLD, ierr )
+             CALL MPI_BCAST( init_3d%tsoil_init, SIZE(init_3d%tsoil_init),     &
+                             MPI_REAL, 0, MPI_COMM_WORLD, ierr )
+#endif
+
+          ENDIF
+!
+!--       Proceed with Level 2 initialization. Information from dynamic input
+!--       is now available on all processes. 
           IF ( init_3d%from_file_msoil )  THEN
 
@@ -3983 +4099 @@
                 ENDDO
              ENDIF
-
           ENDIF
 !