Changeset 4037 for LMDZ6/trunk/libf/dyn3d

Timestamp:

Dec 1, 2021, 8:27:00 AM (3 years ago)

Author:

crisi

Message:

correction dans check_isotopes

File:

• LMDZ6/trunk/libf/dyn3d/check_isotopes.F (modified) (7 diffs)

LMDZ6/trunk/libf/dyn3d/check_isotopes.F

@@ -23 +23 @@
         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
@@ -43 +43 @@
                   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 = 1,ip1jmp1
@@ -67 +67 @@
                   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:
@@ -93 +93 @@
                   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
@@ -118 +118 @@
                   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
@@ -149 +149 @@
                   write(*,*) 'q,iq,i,k,=',q(i,k,iq),iq,i,k
                   write(*,*) 'deltaD=',deltaD
-                  stop
+                  call abort_gcm('check_isotopes_loc','plantage iso',0)
                  endif !if ((deltaD.gt.deltaDmax).or.
                 endif !if (q(i,k,iq).gt.qmin) then
@@ -178 +178 @@
                   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.