Changeset 4171 for LMDZ6/branches/LMDZ-ECRAD/libf/dyn3dmem/vlspltqs_loc.F

Timestamp:

Jun 17, 2022, 4:24:49 PM (2 years ago)

Author:

lguez

Message:

Sync latest trunk changes to branch LMDZ-ECRAD.

Location:

LMDZ6/branches/LMDZ-ECRAD

Files:

• . (modified) (1 prop)
• libf/dyn3dmem/vlspltqs_loc.F (modified) (9 diffs)

LMDZ6/branches/LMDZ-ECRAD/libf/dyn3dmem/vlspltqs_loc.F

@@ -12 +12 @@
 c   --------------------------------------------------------------------
       USE parallel_lmdz
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi                 &
-     &                     qperemin,masseqmin,ratiomin ! MVals et CRisi
+      USE infotrac, ONLY : nqtot,tracers, ! CRisi                 &
+     &                     min_qParent,min_qMass,min_ratio ! MVals et CRisi
       IMPLICIT NONE
 c
@@ -264 +264 @@
 c$OMP END DO NOWAIT
 
-cym ATTENTION ICI en OpenMP reduction pas forcement nécessaire
+cym ATTENTION ICI en OpenMP reduction pas forcement necessaire
 cym      IF(n0.gt.1) THEN
 cym        IF(n0.gt.0) THEN
@@ -335 +335 @@
 c$OMP END DO NOWAIT
 
-! 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)  
-
-      if (nqfils(iq).gt.0) then  
-       do ifils=1,nqdesc(iq)
-         iq2=iqfils(ifils,iq)
+! CRisi: appel recursif de l'advection sur les fils.
+! Il faut faire ca avant d'avoir mis a jour q et masse
+      !write(*,*) 'vlspltqs 336: iq,ijb_x,nqChilds(iq)=',
+!     &     iq,ijb_x,tracers(iq)%nqChilds
+
+      do ifils=1,tracers(iq)%nqDescen
+        iq2=tracers(iq)%iqDescen(ifils)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
-         DO l=1,llm
+        DO l=1,llm
           DO ij=ijb,ije
-           !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
-           if (q(ij,l,iq).gt.qperemin) then ! modif 13 nov 2020
-             Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
-           else
-             Ratio(ij,l,iq2)=ratiomin
-           endif 
+            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+            masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass)
+            if (q(ij,l,iq).gt.min_qParent) then ! modif 13 nov 2020
+              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
+            else
+              Ratio(ij,l,iq2)=min_ratio
+            endif 
           enddo   
-         enddo
-c$OMP END DO NOWAIT
-        enddo !do ifils=1,nqfils(iq)
-        do ifils=1,nqfils(iq)
-         iq2=iqfils(ifils,iq)
-         !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)
-      endif !if (nqfils(iq).gt.0) then
+        enddo
+c$OMP END DO NOWAIT
+      enddo
+      do ifils=1,tracers(iq)%nqChilds
+        iq2=tracers(iq)%iqDescen(ifils)
+        !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
+        call vlx_loc(Ratio,pente_max,masse,u_mq,ijb_x,ije_x,iq2)
+      enddo
 ! end CRisi
 
@@ -372 +370 @@
          DO ij=ijb+1,ije
             !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),masseqmin)
+            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),min_qMass)
             q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+
      &      u_mq(ij-1,l)-u_mq(ij,l))
@@ -389 +387 @@
 
 ! retablir les fils en rapport de melange par rapport a l'air:
-      if (nqfils(iq).gt.0) then  
-       do ifils=1,nqdesc(iq)
-         iq2=iqfils(ifils,iq)  
+      do ifils=1,tracers(iq)%nqDescen
+        iq2=tracers(iq)%iqDescen(ifils)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
-         DO l=1,llm
+        DO l=1,llm
           DO ij=ijb+1,ije
             q(ij,l,iq2)=q(ij,l,iq)*Ratio(ij,l,iq2)            
           enddo
           DO ij=ijb+iip1-1,ije,iip1
-             q(ij-iim,l,iq2)=q(ij,l,iq2)
+            q(ij-iim,l,iq2)=q(ij,l,iq2)
           enddo ! DO ij=ijb+iip1-1,ije,iip1
