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

Timestamp:

Apr 10, 2018 11:01:03 AM (6 years ago)

Author:

Giersch

Message:

Bugfix of the allocation of spectrum_x and spectrum_y in case of restart runs, spectrum_x and spectrum_y are now treat as global restart data not local

File:

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

palm/trunk/SOURCE/read_restart_data_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! spectrum_x and spectrum_y have been moved to global data
+! 
+! 2921 2018-03-22 15:05:23Z Giersch
 ! spinup_time, day_of_year_init and time_utc_init are also read now
 ! 
@@ -117 +120 @@
 
        USE spectra_mod,                                                        &
-           ONLY:  average_count_sp
+           ONLY:  average_count_sp, spectrum_x, spectrum_y
 
        USE statistics,                                                         &
@@ -569 +572 @@
              CASE ( 'simulated_time' )
                 READ ( 13 )  simulated_time
+             CASE ( 'spectrum_x' )
+                IF ( .NOT. ALLOCATED( spectrum_x ) )  THEN
+                   ALLOCATE( spectrum_x( 1:nx/2, 1:100, 1:10 ) )
+                ENDIF
+                READ ( 13 )  spectrum_x
+             CASE ( 'spectrum_y' )
+                IF ( .NOT. ALLOCATED( spectrum_y ) )  THEN
+                   ALLOCATE( spectrum_y( 1:ny/2, 1:100, 1:10 ) )
+                ENDIF
+                READ ( 13 )  spectrum_y
              CASE ( 'spinup_time' )
                 READ ( 13 )  spinup_time
@@ -1052 +1065 @@
         ONLY:  id_random_array, seq_random_array
 
-    USE spectra_mod,                                                           &
-        ONLY:  spectrum_x, spectrum_y
-
     USE surface_mod,                                                           &
         ONLY :  surface_rrd_local
@@ -1690 +1700 @@
                    shf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  =          &
                       tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
-
-                CASE ( 'spectrum_x' )
-                   IF ( k == 1 )  THEN
-                      IF ( nx_on_file /= nx )  THEN
-                         message_string = 'rrd_local: spectrum_x ' //          &
-                                     'on restart file ignored because' //      &
-                                     '&total numbers of grid points (nx) ' //  &
-                                     'do not match'
-                         CALL message( 'rrd_local', 'PA0293',                  &
-                                       0, 1, 0, 6, 0 )
-                         READ ( 13 )  rdummy
-                      ELSE
-                         READ ( 13 )  spectrum_x
-                      ENDIF
-                   ENDIF
-
-                CASE ( 'spectrum_y' )
-                   IF ( k == 1 )  THEN
-                      IF ( ny_on_file /= ny )  THEN
-                         message_string = 'rrd_local: spectrum_y ' //          &
-                                     'on restart file ignored because' //      &
-                                     '&total numbers of grid points (ny) '//   &
-                                     'do not match'
-                         CALL message( 'rrd_local', 'PA0294',                  &
-                                       0, 1, 0, 6, 0 )
-                         READ ( 13 )  rdummy
-                      ELSE
-                         READ ( 13 )  spectrum_y
-                      ENDIF
-                   ENDIF
                    
                 CASE ( 'ssws_av' )