Changeset 2306 for LMDZ5/trunk/libf/phylmd/concvl.F90

Timestamp:

Jun 23, 2015, 11:33:09 AM (9 years ago)

Author:

jyg

Message:

Improved diagnostics: pmflxr and pmflxs are now
the true vertical profiles of liquid and solid
convective precipitation (previously they where
merely diagnosed through a test on temperature).
The convective scheme internal variable for ice
precipitation is Vprecipi.

File:

• LMDZ5/trunk/libf/phylmd/concvl.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/concvl.F90

@@ -60 +60 @@
 ! Ma---------output-R-adiabatic ascent mass flux (kg/m2/s)
 ! mip--------output-R-mass flux shed by adiabatic ascent (kg/m2/s)
-! Vprecip----output-R-vertical profile of precipitations (kg/m2/s)
+! Vprecip----output-R-vertical profile of total precipitation (kg/m2/s)
 ! Tconv------output-R-environment temperature seen by convective scheme (K)
 ! Cape-------output-R-CAPE (J/kg)
@@ -109 +109 @@
 !!       REAL Ma(klon,klev), mip(klon,klev),Vprecip(klon,klev)     !jyg
   REAL Ma(klon, klev), mip(klon, klev), Vprecip(klon, klev+1)      !jyg
+  REAL Vprecipi(klon, klev+1)                                      !jyg
   REAL wght(klon, klev)                                            !RL
 
@@ -415 +416 @@
                     iflag, d_t, d_q, d_u, d_v, d_tra, rain, kbas, ktop, &
                     cbmf, plcl, plfc, wbeff, sig1, w01, ptop2, sigd, &
-                    Ma, mip, Vprecip, upwd, dnwd, dnwdbis, qcondc, wd, &
+                    Ma, mip, Vprecip, Vprecipi, upwd, dnwd, dnwdbis, qcondc, wd, &
                     cape, cin, tvp, &
                     dd_t, dd_q, plim1, plim2, asupmax, supmax0, &
@@ -476 +477 @@
     END DO
   END DO
-!                                                   !jyg
-! --Separation neige/pluie (pour diagnostics)       !jyg
-  DO k = 1, klev                                    !jyg
-    DO i = 1, klon                                  !jyg
-      IF (t1(i,k)<rtt) THEN                         !jyg
-        pmflxs(i, k) = Vprecip(i, k)                !jyg
-      ELSE                                          !jyg
-        pmflxr(i, k) = Vprecip(i, k)                !jyg
-      END IF                                        !jyg
-    END DO                                          !jyg
-  END DO                                            !jyg
+!                                                     !jyg
+  IF (iflag_con == 30 .OR. iflag_ice_thermo ==0) THEN
+! --Separation neige/pluie (pour diagnostics)         !jyg
+    DO k = 1, klev                                    !jyg
+      DO i = 1, klon                                  !jyg
+        IF (t1(i,k)<rtt) THEN                         !jyg
+          pmflxs(i, k) = Vprecip(i, k)                !jyg
+        ELSE                                          !jyg
+          pmflxr(i, k) = Vprecip(i, k)                !jyg
+        END IF                                        !jyg
+      END DO                                          !jyg
+    END DO                                            !jyg
+  ELSE
+    DO k = 1, klev                                    !jyg
+      DO i = 1, klon                                  !jyg
+        pmflxs(i, k) = Vprecipi(i, k)                 !jyg
+        pmflxr(i, k) = Vprecip(i, k)-Vprecipi(i, k)   !jyg
+      END DO                                          !jyg
+    END DO                                            !jyg
+  ENDIF
 
 ! c      IF (if_ebil.ge.2) THEN