Changeset 5968 for LMDZ6/trunk/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Dec 19, 2025, 11:46:41 AM (2 weeks ago)

Author:

yann meurdesoif

Message:

GPU port of surf_ocean
Merge of commit r5896
YM

File:

• LMDZ6/trunk/libf/phylmd/surf_ocean_mod.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -30 +30 @@
 #endif               
         &       )
+!$gpum horizontal knon klon
 
     use albedo, only: alboc, alboc_cd
@@ -39 +40 @@
     USE ocean_cpl_mod, ONLY    : ocean_cpl_noice
     USE indice_sol_mod, ONLY : nbsrf, is_oce
+    USE ocean_albedo_mod, ONLY : ocean_albedo 
 #ifdef ISO
     USE infotrac_phy, ONLY : ntraciso=>ntiso,niso
@@ -175 +177 @@
     INTEGER               :: i, k
     REAL                  :: tmp
-    REAL, PARAMETER       :: cepdu2=(0.1)**2
+    REAL, PARAMETER       :: cepdu2=(0.1)*(0.1)
     REAL, DIMENSION(knon) :: alb_eau, z0_lim
     REAL, DIMENSION(knon) :: radsol
@@ -188 +190 @@
     REAL, DIMENSION(knon) :: yrlat ! compressed latitude
     !**************************************************************************
+!$gpum nocall checksum
 
 #ifdef ISO
@@ -245 +248 @@
     SELECT CASE(type_ocean)
     CASE('couple')
+!$gpum nocall
        CALL ocean_cpl_noice( &
             swnet, lwnet, alb1, &
@@ -259 +263 @@
 
     CASE('slab')
+!$gpum nocall
        CALL ocean_slab_noice( &
             itime, dtime, jour, knon, knindex, &
@@ -403 +408 @@
 !ym not porting to GPU for now      
        if (type_ocean /= 'couple') sss(:knon) = 35.
+!$gpum nocall
        call bulk_flux(tkt, tks, taur, dter, dser, t_int, s_int, ds_ns, dt_ns, &
             u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = sss(:knon), &