Changeset 2105 for LMDZ5/trunk/libf

Timestamp:

Aug 7, 2014, 9:14:27 AM (10 years ago)

Author:

fhourdin

Message:

Modification des options pour inclure un aspect grande échelle
dans le calcul de Alp et pour tenir compte des précpitatios
sratiformes pour alimenter les poches froides.
Contrôlés respectivement par les flags alp_ofet=X, avec X<0.
(la valeur par défaut est 0) et iflag_wake=3.

Pour le calcul de Alp :
! Ajout d'une composante 3 * A * w w'2 a w'3 avec w=www : w max sous pbase
! ou A est la fraction couverte par les ascendances w'
! on utilise le fait que A * w'3 = ALP
! et donc A * w'2 ~ ALP / sqrt(ALE) (on ajoute 0.1 pour les
! singularites)

Modified options to account for large scale velocity in ALP
and for feeding of wakes by large scale rainfall.
Controled respectively by alp_offset=X, with X<0. (default being 0) and iflag_wake=3.

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -388 +388 @@
   REAL dt_a(klon,klev)
   REAL dq_a(klon,klev)
+  REAL, dimension(klon) :: www
   REAL, SAVE :: alp_offset
   !$OMP THREADPRIVATE(alp_offset)
@@ -2010 +2011 @@
         ! proportionnel ??a w si <0
         ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Estimation d'une vitesse verticale effective pour ALP
+        www(1:klon)=0.
+        do k=2,klev-1
+           do i=1,klon
+              www(i)=max(www(i),-omega(i,k)*RD*t_seri(i,k)/(RG*paprs(i,k)) &
+&                    *zw2(i,k)*zw2(i,k))
+!             if (paprs(i,k)>pbase(i)) then
+! calcul approche de la vitesse verticale en m/s
+!                www(i)=max(www(i),-omega(i,k)*RD*temp(i,k)/(RG*paprs(i,k))
+!             endif
+!   Le 0.1 est en gros H / ps = 1e5 / 1e4
+           enddo
+        enddo
+        do i=1,klon
+           if (www(i)>0. .and. ale_bl(i)>0. ) www(i)=www(i)/ale_bl(i)
+        enddo
+
+
         do i = 1,klon
            ALE(i) = max(ale_wake(i),Ale_bl(i))
@@ -2020 +2039 @@
               ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
            else
-              ALP(i)=alp_wake(i)+Alp_bl(i)+alp_offset*min(omega(i,6),0.)
-              if (alp(i)<0.) then
-                 print*,'ALP ',alp(i),alp_wake(i) &
-                      ,Alp_bl(i),alp_offset*min(omega(i,6),0.)
-              endif
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!                                _                  _
+! Ajout d'une composante 3 * A * w w'2 a  w'3  avec w=www : w max sous pbase
+!       ou A est la fraction couverte par les ascendances w'
+!       on utilise le fait que A * w'3 = ALP
+!       et donc A * w'2 ~ ALP / sqrt(ALE)  (on ajoute 0.1 pour les
+!       singularites)
+             ALP(i)=alp_wake(i)*(1.+3.*www(i)/( sqrt(ale_wake(i))+0.1) ) &
+  &                +alp_bl(i)  *(1.+3.*www(i)/( sqrt(ale_bl(i))  +0.1) )
+!             ALP(i)=alp_wake(i)+Alp_bl(i)+alp_offset*min(omega(i,6),0.)
+!             if (alp(i)<0.) then
+!                print*,'ALP ',alp(i),alp_wake(i) &
+!                     ,Alp_bl(i),alp_offset*min(omega(i,6),0.)
+!             endif
            endif
         enddo
-        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
      endif
@@ -2258 +2286 @@
      ENDDO
 
-     if (iflag_wake==2) then
+     IF (iflag_wake==2) THEN
         ok_wk_lsp(:)=max(sign(1.,wake_s(:)-wake_s_min_lsp),0.)
         DO k = 1,klev
@@ -2266 +2294 @@
                 ok_wk_lsp(:)*(d_q_eva(:,k)+d_q_lsc(:,k))/dtime
         ENDDO
-     endif
+     ELSEIF (iflag_wake==3) THEN
+        ok_wk_lsp(:)=max(sign(1.,wake_s(:)-wake_s_min_lsp),0.)
+        DO k = 1,klev
+           DO i=1,klon
+              IF (rneb(i,k)==0.) THEN
+! On ne tient compte des tendances qu'en dehors des nuages (c'est �|  dire
+! a priri dans une region ou l'eau se reevapore).
+                dt_dwn(i,k)= dt_dwn(i,k)+ &
+                ok_wk_lsp(i)*d_t_lsc(i,k)/dtime
+                dq_dwn(i,k)= dq_dwn(i,k)+ &
+                ok_wk_lsp(i)*d_q_lsc(i,k)/dtime
+              ENDIF
+           ENDDO
+        ENDDO
+     ENDIF
+
      !
      !calcul caracteristiques de la poche froide