Changeset 4361 for palm/trunk

Timestamp:

Jan 7, 2020 12:22:38 PM (5 years ago)

Author:

suehring

Message:

Plant-canopy: Avoid an unnecessary exchange of ghost points; remove unused arrays in pcm_rrd_local

Location:

palm/trunk/SOURCE

Files:

• module_interface.f90 (modified) (2 diffs)
• plant_canopy_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/module_interface.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove unused arrays in pmc_rrd_local
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Add pcm_rrd_local and pcm_wrd_local
 ! 
@@ -1789 +1792 @@
                                nynf, nync, nyn_on_file,                        &
                                nysf, nysc, nys_on_file,                        &
-                               tmp_2d, tmp_3d, found                           &
+                               found                                           &
                                            )
 

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! - Remove unused arrays in pmc_rrd_local
+! - Remove one exchange of ghost points
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! - Bugfix, read restart data for time-averaged pcm output quantities
 ! - Output of plant-canopy quantities will fill values
@@ -1152 +1156 @@
        pch_index_ji = 0
        
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
+       DO  i = nxlg, nxrg
+          DO  j = nysg, nyng
              DO  k = 0, pch_index
                 IF ( lad_s(k,j,i) /= 0 )  pch_index_ji(j,i) = k
@@ -1168 +1172 @@
           ENDDO
        ENDDO
-
-       CALL exchange_horiz_2d_int( pch_index_ji, nys, nyn, nxl, nxr, nbgp )
 !
 !--    Calculate global pch_index value (index of top of plant canopy from ground)
        pch_index = MAXVAL( pch_index_ji )
-       
-       
 !
 !--    Exchange pch_index from all processors
@@ -1475 +1475 @@
     SUBROUTINE pcm_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc,          &
                               nxr_on_file, nynf, nync, nyn_on_file, nysf,      &
-                              nysc, nys_on_file, tmp_2d, tmp_3d, found )
+                              nysc, nys_on_file, found )
 
        INTEGER(iwp) ::  k               !< 
@@ -1492 +1492 @@
 
        LOGICAL, INTENT(OUT)  :: found
-
-       REAL(wp), DIMENSION(nys_on_file-nbgp:nyn_on_file+nbgp,                  &
-                           nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_2d   !< temporary 2D array
-
-       REAL(wp), DIMENSION(nzb:nzt+1,                                          &
-                           nys_on_file-nbgp:nyn_on_file+nbgp,                  &
-                           nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d   !< temporary 3D array with entire vertical dimension
 
        REAL(wp), DIMENSION(0:pch_index,                                        &