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

Timestamp:

May 20, 2015, 3:27:07 PM (9 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/vlspltqs_loc.F (modified) (9 diffs)

LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F

@@ -47 +47 @@
       INTEGER ijb,ije,ijb_x,ije_x
       
-      write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
-     &   iq,ijb_x
+      !write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
+!     &   iq,ijb_x
 
 c   calcul de la pente a droite et a gauche de la maille
@@ -280 +280 @@
                ENDDO
                niju=iju
-               PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
+               !PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
 
 c  traitement des mailles
@@ -336 +336 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=',
-     &     iq,ijb_x,nqfils(iq)  
+      !write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=',
+!     &     iq,ijb_x,nqfils(iq)  
 
       if (nqfils(iq).gt.0) then  
@@ -353 +353 @@
         do ifils=1,nqfils(iq)
          iq2=iqfils(ifils,iq)
-         write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
+         !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
          call vlx_loc(Ratio,pente_max,masse,u_mq,ijb_x,ije_x,iq2)
         enddo !do ifils=1,nqfils(iq)
@@ -359 +359 @@
 ! end CRisi
 
-      write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x   
+      !write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x   
 
 c   calcul des tendances
@@ -379 +379 @@
 c$OMP END DO NOWAIT
 
-      write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -398 +398 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
 
 c     CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1)
@@ -477 +477 @@
       l=3
       if ((ij.ge.ijb).and.(ij.le.ije)) then
-        write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
-     &             ij,l,iq,ijb,q(ij,l,:)
+        !write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
+!     &             ij,l,iq,ijb,q(ij,l,:)
       endif  
 
@@ -725 +725 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq)
+      !write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq)
       
       ijb=ij_begin-2*iip1