-         enddo
-c$OMP END DO NOWAIT
-        enddo !do ifils=1,nqdesc(iq)
-      endif !if (nqfils(iq).gt.0) then
+        enddo
+c$OMP END DO NOWAIT
+      enddo
 
       !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
@@ -426 +422 @@
 c   --------------------------------------------------------------------
       USE parallel_lmdz
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi                 &
-     &                     qperemin,masseqmin,ratiomin ! MVals et CRisi
+      USE infotrac, ONLY : nqtot,tracers, ! CRisi                 &
+     &                     min_qParent,min_qMass,min_ratio ! MVals et CRisi
       USE comconst_mod, ONLY: pi
       IMPLICIT NONE
@@ -731 +727 @@
 c$OMP END DO NOWAIT
 
-! 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)
+! CRisi: appel recursif de l'advection sur les fils.
+! Il faut faire ca avant d'avoir mis a jour q et masse
+      !write(*,*) 'vlyqs 689: iq,nqChilds(iq)=',iq,tracers(iq)%nqChilds
       
       ijb=ij_begin-2*iip1
@@ -747 +743 @@
       !write(lunout,*) 'ij_begin,ij_end=',ij_begin,ij_end
       !write(lunout,*) 'pole_nord,pole_sud=',pole_nord,pole_sud
-      if (nqfils(iq).gt.0) then  
-       do ifils=1,nqdesc(iq)
-         iq2=iqfils(ifils,iq)
+      do ifils=1,tracers(iq)%nqDescen
+        iq2=tracers(iq)%iqDescen(ifils)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
-         DO l=1,llm
+        DO l=1,llm
           ! modif des bornes: CRisi 16 nov 2020
-          ! d'abord masse avec bornes corrigées
+          ! d'abord masse avec bornes corrigees
           DO ij=ijbm,ijem 
-           !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+            masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass)
           enddo !DO ij=ijbm,ijem 
 
           ! ensuite Ratio avec anciennes bornes
           DO ij=ijb,ije
-           !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           !write(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq)
-           if (q(ij,l,iq).gt.qperemin) then ! modif 13 nov 2020
-             Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
-           else
-             Ratio(ij,l,iq2)=ratiomin   
-           endif
+            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+            !write(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq)
+            if (q(ij,l,iq).gt.min_qParent) then ! modif 13 nov 2020
+              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
+            else
+              Ratio(ij,l,iq2)=min_ratio   
+            endif
           enddo !DO ij=ijbm,ijem
-         enddo !DO l=1,llm
-c$OMP END DO NOWAIT
-        enddo !do ifils=1,nqdesc(iq)
-        do ifils=1,nqfils(iq)
-         iq2=iqfils(ifils,iq)
-         !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2
-         call vly_loc(Ratio,pente_max,masse,qbyv,iq2)
-        enddo !do ifils=1,nqfils(iq)
-      endif !if (nqfils(iq).gt.0) then
+        enddo !DO l=1,llm
+c$OMP END DO NOWAIT
+      enddo
+      do ifils=1,tracers(iq)%nqChilds
+        iq2=tracers(iq)%iqDescen(ifils)
+        !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2
+        call vly_loc(Ratio,pente_max,masse,qbyv,iq2)
+      enddo
 
         
@@ -856 +850 @@
 !      if (pole_sud)  ije=ij_end-iip1
   
-      if (nqfils(iq).gt.0) then  
-       do ifils=1,nqdesc(iq)
-         iq2=iqfils(ifils,iq)  
+      do ifils=1,tracers(iq)%nqDescen
+        iq2=tracers(iq)%iqDescen(ifils)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
-         DO l=1,llm
+        DO l=1,llm
           DO ij=ijb,ije
             q(ij,l,iq2)=q(ij,l,iq)*Ratio(ij,l,iq2)            
           enddo
-         enddo
-c$OMP END DO NOWAIT
-        enddo !do ifils=1,nqdesc(iq)
-      endif !if (nqfils(iq).gt.0) then
+        enddo
+c$OMP END DO NOWAIT
+      enddo