Changeset 5445 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90
- Timestamp:
- Dec 21, 2024, 11:03:43 PM (15 hours ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r5425 r5445 80 80 USE lmdz_call_blowing_snow, ONLY : call_blowing_snow_sublim_sedim 81 81 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 82 84 USE yamada_ini_mod, ONLY : yamada_ini 83 85 USE lmdz_atke_turbulence_ini, ONLY : atke_ini … … 1263 1265 CHARACTER(len=512) :: namelist_ecrad_file 1264 1266 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 1265 1276 !======================================================================! 1266 1277 ! Bifurcation vers un nouveau moniteur physique pour experimenter ! … … 1818 1829 CALL iniradia(klon,klev,paprs(1,1:klev+1)) 1819 1830 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 1820 1836 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1821 1837 CALL wake_ini(rg,rd,rv,prt_level) … … 3802 3818 ! 3803 3819 !=================================================================== 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 !=================================================================== 3804 3824 ! Computation of ratqs, the width (normalized) of the subrid scale 3805 3825 ! water distribution
Note: See TracChangeset
for help on using the changeset viewer.