Index: LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F
===================================================================
--- LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F	(revision 3995)
+++ LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F	(revision 4038)
@@ -156,9 +156,4 @@
      *                        pk,teta )
 
-          !write(*,*) 'advtrac 162: apres appel vlspltgen_loc'
-      if (ok_iso_verif) then
-           call check_isotopes(q,ijb_u,ije_u,'advtrac 162')
-      endif !if (ok_iso_verif) then
-
 #ifdef DEBUG_IO     
           do j=1,nqtot
Index: LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F
===================================================================
--- LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F	(revision 3995)
+++ LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F	(revision 4038)
@@ -24,5 +24,5 @@
         parameter (errmaxrel=1e-3)
         parameter (qmin=1e-11)
-        parameter (deltaDmax=200.0,deltaDmin=-999.9)
+        parameter (deltaDmax=1000.0,deltaDmin=-999.0)
         parameter (ridicule=1e-12)
         real deltaD
@@ -30,5 +30,5 @@
         if (ok_isotopes) then
 
-        !write(*,*) 'check_isotopes 31: err_msg=',err_msg
+!        write(*,*) 'check_isotopes 31: err_msg=',err_msg
         ! verifier que rien n'est NaN
         do ixt=1,ntraciso
@@ -45,5 +45,5 @@
                   write(*,*) 'q,i,k,iq=',q(i,k,iq),i,k,iq
                   write(*,*) 'borne=',borne
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
                 endif  !if ((x(ixt,i,j).gt.-borne).and.
               enddo !DO i = ijb,ije
@@ -53,5 +53,5 @@
         enddo !do ixt=1,ntraciso
 
-        !write(*,*) 'check_isotopes 52'
+!        write(*,*) 'check_isotopes 52'
         ! verifier que l'eau normale est OK
         if (use_iso(1)) then
@@ -71,5 +71,5 @@
                   write(*,*) 'q,iq,i,k=',q(i,k,iq),iq,i,k
                   write(*,*) 'q(i,k,phase)=',q(i,k,phase)
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
               endif !if ((abs(q(i,k,phase)-q(i,k,iq)).gt.errmax).and.
               ! bidouille pour éviter divergence:
@@ -81,5 +81,5 @@
         endif !if (use_iso(1)) then
         
-        !write(*,*) 'check_isotopes 78'
+!        write(*,*) 'check_isotopes 78'
         ! verifier que HDO est raisonable
         if (use_iso(2)) then
@@ -99,5 +99,5 @@
                   write(*,*) 'q=',q(i,k,:)
                   write(*,*) 'deltaD=',deltaD
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
              endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
             endif !if (q(i,k,iq).gt.qmin) then
@@ -108,5 +108,5 @@
         endif !if (use_iso(2)) then
 
-        !write(*,*) 'check_isotopes 103'
+!        write(*,*) 'check_isotopes 103'
         ! verifier que O18 est raisonable
         if (use_iso(3)) then
@@ -126,5 +126,5 @@
                   write(*,*) 'xt=',q(i,k,:)
                   write(*,*) 'deltaO18=',deltaD
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
              endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
             endif !if (q(i,k,iq).gt.qmin) then
@@ -136,5 +136,5 @@
 
 
-        !write(*,*) 'check_isotopes 129'
+!        write(*,*) 'check_isotopes 129'
         if (ok_isotrac) then
 
@@ -158,6 +158,6 @@
                   write(*,*) 'ixt,ieau=',ixt,ieau
                   write(*,*) 'q,iq,i,k,=',q(i,k,iq),iq,i,k
-                  write(*,*) 'deltaD=',deltaD
-                  stop
+                  write(*,*) 'deltaD=',deltaD                  
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
                  endif !if ((deltaD.gt.deltaDmax).or.
                 endif !if (q(i,k,iq).gt.qmin) then
@@ -190,5 +190,5 @@
                   write(*,*) 'i,k,=',i,k
                   write(*,*) 'q(i,k,:)=',q(i,k,:)
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
                  endif !if ((abs(q(i,k,phase)-q(i,k,iq)).gt.errmax).and.
                   
@@ -209,5 +209,5 @@
 
         endif ! if (ok_isotopes)
-        !write(*,*) 'check_isotopes 198'
+!        write(*,*) 'check_isotopes 198'
         
         end