Ignore:
Timestamp:
Nov 20, 2018 5:04:13 PM (2 years ago)
Author:
suehring
Message:

Revise ghost point exchange in netcdf-data input; new routine for ghost point exchange of 1-Byte Integer; Remove tabs in chemistry model which prevent compilation with gfortran and debug options

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_pegrid.f90

    r3347 r3542  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Introduce new MPI-datatype for ghostpoint exchange of 2D 8-bit Integer arrays
    2323!
    2424! Former revisions:
     
    11101110!
    11111111!-- Define new MPI derived datatypes for the exchange of ghost points in
    1112 !-- x- and y-direction for 2D-INTEGER arrays (line) - on normal grid
     1112!-- x- and y-direction for 2D-INTEGER arrays (line) - on normal grid.
     1113!-- Define types for 32-bit and 8-bit Integer. The 8-bit Integer are only
     1114!-- required on normal grid, while 32-bit Integer may be also required on
     1115!-- coarser grid level in case of multigrid solver.
     1116!
     1117!-- 8-bit Integer
     1118    CALL MPI_TYPE_VECTOR( nxr-nxl+1+2*nbgp, nbgp, ngp_y, MPI_BYTE,             &
     1119                          type_x_byte, ierr )
     1120    CALL MPI_TYPE_COMMIT( type_x_byte, ierr )
     1121
     1122    CALL MPI_TYPE_VECTOR( nbgp, ngp_y, ngp_y, MPI_BYTE,                        &
     1123                          type_y_byte, ierr )
     1124    CALL MPI_TYPE_COMMIT( type_y_byte, ierr )
     1125!
     1126!-- 32-bit Integer
    11131127    ALLOCATE( type_x_int(0:maximum_grid_level),                                &
    11141128              type_y_int(0:maximum_grid_level) )
    1115 
     1129             
    11161130    CALL MPI_TYPE_VECTOR( nxr-nxl+1+2*nbgp, nbgp, ngp_y, MPI_INTEGER,          &
    11171131                          type_x_int(0), ierr )
     
    12111225
    12121226!--       For 2D-exchange of INTEGER arrays on coarser grid level, where 2 ghost
    1213 !--       points need to be exchanged.
    1214           CALL MPI_TYPE_VECTOR( nxr_l-nxl_l+5, 2, nyn_l-nys_l+5, MPI_INTEGER,          &
     1227!--       points need to be exchanged. Only required for 32-bit Integer arrays.
     1228          CALL MPI_TYPE_VECTOR( nxr_l-nxl_l+5, 2, nyn_l-nys_l+5, MPI_INTEGER,  &
    12151229                                type_x_int(i), ierr )
    12161230          CALL MPI_TYPE_COMMIT( type_x_int(i), ierr )
    12171231
    12181232
    1219           CALL MPI_TYPE_VECTOR( 2, nyn_l-nys_l+5, nyn_l-nys_l+5, MPI_INTEGER,          &
     1233          CALL MPI_TYPE_VECTOR( 2, nyn_l-nys_l+5, nyn_l-nys_l+5, MPI_INTEGER,  &
    12201234                                type_y_int(i), ierr )
    12211235          CALL MPI_TYPE_COMMIT( type_y_int(i), ierr )
Note: See TracChangeset for help on using the changeset viewer.