Changeset 5831 for LMDZ6

Timestamp:

Sep 24, 2025, 12:07:20 PM (2 months ago)

Author:

rkazeroni

Message:

Finish GPU porting work on ratqs_main:

• Computation of resolmax takes parallelization into account using reduce_max (variable now defined as resolmax_glo)

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_ratqs_main.f90

@@ -3 +3 @@
   PRIVATE
 
-  LOGICAL, SAVE :: first=.TRUE.  ! fisrt call to ratqs_main
+  LOGICAL, SAVE :: first=.TRUE.  ! first call to ratqs_main
   !$OMP THREADPRIVATE(first)
 
-  REAL, SAVE :: resolmax
-  !$OMP THREADPRIVATE(resolmax)
+  REAL, SAVE :: resolmax_glo
+  !$OMP THREADPRIVATE(resolmax_glo)
 
   PUBLIC ratqs_main_first, ratqs_main
@@ -14 +14 @@
 
 SUBROUTINE ratqs_main_first(klon, cell_area)
-
+  USE mod_phys_lmdz_para
   IMPLICIT NONE
   INTEGER, INTENT(in) :: klon
   REAL, DIMENSION(klon), INTENT(in) :: cell_area
+  REAL :: resolmax
 
   IF (first) THEN
      resolmax=sqrt(maxval(cell_area))
+     CALL reduce_max(resolmax, resolmax_glo)
+     CALL bcast(resolmax_glo)
      first = .FALSE.
   END IF
@@ -190 +193 @@
             do i=1,klon
               resol=sqrt(cell_area(i))
-              fact = sqrt(resol/resolmax)
+              fact = sqrt(resol/resolmax_glo)
               ratqss(i,k)=ratqsbas*fact+0.5*(ratqshaut-ratqsbas)*fact &
               *( tanh( (ratqsp0-pplay(i,k))/ratqsdp) + 1.)