Changeset 3331 for LMDZ5/branches/IPSLCM5A2.1_ISO/libf/dyn3dmem/check_isotopes_loc.F

Timestamp:

May 22, 2018, 4:30:21 PM (7 years ago)

Author:

acozic

Message:

Add modification for isotopes

File:

• LMDZ5/branches/IPSLCM5A2.1_ISO/libf/dyn3dmem/check_isotopes_loc.F (modified) (3 diffs)

LMDZ5/branches/IPSLCM5A2.1_ISO/libf/dyn3dmem/check_isotopes_loc.F

@@ -12 +12 @@
 
         ! locals
-        integer ixt,phase,k,i,iq,iiso,izone,ieau,iqeau
+        integer ixt,phase,k,i,iq,iiso,izone,ieau,iqeau,ixt2,iq2
         real xtractot,xiiso
         real borne
@@ -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
+        real dexcessmax,dexcessmin
+        parameter (dexcessmax=6000.0,dexcessmin=-100.0)
 
         if (ok_isotopes) then
@@ -134 +136 @@
           enddo ! do phase=1,nqo 
         endif !if (use_iso(2)) then
+
+        !write(*,*) 'check_isotopes 103'
+        ! verifier que dexcess est raisonable
+        if (use_iso(3).and.use_iso(2)) then
+          ixt=indnum_fn_num(3)
+          ixt2=indnum_fn_num(2)
+          do phase=1,nqo
+            iq=iqiso(ixt,phase)
+            iq2=iqiso(ixt2,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+            do k=1,llm
+            DO i = ijb,ije
+            if (q(i,k,iq).gt.qmin) then
+             deltaD=(q(i,k,iq2)/q(i,k,phase)/tnat(2)-1)*1000
+     &            -8.0*(q(i,k,iq)/q(i,k,phase)/tnat(3)-1)*1000
+             if ((deltaD.gt.dexcessmax).or.(deltaD.lt.dexcessmin)) then
+                  write(*,*) 'erreur detectee iso_verif_aberrant O18:'
+                  write(*,*) err_msg
+                  write(*,*) 'ixt,ixt2,phase=',ixt,ixt2,phase
+                  write(*,*) 'q,iq,i,k,=',q(i,k,phase),iq,i,k
+                  write(*,*) 'xt=',q(i,k,:)
+                  write(*,*) 'dexcess=',deltaD
+                  stop
+             endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
+            endif !if (q(i,k,iq).gt.qmin) then
+            enddo !DO i = ijb,ije
+            enddo !do k=1,llm
+c$OMP END DO NOWAIT
+          enddo ! do phase=1,nqo 
+        endif !if (use_iso(2)) then
+