Changeset 199 for LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F

Timestamp:

Apr 2, 2001, 11:12:13 AM (23 years ago)

Author:

lmdz

Message:

Introduction Kerry Emmanuel
ratqs interactif
mise a jour offline
FH
LF

File:

• LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F (modified) (3 diffs)

LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F

@@ -1 @@
-      SUBROUTINE fisrtilp_tr(dtime,paprs,pplay,t,q,
+c $Header
+c
+      SUBROUTINE fisrtilp_tr(dtime,paprs,pplay,t,q,ratqs,
      s                   d_t, d_q, d_ql, rneb, radliq, rain, snow,
      s                   pfrac_impa, pfrac_nucl, pfrac_1nucl,
@@ -66 +68 @@
       PARAMETER (coef_eva=2.0E-05)
       LOGICAL calcrat ! calculer ratqs au lieu de fixer sa valeur
-      REAL ratqs ! determine la largeur de distribution de vapeur
+      REAL ratqs(klon,klev) ! determine la largeur de distribution de vapeur
       PARAMETER (calcrat=.TRUE.)
       REAL zx_min, rat_max
@@ -281 +283 @@
          DO i = 1, klon
 c
-            zx = pplay(i,k)/paprs(i,1)
-            zx = (zx_max-zx)/(zx_max-zx_min)
-            zx = MIN(MAX(zx,0.0),1.0)
-            zx = zx * zx * zx
-            ratqs = zx * (rat_max-rat_min) + rat_min
-            IF (.NOT.calcrat) ratqs=0.05
-c
-            zdelq = ratqs * zq(i)
+            zdelq = ratqs(i,k) * zq(i)
             rneb(i,k) = (zq(i)+zdelq-zqs(i)) / (2.0*zdelq)
             zqn(i) = (zq(i)+zdelq+zqs(i))/2.0