Changeset 5969 for LMDZ6/trunk/libf/phylmd/fonte_neige_mod.F90

Timestamp:

Dec 19, 2025, 11:55:37 AM (2 weeks ago)

Author:

yann meurdesoif

Message:

Bug Fix + gpu port of surf_landice

YM

File:

• LMDZ6/trunk/libf/phylmd/fonte_neige_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/fonte_neige_mod.F90

@@ -237 +237 @@
 #endif
      &   )
-!$gpum horizontal knon
+!$gpum horizontal knon klon
     USE indice_sol_mod
 #ifdef ISO
@@ -446 +446 @@
           j = knindex(i)
           !--temporal filtering
-          run_off_lic(i)   = coeff_rel*fqcalving(i) + (1.-coeff_rel)*run_off_lic_0(j)
-          run_off_lic_0(j) = run_off_lic(i)
+          run_off_lic(j)   = coeff_rel*fqcalving(i) + (1.-coeff_rel)*run_off_lic_0(j)
+          run_off_lic_0(j) = run_off_lic(j)
           !--add melting snow and liquid precip to runoff of ice cap
-          run_off_lic(i)   = run_off_lic(i) + fqfonte(i) + precip_rain(i)
+          run_off_lic(j)   = run_off_lic(j) + fqfonte(i) + precip_rain(i)
        ENDDO
     ENDIF
@@ -455 +455 @@
 #ifdef ISO
     DO i = 1, knon    
-      run_off_lic_diag(i) = run_off_lic(i)
+      run_off_lic_diag(i) = run_off_lic(knindex(i))
     ENDDO ! DO i = 1, knon    
 #endif
@@ -472 +472 @@
     IF (nisurf == is_lic) THEN
     DO i = 1, knon
-       runofflic_global(knindex(i)) = run_off_lic(i)
+       runofflic_global(knindex(i)) = run_off_lic(knindex(i))
     ENDDO
     ENDIF