Changeset 5826 for LMDZ6

Timestamp:

Sep 22, 2025, 3:22:58 PM (2 months ago)

Author:

rkazeroni

Message:

For GPU porting of ratqs_main routine:

• Simplify refactoring of rev 5820
• Add forgotten THREADPRIVATE

File:

• LMDZ6/trunk/libf/phylmd/lmdz_ratqs_main.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_ratqs_main.f90

@@ -7 +7 @@
 
   REAL, SAVE :: resolmax
-
-  REAL, SAVE, ALLOCATABLE :: resol(:)
-  !$OMP THREADPRIVATE(resol)
+  !$OMP THREADPRIVATE(resolmax)
 
   PUBLIC ratqs_main_first, ratqs_main
@@ -20 +18 @@
   INTEGER, INTENT(in) :: klon
   REAL, DIMENSION(klon), INTENT(in) :: cell_area
-  INTEGER :: i
 
   IF (first) THEN
-     ALLOCATE(resol(klon))
      resolmax=sqrt(maxval(cell_area))
-     DO i=1,klon
-        resol(i) = sqrt(cell_area(i))
-     END DO
      first = .FALSE.
   END IF
@@ -93 +86 @@
 real facteur,zfratqs1,zfratqs2
 real, dimension(klon,klev) :: ratqs_hetero_,ratqs_oro_,ratqs_tke_
-real, dimension(klon) :: fact
+real :: resol, fact
 
 !-------------------------------------------------------------------------
@@ -196 +189 @@
          do k=1,klev
             do i=1,klon
-              fact(i) = sqrt(resol(i)/resolmax)
-              ratqss(i,k)=ratqsbas*fact(i)+0.5*(ratqshaut-ratqsbas)*fact(i) &
+              resol=sqrt(cell_area(i))
+              fact = sqrt(resol/resolmax)
+              ratqss(i,k)=ratqsbas*fact+0.5*(ratqshaut-ratqsbas)*fact &
               *( tanh( (ratqsp0-pplay(i,k))/ratqsdp) + 1.)
            enddo