Ignore:
Timestamp:
Jan 17, 2017 4:38:49 PM (5 years ago)
Author:
raasch
Message:

all OpenACC directives and related parts removed from the code

File:
1 edited

Legend:

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

    r2101 r2118  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! OpenACC directives removed
    2323!
    2424! Former revisions:
     
    195195
    196196          REAL(wp)    ::  ll(nxl_z:nxr_z,nys_z:nyn_z) !<
    197           !$acc declare create( ll )
    198197
    199198
     
    201200          nnyh = ( ny + 1 ) / 2
    202201
    203           !$acc kernels present( tric )
    204202          DO  j = nys_z, nyn_z
    205203             DO  i = nxl_z, nxr_z
     
    239237             ENDDO
    240238          ENDDO
    241           !$acc end kernels
    242239
    243240          IF ( ibc_p_b == 1 )  THEN
    244              !$acc kernels present( tric )
    245241             DO  j = nys_z, nyn_z
    246242                DO  i = nxl_z, nxr_z
     
    248244                ENDDO
    249245             ENDDO
    250              !$acc end kernels
    251246          ENDIF
    252247          IF ( ibc_p_t == 1 )  THEN
    253              !$acc kernels present( tric )
    254248             DO  j = nys_z, nyn_z
    255249                DO  i = nxl_z, nxr_z
     
    257251                ENDDO
    258252             ENDDO
    259              !$acc end kernels
    260253          ENDIF
    261254
     
    288281
    289282          REAL(wp), DIMENSION(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1)   ::  ar1 !<
    290           !$acc declare create( ar1 )
    291283
    292284!
    293285!--       Forward substitution
    294286          DO  k = 0, nz - 1
    295              !$acc kernels present( ar, tri )
    296287             DO  j = nys_z, nyn_z
    297288                DO  i = nxl_z, nxr_z
     
    305296                ENDDO
    306297             ENDDO
    307              !$acc end kernels
    308298          ENDDO
    309299
     
    314304!--       the model domain.
    315305          DO  k = nz-1, 0, -1
    316              !$acc kernels present( ar, tri )
    317306             DO  j = nys_z, nyn_z
    318307                DO  i = nxl_z, nxr_z
     
    326315                ENDDO
    327316             ENDDO
    328              !$acc end kernels
    329317          ENDDO
    330318
     
    335323          IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1 )  THEN
    336324             IF ( nys_z == 0  .AND.  nxl_z == 0 )  THEN
    337                 !$acc kernels loop present( ar )
    338325                DO  k = 1, nz
    339326                   ar(nxl_z,nys_z,k) = 0.0_wp
    340327                ENDDO
    341                 !$acc end kernels loop
    342328             ENDIF
    343329          ENDIF
     
    372358
    373359          REAL(wp), DIMENSION(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1) ::  ar1 !<
    374           !$acc declare create( ar1 )
    375360
    376361!
    377362!--       Forward substitution
    378363          DO  k = 0, nz - 1
    379              !$acc kernels present( ar, tri )
    380              !$acc loop
    381364             DO  j = nys_z, nyn_z
    382365                DO  i = nxl_z, nxr_z
     
    390373                ENDDO
    391374             ENDDO
    392              !$acc end kernels
    393375          ENDDO
    394376
     
    399381!--       the model domain.
    400382          DO  k = nz-1, 0, -1
    401              !$acc kernels present( ar, tri )
    402              !$acc loop
    403383             DO  j = nys_z, nyn_z
    404384                DO  i = nxl_z, nxr_z
     
    412392                ENDDO
    413393             ENDDO
    414              !$acc end kernels
    415394          ENDDO
    416395
     
    421400          IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1 )  THEN
    422401             IF ( nys_z == 0  .AND.  nxl_z == 0 )  THEN
    423                 !$acc kernels loop present( ar )
    424402                DO  k = 1, nz
    425403                   ar(nxl_z,nys_z,k) = 0.0_wp
     
    451429!
    452430!--       Splitting
    453           !$acc kernels present( tri, tric )
    454           !$acc loop
    455431          DO  j = nys_z, nyn_z
    456              !$acc loop vector( 32 )
    457432             DO  i = nxl_z, nxr_z
    458433                tri(i,j,0,1) = tric(i,j,0)
    459434             ENDDO
    460435          ENDDO
    461           !$acc end kernels
    462436
    463437          DO  k = 1, nz-1
    464              !$acc kernels present( tri, tric )
    465              !$acc loop
    466              DO  j = nys_z, nyn_z
    467                 !$acc loop vector( 32 )
     438             DO  j = nys_z, nyn_z
    468439                DO  i = nxl_z, nxr_z
    469440                   tri(i,j,k,2) = ddzuw(k,1) / tri(i,j,k-1,1)
     
    471442                ENDDO
    472443             ENDDO
    473              !$acc end kernels
    474444          ENDDO
    475445
Note: See TracChangeset for help on using the changeset viewer.