Changeset 4052 for LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

Timestamp:

Dec 26, 2021, 10:27:36 PM (2 years ago)

Author:

dcugnet

Message:

Fixes for previous commit:

• few tracers(:)%nqChilds were errouneously replaced with tarcers(:)%nqDesc
• minor changes

File:

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

LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

@@ -174 +174 @@
    !           Rayleigh distillation las.
     iName = tracers(iq)%iso_iName
+    IF(.NOT.ok_isotopes .OR. iName <= 0) CYCLE
     iZone = tracers(iq)%iso_iZone
     iPhase= tracers(iq)%iso_iPhase
     iqParent = tracers(iq)%iqParent
-    IF(ok_isotopes .AND. iName>0) THEN
-      IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName)    &
-     &           *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.)
-      IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase))
-    END IF
+    IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName)    &
+     &         *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.)
+    IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase))
   END DO
   DEALLOCATE(q_glo)