Changeset 2497

Timestamp:

Apr 18, 2016, 9:57:11 PM (8 years ago)

Author:

fhourdin

Message:

Declaration de ql_ancien et qs_ancien.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phyetat0.F90 (modified) (2 diffs)
• phyredem.F90 (modified) (2 diffs)
• phys_state_var_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -13 +13 @@
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
        falb_dir, falb_dif, &
-       ftsol, pbl_tke, pctsrf, q_ancien, radpas, radsol, rain_fall, ratqs, &
+       ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, radpas, radsol, rain_fall, ratqs, &
        rnebcon, rugoro, sig1, snow_fall, solaire_etat0, sollw, sollwdown, &
        solsw, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
@@ -322 +322 @@
 
   ancien_ok=.true.
+  ql_ancien=0.;qs_ancien=0.
   ancien_ok=ancien_ok.AND.phyetat0_get(klev,t_ancien,"TANCIEN","TANCIEN",0.)
   ancien_ok=ancien_ok.AND.phyetat0_get(klev,q_ancien,"QANCIEN","QANCIEN",0.)
+  ancien_ok=ancien_ok.AND.phyetat0_get(klev,ql_ancien,"QLANCIEN","QLANCIEN",0.)
+  ancien_ok=ancien_ok.AND.phyetat0_get(klev,qs_ancien,"QSANCIEN","QSANCIEN",0.)
   ancien_ok=ancien_ok.AND.phyetat0_get(klev,u_ancien,"UANCIEN","UANCIEN",0.)
   ancien_ok=ancien_ok.AND.phyetat0_get(klev,v_ancien,"VANCIEN","VANCIEN",0.)

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -13 +13 @@
                                 sollwdown, rain_fall, snow_fall, z0m, z0h, &
                                 agesno, zmea, zstd, zsig, zgam, zthe, zpic, &
-                                zval, rugoro, t_ancien, q_ancien, u_ancien, &
+                                zval, rugoro, t_ancien, q_ancien, ql_ancien, qs_ancien,  u_ancien, &
                                 v_ancien, clwcon, rnebcon, ratqs, pbl_tke, &
                                 wake_delta_pbl_tke, zmax0, f0, sig1, w01, &
@@ -200 +200 @@
 
   CALL put_field("QANCIEN", "QANCIEN", q_ancien)
+
+  CALL put_field("QLANCIEN", "QLANCIEN", ql_ancien)
+
+  CALL put_field("QSANCIEN", "QSANCIEN", qs_ancien)
 
   CALL put_field("UANCIEN", "", u_ancien)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -64 +64 @@
       REAL, ALLOCATABLE, SAVE :: t_ancien(:,:), q_ancien(:,:)
 !$OMP THREADPRIVATE(t_ancien, q_ancien)
+      REAL, ALLOCATABLE, SAVE :: ql_ancien(:,:), qs_ancien(:,:)
+!$OMP THREADPRIVATE(ql_ancien, qs_ancien)
       REAL, ALLOCATABLE, SAVE :: u_ancien(:,:), v_ancien(:,:)
 !$OMP THREADPRIVATE(u_ancien, v_ancien)
@@ -441 +443 @@
       ALLOCATE(rugoro(klon))
       ALLOCATE(t_ancien(klon,klev), q_ancien(klon,klev))
+      ALLOCATE(ql_ancien(klon,klev), qs_ancien(klon,klev))
       ALLOCATE(u_ancien(klon,klev), v_ancien(klon,klev))
 !!! Rom P >>>
@@ -601 +604 @@
       deallocate(zthe, zpic, zval)
       deallocate(rugoro, t_ancien, q_ancien, clwcon, rnebcon)
+      deallocate(qs_ancien, ql_ancien)
       deallocate(qtc_cv,sigt_cv)
       deallocate(        u_ancien, v_ancien                 )