Changeset 5668 for LMDZ6

Timestamp:

May 21, 2025, 11:54:30 AM (3 weeks ago)

Author:

Laurent Fairhead

Message:

Correction of bug for revision 5662 for the 1D

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• old_lmdz1d.f90 (modified) (4 diffs)
• scm.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.f90

@@ -5 +5 @@
       SUBROUTINE old_lmdz1d
 
-USE flux_arp_mod_h
+   USE flux_arp_mod_h
       USE compbl_mod_h
          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, &
@@ -24 +25 @@
        ql_ancien, qs_ancien, qbs_ancien, cf_ancien, rvc_ancien, &
        prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
-       u10m,v10m,ale_wake,ale_bl_stat
+       u10m,v10m,ale_wake,ale_bl_stat, &
+       hice, tice, bilg_cumul
 
 
@@ -266 +268 @@
 !                           <> 0, tendencies of forcing are not added
       INTEGER :: flag_inhib_forcing = 0
+      LOGICAL :: found
 
 !=====================================================================
@@ -820 +823 @@
 !       tsoil(1,11,1)=308.00
 !-----------------------------------------------------------------------
-        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 -----------------

LMDZ6/trunk/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, &
@@ -24 +25 @@
        ql_ancien, qs_ancien, qbs_ancien, cf_ancien, rvc_ancien, &
        prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
-       u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_
+       u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_, &
+       hice, tice, bilg_cumul
 
 
@@ -231 +233 @@
       CHARACTER(len=80) :: abort_message
       CHARACTER(len=20) :: modname = 'scm'
+      LOGICAL :: found
 
       print*,'VOUS ENTREZ DANS LE 1D FORMAT STANDARD'
@@ -613 +616 @@
 
 !-----------------------------------------------------------------------
-        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 -----------------