Changeset 5445

Timestamp:

Dec 21, 2024, 11:03:43 PM (9 hours ago)

Author:

fhourdin

Message:

Adding wind subrid scale distribution

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_surf_wind.f90 (added)
• lmdz_surf_wind_ini.f90 (added)
• physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -80 +80 @@
     USE lmdz_call_blowing_snow, ONLY : call_blowing_snow_sublim_sedim
     USE lmdz_wake_ini, ONLY : wake_ini
+    USE lmdz_surf_wind_ini, ONLY : surf_wind_ini, iflag_surf_wind
+    USE lmdz_surf_wind, ONLY : surf_wind
     USE yamada_ini_mod, ONLY : yamada_ini
     USE lmdz_atke_turbulence_ini, ONLY : atke_ini
@@ -1263 +1265 @@
     CHARACTER(len=512) :: namelist_ecrad_file
 
+
+    ! Subgrid scale wind :
+    ! Need to be allocatable/save because the number of bin is not known (provided by surf_wind_ini)
+    integer, save :: nsrfwnd=1
+    real, dimension(:,:), allocatable, save :: surf_wind_value, surf_wind_proba ! module and probability of sugrdi wind wind sample
+    !$OMP THREADPRIVATE(nsrfwnd,surf_wind_value, surf_wind_proba)
+    
+
+
     !======================================================================!
     ! Bifurcation vers un nouveau moniteur physique pour experimenter      !
@@ -1818 +1829 @@
        CALL iniradia(klon,klev,paprs(1,1:klev+1))
 
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+       CALL surf_wind_ini(klon,lunout)
+       CALL getin_p('nsrfwnd',nsrfwnd)
+       allocate(surf_wind_value(klon,nsrfwnd),surf_wind_proba(klon,nsrfwnd))
+    
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        CALL wake_ini(rg,rd,rv,prt_level)
@@ -3802 +3818 @@
     !
     !===================================================================
+    ! Computation of subrgid scale near-surface wind distribution
+    call surf_wind(klon,nsrfwnd,u10m,v10m,wake_s,wake_Cstar,ustar,wstar,surf_wind_value,surf_wind_proba)
+
+    !===================================================================
     ! Computation of ratqs, the width (normalized) of the subrid scale
     ! water distribution