Changeset 5489 for LMDZ6/branches/contrails/libf/phylmd/fonte_neige_mod.F90

Timestamp:

Jan 17, 2025, 6:16:25 PM (13 days ago)

Author:

aborella

Message:

Merge with trunk

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/fonte_neige_mod.F90 (modified) (5 diffs)

LMDZ6/branches/contrails/libf/phylmd/fonte_neige_mod.F90

@@ -231 +231 @@
   SUBROUTINE fonte_neige( knon, nisurf, knindex, dtime, &
        tsurf, precip_rain, precip_snow, &
-       snow, qsol, tsurf_new, evap &
+       snow, qsol, tsurf_new, evap, ice_sub &
 #ifdef ISO    
      & ,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag   &
@@ -288 +288 @@
     REAL, DIMENSION(klon), INTENT(INOUT) :: evap
 
+
+    REAL, DIMENSION(klon), INTENT(OUT)   :: ice_sub
 #ifdef ISO    
         ! sortie de quelques diagnostiques
@@ -297 +299 @@
     REAL, DIMENSION(klon), INTENT(OUT) ::  runoff_diag   
     REAL, DIMENSION(klon), INTENT(OUT) :: run_off_lic_diag  
-    REAL,                  INTENT(OUT) :: coeff_rel_diag
-#endif
+    REAL,                  INTENT(OUT) :: coeff_rel_diag    
+#endif
+
 
 ! Local variables
@@ -345 +348 @@
 
     snow_evap = 0.
+    ice_sub(:) = 0.
   
     IF (.NOT. ok_lic_cond) THEN
@@ -363 +367 @@
     
     bil_eau_s(:) = (precip_rain(:) * dtime) - (evap(:) - snow_evap(:)) * dtime
+
+    IF (nisurf==is_lic) THEN
+       DO i=1,knon
+          ice_sub(i)=evap(i)-snow_evap(i)
+       ENDDO
+    ENDIF
+
 #ifdef ISO
     snow_evap_diag(:) = snow_evap(:)