Changeset 2286 for LMDZ5/trunk/libf/dyn3d/vlsplt.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/dyn3d/vlsplt.F (modified) (10 diffs)

LMDZ5/trunk/libf/dyn3d/vlsplt.F

@@ -454 +454 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlsplt 326: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vlsplt 326: iq,nqfils(iq)=',iq,nqfils(iq) 
       
       if (nqdesc(iq).gt.0) then  
@@ -579 +579 @@
       DATA temps0,temps1,temps2,temps3,temps4,temps5/0.,0.,0.,0.,0.,0./
 
-      write(*,*) 'vly 578: entree, iq=',iq
+      !write(*,*) 'vly 578: entree, iq=',iq
 
       IF(first) THEN
@@ -758 +758 @@
       ENDDO
 
-      write(*,*) 'vly 756'
+      !write(*,*) 'vly 756'
       DO l=1,llm
        DO ij=1,ip1jm
@@ -776 +776 @@
 ! 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) 
    
       if (nqfils(iq).gt.0) then  
@@ -874 +874 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vly 853: sortie'
+      !write(*,*) 'vly 853: sortie'
 
       RETURN
@@ -933 +933 @@
 c    sens de W
 
-      write(*,*) 'vlz 923: entree'
+      !write(*,*) 'vlz 923: entree'
 
 #ifdef BIDON
@@ -964 +964 @@
       ENDDO
 
-      write(*,*) 'vlz 954'
+      !write(*,*) 'vlz 954'
       DO ij=1,ip1jmp1
          dzq(ij,1)=0.
@@ -981 +981 @@
 c calcul de  - d( q   * w )/ d(sigma)    qu'on ajoute a  dq pour calculer dq
 
-       write(*,*) 'vlz 969'
+       !write(*,*) 'vlz 969'
        DO l = 1,llm-1
          do  ij = 1,ip1jmp1
@@ -1002 +1002 @@
 ! 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)
@@ -1041 +1041 @@
         enddo !do ifils=1,nqdesc(iq)
       endif !if (nqfils(iq).gt.0) then
-      write(*,*) 'vlsplt 1032'
+      !write(*,*) 'vlsplt 1032'
 
       RETURN