Changeset 4490

Timestamp:

Apr 3, 2023, 2:56:52 AM (14 months ago)

Author:

crisi

Message:

correction initialization water tagging in dyn3dmem

File:

• LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

@@ -183 +183 @@
          CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized with a simplified Rayleigh distillation law.', modname)
          q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat*(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal-1.)
+         ! Camille 9 mars 2023: point de vigilence: initialisation incohérente
+         ! avec celle de xt_ancien dans la physiq.
       ELSE
          CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized to its parent isotope concentration.', modname)
-         q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqIsoPha(iName,iPhase))
+         ! Camille 9 mars 2023: attention!! seuls les tags qui correspondent à
+         ! izone=izone_init (définie dans isotrac_mod) sont initialisés comme
+         ! les parents. Sinon, c'est nul.
+         ! j'ai fait ça en attendant, mais il faudrait initialiser proprement en
+         ! remplacant 1 par izone_init dans la ligne qui suit.
+         IF(tracers(iq)%iso_iZone == 1) THEN
+          q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqIsoPha(iName,iPhase))
+         ELSE
+           q(ijb_u:ije_u,:,iq) =  0.
+         ENDIF
       END IF
     !--------------------------------------------------------------------------------------------------------------------------