Changeset 2240 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Mar 23, 2015, 5:07:43 PM (9 years ago)

Author:

fhourdin

Message:

Revisite de la formule des flux de surface
(en priorité sur l'océan) en tenant compte des bourrasques de
vent et de la différence entre les hauteurs de rugosités pour
la quantité de mouvement, l'enthalpie et éventuellement l'humidité.

Etape 1 :
Introduction d'un calcul de gustiness dans la physique
gustiness(:)=f_gust_bl * ale_bl + f_gust_wk * ame_wk
Cette variable est passée ensuite jusqu'au fin fond de la couche limite.
L'étape 1 est prête à commettre, ne nécessite pas de nouvelles
variables dans les startphy et assure la convergence numérique.

Introduction of gustiness in the surface flux computation.
Gustiness is computed from as
gustiness(:)=f_gust_bl * ale_bl + f_gust_wk * ame_wk
and pass through pbl_surface down to the routines that compute
surface fluxes.

File:

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

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1841 +1841 @@
 !>jyg+nrlmd
 !
+!-------gustiness calculation-------!
+     do i = 1, klon
+        gustiness(i)=f_gust_bl*ale_bl(i)+f_gust_wk*ale_wake(i)
+     enddo
+
      CALL pbl_surface(  &
           dtime,     date0,     itap,    days_elapsed+1, &
@@ -1847 +1852 @@
           zsig,      sollwdown, pphi,    cldt,      &
           rain_fall, snow_fall, solsw,   sollw,     &
+          gustiness,                                &
           t_seri,    q_seri,    u_seri,  v_seri,    &
 !nrlmd+jyg<
@@ -2151 +2157 @@
         ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Estimation d'une vitesse verticale effective pour ALP
+        if (1==0) THEN
         www(1:klon)=0.
         do k=2,klev-1
@@ -2166 +2173 @@
            if (www(i)>0. .and. ale_bl(i)>0. ) www(i)=www(i)/ale_bl(i)
         enddo
+        ENDIF
 
 
@@ -2178 +2186 @@
               ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
            else
+       abort_message ='Ne pas passer la car www non calcule'
+       CALL abort_gcm (modname,abort_message,1)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -4313 +4323 @@
   ENDDO
   RETURN
-END SUBROUTINE gr_fi_ecrit
+  END SUBROUTINE gr_fi_ecrit
+