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

Timestamp:

Dec 18, 2018 12:31:28 PM (6 years ago)

Author:

knoop

Message:

OpenACC port for SPEC

File:

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

palm/trunk/SOURCE/boundary_conds.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! OpenACC port for SPEC
+! 
+! 3589 2018-11-30 15:09:51Z suehring
 ! Move the control parameter "salsa" from salsa_mod to control_parameters 
 ! (M. Kurppa)
@@ -289 +292 @@
        kb = MERGE( -1, 1, l == 0 )
        !$OMP PARALLEL DO PRIVATE( i, j, k )
+       !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+       !$ACC PRESENT(bc_h, w_p)
        DO  m = 1, bc_h(l)%ns
           i = bc_h(l)%i(m)            
@@ -300 +305 @@
 !-- Top boundary. A nested domain ( ibc_uv_t = 3 ) does not require settings.
     IF ( ibc_uv_t == 0 )  THEN
+        !$ACC KERNELS PRESENT(u_p, v_p, u_init, v_init)
         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
         u_p(nzt+1,:,:) = u_p(nzt,:,:)
@@ -311 +318 @@
     IF (  .NOT.  child_domain  .AND.  .NOT.  nesting_offline  .AND.            &
                  TRIM(coupling_mode) /= 'vnested_fine' )  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
 
@@ -342 +351 @@
           kb = MERGE( -1, 1, l == 0 )
           !$OMP PARALLEL DO PRIVATE( i, j, k )
+          !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+          !$ACC PRESENT(bc_h, pt_p)
           DO  m = 1, bc_h(l)%ns
              i = bc_h(l)%i(m)            
@@ -364 +375 @@
         pt_p(nzt+1,:,:) = pt_p(nzt,:,:)
     ELSEIF ( ibc_pt_t == 2 )  THEN
+        !$ACC KERNELS PRESENT(pt_p, dzu)
         pt_p(nzt+1,:,:) = pt_p(nzt,:,:) + bc_pt_t_val * dzu(nzt+1)
+        !$ACC END KERNELS
     ENDIF
 
@@ -379 +392 @@
              kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
+             !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+             !$ACC PRESENT(bc_h, e_p)
              DO  m = 1, bc_h(l)%ns
                 i = bc_h(l)%i(m)            
@@ -444 +459 @@
 
        IF ( .NOT. child_domain )  THEN
+          !$ACC KERNELS PRESENT(e_p)
           e_p(nzt+1,:,:) = e_p(nzt,:,:)
+          !$ACC END KERNELS
 !
 !--    Nesting case: if parent operates in RANS mode and child in LES mode,