Changeset 5776 for LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

Timestamp:

Jul 15, 2025, 2:13:19 PM (2 weeks ago)

Author:

evignon

Message:

ajout de l'advection horizontale de la TKE. Travaux dans le cadre de la these de Valentin Wiener

File:

• LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

@@ -26 +26 @@
        falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien, &
        ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, &
-       cf_ancien, rvc_ancien, radpas, radsol, rain_fall, ratqs, &
+       cf_ancien, rvc_ancien, tke_ancien, radpas, radsol, rain_fall, ratqs, &
        rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &
        solsw, solswfdiff, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
@@ -518 +518 @@
   ENDIF
 
+
   found=phyetat0_get(clwcon,"CLWCON","CLWCON",0.)
   found=phyetat0_get(rnebcon,"RNEBCON","RNEBCON",0.)
@@ -528 +529 @@
 !==================================
 !
-  IF (iflag_pbl>1) then
+  ! cas specifique de l'advection de TKE
+  IF (ok_advtke) THEN
+       ancien_ok=ancien_ok.AND.phyetat0_get(tke_ancien,"TKEANCIEN","TKEANCIEN",0.)
+  ELSE
+    tke_ancien(:,:)=0.
+  ENDIF
+
+  IF (ok_advtke) THEN
+    IF ( (maxval(tke_ancien).EQ.minval(tke_ancien))) THEN
+       ancien_ok=.false.
+    ENDIF
+  ENDIF
+
+  IF ((iflag_pbl>1)) then
      found=phyetat0_srf(pbl_tke,"TKE","Turb. Kinetic. Energ. ",1.e-8)
   ENDIF