Changeset 4038 for LMDZ6/trunk/libf

Timestamp:

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

Author:

crisi

Message:

correction dans check_isotopes + suppression d'un call check_isotopes problematique dans advtrac

Location:

LMDZ6/trunk/libf/dyn3dmem

Files:

• advtrac_loc.F (modified) (1 diff)
• check_isotopes_loc.F (modified) (13 diffs)

LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F

@@ -156 +156 @@
      *                        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

LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F

@@ -24 +24 @@
         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 +30 @@
         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 +45 @@
                   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 +53 @@
         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 +71 @@
                   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 +81 @@
         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 +99 @@
                   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 +108 @@
         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 +126 @@
                   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 +136 @@
 
 
-        !write(*,*) 'check_isotopes 129'
+!        write(*,*) 'check_isotopes 129'
         if (ok_isotrac) then
 
@@ -158 +158 @@
                   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 +190 @@
                   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 +209 @@
 
         endif ! if (ok_isotopes)
-        !write(*,*) 'check_isotopes 198'
+!        write(*,*) 'check_isotopes 198'
         
         end