Changeset 1128 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Apr 12, 2013 6:19:32 AM (11 years ago)

Author:

raasch

Message:

asynchronous transfer of ghost point data for acc-optimized version

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (16 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! those parts requiring global communication moved to time_integration,
+! loop index bounds in accelerator version replaced by i_left, i_right, j_south,
+! j_north
 !
 ! Former revisions:
@@ -223 +225 @@
     CALL cpu_log( log_point(32), 'all progn.equations', 'start' )
 
-
-!
-!-- Calculate those variables needed in the tendency terms which need
-!-- global communication
-    IF ( .NOT. neutral )  CALL calc_mean_profile( pt, 4 )
-    IF ( ocean         )  CALL calc_mean_profile( rho, 64 )
-    IF ( humidity      )  CALL calc_mean_profile( vpt, 44 )
-    IF ( .NOT. constant_diffusion )  CALL production_e_init
-    IF ( ( ws_scheme_mom .OR. ws_scheme_sca )  .AND.  &
-         intermediate_timestep_count == 1 )  CALL ws_statistics
 !
 !-- Loop over all prognostic equations
@@ -724 +716 @@
              ENDIF
 
-         ENDIF
+          ENDIF
+
 !
 !--       If required, compute prognostic equation for turbulent kinetic 
@@ -811 +804 @@
     REAL    ::  sbt
 
-!
-!-- Calculate those variables needed in the tendency terms which need
-!-- global communication
-    IF ( .NOT. neutral )  CALL calc_mean_profile( pt, 4 )
-    IF ( ocean         )  CALL calc_mean_profile( rho, 64 )
-    IF ( humidity      )  CALL calc_mean_profile( vpt, 44 )
-    IF ( ( ws_scheme_mom .OR. ws_scheme_sca )  .AND.  &
-         intermediate_timestep_count == 1 )  CALL ws_statistics
 
 !
@@ -1337 +1322 @@
 
        CALL cpu_log( log_point(16), 'tke-equation', 'start' )
-
-!
-!--    TKE-tendency terms with communication
-       CALL production_e_init
 
        sbt = tsc(2)
@@ -1464 +1445 @@
 
 !
-!-- Calculate those variables needed in the tendency terms which need
-!-- global communication
-    IF ( .NOT. neutral )  CALL calc_mean_profile( pt, 4 )
-    IF ( ocean         )  CALL calc_mean_profile( rho, 64 )
-    IF ( humidity      )  CALL calc_mean_profile( vpt, 44 )
-    IF ( ( ws_scheme_mom .OR. ws_scheme_sca )  .AND.  &
-         intermediate_timestep_count == 1 )  CALL ws_statistics
-
-!
 !-- u-velocity component
 !++ Statistics still not ported to accelerators
@@ -1501 +1473 @@
 !-- External pressure gradient
     IF ( dp_external )  THEN
-       DO  i = nxlu, nxr
-          DO  j = nys, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              DO  k = dp_level_ind_b+1, nzt
                 tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
@@ -1516 +1488 @@
     !$acc kernels present( nzb_u_inner, rdf, tend, tu_m, u, ug, u_p )
     !$acc loop
-    DO  i = nxlu, nxr
-       DO  j = nys, nyn
+    DO  i = i_left, i_right
+       DO  j = j_south, j_north
           !$acc loop vector( 32 )
           DO  k = 1, nzt
@@ -1563 +1535 @@
 !-- External pressure gradient
     IF ( dp_external )  THEN
-       DO  i = nxl, nxr
-          DO  j = nysv, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              DO  k = dp_level_ind_b+1, nzt
                 tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
@@ -1578 +1550 @@
     !$acc kernels present( nzb_v_inner, rdf, tend, tv_m, v, vg, v_p )
     !$acc loop
-    DO  i = nxl, nxr
-       DO  j = nysv, nyn
+    DO  i = i_left, i_right
+       DO  j = j_south, j_north
           !$acc loop vector( 32 )
           DO  k = 1, nzt
@@ -1640 +1612 @@
     !$acc kernels present( nzb_w_inner, rdf, tend, tw_m, w, w_p )
     !$acc loop
-    DO  i = nxl, nxr
-       DO  j = nys, nyn
+    DO  i = i_left, i_right
+       DO  j = j_south, j_north
           !$acc loop vector( 32 )
           DO  k = 1, nzt-1
@@ -1734 +1706 @@
        !$acc         present( tend, tpt_m, pt, pt_p )
        !$acc loop
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              !$acc loop vector( 32 )
              DO  k = 1, nzt
@@ -1802 +1774 @@
 !
 !--    Prognostic equation for salinity
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              DO  k = nzb_s_inner(j,i)+1, nzt
                 sa_p(k,j,i) = sa(k,j,i) + dt_3d * ( sbt * tend(k,j,i) +        &
@@ -1890 +1862 @@
 !
 !--    Prognostic equation for total water content / scalar
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              DO  k = nzb_s_inner(j,i)+1, nzt
                 q_p(k,j,i) = q(k,j,i) + dt_3d * ( sbt * tend(k,j,i) +          &
@@ -1919 +1891 @@
 
        CALL cpu_log( log_point(16), 'tke-equation', 'start' )
-
-!
-!--    TKE-tendency terms with communication
-       CALL production_e_init
 
        sbt = tsc(2)
@@ -1984 +1952 @@
        !$acc kernels present( e, e_p, nzb_s_inner, tend, te_m )
        !$acc loop
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
+       DO  i = i_left, i_right
+          DO  j = j_south, j_north
              !$acc loop vector( 32 )
              DO  k = 1, nzt