Changeset 5828 for LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop.f90

Timestamp:

Sep 23, 2025, 4:32:02 PM (2 months ago)

Author:

rkazeroni

Message:

For GPU porting of call_cloud_optics_prop routine:

• Add "horizontal" comment to specify possible names of horizontal variables
• Put routine into module (speeds up source-to-source transformation)
• Move declaration of variable with SAVE attribute outside of the compute routine to the module
• Record event with a 2D "first" array instead of a scalar to enable GPU porting
• Perform reduction on this "first" array and print (once) outside of the compute routine since this cannot be done on GPU in the current form

File:

• LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop.f90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop.f90

@@ -1 +1 @@
 ! $Id$
 MODULE lmdz_cloud_optics_prop
+  PRIVATE
+
+  LOGICAL, SAVE :: first_first = .TRUE.
+  !$OMP THREADPRIVATE(first_first)
+
+  PUBLIC cloud_optics_prop, cloud_optics_prop_post
 
 CONTAINS
+
+SUBROUTINE cloud_optics_prop_post()
+  USE lmdz_cloud_optics_prop_ini, ONLY: novlp
+  USE lmdz_cloud_optics_prop_ini, ONLY: first
+  IMPLICIT NONE
+
+  IF (first_first) THEN
+    IF (ANY(first)) THEN
+      IF (novlp==1) THEN
+        WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM_ &
+                &                                             &
+                &                                          RANDOM'
+        first_first = .FALSE.
+      ELSEIF (novlp==2) THEN
+        WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM'
+        first_first = .FALSE.
+      ELSEIF (novlp==3) THEN
+        WRITE (*, *) 'Hypothese de recouvrement: RANDOM'
+        first_first = .FALSE.
+      ENDIF
+    ENDIF
+  ENDIF
+
+END SUBROUTINE cloud_optics_prop_post
 
 SUBROUTINE cloud_optics_prop(klon, klev, paprs, pplay, temp, radocond, picefra, pclc, &
@@ -29 +59 @@
   USE lmdz_cloud_optics_prop_ini , ONLY : rei_coef, rei_min_temp
   USE lmdz_cloud_optics_prop_ini , ONLY : zepsec, novlp, iflag_ice_thermo, ok_new_lscp
+  USE lmdz_cloud_optics_prop_ini , ONLY : first
   
 
@@ -115 +146 @@
   ! Local variables
   !----------------
-
-  LOGICAL, SAVE :: first = .TRUE.
-  !$OMP THREADPRIVATE(first)
   INTEGER flag_max
 
@@ -652 +680 @@
 
           IF (novlp.EQ.2) THEN
-            IF (first) THEN
-              WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM'
-              first = .FALSE.
+            IF (first_first) THEN
+              first(i,k) = .TRUE.
             ENDIF
             flag_max = -1.
@@ -661 +688 @@
 
           IF (novlp.EQ.3) THEN
-            IF (first) THEN
-              WRITE (*, *) 'Hypothese de recouvrement: RANDOM'
-              first = .FALSE.
+            IF (first_first) THEN
+              first(i,k) = .TRUE.
             ENDIF
             flag_max = 1.
@@ -670 +696 @@
 
           IF (novlp.EQ.1) THEN
-            IF (first) THEN
-              WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM_ &
-                &                                             &
-                &                                          RANDOM'
-              first = .FALSE.
+            IF (first_first) THEN
+              first(i,k) = .TRUE.
             ENDIF
             flag_max = 1.
@@ -767 +790 @@
   ENDIF !ok_cdnc
 
-  first=.false. !to be sure
-
   RETURN