Changeset 5883 for LMDZ6/branches/PBLSURF_GPUPORT

Timestamp:

Nov 24, 2025, 2:51:09 PM (5 weeks ago)

Author:

yann meurdesoif

Message:

GPU port of freinage subroutine
untested (need orchidee)
YM

Location:

LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd

Files:

• freinage.f90 (modified) (3 diffs)
• pbl_surface_mod.F90 (modified) (1 diff)

LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/freinage.f90

@@ +1 @@
+MODULE freinage_mod
 !
 ! $Id$
 !
+
+CONTAINS
+
   SUBROUTINE freinage(klon, knon, uu, vv,  &
        tt,veget,lai, height,ypaprs,ypplay,drag_pro,d_u,d_v)
+!$gpum horizontal knon klon
 
     !ONLINE:
@@ -120 +125 @@
       enddo
      enddo 
+     
       do k=1,klev
-        where (sumveg(1:knon,k) > 0.05 ) 
-!        drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k)/sumveg(1:knon,k)
-        drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k)
-        elsewhere
-        drag_pro(1:knon,k)=0.0
-       endwhere
+
+!ym where are not correctly supported by gpumorphosis switch to LOOP, IF/ENDIF
+!ym         where (sumveg(1:knon,k) > 0.05 ) 
+!ym !        drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k)/sumveg(1:knon,k)
+!ym         drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k)
+!ym         elsewhere
+!ym         drag_pro(1:knon,k)=0.0
+!ym        endwhere
+
+        DO i=1,knon
+          IF (sumveg(i,k) > 0.05) THEN
+            drag_pro(i,k)=Cd_frein*drag_pro(i,k)
+          ELSE
+           drag_pro(i,k)=0.0
+          ENDIF
+        ENDDO
         d_u(1:knon,k) =(-1)*drag_pro(1:knon,k)*uu(1:knon,k)*wind(1:knon,k)
         d_v(1:knon,k) =(-1)*drag_pro(1:knon,k)*vv(1:knon,k)*wind(1:knon,k)
@@ -134 +150 @@
  END SUBROUTINE freinage
 
+END MODULE freinage_mod
+

LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/pbl_surface_mod.F90

@@ -2450 +2450 @@
     USE mod_phys_lmdz_para, ONLY : is_master
     USE cdrag_mod, ONLY : cdrag
+    USE freinage_mod, ONLY : freinage
 
 IMPLICIT NONE