Changeset 4052 for LMDZ6/trunk/libf/dyn3d/dynetat0.f90

Timestamp:

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

Author:

dcugnet

Message:

Fixes for previous commit:

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

File:

• LMDZ6/trunk/libf/dyn3d/dynetat0.f90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3d/dynetat0.f90

@@ -137 +137 @@
    !           Rayleigh distillation law.
     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