Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/dyn1d/scm.f90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/dyn1d/scm.f90 (modified) (4 diffs)

LMDZ6/branches/contrails/libf/phylmd/dyn1d/scm.f90

@@ -5 +5 @@
    USE clesphys_mod_h
    USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin
+   USE phyetat0_get_mod, ONLY : phyetat0_get, phyetat0_srf
    USE phys_state_var_mod, ONLY : phys_state_var_init, phys_state_var_end, &
        clwcon, detr_therm, &
@@ -182 +183 @@
       real :: run_off_lic_0(1)
       real :: fder(1),snsrf(1,nbsrf),qsurfsrf(1,nbsrf)
+      real :: hice(1), tice(1), bilg_cumul(1)
       real :: tsoil(1,nsoilmx,nbsrf)
       ! AM
@@ -232 +234 @@
       CHARACTER(len=80) :: abort_message
       CHARACTER(len=20) :: modname = 'scm'
+      LOGICAL :: found
 
       print*,'VOUS ENTREZ DANS LE 1D FORMAT STANDARD'
@@ -614 +617 @@
 
 !-----------------------------------------------------------------------
-        call pbl_surface_init(fder, snsrf, qsurfsrf, tsoil)
+  !GG
+  ! Sea ice
+  !IF (iflag_seaice == 2) THEN
+
+  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
+  END IF
+  found=phyetat0_get(tice,"tice","Sea Ice temperature",0.)
+  IF (.NOT. found) THEN
+       PRINT*, "phyetat0: Le champ <tice> est absent"
+       PRINT*, "Initialisation a tsol_sic"
+       tice(:)=ftsol(:,is_sic)
+  END IF
+  found=phyetat0_get(bilg_cumul,"bilg_cumul","Flux conductivite + transmit sea-ice",0.)
+  IF (.NOT. found) THEN
+       PRINT*, "phyetat0: Le champ <bilg_cumul> est absent"
+       PRINT*, "Initialisation a zero"
+       bilg_cumul(:)=0.0
+  END IF
+
+  !END IF
+  !GG
+        call pbl_surface_init(fder, snsrf, qsurfsrf, tsoil, hice, tice, bilg_cumul)
 
 !------------------ prepare limit conditions for limit.nc -----------------