Changeset 4613 for LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jul 7, 2023, 4:01:37 AM (12 months ago)

Author:

fhourdin

Message:

Integration/replay_isation travail Louis (ratqs)

File:

• LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -102 +102 @@
       REAL, ALLOCATABLE, SAVE :: rneb_ancien(:,:)
 !$OMP THREADPRIVATE(rneb_ancien)
-      REAL, ALLOCATABLE, SAVE :: qtc_cv(:,:),sigt_cv(:,:)
-!$OMP THREADPRIVATE(qtc_cv,sigt_cv)
+      REAL, ALLOCATABLE, SAVE :: qtc_cv(:,:),sigt_cv(:,:),detrain_cv(:,:),fm_cv(:,:)
+!$OMP THREADPRIVATE(qtc_cv,sigt_cv,detrain_cv,fm_cv)
       REAL, ALLOCATABLE, SAVE :: ratqs(:,:)
 !$OMP THREADPRIVATE(ratqs)
@@ -229 +229 @@
 !$OMP THREADPRIVATE(ftd)
 ! fqd : convective moistening due to unsaturated downdraughts
-      REAL,ALLOCATABLE,SAVE :: fqd(:,:)     
-!$OMP THREADPRIVATE(fqd)
+      REAL,ALLOCATABLE,SAVE :: fqd(:,:),fqcomp(:,:)     
+!$OMP THREADPRIVATE(fqd,fqcomp)
 #ifdef ISO
       REAL, ALLOCATABLE, SAVE :: fxtd(:,:,:)
@@ -454 +454 @@
 !$OMP THREADPRIVATE(ale_bl_trig)
 
-      REAL,SAVE,ALLOCATABLE :: ratqs_inter(:,:)
-!$OMP THREADPRIVATE(ratqs_inter)
+      REAL,SAVE,ALLOCATABLE :: ratqs_inter_(:,:)
+!$OMP THREADPRIVATE(ratqs_inter_)
 
 #ifdef ISO
@@ -567 +567 @@
       ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev))
       ALLOCATE(rneb_ancien(klon,klev))
-      ALLOCATE(qtc_cv(klon,klev),sigt_cv(klon,klev))
+      ALLOCATE(qtc_cv(klon,klev),sigt_cv(klon,klev),detrain_cv(klon,klev),fm_cv(klon,klev))
       ALLOCATE(ratqs(klon,klev))
       ALLOCATE(pbl_tke(klon,klev+1,nbsrf+1))
@@ -638 +638 @@
       ALLOCATE(sigd(klon))
       ALLOCATE(cin(klon), ALE(klon), ALP(klon))
-      ALLOCATE(ftd(klon,klev), fqd(klon,klev))
+      ALLOCATE(ftd(klon,klev), fqd(klon,klev),fqcomp(klon,klev))
       ALLOCATE(Ale_bl(klon))
       ALLOCATE(ale_wake(klon))
@@ -733 +733 @@
 
       ALLOCATE(ale_bl_trig(klon))
-      ALLOCATE(ratqs_inter(klon,klev))
+      ALLOCATE(ratqs_inter_(klon,klev))
       IF (ok_gwd_rando) THEN
         ALLOCATE(du_gwd_rando(klon, klev))
@@ -773 +773 @@
       DEALLOCATE(qs_ancien, ql_ancien, qbs_ancien, rneb_ancien)
       DEALLOCATE(prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien)
-      DEALLOCATE(qtc_cv,sigt_cv)
+      DEALLOCATE(qtc_cv,sigt_cv,detrain_cv,fm_cv)
       DEALLOCATE(u_ancien, v_ancien)
       DEALLOCATE(tr_ancien)                           !RomP
@@ -821 +821 @@
       DEALLOCATE(wd, sigd)
       DEALLOCATE(cin, ALE, ALP)
-      DEALLOCATE(ftd, fqd)
+      DEALLOCATE(ftd, fqd, fqcomp)
       DEALLOCATE(Ale_bl, Alp_bl)
       DEALLOCATE(ale_wake)
@@ -886 +886 @@
       if (.not. ok_hines .and. ok_gwd_rando) DEALLOCATE(du_gwd_front)
       DEALLOCATE(ale_bl_trig)
-      DEALLOCATE(ratqs_inter)
+      DEALLOCATE(ratqs_inter_)
 
       if (activate_ocean_skin >= 1) then