Ignore:
Timestamp:
Mar 12, 2014 10:29:42 AM (7 years ago)
Author:
raasch
Message:

openmp bugfix + bugfix for single core MPI runs
ulimit option in mrun changed from -Ss to -s

File:
1 edited

Legend:

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

    r1213 r1304  
    356356!-- 1. transposition  z --> x
    357357!-- This transposition is not neccessary in case of a 1d-decomposition along x
     358    nys_x = nys
     359    nyn_x = nyn
     360    nny_x = nny
     361    nnz_x = nz / pdims(1)
     362    nzb_x = 1 + myidx * nnz_x
     363    nzt_x = ( myidx + 1 ) * nnz_x
     364    sendrecvcount_zx = nnx * nny * nnz_x
     365
    358366    IF ( pdims(2) /= 1 )  THEN
    359 
    360        nys_x = nys
    361        nyn_x = nyn
    362        nny_x = nny
    363367       IF ( MOD( nz , pdims(1) ) /= 0 )  THEN
    364368          WRITE( message_string, * ) 'transposition z --> x:',                &
     
    367371          CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 )
    368372       ENDIF
    369        nnz_x = nz / pdims(1)
    370        nzb_x = 1 + myidx * nnz_x
    371        nzt_x = ( myidx + 1 ) * nnz_x
    372        sendrecvcount_zx = nnx * nny * nnz_x
    373 
    374     ELSE
    375 !
    376 !---   Setting of dummy values because otherwise variables are undefined in
    377 !---   the next step  x --> y
    378 !---   WARNING: This case has still to be clarified!!!!!!!!!!!!
    379        nnz_x = 1
    380        nzb_x = 1
    381        nzt_x = 1
    382        nny_x = nny
    383 
    384373    ENDIF
    385374
     
    403392!-- 3. transposition  y --> z  (ELSE:  x --> y  in case of 1D-decomposition
    404393!-- along x)
     394    nnx_z = nnx_y
     395    nxl_z = nxl_y
     396    nxr_z = nxr_y
     397    nny_z = (ny+1) / pdims(1)
     398    nys_z = myidx * nny_z
     399    nyn_z = ( myidx + 1 ) * nny_z - 1
     400    sendrecvcount_yz = nnx_y * nny_z * nnz_y
     401
    405402    IF ( pdims(2) /= 1 )  THEN
    406403!
     
    408405!--    This transposition is not neccessary in case of a 1d-decomposition
    409406!--    along x, except that the uptream-spline method is switched on
    410        nnx_z = nnx_y
    411        nxl_z = nxl_y
    412        nxr_z = nxr_y
    413407       IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
    414408          WRITE( message_string, * ) 'transposition y --> z:',            &
     
    417411          CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 )
    418412       ENDIF
    419        nny_z = (ny+1) / pdims(1)
    420        nys_z = myidx * nny_z
    421        nyn_z = ( myidx + 1 ) * nny_z - 1
    422        sendrecvcount_yz = nnx_y * nny_z * nnz_y
    423413
    424414    ELSE
Note: See TracChangeset for help on using the changeset viewer.