Changeset 4052 for LMDZ6/trunk/libf/dyn3dmem

Timestamp:

Dec 26, 2021, 10:27:36 PM (4 years ago)

Author:

dcugnet

Message:

Fixes for previous commit:

• few tracers(:)%nqChilds were errouneously replaced with tarcers(:)%nqDesc
• minor changes

Location:

LMDZ6/trunk/libf/dyn3dmem

Files:

• dynetat0_loc.F90 (modified) (1 diff)
• iniacademic_loc.F90 (modified) (2 diffs)
• vlsplt_loc.F (modified) (1 diff)
• vlspltgen_loc.F (modified) (1 diff)
• vlspltqs_loc.F (modified) (6 diffs)

LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

@@ -174 +174 @@
    !           Rayleigh distillation las.
     iName = tracers(iq)%iso_iName
+    IF(.NOT.ok_isotopes .OR. iName <= 0) CYCLE
     iZone = tracers(iq)%iso_iZone
     iPhase= tracers(iq)%iso_iPhase
     iqParent = tracers(iq)%iqParent
-    IF(ok_isotopes .AND. iName>0) THEN
-      IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName)    &
-     &           *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.)
-      IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase))
-    END IF
+    IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName)    &
+     &         *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.)
+    IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase))
   END DO
   DEALLOCATE(q_glo)

LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90

@@ -284 +284 @@
            do iq=1,nqtot
               q(ijb_u:ije_u,:,iq)=0.
+!              IF(tracers(iq)%name == 'H2O'//Phases_sep//'g') q(ijb_u:ije_u,:,iq)=1.e-10
+!              IF(tracers(iq)%name == 'H2O'//phases_sep//'l') q(ijb_u:ije_u,:,iq)=1.e-15
               if (tracers(iq)%name == 'H2Ov') q(ijb_u:ije_u,:,iq)=1.e-10
               if (tracers(iq)%name == 'H2Ol') q(ijb_u:ije_u,:,iq)=1.e-15
@@ -290 +292 @@
               ! distill de Rayleigh très simplifiée
               iName = tracers(iq)%iso_iName
+              if (.NOT.ok_isotopes .OR. iName <= 0) CYCLE
               iZone = tracers(iq)%iso_iZone
               iPhase= tracers(iq)%iso_iPhase
               iqParent = tracers(iq)%iqParent
-              if (ok_isotopes .AND. iName > 0) then
-                if (iZone == 0) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat(iName) &
-                                                    *(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1)
-                if (iZone == 1) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqiso(iso_indnum(iq),iPhase))
-              endif !if (ok_isotopes) then
-
+              if (iZone == 0) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat(iName) &
+                                                  *(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1)
+              if (iZone == 1) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqiso(iso_indnum(iq),iPhase))
            enddo
         else

LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F

@@ -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,nqDesc(iq)=',iq,tracers(iq)%nqDescen
+      !write(*,*) 'vly 689: iq,nqChilds(iq)=',iq,tracers(iq)%nqChilds
 
       ijb=ij_begin-2*iip1

LMDZ6/trunk/libf/dyn3dmem/vlspltgen_loc.F

@@ -346 +346 @@
         SELECT CASE(tracers(iq)%iadv)
           CASE(0); CYCLE
-          CASE(10); call vly_loc(zq,pente_max,zm,mv,iq)
+          CASE(10); call   vly_loc(zq,pente_max,zm,mv,     iq)
           CASE(14); call vlyqs_loc(zq,pente_max,zm,mv,qsat,iq)
           CASE DEFAULT

LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F

@@ -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,iqDescen(iq)=',
-!     &     iq,ijb_x,tracers(iq)%iqDescen  
+! 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
@@ -356 +356 @@
 c$OMP END DO NOWAIT
       enddo
-      do ifils=1,tracers(iq)%nqDescen
+      do ifils=1,tracers(iq)%nqChilds
         iq2=tracers(iq)%iqDescen(ifils)
         !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
@@ -727 +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,iqDescen(iq)=',iq,tracers(iq)%iqDescen
+! 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
@@ -748 +748 @@
         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
@@ -767 +767 @@
 c$OMP END DO NOWAIT
       enddo
-      do ifils=1,tracers(iq)%nqDescen
+      do ifils=1,tracers(iq)%nqChilds
         iq2=tracers(iq)%iqDescen(ifils)
         !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2