Changeset 6070 for LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

Timestamp:

Feb 6, 2026, 10:58:17 AM (2 months ago)

Author:

evignon

Message:

reorganisation de la routine ocean_forced_mod pour gerer les nouvelles
options pour le traitement de la banquise en force.
Travail issu des echanges avec G. Gastineau en vu des prochains stages de master

File:

• LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

@@ -18 +18 @@
 !GG  USE surface_data,     ONLY : type_ocean, version_ocean
   USE surface_data,     ONLY : type_ocean, version_ocean, iflag_seaice, &
-                               iflag_seaice_alb, iflag_leads
+                               iflag_seaice_alb, iflag_leads, sic_hice_fixed
 !GG
   USE phyetat0_get_mod, ONLY : phyetat0_get, phyetat0_srf
@@ -725 +725 @@
   !IF (iflag_seaice == 2) THEN
 
+  ! initial ice thickness (if not read in the startphy) could be included in the .def files
   found=phyetat0_get(hice,"hice","Ice thickness",0.)
   IF (.NOT. found) THEN
        PRINT*, "phyetat0: Le champ <hice> est absent"
        PRINT*, "Initialisation a hice=1m "
-       hice(:)=1.0
+       hice(:)=sic_hice_fixed
   END IF
   found=phyetat0_get(tice,"tice","Sea Ice temperature",0.)