Changeset 191 for palm/trunk/UTIL


Ignore:
Timestamp:
Aug 27, 2008 2:34:51 PM (16 years ago)
Author:
raasch
Message:

further adjustments and bugfixes for SGI system

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/UTIL/combine_plot_fields.f90

    r179 r191  
    44! Actual revisions:
    55! -----------------
     6! Subdomain data are read into temporary arrays pf_tmp/pf3d_tmp in order to
     7! avoid INTEL compiler warnings about (automatic) creation of temporary arrays
    68! Bugfix: three misplaced #endif directives
    79!
     
    109111    REAL ::  dx, simulated_time
    110112    REAL, DIMENSION(:),   ALLOCATABLE   ::  eta, ho, hu
    111     REAL, DIMENSION(:,:), ALLOCATABLE   ::  pf
    112     REAL(spk), DIMENSION(:,:,:), ALLOCATABLE ::  pf3d
     113    REAL, DIMENSION(:,:), ALLOCATABLE   ::  pf, pf_tmp
     114    REAL(spk), DIMENSION(:,:,:), ALLOCATABLE ::  pf3d, pf3d_tmp
    113115
    114116    PRINT*, ''
     
    357359!
    358360!--                Read the subdomain grid-point values
    359                    READ ( id+110 )  pf(xa:xe,ya:ye)
     361                   ALLOCATE( pf_tmp(xa:xe,ya:ye) )
     362                   READ ( id+110 )  pf_tmp
     363                   pf(xa:xe,ya:ye) = pf_tmp
     364                   DEALLOCATE( pf_tmp )
    360365                ENDIF
    361366                IF ( id == 0 )  fanz(av) = fanz(av) + 1
     
    654659!--          Read subdomain indices and grid point values
    655660             READ ( id+110, END=999 )  xa, xe, ya, ye, za, ze
    656              READ ( id+110 )  pf3d(xa:xe,ya:ye,za:ze)
     661             ALLOCATE( pf3d_tmp(xa:xe,ya:ye,za:ze) )
     662             READ ( id+110 )  pf3d_tmp
     663             pf3d(xa:xe,ya:ye,za:ze) = pf3d_tmp
     664             DEALLOCATE( pf3d_tmp )
    657665             IF ( id == 0 )  fanz(av) = fanz(av) + 1
    658666
Note: See TracChangeset for help on using the changeset viewer.