Changeset 4150 for LMDZ6

Timestamp:

May 16, 2022, 10:43:24 AM (3 years ago)

Author:

Ehouarn Millour

Message:

Fix previously introduced bug in the code reorganization for OpenACC.
EM

File:

• LMDZ6/branches/Portage_acc/libf/phylmd/physiq_mod.F90 (modified) (1 diff)

LMDZ6/branches/Portage_acc/libf/phylmd/physiq_mod.F90

@@ -2999 +2999 @@
              !$acc update self(d_t_adjwk, d_q_adjwk)
              !
-             IF (wake_s(i) .GT. 1.e-3) THEN
-                !$acc kernels default(none) async
+                 !$acc kernels default(none) async
                 DO k=1,klev
                    DO i=1,klon
+                      IF (wake_s(i) .GT. 1.e-3) THEN
                       t_w(i,k) = t_w(i,k) + d_t_adjwk(i,k)
                       q_w(i,k) = q_w(i,k) + d_q_adjwk(i,k)
                       d_deltat_ajs_cv(i,k) = d_t_adjwk(i,k)
                       d_deltaq_ajs_cv(i,k) = d_q_adjwk(i,k)
+                      ELSE
+                      d_deltat_ajs_cv(i,k) = 0.
+                      d_deltaq_ajs_cv(i,k) = 0.
+                      ENDIF
                    ENDDO
                 ENDDO
                 !$acc end kernels
-             ELSE
-                !$acc kernels default(none)
-                d_deltat_ajs_cv(i,k) = 0.
-                d_deltaq_ajs_cv(i,k) = 0.
-                !$acc end kernels
-             ENDIF
+
              !$acc update self(d_deltat_ajs_cv,d_deltaq_ajs_cv)               
              IF (iflag_adjwk == 2) THEN