Changeset 2118 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Jan 17, 2017 4:38:49 PM (8 years ago)

Author:

raasch

Message:

all OpenACC directives and related parts removed from the code

File:

• palm/trunk/SOURCE/boundary_conds.f90 (modified) (8 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! OpenACC directives removed
 ! 
 ! Former revisions:
@@ -188 +188 @@
 !-- Bottom boundary 
     IF ( ibc_uv_b == 1 )  THEN
-       !$acc kernels present( u_p, v_p )
        u_p(nzb,:,:) = u_p(nzb+1,:,:)
        v_p(nzb,:,:) = v_p(nzb+1,:,:)
-       !$acc end kernels
-    ENDIF
-
-    !$acc kernels present( nzb_w_inner, w_p )
+    ENDIF
+
     DO  i = nxlg, nxrg
        DO  j = nysg, nyng
@@ -200 +197 @@
        ENDDO
     ENDDO
-    !$acc end kernels
 
 !
 !-- Top boundary. A nested domain ( ibc_uv_t = 3 ) does not require settings.
     IF ( ibc_uv_t == 0 )  THEN
-       !$acc kernels present( u_init, u_p, v_init, v_p )
         u_p(nzt+1,:,:) = u_init(nzt+1)
         v_p(nzt+1,:,:) = v_init(nzt+1)
-       !$acc end kernels
     ELSEIF ( ibc_uv_t == 1 )  THEN
-       !$acc kernels present( u_p, v_p )
         u_p(nzt+1,:,:) = u_p(nzt,:,:)
         v_p(nzt+1,:,:) = v_p(nzt,:,:)
-       !$acc end kernels
     ENDIF
 
     IF ( .NOT. nest_domain )  THEN
-       !$acc kernels present( w_p )
        w_p(nzt:nzt+1,:,:) = 0.0_wp  ! nzt is not a prognostic level (but cf. pres)
-       !$acc end kernels
     ENDIF
 
@@ -227 +217 @@
 !-- the sea surface temperature of the coupled ocean model.
     IF ( ibc_pt_b == 0 )  THEN
-       !$acc kernels present( nzb_s_inner, pt, pt_p )
-       !$acc loop independent
        DO  i = nxlg, nxrg
-          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt(nzb_s_inner(j,i),j,i)
           ENDDO
        ENDDO
-       !$acc end kernels
     ELSEIF ( ibc_pt_b == 1 )  THEN
-       !$acc kernels present( nzb_s_inner, pt_p )
-       !$acc loop independent
        DO  i = nxlg, nxrg
-          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt_p(nzb_s_inner(j,i)+1,j,i)
           ENDDO
        ENDDO
-      !$acc end kernels
     ENDIF
 
@@ -251 +233 @@
 !-- Temperature at top boundary
     IF ( ibc_pt_t == 0 )  THEN
-       !$acc kernels present( pt, pt_p )
         pt_p(nzt+1,:,:) = pt(nzt+1,:,:)
 !
@@ -259 +240 @@
            pt_p(nzt+1,:,:) = pt_init(nzt+1)
         ENDIF
-       !$acc end kernels
     ELSEIF ( ibc_pt_t == 1 )  THEN
-       !$acc kernels present( pt_p )
         pt_p(nzt+1,:,:) = pt_p(nzt,:,:)
-       !$acc end kernels
     ELSEIF ( ibc_pt_t == 2 )  THEN
-       !$acc kernels present( dzu, pt_p )
         pt_p(nzt+1,:,:) = pt_p(nzt,:,:) + bc_pt_t_val * dzu(nzt+1)
-       !$acc end kernels
     ENDIF
 
@@ -274 +250 @@
 !-- Generally Neumann conditions with de/dz=0 are assumed
     IF ( .NOT. constant_diffusion )  THEN
-       !$acc kernels present( e_p, nzb_s_inner )
-       !$acc loop independent
        DO  i = nxlg, nxrg
-          !$acc loop independent
           DO  j = nysg, nyng
              e_p(nzb_s_inner(j,i),j,i) = e_p(nzb_s_inner(j,i)+1,j,i)
@@ -285 +258 @@
           e_p(nzt+1,:,:) = e_p(nzt,:,:)
        ENDIF
-       !$acc end kernels
     ENDIF