Changeset 2281 for LMDZ5/trunk/libf/dyn3dmem/check_isotopes_loc.F

Timestamp:

May 18, 2015, 11:09:37 AM (9 years ago)

Author:

crisi

Message:

Camille Risi: corrections of bugs for the isotopic part

File:

• LMDZ5/trunk/libf/dyn3dmem/check_isotopes_loc.F (modified) (14 diffs)

LMDZ5/trunk/libf/dyn3dmem/check_isotopes_loc.F

@@ -30 +30 @@
         if (ok_isotopes) then
 
+        !write(*,*) 'check_isotopes 31: err_msg=',err_msg
         ! verifier que rien n'est NaN
         do ixt=1,ntraciso
           do phase=1,nqo
             iq=iqiso(ixt,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
             do k=1,llm
               DO i = ijb,ije
@@ -47 +49 @@
               enddo !DO i = ijb,ije
             enddo !do k=1,llm
+c$OMP END DO NOWAIT
           enddo !do phase=1,nqo
         enddo !do ixt=1,ntraciso
 
+        !write(*,*) 'check_isotopes 52'
         ! verifier que l'eau normale est OK
         if (use_iso(1)) then
@@ -55 +59 @@
           do phase=1,nqo
             iq=iqiso(ixt,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
             do k=1,llm
             DO i = ijb,ije  
@@ -63 +68 @@
                   write(*,*) 'erreur detectee par iso_verif_egalite:'
                   write(*,*) err_msg
-                  write(*,*) 'ixt,phase=',ixt,phase
+                  write(*,*) 'ixt,phase,ijb=',ixt,phase,ijb
                   write(*,*) 'q,iq,i,k=',q(i,k,iq),iq,i,k
                   write(*,*) 'q(i,k,phase)=',q(i,k,phase)
@@ -72 +77 @@
             enddo ! DO i = ijb,ije
             enddo !do k=1,llm
+c$OMP END DO NOWAIT
           enddo ! do phase=1,nqo 
         endif !if (use_iso(1)) then
         
+        !write(*,*) 'check_isotopes 78'
         ! verifier que HDO est raisonable
         if (use_iso(2)) then
@@ -80 +87 @@
           do phase=1,nqo
             iq=iqiso(ixt,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
             do k=1,llm
             DO i = ijb,ije
@@ -96 +104 @@
             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
 
+        !write(*,*) 'check_isotopes 103'
         ! verifier que O18 est raisonable
         if (use_iso(3)) then
@@ -104 +114 @@
           do phase=1,nqo
             iq=iqiso(ixt,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
             do k=1,llm
             DO i = ijb,ije
@@ -120 +131 @@
             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
 
 
+        !write(*,*) 'check_isotopes 129'
         if (ok_isotrac) then
 
@@ -133 +146 @@
                iq=iqiso(ixt,phase)
                iqeau=iqiso(ieau,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
                do k=1,llm
                 DO i = ijb,ije
@@ -150 +164 @@
                 enddo !DO i = ijb,ije
                 enddo  ! do k=1,llm
+c$OMP END DO NOWAIT
               enddo ! do phase=1,nqo    
             enddo !do izone=1,ntraceurs_zone
@@ -157 +172 @@
            do phase=1,nqo
               iq=iqiso(iiso,phase)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
               do k=1,llm
                 DO i = ijb,ije
@@ -186 +202 @@
                 enddo !DO i = ijb,ije
               enddo !do k=1,llm
+c$OMP END DO NOWAIT
            enddo !do phase=1,nqo
           enddo !do iiso=1,niso
@@ -192 +209 @@
 
         endif ! if (ok_isotopes)
+        !write(*,*) 'check_isotopes 198'
         
         end