Changeset 2286 for LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F

Timestamp:

May 20, 2015, 3:27:07 PM (10 years ago)

Author:

Ehouarn Millour

Message:

• fix in dyn3d the array out of bound issue (in qminimum) that was only corrected in dyn3dmem in rev 2285.
• comment out many invasive debug writes.

EM

File:

• LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F (modified) (12 diffs)

LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F

@@ -56 +56 @@
       INTEGER ijb,ije,ijb_x,ije_x
       
-      write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
-     &   iq,ijb_x
+      !write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
+!     &   iq,ijb_x
 c   calcul de la pente a droite et a gauche de la maille
 
@@ -154 +154 @@
       ENDIF ! (pente_max.lt.-1.e-5)
 
-      write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
 
 c   bouclage de la pente en iip1:
@@ -273 +273 @@
                ENDDO
                niju=iju
-               PRINT*,'vlx 278, niju,nl',niju,nl(l)
+               !PRINT*,'vlx 278, niju,nl',niju,nl(l)
 
 c  traitement des mailles
@@ -331 +331 @@
 ! Il faut faire ça avant d'avoir mis à jour q et masse
 
-      write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq) 
+      !write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq) 
 
       if (nqfils(iq).gt.0) then  
@@ -354 +354 @@
 ! end CRisi
 
-      write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x
 
 c   calcul des tENDances
@@ -373 +373 @@
       ENDDO
 c$OMP END DO NOWAIT
-      write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -394 +394 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x
 c     CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1)
 c     CALL SCOPY((jjm-1)*llm,masse(iip1+iip1,1),iip1,masse(iip2,1),iip1)
@@ -726 +726 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
 
       ijb=ij_begin-2*iip1
@@ -926 +926 @@
 c    sens de W
 
-      write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
+      !write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
 #ifdef BIDON
       IF(testcpu) THEN
@@ -982 +982 @@
 c calcul de  - d( q   * w )/ d(sigma)    qu'on ajoute a  dq pour calculer dq
 
-       write(*,*) 'vlz 982,ijb,ije=',ijb,ije
+       !write(*,*) 'vlz 982,ijb,ije=',ijb,ije
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
        DO l = 1,llm-1
@@ -998 +998 @@
        ENDDO
 c$OMP END DO NOWAIT    
-       write(*,*) 'vlz 1001'    
+       !write(*,*) 'vlz 1001'    
 
 c$OMP MASTER
@@ -1010 +1010 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
+      !write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
       if (nqfils(iq).gt.0) then  
        do ifils=1,nqdesc(iq)