Changeset 2240 for LMDZ5/trunk/libf/phylmd/surf_land_mod.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/surf_land_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/surf_land_mod.F90

@@ -14 +14 @@
        AcoefH, AcoefQ, BcoefH, BcoefQ, &
        AcoefU, AcoefV, BcoefU, BcoefV, & 
-       pref, u1, v1, rugoro, pctsrf, &
+       pref, u1, v1, gustiness, rugoro, pctsrf, &
        lwdown_m, q2m, t2m, &
        snow, qsol, agesno, tsoil, &
@@ -61 +61 @@
     REAL, DIMENSION(klon), INTENT(IN)       :: AcoefU, AcoefV, BcoefU, BcoefV
     REAL, DIMENSION(klon), INTENT(IN)       :: pref   ! pressure reference
-    REAL, DIMENSION(klon), INTENT(IN)       :: u1, v1
+    REAL, DIMENSION(klon), INTENT(IN)       :: u1, v1, gustiness
     REAL, DIMENSION(klon), INTENT(IN)       :: rugoro
     REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
@@ -137 +137 @@
             knindex, rlon, rlat, pctsrf, &
             debut, lafin, &
-            zlev,  u1, v1, temp_air, spechum, epot_air, ccanopy, & 
+            zlev,  u1, v1, gustiness, temp_air, spechum, epot_air, ccanopy, & 
             cdragh, AcoefH, AcoefQ, BcoefH, BcoefQ, &
             precip_rain, precip_snow, lwdown_m, swnet, swdown, &
@@ -160 +160 @@
             tsurf, p1lay, cdragh, precip_rain, precip_snow, temp_air, &
             spechum, AcoefH, AcoefQ, BcoefH, BcoefQ, pref, &
-            u1, v1, rugoro, swnet, lwnet, &
+            u1, v1, gustiness, rugoro, swnet, lwnet, &
             snow, qsol, agesno, tsoil, &
             qsurf, z0_new, alb1_new, alb2_new, evap, &
@@ -175 +175 @@
     v0(:)=0.0
     CALL calcul_flux_wind(knon, dtime, &
-         u0, v0, u1, v1, cdragm, &
+         u0, v0, u1, v1, gustiness, cdragm, &
          AcoefU, AcoefV, BcoefU, BcoefV, &
          p1lay, temp_air, &