Changeset 1320 for palm/trunk/SOURCE/fft_xy.f90
- Timestamp:
- Mar 20, 2014 8:40:49 AM (10 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/fft_xy.f90 (modified) (26 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/fft_xy.f90
r1310 r1320 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! ONLY-attribute added to USE-statements, 23 ! kind-parameters added to all INTEGER and REAL declaration statements, 24 ! kinds are defined in new module kinds, 25 ! old module precision_kind is removed, 26 ! revision history before 2012 removed, 27 ! comment fields (!:) to be used for variable explanations added to 28 ! all variable declaration statements 23 29 ! 24 30 ! Former revisions: … … 65 71 ! 1036 2012-10-22 13:43:42Z raasch 66 72 ! code put under GPL (PALM 3.9) 67 !68 ! 274 2009-03-26 15:11:21Z heinze69 ! Output of messages replaced by message handling routine.70 !71 ! Feb. 200772 ! RCS Log replace by Id keyword, revision history cleaned up73 !74 ! Revision 1.4 2006/03/28 12:27:09 raasch75 ! Stop when system-specific fft is selected on NEC. For unknown reasons this76 ! causes a program abort during first allocation in init_grid.77 !78 ! Revision 1.2 2004/04/30 11:44:27 raasch79 ! Module renamed from fft_for_1d_decomp to fft_xy, 1d-routines renamed to80 ! fft_x and fft_y,81 ! function FFT replaced by subroutine FFTN due to problems with 64-bit82 ! mode on ibm,83 ! shape of array cwork is explicitly stored in ishape/jshape and handled84 ! to routine FFTN instead of shape-function (due to compiler error on85 ! decalpha),86 ! non vectorized FFT for nec included87 73 ! 88 74 ! Revision 1.1 2002/06/11 13:00:49 raasch … … 96 82 !------------------------------------------------------------------------------! 97 83 98 USE control_parameters 99 USE indices 84 USE control_parameters, & 85 ONLY: fft_method, message_string 86 87 USE indices, & 88 ONLY: nx, ny, nz 89 100 90 #if defined( __cuda_fft ) 101 91 USE ISO_C_BINDING … … 103 93 USE, INTRINSIC :: ISO_C_BINDING 104 94 #endif 105 USE precision_kind 106 USE singleton 95 96 USE kinds 97 98 USE singleton, & 99 ONLY: fftn 100 107 101 USE temperton_fft 108 USE transpose_indices 102 103 USE transpose_indices, & 104 ONLY: nyn_x, nys_x, nzb_x, nzb_y, nzt_x, nzt_y 109 105 110 106 IMPLICIT NONE … … 113 109 PUBLIC fft_x, fft_x_1d, fft_y, fft_y_1d, fft_init, fft_x_m, fft_y_m 114 110 115 INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: ifax_x, ifax_y 116 117 LOGICAL, SAVE :: init_fft = .FALSE. 118 119 REAL, SAVE :: dnx, dny, sqr_dnx, sqr_dny 120 REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trigs_x, trigs_y 111 INTEGER(iwp), DIMENSION(:), ALLOCATABLE, SAVE :: ifax_x !: 112 INTEGER(iwp), DIMENSION(:), ALLOCATABLE, SAVE :: ifax_y !: 113 114 LOGICAL, SAVE :: init_fft = .FALSE. !: 115 116 REAL(wp), SAVE :: dnx !: 117 REAL(wp), SAVE :: dny !: 118 REAL(wp), SAVE :: sqr_dnx !: 119 REAL(wp), SAVE :: sqr_dny !: 120 121 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trigs_x !: 122 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trigs_y !: 121 123 122 124 #if defined( __ibm ) 123 INTEGER, PARAMETER :: nau1 = 20000, nau2 = 22000 125 INTEGER(iwp), PARAMETER :: nau1 = 20000 !: 126 INTEGER(iwp), PARAMETER :: nau2 = 22000 !: 124 127 ! 125 128 !-- The following working arrays contain tables and have to be "save" and 126 129 !-- shared in OpenMP sense 127 REAL, DIMENSION(nau1), SAVE :: aux1, auy1, aux3, auy3 130 REAL(wp), DIMENSION(nau1), SAVE :: aux1 !: 131 REAL(wp), DIMENSION(nau1), SAVE :: auy1 !: 132 REAL(wp), DIMENSION(nau1), SAVE :: aux3 !: 133 REAL(wp), DIMENSION(nau1), SAVE :: auy3 !: 134 128 135 #elif defined( __nec ) 129 INTEGER, SAVE :: nz1 130 REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trig_xb, trig_xf, trig_yb, & 131 trig_yf 136 INTEGER(iwp), SAVE :: nz1 !: 137 138 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trig_xb !: 139 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trig_xf !: 140 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trig_yb !: 141 REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE :: trig_yf !: 142 132 143 #elif defined( __cuda_fft ) 133 INTEGER(C_INT), SAVE :: plan_xf, plan_xi, plan_yf, plan_yi 134 INTEGER, SAVE :: total_points_x_transpo, total_points_y_transpo 144 INTEGER(C_INT), SAVE :: plan_xf !: 145 INTEGER(C_INT), SAVE :: plan_xi !: 146 INTEGER(C_INT), SAVE :: plan_yf !: 147 INTEGER(C_INT), SAVE :: plan_yi !: 148 149 INTEGER(iwp), SAVE :: total_points_x_transpo !: 150 INTEGER(iwp), SAVE :: total_points_y_transpo !: 135 151 #endif 136 152 137 153 #if defined( __fftw ) 138 154 INCLUDE 'fftw3.f03' 139 INTEGER(KIND=C_INT) :: nx_c, ny_c 140 COMPLEX(KIND=C_DOUBLE_COMPLEX), DIMENSION(:), ALLOCATABLE, SAVE :: x_out, y_out 141 REAL(KIND=C_DOUBLE), DIMENSION(:), ALLOCATABLE, SAVE :: x_in, y_in 155 INTEGER(KIND=C_INT) :: nx_c !: 156 INTEGER(KIND=C_INT) :: ny_c !: 157 158 COMPLEX(KIND=C_DOUBLE_COMPLEX), DIMENSION(:), ALLOCATABLE, SAVE :: & 159 x_out !: 160 COMPLEX(KIND=C_DOUBLE_COMPLEX), DIMENSION(:), ALLOCATABLE, SAVE :: & 161 y_out !: 162 163 REAL(KIND=C_DOUBLE), DIMENSION(:), ALLOCATABLE, SAVE :: & 164 x_in !: 165 REAL(KIND=C_DOUBLE), DIMENSION(:), ALLOCATABLE, SAVE :: & 166 y_in !: 167 168 142 169 TYPE(C_PTR), SAVE :: plan_xf, plan_xi, plan_yf, plan_yi 143 170 #endif … … 186 213 !-- in OpenMP sense 187 214 #if defined( __ibm ) 188 REAL, DIMENSION(0:nx+2) :: workx 189 REAL, DIMENSION(0:ny+2) :: worky 190 REAL, DIMENSION(nau2) :: aux2, auy2, aux4, auy4 215 REAL(wp), DIMENSION(0:nx+2) :: workx !: 216 REAL(wp), DIMENSION(0:ny+2) :: worky !: 217 REAL(wp), DIMENSION(nau2) :: aux2 !: 218 REAL(wp), DIMENSION(nau2) :: auy2 !: 219 REAL(wp), DIMENSION(nau2) :: aux4 !: 220 REAL(wp), DIMENSION(nau2) :: auy4 !: 191 221 #elif defined( __nec ) 192 REAL , DIMENSION(0:nx+3,nz+1) :: work_x193 REAL , DIMENSION(0:ny+3,nz+1) :: work_y194 REAL , DIMENSION(6*(nx+3),nz+1) :: workx195 REAL , DIMENSION(6*(ny+3),nz+1) :: worky222 REAL(wp), DIMENSION(0:nx+3,nz+1) :: work_x !: 223 REAL(wp), DIMENSION(0:ny+3,nz+1) :: work_y !: 224 REAL(wp), DIMENSION(6*(nx+3),nz+1) :: workx !: 225 REAL(wp), DIMENSION(6*(ny+3),nz+1) :: worky !: 196 226 #endif 197 227 … … 228 258 CALL message( 'fft_init', 'PA0187', 1, 2, 0, 6, 0 ) 229 259 230 ALLOCATE( trig_xb(2*(nx+1)), trig_xf(2*(nx+1)), &260 ALLOCATE( trig_xb(2*(nx+1)), trig_xf(2*(nx+1)), & 231 261 trig_yb(2*(ny+1)), trig_yf(2*(ny+1)) ) 232 262 … … 240 270 CALL DZFFT( 0, nx+1, sqr_dnx, work_x, work_x, trig_xf, workx, 0 ) 241 271 CALL ZDFFT( 0, nx+1, sqr_dnx, work_x, work_x, trig_xb, workx, 0 ) 242 CALL DZFFTM( 0, nx+1, nz1, sqr_dnx, work_x, nx+4, work_x, nx+4, &272 CALL DZFFTM( 0, nx+1, nz1, sqr_dnx, work_x, nx+4, work_x, nx+4, & 243 273 trig_xf, workx, 0 ) 244 CALL ZDFFTM( 0, nx+1, nz1, sqr_dnx, work_x, nx+4, work_x, nx+4, &274 CALL ZDFFTM( 0, nx+1, nz1, sqr_dnx, work_x, nx+4, work_x, nx+4, & 245 275 trig_xb, workx, 0 ) 246 276 ! … … 248 278 CALL DZFFT( 0, ny+1, sqr_dny, work_y, work_y, trig_yf, worky, 0 ) 249 279 CALL ZDFFT( 0, ny+1, sqr_dny, work_y, work_y, trig_yb, worky, 0 ) 250 CALL DZFFTM( 0, ny+1, nz1, sqr_dny, work_y, ny+4, work_y, ny+4, &280 CALL DZFFTM( 0, ny+1, nz1, sqr_dny, work_y, ny+4, work_y, ny+4, & 251 281 trig_yf, worky, 0 ) 252 CALL ZDFFTM( 0, ny+1, nz1, sqr_dny, work_y, ny+4, work_y, ny+4, &282 CALL ZDFFTM( 0, ny+1, nz1, sqr_dny, work_y, ny+4, work_y, ny+4, & 253 283 trig_yb, worky, 0 ) 254 284 #elif defined( __cuda_fft ) … … 278 308 nx_c = nx+1 279 309 ny_c = ny+1 280 ALLOCATE( x_in(0:nx+2), y_in(0:ny+2), x_out(0:(nx+1)/2), &310 ALLOCATE( x_in(0:nx+2), y_in(0:ny+2), x_out(0:(nx+1)/2), & 281 311 y_out(0:(ny+1)/2) ) 282 312 plan_xf = FFTW_PLAN_DFT_R2C_1D( nx_c, x_in, x_out, FFTW_ESTIMATE ) … … 322 352 IMPLICIT NONE 323 353 324 CHARACTER (LEN=*) :: direction 325 INTEGER :: i, ishape(1), j, k 326 327 LOGICAL :: forward_fft 328 329 REAL, DIMENSION(0:nx+2) :: work 330 REAL, DIMENSION(nx+2) :: work1 331 COMPLEX, DIMENSION(:), ALLOCATABLE :: cwork 354 CHARACTER (LEN=*) :: direction !: 355 356 COMPLEX(wp), DIMENSION(:), ALLOCATABLE :: cwork !: 357 358 INTEGER(iwp) :: i !: 359 INTEGER(iwp) :: ishape(1) !: 360 INTEGER(iwp) :: j !: 361 INTEGER(iwp) :: k !: 362 363 LOGICAL :: forward_fft !: 364 365 REAL(wp), DIMENSION(0:nx+2) :: work !: 366 REAL(wp), DIMENSION(nx+2) :: work1 !: 367 332 368 #if defined( __ibm ) 333 REAL, DIMENSION(nau2) :: aux2, aux4 369 REAL(wp), DIMENSION(nau2) :: aux2 !: 370 REAL(wp), DIMENSION(nau2) :: aux4 !: 334 371 #elif defined( __nec ) 335 REAL , DIMENSION(6*(nx+1)) :: work2372 REAL(wp), DIMENSION(6*(nx+1)) :: work2 !: 336 373 #elif defined( __cuda_fft ) 337 COMPLEX(dpk), DIMENSION(0:(nx+1)/2,nys_x:nyn_x,nzb_x:nzt_x) :: ar_tmp 374 COMPLEX(dpk), DIMENSION(0:(nx+1)/2,nys_x:nyn_x,nzb_x:nzt_x) :: & 375 ar_tmp !: 338 376 !$acc declare create( ar_tmp ) 339 377 #endif 340 REAL, DIMENSION(0:nx,nys_x:nyn_x), OPTIONAL :: ar_2d 341 REAL, DIMENSION(0:nx,nys_x:nyn_x,nzb_x:nzt_x) :: ar 378 379 REAL(wp), DIMENSION(0:nx,nys_x:nyn_x), OPTIONAL :: & 380 ar_2d !: 381 REAL(wp), DIMENSION(0:nx,nys_x:nyn_x,nzb_x:nzt_x) :: & 382 ar !: 342 383 343 384 IF ( direction == 'forward' ) THEN … … 540 581 DO j = nys_x, nyn_x 541 582 542 CALL DRCFT( 0, ar, 1, work, 1, nx+1, 1, 1, sqr_dnx, aux1, nau1,&543 aux2, nau2 )583 CALL DRCFT( 0, ar, 1, work, 1, nx+1, 1, 1, sqr_dnx, aux1, & 584 nau1, aux2, nau2 ) 544 585 545 586 DO i = 0, (nx+1)/2 … … 570 611 work(nx+2) = 0.0 571 612 572 CALL DCRFT( 0, work, 1, work, 1, nx+1, 1, -1, sqr_dnx, aux3, nau1, &573 aux 4, nau2 )613 CALL DCRFT( 0, work, 1, work, 1, nx+1, 1, -1, sqr_dnx, & 614 aux3, nau1, aux4, nau2 ) 574 615 575 616 DO i = 0, nx … … 709 750 IMPLICIT NONE 710 751 711 CHARACTER (LEN=*) :: direction 712 INTEGER :: i, ishape(1) 713 714 LOGICAL :: forward_fft 715 716 REAL, DIMENSION(0:nx) :: ar 717 REAL, DIMENSION(0:nx+2) :: work 718 REAL, DIMENSION(nx+2) :: work1 719 COMPLEX, DIMENSION(:), ALLOCATABLE :: cwork 752 CHARACTER (LEN=*) :: direction !: 753 754 INTEGER(iwp) :: i !: 755 INTEGER(iwp) :: ishape(1) !: 756 757 LOGICAL :: forward_fft !: 758 759 REAL(wp), DIMENSION(0:nx) :: ar !: 760 REAL(wp), DIMENSION(0:nx+2) :: work !: 761 REAL(wp), DIMENSION(nx+2) :: work1 !: 762 763 COMPLEX(wp), DIMENSION(:), ALLOCATABLE :: cwork !: 764 720 765 #if defined( __ibm ) 721 REAL, DIMENSION(nau2) :: aux2, aux4 766 REAL(wp), DIMENSION(nau2) :: aux2 !: 767 REAL(wp), DIMENSION(nau2) :: aux4 !: 722 768 #elif defined( __nec ) 723 REAL , DIMENSION(6*(nx+1)) :: work2769 REAL(wp), DIMENSION(6*(nx+1)) :: work2 !: 724 770 #endif 725 771 … … 838 884 IF ( forward_fft ) THEN 839 885 840 CALL DRCFT( 0, ar, 1, work, 1, nx+1, 1, 1, sqr_dnx, aux1, nau1, &886 CALL DRCFT( 0, ar, 1, work, 1, nx+1, 1, 1, sqr_dnx, aux1, nau1, & 841 887 aux2, nau2 ) 842 888 … … 945 991 IMPLICIT NONE 946 992 947 CHARACTER (LEN=*) :: direction 948 INTEGER :: i, j, jshape(1), k 949 INTEGER :: nxl_y_bound, nxl_y_l, nxr_y_bound, nxr_y_l 950 951 LOGICAL :: forward_fft 952 953 REAL, DIMENSION(0:ny+2) :: work 954 REAL, DIMENSION(ny+2) :: work1 955 COMPLEX, DIMENSION(:), ALLOCATABLE :: cwork 993 CHARACTER (LEN=*) :: direction !: 994 995 INTEGER(iwp) :: i !: 996 INTEGER(iwp) :: j !: 997 INTEGER(iwp) :: jshape(1) !: 998 INTEGER(iwp) :: k !: 999 INTEGER(iwp) :: nxl_y_bound !: 1000 INTEGER(iwp) :: nxl_y_l !: 1001 INTEGER(iwp) :: nxr_y_bound !: 1002 INTEGER(iwp) :: nxr_y_l !: 1003 1004 LOGICAL :: forward_fft !: 1005 1006 REAL(wp), DIMENSION(0:ny+2) :: work !: 1007 REAL(wp), DIMENSION(ny+2) :: work1 !: 1008 1009 COMPLEX(wp), DIMENSION(:), ALLOCATABLE :: cwork !: 1010 956 1011 #if defined( __ibm ) 957 REAL, DIMENSION(nau2) :: auy2, auy4 1012 REAL(wp), DIMENSION(nau2) :: auy2 !: 1013 REAL(wp), DIMENSION(nau2) :: auy4 !: 958 1014 #elif defined( __nec ) 959 REAL , DIMENSION(6*(ny+1)) :: work21015 REAL(wp), DIMENSION(6*(ny+1)) :: work2 !: 960 1016 #elif defined( __cuda_fft ) 961 COMPLEX(dpk), DIMENSION(0:(ny+1)/2,nxl_y:nxr_y,nzb_y:nzt_y) :: ar_tmp 1017 COMPLEX(dpk), DIMENSION(0:(ny+1)/2,nxl_y:nxr_y,nzb_y:nzt_y) :: & 1018 ar_tmp !: 962 1019 !$acc declare create( ar_tmp ) 963 1020 #endif 964 REAL, DIMENSION(0:ny,nxl_y_l:nxr_y_l,nzb_y:nzt_y) :: ar 965 REAL, DIMENSION(0:ny,nxl_y_bound:nxr_y_bound,nzb_y:nzt_y) :: ar_tr 1021 1022 REAL(wp), DIMENSION(0:ny,nxl_y_l:nxr_y_l,nzb_y:nzt_y) :: & 1023 ar !: 1024 REAL(wp), DIMENSION(0:ny,nxl_y_bound:nxr_y_bound,nzb_y:nzt_y) :: & 1025 ar_tr !: 966 1026 967 1027 IF ( direction == 'forward' ) THEN … … 1140 1200 DO i = nxl_y_l, nxr_y_l 1141 1201 1142 CALL DRCFT( 0, ar, 1, work, 1, ny+1, 1, 1, sqr_dny, auy1, nau1, &1143 auy2, nau2 )1202 CALL DRCFT( 0, ar, 1, work, 1, ny+1, 1, 1, sqr_dny, auy1, & 1203 nau1, auy2, nau2 ) 1144 1204 1145 1205 DO j = 0, (ny+1)/2 … … 1170 1230 work(ny+2) = 0.0 1171 1231 1172 CALL DCRFT( 0, work, 1, work, 1, ny+1, 1, -1, sqr_dny, auy3, nau1,&1173 auy 4, nau2 )1232 CALL DCRFT( 0, work, 1, work, 1, ny+1, 1, -1, sqr_dny, & 1233 auy3, nau1, auy4, nau2 ) 1174 1234 1175 1235 DO j = 0, ny … … 1307 1367 1308 1368 CHARACTER (LEN=*) :: direction 1309 INTEGER :: j, jshape(1) 1310 1311 LOGICAL :: forward_fft 1312 1313 REAL, DIMENSION(0:ny) :: ar 1314 REAL, DIMENSION(0:ny+2) :: work 1315 REAL, DIMENSION(ny+2) :: work1 1316 COMPLEX, DIMENSION(:), ALLOCATABLE :: cwork 1369 1370 INTEGER(iwp) :: j !: 1371 INTEGER(iwp) :: jshape(1) !: 1372 1373 LOGICAL :: forward_fft !: 1374 1375 REAL(wp), DIMENSION(0:ny) :: ar !: 1376 REAL(wp), DIMENSION(0:ny+2) :: work !: 1377 REAL(wp), DIMENSION(ny+2) :: work1 !: 1378 1379 COMPLEX(wp), DIMENSION(:), ALLOCATABLE :: cwork !: 1380 1317 1381 #if defined( __ibm ) 1318 REAL, DIMENSION(nau2) :: auy2, auy4 1382 REAL(wp), DIMENSION(nau2) :: auy2 !: 1383 REAL(wp), DIMENSION(nau2) :: auy4 !: 1319 1384 #elif defined( __nec ) 1320 REAL , DIMENSION(6*(ny+1)) :: work21385 REAL(wp), DIMENSION(6*(ny+1)) :: work2 !: 1321 1386 #endif 1322 1387 … … 1437 1502 IF ( forward_fft ) THEN 1438 1503 1439 CALL DRCFT( 0, ar, 1, work, 1, ny+1, 1, 1, sqr_dny, auy1, nau1, &1504 CALL DRCFT( 0, ar, 1, work, 1, ny+1, 1, 1, sqr_dny, auy1, nau1, & 1440 1505 auy2, nau2 ) 1441 1506 … … 1458 1523 work(ny+2) = 0.0 1459 1524 1460 CALL DCRFT( 0, work, 1, work, 1, ny+1, 1, -1, sqr_dny, auy3, nau1,&1461 auy4, nau2 )1525 CALL DCRFT( 0, work, 1, work, 1, ny+1, 1, -1, sqr_dny, auy3, & 1526 nau1, auy4, nau2 ) 1462 1527 1463 1528 DO j = 0, ny … … 1527 1592 IMPLICIT NONE 1528 1593 1529 CHARACTER (LEN=*) :: direction 1530 INTEGER :: i, k, siza 1531 1532 REAL, DIMENSION(0:nx,nz) :: ar 1533 REAL, DIMENSION(0:nx+3,nz+1) :: ai 1534 REAL, DIMENSION(6*(nx+4),nz+1) :: work1 1594 CHARACTER (LEN=*) :: direction !: 1595 1596 INTEGER(iwp) :: i !: 1597 INTEGER(iwp) :: k !: 1598 INTEGER(iwp) :: siza !: 1599 1600 REAL(wp), DIMENSION(0:nx,nz) :: ar !: 1601 REAL(wp), DIMENSION(0:nx+3,nz+1) :: ai !: 1602 REAL(wp), DIMENSION(6*(nx+4),nz+1) :: work1 !: 1603 1535 1604 #if defined( __nec ) 1536 INTEGER :: sizw 1537 COMPLEX, DIMENSION((nx+4)/2+1,nz+1) :: work 1605 INTEGER(iwp) :: sizw !: 1606 1607 COMPLEX(wp), DIMENSION((nx+4)/2+1,nz+1) :: work !: 1538 1608 #endif 1539 1609 … … 1588 1658 !-- Tables are initialized once more. This call should not be 1589 1659 !-- necessary, but otherwise program aborts in asymmetric case 1590 CALL DZFFTM( 0, nx+1, nz1, sqr_dnx, work, nx+4, work, nx+4, &1660 CALL DZFFTM( 0, nx+1, nz1, sqr_dnx, work, nx+4, work, nx+4, & 1591 1661 trig_xf, work1, 0 ) 1592 1662 … … 1596 1666 ENDIF 1597 1667 1598 CALL DZFFTM( 1, nx+1, nz1, sqr_dnx, ai, siza, work, sizw, &1668 CALL DZFFTM( 1, nx+1, nz1, sqr_dnx, ai, siza, work, sizw, & 1599 1669 trig_xf, work1, 0 ) 1600 1670 … … 1613 1683 !-- Tables are initialized once more. This call should not be 1614 1684 !-- necessary, but otherwise program aborts in asymmetric case 1615 CALL ZDFFTM( 0, nx+1, nz1, sqr_dnx, work, nx+4, work, nx+4, &1685 CALL ZDFFTM( 0, nx+1, nz1, sqr_dnx, work, nx+4, work, nx+4, & 1616 1686 trig_xb, work1, 0 ) 1617 1687 … … 1664 1734 IMPLICIT NONE 1665 1735 1666 CHARACTER (LEN=*) :: direction 1667 INTEGER :: j, k, ny1, siza 1668 1669 REAL, DIMENSION(0:ny1,nz) :: ar 1670 REAL, DIMENSION(0:ny+3,nz+1) :: ai 1671 REAL, DIMENSION(6*(ny+4),nz+1) :: work1 1736 CHARACTER (LEN=*) :: direction !: 1737 1738 INTEGER(iwp) :: j !: 1739 INTEGER(iwp) :: k !: 1740 INTEGER(iwp) :: ny1 !: 1741 INTEGER(iwp) :: siza !: 1742 1743 REAL(wp), DIMENSION(0:ny1,nz) :: ar !: 1744 REAL(wp), DIMENSION(0:ny+3,nz+1) :: ai !: 1745 REAL(wp), DIMENSION(6*(ny+4),nz+1) :: work1 !: 1746 1672 1747 #if defined( __nec ) 1673 INTEGER :: sizw 1674 COMPLEX, DIMENSION((ny+4)/2+1,nz+1) :: work 1748 INTEGER(iwp) :: sizw !: 1749 1750 COMPLEX(wp), DIMENSION((ny+4)/2+1,nz+1) :: work !: 1675 1751 #endif 1676 1752
Note: See TracChangeset
for help on using the changeset viewer.
