Ignore:
Timestamp:
Jun 6, 2019 2:00:35 PM (5 years ago)
Author:
eckhard
Message:

inifor: bugfix: Fix wrong z_top height with gfortran -O2/-O3

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/UTIL/inifor/src/inifor_grid.f90

    r3997 r4019  
    2626! -----------------
    2727! $Id$
     28! bugfix: Replaced z_top array literal with array, fixes wrong z_top height with gfortran -O2/-O3
     29!
     30!
     31! 3997 2019-05-23 12:35:57Z eckhard
    2832! Read origin_z from static driver if given on the command line
    2933! Warn if origin_z is specified twice (by static driver and --elevation)
     
    206210    REAL(wp) ::  y0                = 0.0_wp       !< y coordinate of PALM-4U Earth tangent [m]
    207211    REAL(wp) ::  z0                = 0.0_wp       !< Elevation of the PALM-4U domain above sea level [m]
    208     REAL(wp) ::  z_top             = 0.0_wp       !< height of the scalar top boundary [m]
    209     REAL(wp) ::  zw_top            = 0.0_wp       !< height of the vertical velocity top boundary [m]
     212    REAL(wp) ::  z_top(1)          = 0.0_wp       !< height of the scalar top boundary [m]
     213    REAL(wp) ::  zw_top(1)         = 0.0_wp       !< height of the vertical velocity top boundary [m]
    210214    REAL(wp) ::  lonmin_cosmo      = 0.0_wp       !< Minimunm longitude of COSMO-DE's rotated-pole grid [COSMO rotated-pole rad]
    211215    REAL(wp) ::  lonmax_cosmo      = 0.0_wp       !< Maximum longitude of COSMO-DE's rotated-pole grid [COSMO rotated-pole rad]
     
    750754                     dz_stretch_level_end=dz_stretch_level_end,             &
    751755                     dz_stretch_factor_array=dz_stretch_factor_array)
    752     z(1:nz) = z_column(1:nz)
    753     z_top = z_column(nz+1)
     756    z(1:nz)  = z_column(1:nz)
     757    z_top(:) = z_column(nz+1)
    754758
    755759!
     
    761765    CALL midpoints(z_column, zw_column)
    762766    zw(1:nz-1) = zw_column(1:nz-1)
    763     zw_top     = zw_column(nz)
     767    zw_top(:)  = zw_column(nz)
    764768
    765769
     
    864868               ymin =  0.5_wp * dy, ymax = ly - 0.5_wp * dy,                   &
    865869               x0=x0, y0=y0, z0 = z0,                                          &
    866                nx = nx, ny = ny, nz = 1, z=(/z_top/))
     870               nx = nx, ny = ny, nz = 1, z=z_top)
    867871
    868872       CALL init_grid_definition('boundary', grid=u_east_grid,                 &
     
    894898               ymin = 0.5_wp * dy, ymax = ly - 0.5_wp * dy,                    &
    895899               x0=x0, y0=y0, z0 = z0,                                          &
    896                nx = nx-1, ny = ny, nz = 1, z=(/z_top/))
     900               nx = nx-1, ny = ny, nz = 1, z=z_top)
    897901
    898902       CALL init_grid_definition('boundary', grid=v_east_grid,                 &
     
    954958               ymin = 0.5_wp * dy, ymax = ly - 0.5_wp * dy,                    &
    955959               x0=x0, y0=y0, z0 = z0,                                          &
    956                nx = nx, ny = ny, nz = 1, z=(/zw_top/))
     960               nx = nx, ny = ny, nz = 1, z=zw_top)
    957961
    958962       CALL init_grid_definition('boundary intermediate', grid=scalars_east_intermediate,   &
Note: See TracChangeset for help on using the changeset viewer.