Changeset 4661 for LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/physiq_mod.F90

Timestamp:

Aug 31, 2023, 2:58:34 PM (15 months ago)

Author:

Laurent Fairhead

Message:

Inclusion and merge of Olivier's modifications

File:

• LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/physiq_mod.F90

@@ -1212 +1212 @@
     REAL qql1(klon),qql2(klon),corrqql
 
+
+    REAL, dimension(klon) :: pr_et
+    REAL, dimension(klon) :: w_et, jlr_g_c, jlr_g_s
+
     REAL pi
 
@@ -2667 +2671 @@
        ELSE IF (iflag_gusts==2) THEN
           gustiness(1:klon)=f_gust_bl*ale_bl_stat(1:klon)+f_gust_wk*ale_wake(1:klon)
+       !!!! modif olivier torres
+       ELSE IF (iflag_gusts==3) THEN
+          w_et=wstar(1,3)
+          jlr_g_s=(0.65*w_et)**2
+          pr_et=rain_con*8640
+          jlr_g_c = (((19.8*(pr_et(1:klon)**2))/(1.5+pr_et(1:klon)+pr_et(1:klon)**2))**(0.4))**2
+          gustiness(1:klon)=jlr_g_c+jlr_g_s
+!!       write(*,*) "rain ",pr_et
+!!       write(*,*) "jlr_g_c",jlr_g_c
+!!       write(*,*) "wstar",wstar(1,3)
+!!       write(*,*) "jlr_g_s",jlr_g_s
           ! ELSE IF (iflag_gusts==2) THEN
           !    do i = 1, klon