Changeset 1815 for palm/trunk/SOURCE/tridia_solver.f90
- Timestamp:
- Apr 6, 2016 1:49:59 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/tridia_solver.f90
r1809 r1815 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! cpp-switch intel11 removed 22 22 ! 23 23 ! Former revisions: … … 521 521 522 522 IF ( j <= nnyh ) THEN 523 #if defined( __intel11 )524 CALL maketri_1dd( j, tri_for_1d )525 #else526 523 CALL maketri_1dd( j ) 527 #endif528 524 ELSE 529 #if defined( __intel11 )530 CALL maketri_1dd( ny+1-j, tri_for_1d )531 #else532 525 CALL maketri_1dd( ny+1-j ) 533 #endif534 526 ENDIF 535 #if defined( __intel11 ) 536 CALL split_1dd( tri_for_1d ) 537 #else 527 538 528 CALL split_1dd 539 #endif540 529 CALL substi_1dd( ar, tri_for_1d ) 541 530 … … 548 537 !> computes the i- and j-dependent component of the matrix 549 538 !------------------------------------------------------------------------------! 550 #if defined( __intel11 )551 SUBROUTINE maketri_1dd( j, tri_for_1d )552 #else553 539 SUBROUTINE maketri_1dd( j ) 554 #endif555 540 556 541 USE constants, & … … 570 555 571 556 REAL(wp), DIMENSION(0:nx) :: l !< 572 573 #if defined( __intel11 )574 REAL(wp), DIMENSION(5,0:nx,0:nz-1) :: tri_for_1d !<575 #endif576 557 577 558 … … 623 604 !> Splitting of the tridiagonal matrix (Thomas algorithm) 624 605 !------------------------------------------------------------------------------! 625 #if defined( __intel11 )626 SUBROUTINE split_1dd( tri_for_1d )627 #else628 606 SUBROUTINE split_1dd 629 #endif630 631 607 632 608 IMPLICIT NONE … … 634 610 INTEGER(iwp) :: i !< 635 611 INTEGER(iwp) :: k !< 636 637 #if defined( __intel11 )638 REAL(wp), DIMENSION(5,0:nx,0:nz-1) :: tri_for_1d !<639 #endif640 612 641 613
Note: See TracChangeset
for help on using the changeset viewer.