Changeset 3856 for LMDZ6/trunk/libf/phylmd/calcratqs.F90

Timestamp:

Feb 25, 2021, 7:02:16 PM (3 years ago)

Author:

fhourdin

Message:

Introduction d'un ratqs (largeur relative de la distribution sous maille de l'eau)
interactif (travail de Louis D'Alençon, Frédéric Hourdin et Catherine Rio).
Contrôlé par le flag iflag_ratqs=10 (précédemment 4).
N'affecte pas les calculs si iflag_ratqs=4.
Fredho

File:

• LMDZ6/trunk/libf/phylmd/calcratqs.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/calcratqs.F90

@@ -2 +2 @@
            iflag_ratqs,iflag_con,iflag_cld_th,pdtphys, &
            ratqsbas,ratqshaut,ratqsp0,ratqsdp, &
-           tau_ratqs,fact_cldcon,   &
+           tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
            ptconv,ptconvth,clwcon0th, rnebcon0th,      &
            paprs,pplay,q_seri,zqsat,fm_therm,          &
-           ratqs,ratqsc)
+           ratqs,ratqsc,ratqs_inter)
 
 implicit none
@@ -27 +27 @@
 logical, dimension(klon,klev),intent(in) :: ptconv
 real, dimension(klon,klev),intent(in) :: rnebcon0th,clwcon0th
-
+real, dimension(klon,klev),intent(in) :: wake_deltaq,wake_s
 ! Output
-real, dimension(klon,klev),intent(inout) :: ratqs,ratqsc
+real, dimension(klon,klev),intent(inout) :: ratqs,ratqsc,ratqs_inter
 logical, dimension(klon,klev),intent(inout) :: ptconvth
 
@@ -124 +124 @@
          enddo
 
-      else if (iflag_ratqs==4) then
+      else if (iflag_ratqs==4) then 
          do k=1,klev
            ratqss(:,k)=ratqsbas+0.5*(ratqshaut-ratqsbas) &
@@ -131 +131 @@
          enddo
 
+      else if (iflag_ratqs==10) then ! ratqs interactif dépendant de la présence de poches froides
+         call calcratqs_inter(klon,klev,pdtphys,ratqsbas,wake_deltaq,wake_s,q_seri,ratqs_inter)
+         do k=1,klev
+            do i=1, klon
+                ratqss(i,k)=ratqs_inter(i,k)+0.5*(ratqshaut-ratqs_inter(i,k)) &
+                *( tanh( (ratqsp0-pplay(i,k))/ratqsdp) + 1.)
+            enddo
+         enddo
+
       endif
-
-