Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2209 r2240  
    3232         o_t2m_min_mon, o_t2m_max_mon, &
    3333         o_q2m, o_ustar, o_u10m, o_v10m, &
    34          o_wind10m, o_wind10max, o_sicf, &
     34         o_wind10m, o_wind10max, o_gusts, o_sicf, &
    3535         o_psol, o_mass, o_qsurf, o_qsol, &
    3636         o_precip, o_ndayrain, o_plul, o_pluc, &
     
    159159         SWdn200, SWdn200clr, LWup200, LWup200clr, &
    160160         LWdn200, LWdn200clr, solsw, solsw0, sollw, &
    161          radsol, sollw0, sollwdown, sollw, &
     161         radsol, sollw0, sollwdown, sollw, gustiness, &
    162162         sollwdownclr, lwdn0, ftsol, ustar, u10m, &
    163163         v10m, pbl_tke, wake_delta_pbl_TKE, &
     
    354354       ENDIF
    355355       CALL histwrite_phy(o_wind10max, zx_tmp_fi2d)
     356
     357       CALL histwrite_phy(o_gusts, gustiness)
    356358
    357359       IF (vars_defined) THEN
Note: See TracChangeset for help on using the changeset viewer.