Changeset 4285 for LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90

Timestamp:

Sep 27, 2022, 10:39:25 AM (2 years ago)

Author:

jghattas

Message:

Reported modifications done in rev [4283] into phydlmdiso which now compiles again.

File:

• LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90

@@ -144 +144 @@
     INTEGER                  :: i,j,nt
     REAL, DIMENSION(klon)    :: fqfonte,ffonte
+    REAL, DIMENSION(klon)    :: run_off_lic_frac
 #ifdef ISO       
       real, parameter :: t_coup = 273.15
@@ -484 +485 @@
 ! Send run-off on land-ice to coupler if coupled ocean.
 ! run_off_lic has been calculated in fonte_neige or surf_inlandsis
-!
-!****************************************************************************************
-    IF (type_ocean=='couple') THEN
-       CALL cpl_send_landice_fields(itime, knon, knindex, run_off_lic)
+! If landice_opt>=2, corresponding call is done from surf_land_orchidee
+!****************************************************************************************
+    IF (type_ocean=='couple' .AND. landice_opt .LT. 2) THEN
+       ! Compress fraction where run_off_lic is active (here all pctsrf(is_lic))
+       run_off_lic_frac(:)=0.0
+       DO j = 1, knon
+          i = knindex(j)
+          run_off_lic_frac(j) = pctsrf(i,is_lic)
+       ENDDO
+       
+       CALL cpl_send_landice_fields(itime, knon, knindex, run_off_lic, run_off_lic_frac)
     ENDIF