Changeset 1304 for palm/trunk/SOURCE/init_pegrid.f90
 Timestamp:
 Mar 12, 2014 10:29:42 AM (7 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/init_pegrid.f90
r1213 r1304 356 356 ! 1. transposition z > x 357 357 ! This transposition is not neccessary in case of a 1ddecomposition 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 358 366 IF ( pdims(2) /= 1 ) THEN 359 360 nys_x = nys361 nyn_x = nyn362 nny_x = nny363 367 IF ( MOD( nz , pdims(1) ) /= 0 ) THEN 364 368 WRITE( message_string, * ) 'transposition z > x:', & … … 367 371 CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 ) 368 372 ENDIF 369 nnz_x = nz / pdims(1)370 nzb_x = 1 + myidx * nnz_x371 nzt_x = ( myidx + 1 ) * nnz_x372 sendrecvcount_zx = nnx * nny * nnz_x373 374 ELSE375 !376 ! Setting of dummy values because otherwise variables are undefined in377 ! the next step x > y378 ! WARNING: This case has still to be clarified!!!!!!!!!!!!379 nnz_x = 1380 nzb_x = 1381 nzt_x = 1382 nny_x = nny383 384 373 ENDIF 385 374 … … 403 392 ! 3. transposition y > z (ELSE: x > y in case of 1Ddecomposition 404 393 ! 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 405 402 IF ( pdims(2) /= 1 ) THEN 406 403 ! … … 408 405 ! This transposition is not neccessary in case of a 1ddecomposition 409 406 ! along x, except that the uptreamspline method is switched on 410 nnx_z = nnx_y411 nxl_z = nxl_y412 nxr_z = nxr_y413 407 IF ( MOD( ny+1 , pdims(1) ) /= 0 ) THEN 414 408 WRITE( message_string, * ) 'transposition y > z:', & … … 417 411 CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 ) 418 412 ENDIF 419 nny_z = (ny+1) / pdims(1)420 nys_z = myidx * nny_z421 nyn_z = ( myidx + 1 ) * nny_z  1422 sendrecvcount_yz = nnx_y * nny_z * nnz_y423 413 424 414 ELSE
Note: See TracChangeset
for help on using the changeset viewer.