Changeset 4071 for LMDZ6/trunk/libf/phylmd/phys_output_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/phylmd/phys_output_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_mod.F90

@@ -470 +470 @@
           itr = 0
           DO iq = 1, nqtot 
-            IF(.NOT.tracers(iq)%isAdvected .OR. tracers(iq)%isH2Ofamily) CYCLE
+            IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE
             itr = itr + 1
             dn = 'd'//TRIM(tracers(iq)%name)//'_'
@@ -526 +526 @@
 
 !  DO iq=1,nqtot
-!    IF(.NOT.tracers(iq)%isAdvected .OR. tracers(iq)%isH2Ofamily) CYCLE
+!    IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE
 !    WRITE(*,'(a,i1,a,10i3)')'trac(',iq,')%flag = ',o_trac(iq)%flag
 !    WRITE(*,'(a,i1,a)')'trac(',iq,')%name = '//TRIM(o_trac(iq)%name)