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

LMDZ6/trunk/libf/phylmd/nuage.f90

@@ -1 +1 @@
 ! $Id$
+MODULE nuage_mod
+  PRIVATE
+
+  PUBLIC nuage, diagcld1, diagcld2
+
+  CONTAINS
 
 SUBROUTINE nuage(paprs, pplay, t, pqlwp,picefra, pclc, pcltau, pclemi, pch, pcl, pcm, &
@@ -95 +101 @@
   REAL cldtaupi(klon, klev) ! pre-industrial cloud opt thickness for diag
   REAl dzfice(klon)
+  REAL :: pp_ratio(klon)
   ! jq-end
 
@@ -115 +122 @@
             CALL icefrac_lscp(klon,t(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),zfice(:),dzfice(:))
         ELSE
-            CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:))
+            pp_ratio(1:klon) = pplay(1:klon,k)/paprs(1:klon,1)
+            CALL icefrac_lsc(klon,t(:,k),pp_ratio(:),zfice(:))
 
         ENDIF
@@ -257 +265 @@
   RETURN
 END SUBROUTINE nuage
+
 SUBROUTINE diagcld1(paprs, pplay, rain, snow, kbot, ktop, diafra, dialiq)
   USE dimphy
@@ -340 +349 @@
   RETURN
 END SUBROUTINE diagcld1
+
 SUBROUTINE diagcld2(paprs, pplay, t, q, diafra, dialiq)
   USE dimphy
@@ -431 +441 @@
   RETURN
 END SUBROUTINE diagcld2
+
+END MODULE nuage_mod