Changeset 191 for palm/trunk/UTIL

Timestamp:

Aug 27, 2008 2:34:51 PM (16 years ago)

Author:

raasch

Message:

further adjustments and bugfixes for SGI system

File:

• palm/trunk/UTIL/combine_plot_fields.f90 (modified) (4 diffs)

palm/trunk/UTIL/combine_plot_fields.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Subdomain data are read into temporary arrays pf_tmp/pf3d_tmp in order to
+! avoid INTEL compiler warnings about (automatic) creation of temporary arrays
 ! Bugfix: three misplaced #endif directives
 !
@@ -109 +111 @@
     REAL ::  dx, simulated_time
     REAL, DIMENSION(:),   ALLOCATABLE   ::  eta, ho, hu
-    REAL, DIMENSION(:,:), ALLOCATABLE   ::  pf
-    REAL(spk), DIMENSION(:,:,:), ALLOCATABLE ::  pf3d
+    REAL, DIMENSION(:,:), ALLOCATABLE   ::  pf, pf_tmp
+    REAL(spk), DIMENSION(:,:,:), ALLOCATABLE ::  pf3d, pf3d_tmp
 
     PRINT*, ''
@@ -357 +359 @@
 !
 !--                Read the subdomain grid-point values
-                   READ ( id+110 )  pf(xa:xe,ya:ye)
+                   ALLOCATE( pf_tmp(xa:xe,ya:ye) )
+                   READ ( id+110 )  pf_tmp
+                   pf(xa:xe,ya:ye) = pf_tmp
+                   DEALLOCATE( pf_tmp )
                 ENDIF
                 IF ( id == 0 )  fanz(av) = fanz(av) + 1
@@ -654 +659 @@
 !--          Read subdomain indices and grid point values
              READ ( id+110, END=999 )  xa, xe, ya, ye, za, ze
-             READ ( id+110 )  pf3d(xa:xe,ya:ye,za:ze)
+             ALLOCATE( pf3d_tmp(xa:xe,ya:ye,za:ze) )
+             READ ( id+110 )  pf3d_tmp
+             pf3d(xa:xe,ya:ye,za:ze) = pf3d_tmp
+             DEALLOCATE( pf3d_tmp )
              IF ( id == 0 )  fanz(av) = fanz(av) + 1