Changeset 5777

Timestamp:

Jul 15, 2025, 11:04:32 PM (36 hours ago)

Author:

evignon

Message:

encapsulage d'une ligne propre a l'advection de tke

File:

• LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

@@ -252 +252 @@
 !==============================================================================================================================
    nqtrue = SIZE(tracers)                                                                               !--- "true" tracers
-   nqo    =      COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name)     == 'H2O')     !--- Water phases
+   nqo    =      COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name)     == 'H2O')     !-- Water phases
    nqtke  =      COUNT(tracers(:)%component == 'lmdz' .AND. tracers(:)%name               == 'TKE')     !--- TKE tracers
    nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O')     !--- Passed to phytrac
@@ -370 +370 @@
    !--- Compute indices for TKE when it is advected
    itke = strIdx(tracers(:)%name, 'TKE')
-   !--- For TKE, we force isInPhysics=.False.
-   tracers(itke)%isInPhysics = .FALSE.
-
-
+   IF (nqtke .GE. 1) THEN
+       !--- For TKE, we force isInPhysics=.False. 
+       tracers(itke)%isInPhysics = .FALSE.
+   ENDIF
 
    IF(CPPKEY_STRATAER .AND. type_trac == 'coag') THEN