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/ocean_cpl_mod.F90

    r1907 r2240  
    4949       AcoefH, AcoefQ, BcoefH, BcoefQ, &
    5050       AcoefU, AcoefV, BcoefU, BcoefV, &
    51        ps, u1, v1, &
     51       ps, u1, v1, gustiness, &
    5252       radsol, snow, agesno, &
    5353       qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
    … …  
    8383    REAL, DIMENSION(klon), INTENT(IN)        :: AcoefU, AcoefV, BcoefU, BcoefV
    8484    REAL, DIMENSION(klon), INTENT(IN)        :: ps
    85     REAL, DIMENSION(klon), INTENT(IN)        :: u1, v1
     85    REAL, DIMENSION(klon), INTENT(IN)        :: u1, v1, gustiness
    8686
    8787! In/Output arguments
    … …  
    138138         tsurf_cpl, p1lay, cal, beta, cdragh, ps, &
    139139         precip_rain, precip_snow, snow, qsurf,  &
    140          radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
     140         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
    141141         AcoefH, AcoefQ, BcoefH, BcoefQ, &
    142142         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
    … …  
    144144! - Flux calculation at first modele level for U and V
    145145    CALL calcul_flux_wind(knon, dtime, &
    146          u0_cpl, v0_cpl, u1, v1, cdragm, &
     146         u0_cpl, v0_cpl, u1, v1, gustiness, cdragm, &
    147147         AcoefU, AcoefV, BcoefU, BcoefV, &
    148148         p1lay, temp_air, &
    … …  
    185185       AcoefH, AcoefQ, BcoefH, BcoefQ, &
    186186       AcoefU, AcoefV, BcoefU, BcoefV, &
    187        ps, u1, v1, pctsrf, &
     187       ps, u1, v1, gustiness, pctsrf, &
    188188       radsol, snow, qsurf, &
    189189       alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
    … …  
    219219    REAL, DIMENSION(klon), INTENT(IN)        :: AcoefU, AcoefV, BcoefU, BcoefV
    220220    REAL, DIMENSION(klon), INTENT(IN)        :: ps
    221     REAL, DIMENSION(klon), INTENT(IN)        :: u1, v1
     221    REAL, DIMENSION(klon), INTENT(IN)        :: u1, v1, gustiness
    222222    REAL, DIMENSION(klon,nbsrf), INTENT(IN)  :: pctsrf
    223223
    … …  
    281281         tsurf_cpl, p1lay, cal, beta, cdragh, ps, &
    282282         precip_rain, precip_snow, snow, qsurf,  &
    283          radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
     283         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
    284284         AcoefH, AcoefQ, BcoefH, BcoefQ, &
    285285         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
    … …  
    288288! - Flux calculation at first modele level for U and V
    289289    CALL calcul_flux_wind(knon, dtime, &
    290          u0, v0, u1, v1, cdragm, &
     290         u0, v0, u1, v1, gustiness, cdragm, &
    291291         AcoefU, AcoefV, BcoefU, BcoefV, &
    292292         p1lay, temp_air, &
Note: See TracChangeset for help on using the changeset viewer.