Changeset 4722 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90
- Timestamp:
- Oct 9, 2023, 5:33:07 PM (9 months ago)
- Location:
- LMDZ6/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk
- Property svn:mergeinfo changed
/LMDZ6/branches/LMDZ_cdrag_LSCE (added) merged: 4660-4662,4669,4721
- Property svn:mergeinfo changed
-
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4715 r4722 1259 1259 REAL, dimension(klon,klev) :: t_env,q_env 1260 1260 1261 REAL, dimension(klon) :: pr_et 1262 REAL, dimension(klon) :: w_et, jlr_g_c, jlr_g_s 1263 1261 1264 REAL pi 1262 1265 INTEGER ieru … … 2804 2807 ELSE IF (iflag_gusts==2) THEN 2805 2808 gustiness(1:klon)=f_gust_bl*ale_bl_stat(1:klon)+f_gust_wk*ale_wake(1:klon) 2809 !!!! modif olivier torres 2810 ELSE IF (iflag_gusts==3) THEN 2811 w_et=wstar(1,3) 2812 jlr_g_s=(0.65*w_et)**2 2813 pr_et=rain_con*8640 2814 jlr_g_c = (((19.8*(pr_et(1:klon)**2))/(1.5+pr_et(1:klon)+pr_et(1:klon)**2))**(0.4))**2 2815 gustiness(1:klon)=jlr_g_c+jlr_g_s 2816 !! write(*,*) "rain ",pr_et 2817 !! write(*,*) "jlr_g_c",jlr_g_c 2818 !! write(*,*) "wstar",wstar(1,3) 2819 !! write(*,*) "jlr_g_s",jlr_g_s 2806 2820 ! ELSE IF (iflag_gusts==2) THEN 2807 2821 ! do i = 1, klon
Note: See TracChangeset
for help on using the changeset viewer.