Changeset 5914 for LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd

Timestamp:

Dec 4, 2025, 8:02:43 PM (3 weeks ago)

Author:

yann meurdesoif

Message:

gpu port of surf_land
YM

Location:

LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd

Files:

• calbeta_clim.f90 (modified) (2 diffs)
• surf_land_bucket_hetero_mod.F90 (modified) (1 diff)
• surf_land_bucket_mod.F90 (modified) (5 diffs)
• surf_land_mod.F90 (modified) (6 diffs)
• surface_data.f90 (modified) (1 diff)

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

@@ -3 +3 @@
 ! fairhead Exp $
 !
+MODULE calbeta_clim_mod
+
+
+CONTAINS
 
 SUBROUTINE calbeta_clim(klon,time,lat_radian,beta)
+!$gpum horizontal klon
 
   !======================================================================
@@ -71 +76 @@
 
 end SUBROUTINE calbeta_clim
+
+END MODULE calbeta_clim_mod

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

@@ -43 +43 @@
     USE albsno_mod, ONLY :  albsno
     USE calbeta_mod, ONLY :  calbeta
+    USE calbeta_clim_mod, ONLY :  calbeta_clim
 
 #ifdef ISO

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

@@ -23 +23 @@
 #endif           
             )
-
+!$gpum horizontal knon klon
     USE limit_read_mod
     USE surface_data
@@ -50 +50 @@
     USE soil_mod, ONLY :  soil
     USE calbeta_mod, ONLY :  calbeta
+    USE calbeta_clim_mod, ONLY :  calbeta_clim
+
 
 !****************************************************************************************
@@ -125 +127 @@
     REAL :: coeff_rel_diag
 #endif 
+    REAL :: lat(knon), lon(knon)
+
 !
 !****************************************************************************************
@@ -169 +173 @@
     CALL calbeta(dtime, is_ter, knon, snow, qsol, beta, capsol, dif_grnd)
     if (type_veget=='betaclim') then
-       CALL calbeta_clim(knon,jour,latitude(knindex(1:knon)),beta)
+       lat(1:knon) = latitude(knindex(1:knon))
+       CALL calbeta_clim(knon,jour,lat ,beta)
     endif
        
 ! calculate temperature, heat capacity and conduction flux in soil
     IF (soil_model) THEN
+       lat(1:knon) = latitude(knindex(1:knon))
+       lon(1:knon) = longitude(knindex(1:knon))
        CALL soil(dtime, is_ter, knon, snow, tsurf, qsol,  & 
-        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
+        & lon, lat, tsoil, soilcap, soilflux)
 
        DO i=1, knon
@@ -337 +344 @@
        dummy_riverflow(:)   = 0.0
        dummy_coastalflow(:) = 0.0
+!$gpum nocall
        CALL cpl_send_land_fields(itime, knon, knindex, &
             dummy_riverflow, dummy_coastalflow)

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

@@ -29 +29 @@
 #endif               
                )
+!$gpum horizontal knon klon
 
     USE dimphy
     USE surface_data, ONLY    : ok_veget
-    USE carbon_cycle_mod
-
 
     ! See comments in each module surf_land_orchidee_xxx for compatiblity with ORCHIDEE
@@ -57 +56 @@
 #else
     ! Default version
-    USE surf_land_orchidee_mod
+    USE surf_land_orchidee_mod, ONLY : surf_land_orchidee
 #endif
 #endif
@@ -64 +63 @@
 #endif
     
-    USE surf_land_bucket_mod
-    USE surf_land_bucket_hetero_mod
-    USE calcul_fluxs_mod
-    USE indice_sol_mod
+    USE surf_land_bucket_mod, ONLY : surf_land_bucket
+    USE surf_land_bucket_hetero_mod, ONLY : surf_land_bucket_hetero
+    USE calcul_fluxs_mod, ONLY : calcul_flux_wind
+    USE indice_sol_mod, ONLY: nbsrf
 #ifdef ISO
     use infotrac_phy, ONLY: ntiso,niso
@@ -76 +75 @@
 #endif
 
-USE dimpft_mod_h
-        USE clesphys_mod_h
-    USE yomcst_mod_h
-USE print_control_mod, ONLY: lunout
+    USE dimpft_mod_h, ONLY: nvm_lmdz
+    USE clesphys_mod_h, ONLY : nsw, ok_bs
+    USE yomcst_mod_h, ONLY: RCPD, RKAPPA
+    USE print_control_mod, ONLY: lunout
     USE dimsoil_mod_h, ONLY: nsoilmx
-    USE compbl_mod_h
+    USE compbl_mod_h, ONLY: iflag_hetero_surf
 
 ! Input variables  
@@ -250 +249 @@
 
        ! temporary for keeping same results using lwdown_m instead of lwdown
+!$gpum nocall
        CALL surf_land_orchidee(itime, dtime, date0, knon, &
             knindex, rlon, rlat, yrmu0, pctsrf, &
@@ -306 +306 @@
         IF (klon .EQ. 1) THEN
           !
+!$gpum nocall
           CALL surf_land_bucket_hetero(itime, jour, knon, knindex, debut, dtime,&
               tsurf, p1lay, cdragh, precip_rain, precip_totsnow, temp_air, &

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

@@ -12 +12 @@
   !$OMP THREADPRIVATE(ok_veget)
 
-  CHARACTER(len=10), SAVE :: type_veget   ! orchidee/y/bucket/n/betaclim
+  CHARACTER(len=s2s_str_maxlen), SAVE :: type_veget   ! orchidee/y/bucket/n/betaclim
   !$OMP THREADPRIVATE(type_veget)