Changeset 4036 for LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90

Timestamp:

Dec 1, 2021, 8:23:48 AM (3 years ago)

Author:

crisi

Message:

correction d'un probleme d'initialisation phylmdiso + suite du nettoyage dans phylmdiso

File:

• LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90

@@ -16063 +16063 @@
 !      write(*,*) 'xtsnow(4,8,1)=',xtsnow(4,8,1)
 #ifdef ISOVERIF
-      write(*,*) 'phyiso_etat0 15993 tmp: xtsol(iso_eau,1),qsol(1)=',xtsol(iso_eau,1),qsol(1)
       do i=1,klon
          do ixt=1,niso
@@ -16368 +16367 @@
      &                  *(1-kcin(ixt))/(1.0-kcin(ixt)*h0)
            xt_ancien(ixt,i,k)=q_ancien(i,k)*RMerlivat(ixt) &
-     &                  *(q_ancien(i,k)/q0)**(alpha(ixt)-1.0)
+     &                  *(min(q0,q_ancien(i,k))/q0)**(alpha(ixt)-1.0)
            if (q_ancien(i,k).gt.ridicule) then
            xtl_ancien(ixt,i,k)=ql_ancien(i,k)*alpha(ixt) &
@@ -16468 +16467 @@
 
 #ifdef ISOVERIF
+      write(*,*) 'phyisoetat0 16468: verif init dure'
       do i=1,klon
          do ixt=1,niso
@@ -16499 +16499 @@
            if (iso_eau.gt.0) then
              call iso_verif_egalite(xt_ancien(iso_eau,i,k), &
-     &           q_ancien(i,k),'phyiso_etat0_dur 775')        
+     &           q_ancien(i,k),'phyiso_etat0_dur 775a')        
            endif !if (iso_eau.gt.0) then
+           if (iso_HDO.gt.0) then
+             if (q_ancien(i,k).gt.ridicule) then
+              call iso_verif_aberrant_encadre( &
+     &           xt_ancien(iso_hdo,i,k)/q_ancien(i,k), &
+     &          'phyiso_etat0_dur 775b')
+             endif !if (q_ancien(i,k).gt.ridicule) then
+           endif !if (iso_HDO.gt.0) then
+           if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
+              if (q_ancien(i,k).gt.ridicule) then
+                call iso_verif_O18_aberrant( &
+     &              xt_ancien(iso_hdo,i,k)/q_ancien(i,k), &
+     &              xt_ancien(iso_O18,i,k)/q_ancien(i,k), &
+     &              'phyiso_etat0_dur 775c')
+              endif ! if (q_ancien(i,k).gt.ridicule) then
+           endif ! if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then 
          enddo !do k=1,klev
          do nsrf=1,nbsrf