Changeset 6177

Timestamp:

Apr 19, 2026, 2:16:40 PM (2 weeks ago)

Author:

evignon

Message:

petits ajustements dans les nouvelles params de banquise

Location:

LMDZ6/trunk/libf/phylmd

Files:

• limit_read_mod.f90 (modified) (1 diff)
• ocean_forced_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/limit_read_mod.f90

@@ -462 +462 @@
              IF (ierr /= NF90_NOERR) CALL abort_physic(modname,'Lecture echouee pour <SST>',1)
            !GG 
-             IF (iflag_seaice == 1) THEN
+             IF (iflag_seaice == 2) THEN
                ierr = NF90_INQ_VARID(nid, 'HICE', nvarid)
                IF (ierr /= NF90_NOERR) CALL abort_physic(modname,'Le champ <HICE> est absent',1)

LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90

@@ -467 +467 @@
 
 ! albedo  and radiation parameters
-      INTEGER :: iflag_sic_albedo
-! albedo old or NEMO
       REAL :: alb_sno_dry!=rn_alb_sdry !dry snow albedo
       REAL :: alb_sno_wet!=rn_alb_smlt !wet snow albedo
@@ -511 +509 @@
       IF (check) WRITE (*, *) 'Entering surface_seaice, knon=', knon
 
-      iflag_sic_albedo = iflag_seaice_alb
-
-      IF (iflag_seaice .GE. 1) THEN
+      IF (iflag_seaice .GE. 1 .OR. iflag_seaice_alb .GE. 1) THEN
          sno_den = sisno_den !mean snow density, kg/m3
          ice_cond = sice_cond*ice_den !conductivity of ice
@@ -675 +671 @@
 
          ! if iflag_seaice = 0, we fix the sea-ice thickness to a constant value
-         ! this only serves for albedo calculation if iflag_sic_albedo >=1
+         ! this only serves for albedo calculation if iflag_seaice_alb >=1
          DO i = 1, knon
             ki = knindex(i)
@@ -736 +732 @@
 !#ifdef ISO
 
-      ELSE
+      ELSE ! iflag_seaice >= 1
+
+         icesub(:) = 0.
+         icemelt(:) = 0.
 
          ! if iflag_seaice = 1, we fix the sea-ice thickness to a constant value
+         ! NB another idea could be testing the formula by Krinner et al. 1997 : h = (0.2 + 3.8(f**2))(1 + 2(f- f_min))
+         ! Over the Svalbard region in April 2022, G. Nitti found h = (a + b f**2)(1 + 2(f-fmin)) a=0.166, b=0.911 and fmin=0.351
          IF (iflag_seaice .EQ. 1) THEN
                hice(:) = sic_hice_fixed
@@ -755 +756 @@
          END IF
 
-         ! NB another idea could be testing the formula by Krinner et al. 1997 : h = (0.2 + 3.8(f**2))(1 + 2(f- f_min))
-         ! Over the Svalbard region in April 2022, G. Nitti found h = (a + b f**2)(1 + 2(f-fmin)) a=0.166, b=0.911 and fmin=0.351
-
          DO i = 1, knon
 
@@ -767 +765 @@
             END IF
             ! snow and ice sublimation
-            IF (evap(i) > 0.) THEN
-               snow_evap = MIN(snow(i)/dtime, evap(i))
+               snow_evap = MIN(snow(i)/dtime, evap(i)) 
                snow(i) = snow(i) - snow_evap*dtime
                snow(i) = MAX(0.0, snow(i))
-               IF (iflag_seaice == 2) THEN
-                  hice(ki) = MAX(0.0, hice(ki) - (evap(i) - snow_evap)*dtime/ice_den)
+               icesub(i) = evap(i) - snow_evap ! diagnostics
+               ! in case of interactive calculation of hice
+               IF (iflag_seaice == 3) THEN
+                  hice(ki) = MAX(0.0, hice(ki) - (icesub(i))*dtime/ice_den)
                END IF
-            END IF
             ! Melt / Freeze snow from above if Tsurf>0
             IF (tsurf_new(i) .GT. t_melt) THEN
@@ -785 +783 @@
                   snow(i) = snow(i) - e_melt
                   tsurf_new(i) = t_melt
+                  icemelt(i) = 0.
                ELSE ! all snow is melted
-                  ! add remaining heat flux to ice
+                    ! add remaining heat flux to ice
                   e_melt = e_melt - snow(i)
+                  icemelt(i) = e_melt
                   tice(ki) = tice(ki) + e_melt*ice_lat*2./(ice_cap*hice(ki)*ice_den)
                   tsurf_new(i) = tice(ki)
@@ -940 +940 @@
 !****************************************************************************************
 
-      SELECT CASE (iflag_sic_albedo)
+      SELECT CASE (iflag_seaice_alb)
       CASE (0)
          ! default option. Calculation of albedo at snow (alb_neig) and update the age of snow (agesno)