Changeset 5828 for LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop_ini.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_ini.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

@@ -33 +33 @@
   REAL, PARAMETER :: k_ice0=0.005 ! units=m2/g
   REAL, PARAMETER :: df=1.66 ! diffusivity factor
+  LOGICAL, SAVE, ALLOCATABLE :: first(:,:)  ! Test, if the cloud optical depth exceeds the necessary threshold
+!$OMP THREADPRIVATE(first)
 !$OMP THREADPRIVATE(prt_level, lunout, flag_aerosol, iflag_t_glace)
 !$OMP THREADPRIVATE(iflag_rei, novlp, iflag_ice_thermo)  
@@ -45 +47 @@
 CONTAINS
 
-  SUBROUTINE cloud_optics_prop_ini(klon, prt_level_in, lunout_in, flag_aerosol_in, &
+  SUBROUTINE cloud_optics_prop_ini(klon, klev, prt_level_in, lunout_in, flag_aerosol_in, &
        & ok_cdnc_in, bl95_b0_in, &
        & bl95_b1_in, latitude_deg_in, rpi_in, rg_in, rd_in, zepsec_in, novlp_in, &
@@ -53 +55 @@
 
     IMPLICIT NONE
-    INTEGER, INTENT(IN) :: klon
+    INTEGER, INTENT(IN) :: klon, klev
     INTEGER, INTENT(IN) :: prt_level_in, lunout_in
     INTEGER, INTENT(IN) :: flag_aerosol_in
@@ -64 +66 @@
 
     ALLOCATE(latitude_deg(klon))
+    ALLOCATE(first(klon, klev))
+    first(:,:) = .FALSE.
 
     prt_level = prt_level_in