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

Timestamp:

Mar 5, 2023, 11:18:00 AM (16 months ago)

Author:

crisi

Message:

quelques corrections dans les verifs isotopiques

File:

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

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -2682 +2682 @@
              DO i = 1, klon
               if (q_seri(i,k).gt.3e-3) then
-              call iso_verif_positif(deltaD(xt_seri(iso_eau,i,k) &
+              call iso_verif_positif(deltaD(xt_seri(iso_HDO,i,k) &
      &           /q_seri(i,k))+400.0,'physiq 2045a')
-              call iso_verif_positif(deltaD(xt_ancien(iso_eau,i,k) &
+              call iso_verif_positif(deltaD(xt_ancien(iso_HDO,i,k) &
      &            /q_ancien(i,k))+400.0,'physiq 2045b')
 !              call iso_verif_egalite_choix(d_xt_dyn(iso_hdo,i,k) &
@@ -3602 +3602 @@
            endif ! if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
            if ((iso_O17.gt.0).and.(iso_O18.gt.0)) then
-             if ((q_seri(i,k).gt.ridicule).and.(l.lt.nlevmaxO17)) then
+             if ((q_seri(i,k).gt.ridicule).and.(k.lt.nlevmaxO17)) then
               call iso_verif_aberrant_o17(xt_seri(iso_o17,i,k) &
      &           /q_seri(i,k),xt_seri(iso_o18,i,k) &