Changeset 2286 for LMDZ5/trunk/libf/dyn3d

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

Location:

LMDZ5/trunk/libf/dyn3d

Files:

• advtrac.F90 (modified) (2 diffs)
• caladvtrac.F (modified) (1 diff)
• leapfrog.F (modified) (1 diff)
• qminimum.F (modified) (4 diffs)
• vlsplt.F (modified) (10 diffs)
• vlspltqs.F (modified) (6 diffs)

LMDZ5/trunk/libf/dyn3d/advtrac.F90

@@ -238 +238 @@
            ! CRisi: on fait passer tout q pour avoir acces aux fils
            
-           write(*,*) 'advtrac 239: iq,q(1721,19,:)=',iq,q(1721,19,:)     
+           !write(*,*) 'advtrac 239: iq,q(1721,19,:)=',iq,q(1721,19,:)     
            call vlsplt(q,2.,massem,wg,pbarug,pbarvg,dtvr,iq)
            
@@ -247 +247 @@
         else if(iadv(iq).eq.14) then
            !
-           write(*,*) 'advtrac 248: iq,q(1721,19,:)=',iq,q(1721,19,:)
+           !write(*,*) 'advtrac 248: iq,q(1721,19,:)=',iq,q(1721,19,:)
            CALL vlspltqs( q, 2., massem, wg , &
                 pbarug,pbarvg,dtvr,p,pk,teta,iq)

LMDZ5/trunk/libf/dyn3d/caladvtrac.F

@@ -85 +85 @@
           ENDDO
 
-          write(*,*) 'caladvtrac 87'
+          !write(*,*) 'caladvtrac 87'
           CALL qminimum( q, nqtot, finmasse )
-          write(*,*) 'caladvtrac 89'
+          !write(*,*) 'caladvtrac 89'
 
           CALL SCOPY   ( ip1jmp1*llm, masse, 1, finmasse,       1 )

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -344 +344 @@
      *        p, masse, dq,  teta,
      .        flxw, pk)
-          write(*,*) 'caladvtrac 346'
+          !write(*,*) 'caladvtrac 346'
 
          

LMDZ5/trunk/libf/dyn3d/qminimum.F

@@ -55 +55 @@
 
       zx_defau_diag(:,:,:)=0.0
-      q_follow(:,:,:)=q(:,:,:)  
+      q_follow(:,:,1:2)=q(:,:,1:2)  
       DO 1000 k = 1, llm
         DO 1040 i = 1, ip1jmp1
@@ -110 +110 @@
       ENDIF
 
-      write(*,*) 'qminimum 128'
+      !write(*,*) 'qminimum 128'
       if (ok_isotopes) then
       ! CRisi: traiter de même les traceurs d'eau
@@ -158 +158 @@
       
         ! 3) transfert d'eau de la vapeur au liquide
-        write(*,*) 'qminimum 164'
+        !write(*,*) 'qminimum 164'
         do k=1,llm
         DO i = 1,ip1jmp1
@@ -186 +186 @@
 
       endif !if (ok_isotopes) then
-      write(*,*) 'qminimum 188'
+      !write(*,*) 'qminimum 188'
       
 c

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

LMDZ5/trunk/libf/dyn3d/vlspltqs.F

@@ -175 +175 @@
       enddo !do ifils=1,nqdesc(iq)  
       endif ! if (nqfils(iq).gt.0) then
-      write(*,*) 'vlspltqs 183: fin de la routine'
+      !write(*,*) 'vlspltqs 183: fin de la routine'
 
       RETURN
@@ -486 +486 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlspltqs 326: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vlspltqs 326: iq,nqfils(iq)=',iq,nqfils(iq) 
       
       if (nqfils(iq).gt.0) then  
@@ -797 +797 @@
 ! 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) 
    
       if (nqfils(iq).gt.0) then  
@@ -812 +812 @@
         do ifils=1,nqfils(iq)
          iq2=iqfils(ifils,iq)
-         write(*,*) 'vlyqs 783: appel rec de vly, iq2=',iq2
+         !write(*,*) 'vlyqs 783: appel rec de vly, iq2=',iq2
          call vly(Ratio,pente_max,masseq,qbyv,iq2)
         enddo !do ifils=1,nqfils(iq)
@@ -868 +868 @@
       ENDDO
 
-      write(*,*) 'vly 866'
+      !write(*,*) 'vly 866'
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -881 +881 @@
         enddo !do ifils=1,nqdesc(iq)
       endif !if (nqfils(iq).gt.0) then
-      write(*,*) 'vly 879'
+      !write(*,*) 'vly 879'
 
       RETURN