Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/phys_state_var_mod.F90

Timestamp:

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

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/phys_state_var_mod.F90 (modified) (11 diffs)

LMDZ6/branches/contrails/libf/phylmd/phys_state_var_mod.F90

@@ -15 +15 @@
       INTEGER, PARAMETER :: nlevSTD8=8
       INTEGER, PARAMETER :: nlevSTD3=3
-      INTEGER, PARAMETER :: nout=3
+      INTEGER, PARAMETER :: nout=10
       INTEGER, PARAMETER :: napisccp=1
       INTEGER, SAVE :: radpas  ! radiation is called every "radpas" step
@@ -144 +144 @@
       REAL, ALLOCATABLE, SAVE :: cfc_ancien(:,:), qtc_ancien(:,:), nic_ancien(:,:)
 !$OMP THREADPRIVATE(cfc_ancien, qtc_ancien, nic_ancien)
+       REAL, ALLOCATABLE, SAVE :: tke_ancien(:,:)
+!$OMP THREADPRIVATE(tke_ancien)
 !!! RomP >>>
       REAL, ALLOCATABLE, SAVE :: tr_ancien(:,:,:)
@@ -173 +175 @@
 !$OMP THREADPRIVATE(clwcon0th,rnebcon0th)
 ! radiation outputs
+!FC
+      REAL,ALLOCATABLE,SAVE :: lwtoab(:,:), lwtoa0b(:,:)
+!$OMP THREADPRIVATE(lwtoab,lwtoa0b)
+      REAL,ALLOCATABLE,SAVE :: lwtoabp(:,:), lwtoa0bp(:,:)
+!$OMP THREADPRIVATE(lwtoabp,lwtoa0bp)
+! FC
       REAL,ALLOCATABLE,SAVE :: swdnc0(:,:), swdn0(:,:), swdn(:,:)
 !$OMP THREADPRIVATE(swdnc0,swdn0,swdn)
@@ -222 +230 @@
       REAL,ALLOCATABLE,SAVE :: O3sumSTD(:,:,:), O3daysumSTD(:,:,:)
 !$OMP THREADPRIVATE(O3sumSTD,O3daysumSTD) 
+!GG
+      REAL,ALLOCATABLE,SAVE :: hice(:)
+!$OMP THREADPRIVATE(hice)
+      REAL,ALLOCATABLE,SAVE :: tice(:)
+!$OMP THREADPRIVATE(tice)
+      REAL,ALLOCATABLE,SAVE :: bilg_cumul(:)
+!$OMP THREADPRIVATE(bilg_cumul)
+!GG
 !IM begin
       REAL,ALLOCATABLE,SAVE :: wlevSTD(:,:), ulevSTD(:,:), vlevSTD(:,:)
@@ -471 +487 @@
       REAL,ALLOCATABLE,SAVE :: swdnc0_s2(:,:), swdn0_s2(:,:), swdn_s2(:,:)
       REAL,ALLOCATABLE,SAVE :: swupc0_s2(:,:), swup0_s2(:,:), swup_s2(:,:)
-!$OMP THREADPRIVATE(swdnc0_s2, swdn0_s2, swdn_s2, swupc0_s2, swup0_s2, swup_s2)
+!$OMP THREADPRIVATE(swdnc0_s2, swdn0_s2, swdn_s2 , swupc0_s2, swup0_s2, swup_s2)
+!FC
+      REAL,ALLOCATABLE,SAVE :: lwtoab_s2(:,:), lwtoa0b_s2(:,:)
+!$OMP THREADPRIVATE(lwtoab_s2,lwtoa0b_s2)
+
 
 ! pbase : cloud base pressure
@@ -662 +682 @@
       ALLOCATE(cfc_ancien(klon,klev), qtc_ancien(klon,klev), nic_ancien(klon,klev))
       ALLOCATE(qvcon(klon,klev), qccon(klon,klev))
+      ALLOCATE(tke_ancien(klon,klev+1))
 !!! Rom P >>>
       ALLOCATE(tr_ancien(klon,klev,nbtr))
@@ -685 +706 @@
       ALLOCATE(clwcon0th(klon,klev),rnebcon0th(klon,klev))
 ! radiation outputs
+!FC
+      ALLOCATE(lwtoab(klon,nbands_lw_rrtm), lwtoa0b(klon,nbands_lw_rrtm))
+      lwtoab(:,:)=0.0
+      lwtoa0b(:,:)=0.0
+      ALLOCATE(lwtoab_s2(klon,nbands_lw_rrtm), lwtoa0b_s2(klon,nbands_lw_rrtm))
+      lwtoab_s2(:,:)=0.0
+      lwtoa0b_s2(:,:)=0.0
+      ALLOCATE(lwtoabp(klon,nbands_lw_rrtm), lwtoa0bp(klon,nbands_lw_rrtm))
+      lwtoabp(:,:)=0.0
+      lwtoa0bp(:,:)=0.0
+!FC
       ALLOCATE(swdnc0(klon,klevp1), swdn0(klon,klevp1), swdn(klon,klevp1))
       ALLOCATE(swupc0(klon,klevp1), swup0(klon,klevp1), swup(klon,klevp1))
@@ -726 +758 @@
       ALLOCATE(zuthe(klon),zvthe(klon))
       ALLOCATE(alb_neig(klon))
+!GG
+      ALLOCATE(hice(klon))
+      hice(:) = 0.
+      ALLOCATE(tice(klon))
+      tice(:) = 0.
+      ALLOCATE(bilg_cumul(klon))
+      bilg_cumul(:) = 0.
+!GG
 !cloud base mass flux
       ALLOCATE(ema_cbmf(klon))
@@ -913 +953 @@
       DEALLOCATE(cfc_ancien, qtc_ancien, nic_ancien)
       DEALLOCATE(qvcon, qccon)
+      DEALLOCATE(tke_ancien)
       DEALLOCATE(tr_ancien)                           !RomP
       DEALLOCATE(ratqs, pbl_tke,coefh,coefm)
@@ -921 +962 @@
       DEALLOCATE(clwcon0th, rnebcon0th)
 ! radiation outputs
+!FC 
+      DEALLOCATE(lwtoab, lwtoa0b)
+      DEALLOCATE(lwtoab_s2, lwtoa0b_s2)
+      DEALLOCATE(lwtoabp, lwtoa0bp)
+!FC
       DEALLOCATE(swdnc0, swdn0, swdn)
       DEALLOCATE(swupc0, swup0, swup)
@@ -954 +1000 @@
       DEALLOCATE(zuthe, zvthe)
       DEALLOCATE(alb_neig)
+!GG
+      DEALLOCATE(hice)
+      DEALLOCATE(tice)
+      DEALLOCATE(bilg_cumul)
+!GG
       DEALLOCATE(ema_cbmf)
       DEALLOCATE(ema_pcb, ema_pct)