Changeset 5914 for LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/surf_land_bucket_mod.F90

Timestamp:

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

Author:

yann meurdesoif

Message:

gpu port of surf_land
YM

File:

• LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/surf_land_bucket_mod.F90 (modified) (5 diffs)

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)