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

Timestamp:

May 9, 2019 9:48:32 AM (5 years ago)

Author:

suehring

Message:

In a nested child domain, distinguish between soil moisture and temperature initialization in case the parent domain is initialized via the dynamic input file; in the offline nesting, add a safety factor for the calculation of bulk Richardson number in order to avoid division by zero which can potentially happen if 3D buildings are located directly at the lateral model boundaries

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! In a nested child domain, distinguish between soil moisture and temperature 
+! initialization from parent via dynamic input file. Further, initialize soil
+! moisture/temperature from dynamic input file only when initialization via
+! 'inifor' is desired.
+! 
+! 3943 2019-05-02 09:50:41Z maronga
 ! Removed extra blank character
 ! 
@@ -2636 +2642 @@
        IMPLICIT NONE
 
-       LOGICAL      ::  init_soil_from_parent   !< flag controlling initialization of soil in child domains
+       LOGICAL      ::  init_msoil_from_parent   !< flag controlling initialization of soil moisture in nested child domains
+       LOGICAL      ::  init_tsoil_from_parent   !< flag controlling initialization of soil temperature in nested child domains
 
        INTEGER(iwp) ::  i                       !< running index
@@ -4356 +4363 @@
        IF (  TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
 !
-!--       Read soil properties from dynamic input file
-          CALL netcdf_data_input_init_lsm
+!--       Read soil properties from dynamic input file. 
+          IF ( INDEX( initializing_actions, 'inifor' ) /= 0 )                  &
+             CALL netcdf_data_input_init_lsm
 !
 !--       In case no dynamic input is available for a child domain but root 
@@ -4363 +4371 @@
 !--       properties can lead to significant discrepancies in the atmospheric 
 !--       surface forcing. For this reason, the child domain
-!--       is initialized with mean soil profiles from the root domain.          
-          init_soil_from_parent = .FALSE. 
+!--       is initialized with mean soil profiles from the root domain, even if
+!--       no initialization with inifor is set. 
+          init_tsoil_from_parent = .FALSE. 
+          init_msoil_from_parent = .FALSE. 
           IF ( nested_run )  THEN
 #if defined( __parallel )
-             CALL MPI_ALLREDUCE( init_3d%from_file_tsoil  .OR.                 &
-                                 init_3d%from_file_msoil,                      &
-                                 init_soil_from_parent,                        &
+             CALL MPI_ALLREDUCE( init_3d%from_file_tsoil,                      &
+                                 init_tsoil_from_parent,                       &
+                                 1, MPI_LOGICAL, MPI_LOR, MPI_COMM_WORLD, ierr )
+             CALL MPI_ALLREDUCE( init_3d%from_file_msoil,                      &
+                                 init_msoil_from_parent,                       &
                                  1, MPI_LOGICAL, MPI_LOR, MPI_COMM_WORLD, ierr )
 #endif
@@ -4407 +4419 @@
 !--       transfer soil mean profiles from the root-parent domain onto all 
 !--       child domains.
-          IF ( init_soil_from_parent )  THEN
+          IF ( init_msoil_from_parent )  THEN
 !
 !--          Child domains will be only initialized with horizontally 
@@ -4434 +4446 @@
                 DEALLOCATE( pr_soil_init )
              ENDIF
+          ENDIF
+          IF ( init_tsoil_from_parent )  THEN
              IF ( init_3d%from_file_tsoil  .AND.  init_3d%lod_tsoil == 2 )  THEN
                 ALLOCATE( pr_soil_init(0:init_3d%nzs-1) )
@@ -4455 +4469 @@
 
              ENDIF
-
-#if defined( __parallel )
+          ENDIF
+          IF ( init_msoil_from_parent  .OR.  init_tsoil_from_parent )  THEN
 !
 !--          Distribute soil grid information on file from root to all childs.
 !--          Only process with rank 0 sends the information. 
+#if defined( __parallel )
              CALL MPI_BCAST( init_3d%nzs,    1,                                &
                              MPI_INTEGER, 0, MPI_COMM_WORLD, ierr )
+#endif
 
              IF ( .NOT.  ALLOCATED( init_3d%z_soil ) )                         &
                 ALLOCATE( init_3d%z_soil(1:init_3d%nzs) )
-
+#if defined( __parallel )
              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.
-!--          Note, 1d soil profiles are allocated even though soil information
-!--          is already read from dynamic file in one child domain. 
-!--          This case, however, data is not used for further initialization 
-!--          since the LoD=2. 
+          ENDIF
+!
+!--       ALLOCATE arrays on child domains and set control attributes.
+!--       Note, 1d soil profiles are allocated even though soil information
+!--       is already read from dynamic file in one child domain. 
+!--       This case, however, data is not used for further initialization 
+!--       since the LoD=2. 
+          IF ( init_msoil_from_parent )  THEN
              IF ( .NOT. ALLOCATED( init_3d%msoil_1d ) )  THEN
                 ALLOCATE( init_3d%msoil_1d(0:init_3d%nzs-1) )
@@ -4480 +4498 @@
                 init_3d%from_file_msoil = .TRUE.
              ENDIF
+          ENDIF
+          IF ( init_tsoil_from_parent )  THEN
              IF ( .NOT. ALLOCATED( init_3d%tsoil_1d ) )  THEN
                 ALLOCATE( init_3d%tsoil_1d(0:init_3d%nzs-1) )
@@ -4485 +4505 @@
                 init_3d%from_file_tsoil = .TRUE.
              ENDIF
-
-
+          ENDIF
+!
+!--       Distribute soil profiles from root to all childs
+          IF ( init_msoil_from_parent )  THEN
 #if defined( __parallel )
-!
-!--          Distribute soil profiles from root to all childs
              CALL MPI_BCAST( init_3d%msoil_1d, SIZE(init_3d%msoil_1d),         &
                              MPI_REAL, 0, MPI_COMM_WORLD, ierr )
+#endif
+                    
+          ENDIF
+          IF ( init_tsoil_from_parent )  THEN
+#if defined( __parallel )
              CALL MPI_BCAST( init_3d%tsoil_1d, SIZE(init_3d%tsoil_1d),         &
                              MPI_REAL, 0, MPI_COMM_WORLD, ierr )
 #endif
-
           ENDIF
 !