- Timestamp:
- Oct 9, 2023, 10:38:52 PM (15 months ago)
- Location:
- LMDZ6/trunk/libf/phylmdiso
- Files:
-
- 9 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90
r4687 r4723 438 438 REAL, DIMENSION(klon), INTENT(IN) :: rugoro ! rugosity length 439 439 REAL, DIMENSION(klon), INTENT(IN) :: rmu0 ! cosine of solar zenith angle 440 REAL, DIMENSION(klon), INTENT(IN ):: rain_f ! rain fall440 REAL, DIMENSION(klon), INTENT(INOUT) :: rain_f ! rain fall 441 441 REAL, DIMENSION(klon), INTENT(IN) :: snow_f ! snow fall 442 442 REAL, DIMENSION(klon), INTENT(IN) :: bs_f ! blowing snow fall … … 1286 1286 !albedo SB <<< 1287 1287 yrain_f = 0.0 ; ysnow_f = 0.0 ; ybs_f=0.0 ; yfder = 0.0 ; ysolsw = 0.0 1288 ysollw = 0.0 ; yz0m = 0.0 ; yz0h = 0.0 ; y u1 = 0.01288 ysollw = 0.0 ; yz0m = 0.0 ; yz0h = 0.0 ; yz0h_oupas = 0.0 ; yu1 = 0.0 1289 1289 yv1 = 0.0 ; ypaprs = 0.0 ; ypplay = 0.0 ; yqbs1 = 0.0 1290 1290 ydelp = 0.0 ; yu = 0.0 ; yv = 0.0 ; yt = 0.0 … … 1875 1875 ENDDO 1876 1876 CALL cdrag(knon, nsrf, & 1877 speed, yt(:,1), yq(:,1), zgeo1, ypaprs(:,1), &1877 speed, yt(:,1), yq(:,1), zgeo1, ypaprs(:,1), s_pblh, & 1878 1878 yts, yqsurf, yz0m, yz0h, yri0, 0, & 1879 ycdragm, ycdragh, zri1, pref 1879 ycdragm, ycdragh, zri1, pref, rain_f, zxtsol, ypplay(:,1)) 1880 1880 1881 1881 ! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013 … … 1909 1909 1910 1910 CALL cdrag(knon, nsrf, & 1911 speed_x, yt_x(:,1), yq_x(:,1), zgeo1_x, ypaprs(:,1), &1911 speed_x, yt_x(:,1), yq_x(:,1), zgeo1_x, ypaprs(:,1),s_pblh_x,& 1912 1912 yts_x, yqsurf_x, yz0m, yz0h, yri0, 0, & 1913 ycdragm_x, ycdragh_x, zri1_x, pref_x )1913 ycdragm_x, ycdragh_x, zri1_x, pref_x, rain_f, zxtsol, ypplay(:,1) ) 1914 1914 1915 1915 ! --- special Dice. JYG+MPL 25112013 … … 1936 1936 ENDDO 1937 1937 CALL cdrag(knon, nsrf, & 1938 speed_w, yt_w(:,1), yq_w(:,1), zgeo1_w, ypaprs(:,1), &1938 speed_w, yt_w(:,1), yq_w(:,1), zgeo1_w, ypaprs(:,1),s_pblh_w,& 1939 1939 yts_w, yqsurf_w, yz0m, yz0h, yri0, 0, & 1940 ycdragm_w, ycdragh_w, zri1_w, pref_w ) 1941 ! 1942 zgeo1(:) = wake_s(:)*zgeo1_w(:) + (1.-wake_s(:))*zgeo1_x(:) 1940 ycdragm_w, ycdragh_w, zri1_w, pref_w, rain_f, zxtsol, ypplay(:,1) ) 1941 ! 1942 !!!bug !! zgeo1(:) = wake_s(:)*zgeo1_w(:) + (1.-wake_s(:))*zgeo1_x(:) 1943 zgeo1(1:knon) = wake_s(1:knon)*zgeo1_w(1:knon) + (1.-wake_s(1:knon))*zgeo1_x(1:knon) 1943 1944 1944 1945 ! --- special Dice. JYG+MPL 25112013 puis BOMEX … … 2380 2381 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 2381 2382 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 2382 yt2m, yq2m, yt10m, yq10m, yu10m, yustar )2383 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, ypblh, rain_f, zxtsol) 2383 2384 ENDIF 2384 2385 … … 3608 3609 IF (iflag_split .eq.0) THEN 3609 3610 IF (iflag_new_t2mq2m==1) THEN 3610 CALL stdlevvarn(klon, knon, nsrf, zxli, &3611 CALL stdlevvarn(klon, knon, nsrf, zxli, & 3611 3612 uzon, vmer, tair1, qair1, zgeo1, & 3612 3613 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & … … 3617 3618 uzon, vmer, tair1, qair1, zgeo1, & 3618 3619 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 3619 yt2m, yq2m, yt10m, yq10m, yu10m, yustar )3620 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, ypblh, rain_f, zxtsol) 3620 3621 ENDIF 3621 3622 ELSE !(iflag_split .eq.0) … … 3635 3636 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & 3636 3637 tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, & 3637 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x )3638 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, ypblh_x, rain_f, zxtsol) 3638 3639 CALL stdlevvar(klon, knon, nsrf, zxli, & 3639 3640 uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, & 3640 3641 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 3641 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w )3642 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, ypblh_w, rain_f, zxtsol) 3642 3643 ENDIF 3643 3644 !!! -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4715 r4723 1333 1333 #endif 1334 1334 1335 1336 REAL, dimension(klon,klev) :: t_env,q_env 1337 1338 REAL, dimension(klon) :: pr_et 1339 REAL, dimension(klon) :: w_et, jlr_g_c, jlr_g_s 1340 1335 1341 REAL pi 1336 1342 INTEGER ieru … … 3101 3107 ELSE IF (iflag_gusts==2) THEN 3102 3108 gustiness(1:klon)=f_gust_bl*ale_bl_stat(1:klon)+f_gust_wk*ale_wake(1:klon) 3109 !!!! modif olivier torres 3110 ELSE IF (iflag_gusts==3) THEN 3111 w_et=wstar(1,3) 3112 jlr_g_s=(0.65*w_et)**2 3113 pr_et=rain_con*8640 3114 jlr_g_c = (((19.8*(pr_et(1:klon)**2))/(1.5+pr_et(1:klon)+pr_et(1:klon)**2))**(0.4))**2 3115 gustiness(1:klon)=jlr_g_c+jlr_g_s 3116 !! write(*,*) "rain ",pr_et 3117 !! write(*,*) "jlr_g_c",jlr_g_c 3118 !! write(*,*) "wstar",wstar(1,3) 3119 !! write(*,*) "jlr_g_s",jlr_g_s 3103 3120 ! ELSE IF (iflag_gusts==2) THEN 3104 3121 ! do i = 1, klon
Note: See TracChangeset
for help on using the changeset viewer.