Changeset 5833 for LMDZ6/trunk/libf/phylmd/alpale_th.f90

Timestamp:

Sep 24, 2025, 3:12:42 PM (2 months ago)

Author:

rkazeroni

Message:

For GPU porting of alpale_th and alpale_wk routines:

• Put routine into module (speeds up source-to-source transformation)
• Add "horizontal" comment to specify possible names of horizontal variables
• Move declaration of variables with SAVE attributes from the compute routine to the module
• Move one-time instructions (allocate, getin, print) to a dedicated routine_first

File:

• LMDZ6/trunk/libf/phylmd/alpale_th.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/alpale_th.f90

@@ -2 +2 @@
 ! $Id$
 !
+!$gpum horizontal klon
+MODULE alpale_th_mod
+  PRIVATE
+
+  LOGICAL, SAVE                                              :: first = .TRUE.
+  !$OMP THREADPRIVATE(first)
+  LOGICAL, SAVE                                              :: multiply_proba_notrig = .FALSE.
+  !$OMP THREADPRIVATE(multiply_proba_notrig)
+  REAL, SAVE                                                 :: random_notrig_max=1.
+  !$OMP THREADPRIVATE(random_notrig_max)
+  REAL, SAVE                                                 :: cv_feed_area
+  !$OMP THREADPRIVATE(cv_feed_area)
+
+  PUBLIC alpale_th, alpale_th_first
+
+  CONTAINS
+
+SUBROUTINE alpale_th_first()
+
+  USE alpale_mod, ONLY: iflag_clos_bl
+  USE ioipsl_getin_p_mod, ONLY : getin_p
+
+  IMPLICIT NONE
+
+  IF (first) THEN
+    CALL getin_p('multiply_proba_notrig',multiply_proba_notrig)
+    IF (iflag_clos_bl .LT. 3) THEN
+      random_notrig_max=1.
+      CALL getin_p('random_notrig_max',random_notrig_max)
+    ELSEIF (iflag_clos_bl .EQ. 3) THEN  ! (iflag_clos_bl .LT. 3)
+      cv_feed_area = 1.e10   ! m2
+      CALL getin_p('cv_feed_area', cv_feed_area)
+    ENDIF  !! (iflag_clos_bl .LT. 3)
+    first = .FALSE.
+  ENDIF
+
+END SUBROUTINE alpale_th_first
+
 SUBROUTINE alpale_th ( dtime, lmax_th, t_seri, cell_area,  &
                        cin, s2, n2, strig,  &
@@ -18 +56 @@
 
   USE dimphy
-  USE ioipsl_getin_p_mod, ONLY : getin_p
   USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level
-  USE alpale_mod
+  USE alpale_mod, ONLY: iflag_clos_bl, iflag_coupl, iflag_trig_bl, s_trig, tau_trig_deep, tau_trig_shallow
   IMPLICIT NONE
 
@@ -50 +87 @@
 !----------------
   INTEGER                                                    :: i
-  LOGICAL, SAVE                                              :: first = .TRUE.
-  LOGICAL, SAVE                                              :: multiply_proba_notrig = .FALSE.
-  REAL, SAVE                                                 :: random_notrig_max=1.
-  REAL, SAVE                                                 :: cv_feed_area
   REAL                                                       :: birth_number
   REAL, DIMENSION(klon)                                      :: ale_bl_ref
   REAL, DIMENSION(klon)                                      :: tau_trig
-!
-    !$OMP THREADPRIVATE(multiply_proba_notrig) 
-    !$OMP THREADPRIVATE(random_notrig_max) 
-    !$OMP THREADPRIVATE(cv_feed_area) 
-    !$OMP THREADPRIVATE(first)
-!
+
  REAL umexp  ! expression of (1.-exp(-x))/x valid for all x, especially when x->0
  REAL x
 !
-     CHARACTER (LEN=20) :: modname='alpale_th'
+     CHARACTER (LEN=20), PARAMETER :: modname='alpale_th'
      CHARACTER (LEN=80) :: abort_message
      
@@ -84 +112 @@
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 !
-
-    IF (first) THEN
-      CALL getin_p('multiply_proba_notrig',multiply_proba_notrig)
-      IF (iflag_clos_bl .LT. 3) THEN
-         random_notrig_max=1.
-         CALL getin_p('random_notrig_max',random_notrig_max)
-      ELSEIF (iflag_clos_bl .EQ. 3) THEN  ! (iflag_clos_bl .LT. 3)
-         cv_feed_area = 1.e10   ! m2
-         CALL getin_p('cv_feed_area', cv_feed_area)
-      ENDIF  !! (iflag_clos_bl .LT. 3)
-      first=.FALSE.
-    ENDIF
 
 !!
@@ -344 +360 @@
    END SUBROUTINE alpale_th
 
+END MODULE alpale_th_mod