Changeset 5669 for LMDZ6/trunk/libf/phylmd/dyn1d

Timestamp:

May 21, 2025, 2:48:11 PM (7 weeks ago)

Author:

Laurent Fairhead

Message:

Initialisation problem for 1D

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• old_lmdz1d.f90 (modified) (3 diffs)
• scm.f90 (modified) (2 diffs)

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

@@ -25 +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, &
-       hice, tice, bilg_cumul
+       u10m,v10m,ale_wake,ale_bl_stat
 
 
@@ -220 +219 @@
       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)
 !     real :: agesno(1,nbsrf)
@@ -827 +827 @@
   !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
+     hice(:)=1.0
+     tice(:)=ftsol(:,is_sic)
+     bilg_cumul(:)=0.0
+
 
   !END IF

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

@@ -25 +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_, &
-       hice, tice, bilg_cumul
+       u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_
 
 
@@ -183 +182 @@
       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