Changeset 5889 for LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/ocean_forced_mod.F90

Timestamp:

Nov 27, 2025, 6:44:57 PM (5 weeks ago)

Author:

yann meurdesoif

Message:

GPU port of surf_seaice

YM

File:

• LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/ocean_forced_mod.F90 (modified) (7 diffs)

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

@@ -43 +43 @@
     USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
     use config_ocean_skin_m, only: activate_ocean_skin
+    USE calbeta_mod, ONLY : calbeta
+
 #ifdef ISO
     USE infotrac_phy, ONLY: ntiso,niso
@@ -310 +312 @@
 #endif            
        )
+!$gpum horizontal knon klon
 !
 ! This subroutine treats the ocean where there is ice. 
@@ -329 +332 @@
     USE fonte_neige_mod,  ONLY : fonte_neige
     USE indice_sol_mod
+    USE albsno_mod, ONLY : albsno
+    USE soil_mod, ONLY : soil
+    USE calbeta_mod, ONLY : calbeta
     USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
 #ifdef ISO
@@ -412 +418 @@
 ! Local variables
 !****************************************************************************************
-    LOGICAL                     :: check=.FALSE.
+    LOGICAL,PARAMETER           :: check=.FALSE.
     INTEGER                     :: i, j
     REAL                        :: zfra
@@ -470 +476 @@
 
 ! albedo  and radiation parameters
-    INTEGER, SAVE :: iflag_sic_albedo
+    INTEGER :: iflag_sic_albedo
 ! albedo old or NEMO
     REAL :: alb_sno_dry!=rn_alb_sdry !dry snow albedo
@@ -492 +498 @@
     ! ice (not snow). Should be visible only, not NIR
     REAL :: pen_ext !=si_pen_ext !extinction length of penetrating shortwave (m-1)
+    REAl :: lon(knon), lat(knon)   ! for indexation
 
 ! HF from ocean below ice
@@ -535 +542 @@
     IF (soil_model) THEN 
 ! update tsoil and calculate soilcap and soilflux
+       lon(1:knon) = longitude(knindex(1:knon))
+       lat(1:knon) = latitude(knindex(1:knon))
        CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, qsol, &
-        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil,soilcap, soilflux)
+        & lon, lat, tsoil,soilcap, soilflux)
        cal(1:knon) = RCPD / soilcap(1:knon)
        radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)