Changeset 5876 for LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/coef_diff_turb_mod.f90

Timestamp:

Nov 21, 2025, 3:58:45 PM (5 weeks ago)

Author:

yann meurdesoif

Message:

GPU port of cdrag + call_atke + coef_diff_turb

YM

File:

• LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/coef_diff_turb_mod.f90 (modified) (9 diffs)

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

@@ -6 +6 @@
 ! at surface(cdrag)
 !
-  USE clesphys_mod_h
     IMPLICIT NONE
 
@@ -16 +15 @@
        ypaprs, ypplay, yu, yv, yq, yt, yts, yqsurf, ycdragm, &
        ycoefm, ycoefh ,yq2, yeps, ydrgpro)
-
+!$gpum horizontal knon
     USE dimphy, ONLY : klev
     USE mod_grid_phy_lmdz, ONLY: klon_glo
@@ -24 +23 @@
     USE yoethf_mod_h
     USE compbl_mod_h
+    USE yamada4_mod, ONLY : yamada4
+    USE ustarhb_mod, ONLY : ustarhb
+    USE vdif_kcay_mod, ONLY : vdif_kcay
+    USE coefkzmin_mod, ONLY : coefkzmin
+    USE clesphys_mod_h, ONLY : ksta, ksta_ter, ok_kzmin
+
 !
 ! Calculate coefficients(ycoefm, ycoefh) for turbulent diffusion in the
@@ -206 +211 @@
        qsurf, &
        pcfm, pcfh)
-    
+!$gpum horizontal knon    
     USE yomcst_mod_h
-    USE dimphy
+    USE dimphy, ONLY : klev
     USE indice_sol_mod
     USE print_control_mod, ONLY: prt_level, lunout
@@ -254 +259 @@
 ! Quelques constantes et options:
 !
-    REAL, PARAMETER :: cepdu2=0.1**2
+    REAL, PARAMETER :: cepdu2=0.1*0.1
     REAL, PARAMETER :: CKAP=0.4
     REAL, PARAMETER :: cb=5.0
@@ -295 +300 @@
     REAL gamt(2:klev)
 
-    LOGICAL, SAVE :: appel1er=.TRUE.
-    !$OMP THREADPRIVATE(appel1er)
 !
 ! Fonctions thermodynamiques et fonctions d'instabilite
@@ -306 +309 @@
     isommet=klev
       
-    IF (appel1er) THEN
-       IF (prt_level > 9) THEN
-          WRITE(lunout,*)'coefkz, opt_ec:', opt_ec
-          WRITE(lunout,*)'coefkz, richum:', richum
-          IF (richum) WRITE(lunout,*)'coefkz, ratqs:', ratqs
-          WRITE(lunout,*)'coefkz, isommet:', isommet
-          WRITE(lunout,*)'coefkz, tvirtu:', tvirtu
-          appel1er = .FALSE.
-       ENDIF
-    ENDIF
 !
 ! Initialiser les sorties
@@ -496 +489 @@
   SUBROUTINE coefkz2(nsrf, knon, paprs, pplay,t, &
        pcfm, pcfh)
-
+!$gpum horizontal knon
     USE yomcst_mod_h
-USE dimphy
+    USE dimphy, ONLY : klev
     USE indice_sol_mod
 
@@ -522 +515 @@
     REAL, DIMENSION(knon, klev+1), INTENT(IN) ::  paprs
     REAL, DIMENSION(knon, klev), INTENT(IN)   ::  pplay
-    REAL, DIMENSION(knon, klev), INTENT(IN)   :: t(klon,klev)
+    REAL, DIMENSION(knon, klev), INTENT(IN)   ::  t
 
     REAL, DIMENSION(knon, klev), INTENT(OUT)  :: pcfm, pcfh