Changeset 4356 for palm/trunk/SOURCE/read_restart_data_mod.f90

Timestamp:

Dec 20, 2019 5:09:33 PM (4 years ago)

Author:

suehring

Message:

Bugfix in message calls for local checks; error messages in init_grid slightly revised; bugfix in time_integration (uninitialized emission time index); read_restart_data (change from J.Resler): use dynamic array allocation rather than automatic arrays to avoid problems with stack memory

File:

• palm/trunk/SOURCE/read_restart_data_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/read_restart_data_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Change automatic arrays to allocatable ones in rrd_local, in order to avoid
+! memory problems due to too small stack size for large jobs with intel 
+! compiler. (J.Resler)
+! 
+! 4331 2019-12-10 18:25:02Z suehring
 ! Enable restart data for 2-m potential temperature output
 ! 
@@ -1050 +1055 @@
     INTEGER(iwp), DIMENSION(numprocs_previous_run) ::  overlap_count   !<
 
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  nxlfa      !<
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  nxrfa      !<
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  nynfa      !<
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  nysfa      !<
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  offset_xa  !<
-    INTEGER(iwp), DIMENSION(numprocs_previous_run,1000) ::  offset_ya  !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nxlfa      !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nxrfa      !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nynfa      !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nysfa      !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  offset_xa  !<
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  offset_ya  !<
 
     INTEGER(isp), DIMENSION(:,:),   ALLOCATABLE ::  tmp_2d_id_random   !< temporary array for storing random generator data
@@ -1070 +1075 @@
 !-- Read data from previous model run.
     CALL cpu_log( log_point_s(14), 'rrd_local', 'start' )
+!
+!-- Allocate temporary buffer arrays. In previous versions, there were 
+!-- declared as automated arrays, causing memory problems when these
+!-- were allocate on stack.
+    ALLOCATE( nxlfa(numprocs_previous_run,1000) )
+    ALLOCATE( nxrfa(numprocs_previous_run,1000) )
+    ALLOCATE( nynfa(numprocs_previous_run,1000) )
+    ALLOCATE( nysfa(numprocs_previous_run,1000) )
+    ALLOCATE( offset_xa(numprocs_previous_run,1000) )
+    ALLOCATE( offset_ya(numprocs_previous_run,1000) )
 
 !
@@ -1793 +1808 @@
 
        ENDDO ! dataloop
-
 !
 !--    Close the restart file
@@ -1801 +1815 @@
 
     ENDDO  ! loop over restart files
-
+!
+!-- Deallocate temporary buffer arrays
+    DEALLOCATE( nxlfa )
+    DEALLOCATE( nxrfa )
+    DEALLOCATE( nynfa )
+    DEALLOCATE( nysfa )
+    DEALLOCATE( offset_xa )
+    DEALLOCATE( offset_ya )
 !
 !-- Restore the original filename for the restart file to be written