Changeset 5669 for LMDZ6/trunk

Timestamp:

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

Author:

Laurent Fairhead

Message:

Initialisation problem for 1D

Location:

LMDZ6/trunk/libf/phylmd

Files:

• dyn1d/old_lmdz1d.f90 (modified) (3 diffs)
• dyn1d/scm.f90 (modified) (2 diffs)
• pbl_surface_mod.F90 (modified) (3 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

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -121 +121 @@
 !
 !****************************************************************************************    
+
     ALLOCATE(fder(klon), stat=ierr)
     IF (ierr /= 0) CALL abort_physic('pbl_surface_init', 'pb in allocation',1)
@@ -126 +127 @@
 !GG
     ALLOCATE(hice(klon), stat=ierr)
-    IF (ierr /= 0) CALL abort_physic('pbl_surface_init', 'pb in allocation',1)
+    IF (ierr /= 0) CALL abort_physic('pbl_surface_init hice', 'pb in allocation',1)
 
     ALLOCATE(tice(klon), stat=ierr)
-    IF (ierr /= 0) CALL abort_physic('pbl_surface_init', 'pb in allocation',1)
+    IF (ierr /= 0) CALL abort_physic('pbl_surface_init tice', 'pb in allocation',1)
 
     ALLOCATE(bilg_cumul(klon), stat=ierr)
-    IF (ierr /= 0) CALL abort_physic('pbl_surface_init', 'pb in allocation',1)
+    IF (ierr /= 0) CALL abort_physic('pbl_surface_init bilg', 'pb in allocation',1)
 !GG
 
@@ -4442 +4443 @@
 !   DEALLOCATE(qsol, fder, snow, qsurf, evap, rugos, agesno, ftsoil)
     IF (ALLOCATED(fder)) DEALLOCATE(fder)
+    IF (ALLOCATED(hice)) DEALLOCATE(hice)
+    IF (ALLOCATED(tice)) DEALLOCATE(tice)
+    IF (ALLOCATED(bilg_cumul)) DEALLOCATE(bilg_cumul)
     IF (ALLOCATED(snow)) DEALLOCATE(snow)
     IF (ALLOCATED(qsurf)) DEALLOCATE(qsurf)