Changeset 4771 for palm/trunk/SOURCE/pmc_handle_communicator_mod.f90

Timestamp:

Nov 4, 2020 1:12:46 PM (3 years ago)

Author:

hellstea

Message:

Canopy-restricted anterpolation introduced

File:

• palm/trunk/SOURCE/pmc_handle_communicator_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/pmc_handle_communicator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Canopy-restricted anterpolation introduced. New namelist parameter anterpolation_starting_height 
+! introduced for controlling canopy-restricted anterpolation.
+!
+! 4650 2020-08-25 14:35:50Z raasch
 ! bugfix for r4649
 !
@@ -138 +142 @@
 !--------------------------------------------------------------------------------------------------!
  SUBROUTINE pmc_init_model( comm, nesting_datatransfer_mode, nesting_mode,                         &
-                            anterpolation_buffer_width, pmc_status )
+                            anterpolation_buffer_width, anterpolation_starting_height, pmc_status )
 
     USE control_parameters,                                                                        &
@@ -155 +159 @@
     INTEGER, INTENT(INOUT) ::  pmc_status                  !<
 
+    REAL(wp), INTENT(INOUT) ::  anterpolation_starting_height  !< steering parameter for canopy restricted anterpolation 
+    
     INTEGER ::  childcount     !<
     INTEGER ::  i              !<
@@ -180 +186 @@
 
        CALL read_coupling_layout( nesting_datatransfer_mode, nesting_mode,                         &
-                                  anterpolation_buffer_width, pmc_status )
+            anterpolation_buffer_width, anterpolation_starting_height,                             &
+            pmc_status )
 
        IF ( pmc_status /= pmc_no_namelist_found  .AND.                                             &
@@ -252 +259 @@
                     MPI_COMM_WORLD, istat )
     CALL MPI_BCAST( anterpolation_buffer_width, 1, MPI_INT, 0, MPI_COMM_WORLD, istat )
+    CALL MPI_BCAST( anterpolation_starting_height, 1, MPI_REAL, 0, MPI_COMM_WORLD, istat )
 !
 !-- Assign global MPI processes to individual models by setting the couple id
@@ -422 +430 @@
 !--------------------------------------------------------------------------------------------------!
  SUBROUTINE read_coupling_layout( nesting_datatransfer_mode, nesting_mode,                         &
-                                  anterpolation_buffer_width, pmc_status )
+                                  anterpolation_buffer_width, anterpolation_starting_height,       &
+                                  pmc_status )
 
     IMPLICIT NONE
@@ -431 +440 @@
     INTEGER, INTENT(INOUT)      ::  anterpolation_buffer_width  !< Boundary buffer width for anterpolation
     INTEGER(iwp), INTENT(INOUT) ::  pmc_status                  !<
+
+    REAL(wp), INTENT(INOUT) ::  anterpolation_starting_height   !< steering parameter for canopy restricted anterpolation
 
     INTEGER(iwp) ::  bad_llcorner  !<
@@ -441 +452 @@
                                    nesting_datatransfer_mode,                                      &
                                    nesting_mode,                                                   &
-                                   anterpolation_buffer_width
+                                   anterpolation_buffer_width,                                     &
+                                   anterpolation_starting_height
 
 !