Ignore:
Timestamp:
Dec 21, 2024, 11:03:43 PM (15 hours ago)
Author:
fhourdin
Message:

Adding wind subrid scale distribution

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r5425 r5445  
    8080    USE lmdz_call_blowing_snow, ONLY : call_blowing_snow_sublim_sedim
    8181    USE lmdz_wake_ini, ONLY : wake_ini
     82    USE lmdz_surf_wind_ini, ONLY : surf_wind_ini, iflag_surf_wind
     83    USE lmdz_surf_wind, ONLY : surf_wind
    8284    USE yamada_ini_mod, ONLY : yamada_ini
    8385    USE lmdz_atke_turbulence_ini, ONLY : atke_ini
     
    12631265    CHARACTER(len=512) :: namelist_ecrad_file
    12641266
     1267
     1268    ! Subgrid scale wind :
     1269    ! Need to be allocatable/save because the number of bin is not known (provided by surf_wind_ini)
     1270    integer, save :: nsrfwnd=1
     1271    real, dimension(:,:), allocatable, save :: surf_wind_value, surf_wind_proba ! module and probability of sugrdi wind wind sample
     1272    !$OMP THREADPRIVATE(nsrfwnd,surf_wind_value, surf_wind_proba)
     1273   
     1274
     1275
    12651276    !======================================================================!
    12661277    ! Bifurcation vers un nouveau moniteur physique pour experimenter      !
     
    18181829       CALL iniradia(klon,klev,paprs(1,1:klev+1))
    18191830
     1831!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     1832       CALL surf_wind_ini(klon,lunout)
     1833       CALL getin_p('nsrfwnd',nsrfwnd)
     1834       allocate(surf_wind_value(klon,nsrfwnd),surf_wind_proba(klon,nsrfwnd))
     1835   
    18201836!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    18211837       CALL wake_ini(rg,rd,rv,prt_level)
     
    38023818    !
    38033819    !===================================================================
     3820    ! Computation of subrgid scale near-surface wind distribution
     3821    call surf_wind(klon,nsrfwnd,u10m,v10m,wake_s,wake_Cstar,ustar,wstar,surf_wind_value,surf_wind_proba)
     3822
     3823    !===================================================================
    38043824    ! Computation of ratqs, the width (normalized) of the subrid scale
    38053825    ! water distribution
Note: See TracChangeset for help on using the changeset viewer.