Changeset 4123 for palm/trunk/UTIL

Timestamp:

Jul 26, 2019 1:45:03 PM (5 years ago)

Author:

gronemeier

Message:

bugfix: do not assue that output arrays start with index 0

File:

• palm/trunk/UTIL/binary_to_netcdf.f90 (modified) (6 diffs)

palm/trunk/UTIL/binary_to_netcdf.f90

@@ -781 +781 @@
    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  data_count_per_dimension  !< data count of variable per dimension
    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  bounds_start              !< lower bounds of variable
-   INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  bounds_end                !< upper bounds of variable
    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  bounds_origin             !< lower bounds of dimensions in output file
 
@@ -866 +865 @@
 
             ALLOCATE( bounds_start(1:n_dim) )
-            ALLOCATE( bounds_end(1:n_dim) )
             ALLOCATE( bounds_origin(1:n_dim) )
             ALLOCATE( start_positions(1:n_dim) )
@@ -872 +870 @@
 
             READ( bin_file_unit ) ( bounds_start(i), i = 1, n_dim )
-            READ( bin_file_unit ) ( bounds_end(i), i = 1, n_dim )
+            READ( bin_file_unit ) ( data_count_per_dimension(i), i = 1, n_dim )
             READ( bin_file_unit ) ( bounds_origin(i), i = 1, n_dim )
 
@@ -878 +876 @@
             CALL internal_message( 'debug', routine_name // &
                                    ': bounds_start = ' // TRIM( temp_string ) )
-            WRITE( temp_string, * ) bounds_end
+            WRITE( temp_string, * ) data_count_per_dimension
             CALL internal_message( 'debug', routine_name // &
-                                   ': bounds_end = ' // TRIM( temp_string ) )
+                                   ': data_count_per_dimension = ' // TRIM( temp_string ) )
             WRITE( temp_string, * ) bounds_origin
             CALL internal_message( 'debug', routine_name // &
@@ -888 +886 @@
 
             DO  i = 1, n_dim
-               data_count = data_count * ( bounds_end(i) - bounds_start(i) + 1 )
+               data_count = data_count * data_count_per_dimension(i)
                start_positions(i) = bounds_start(i) - bounds_origin(i) + 1
-               data_count_per_dimension(i) = bounds_end(i) - bounds_start(i) + 1
             ENDDO
 
@@ -981 +978 @@
             DEALLOCATE( data_count_per_dimension )
             DEALLOCATE( bounds_start )
-            DEALLOCATE( bounds_end )
             DEALLOCATE( bounds_origin )