Changeset 4071 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Jan 31, 2022, 9:20:17 PM (2 years ago)

Author:

dcugnet

Message:

• Fix for unadvected tracers (iadv==0)
• The key %isH2Ofamily, from the derived type "trac_type", is replaced with the more general

key %isInPhysics, which is TRUE for tracers both in "qx" and "tr_seri".

Currently, FALSE for tracers descending on H2O (isotopes and tagging tracers included). Could be set to FALSE
for interactive CO2 (type_trac=='inco') or ice supersaturated cloud content (tranfered to "rneb_seri")

• Parser routines taken from version 10 of ​https://svn.lmd.jussieu.fr/tracers-parser

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -2474 +2474 @@
        itr = 0
        DO iq = 1, nqtot
-         IF(tracers(iq)%isH2Ofamily) CYCLE
+         IF(.NOT.tracers(iq)%isInPhysics) CYCLE
          itr = itr+1
 !#ifdef ISOVERIF
@@ -2498 +2498 @@
        itr = 0
        do iq = 1, nqtot
-         IF(tracers(iq)%isH2Ofamily) CYCLE
+         IF(.NOT.tracers(iq)%isInPhysics) CYCLE
          itr = itr+1
          tr_ancien(:,:,itr)=tr_seri(:,:,itr)       
@@ -6531 +6531 @@
     itr = 0
     DO iq = 1, nqtot
-       IF(tracers(iq)%isH2Ofamily) CYCLE
+       IF(.NOT.tracers(iq)%isInPhysics) CYCLE
        itr = itr+1
        DO  k = 1, klev