Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90

Timestamp:

Jul 24, 2024, 1:17:08 PM (4 months ago)

Author:

abarral

Message:

Rename modules in misc from *_mod > lmdz_*
Put cbrt.f90, ch*.f90, pch*.f90 in new lmdz_libmath_pch.f90

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90 (modified) (16 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90

@@ -53 +53 @@
   INTEGER :: ijb, ije, ijb_x, ije_x
 
-  ! !write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
+  !write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
   ! &   iq,ijb_x
   !   calcul de la pente a droite et a gauche de la maille
@@ -68 +68 @@
     !   calcul des pentes avec limitation, Van Leer scheme I:
     !   -----------------------------------------------------
-    ! ! on a besoin de q entre ijb et ije
+    ! on a besoin de q entre ijb et ije
     !   calcul de la pente aux points u
     !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
@@ -145 +145 @@
   ENDIF ! (pente_max.lt.-1.e-5)
 
-  ! !write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
+  !write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
 
   !   bouclage de la pente en iip1:
@@ -169 +169 @@
   ! PRINT*,'Cumule ....'
   !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-  ! ! on a besoin de masse entre ijb et ije
+  ! on a besoin de masse entre ijb et ije
   DO l = 1, llm
     DO ij = ijb, ije - 1
@@ -244 +244 @@
       ENDDO
       niju = iju
-      ! !PRINT*,'vlx 278, niju,nl',niju,nl(l)
+      !PRINT*,'vlx 278, niju,nl',niju,nl(l)
 
       !  traitement des mailles
@@ -302 +302 @@
 
   do ifils = 1, tracers(iq)%nqDescen
-    ! ! attention: comme Ratio est utilisé comme q dans l'appel
-    ! ! recursif, il doit contenir à lui seul tous les indices de tous
-    ! ! les descendants!
+    ! attention: comme Ratio est utilisé comme q dans l'appel
+    ! recursif, il doit contenir à lui seul tous les indices de tous
+    ! les descendants!
     iq2 = tracers(iq)%iqDescen(ifils)
     !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l = 1, llm
       DO ij = ijb, ije
-        ! ! On a besoin de q et masse seulement entre ijb et ije. On ne
-        ! ! les calcule donc que de ijb à ije
-        ! !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+        ! On a besoin de q et masse seulement entre ijb et ije. On ne
+        ! les calcule donc que de 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), min_qMass)
         if (q(ij, l, iq)>min_qParent) then ! modif 13 nov 2020
@@ -333 +333 @@
   DO l = 1, llm
     DO ij = ijb + 1, ije
-      ! !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+      !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), min_qMass)
       q(ij, l, iq) = (q(ij, l, iq) * masse(ij, l, iq) + &
@@ -349 +349 @@
 
   ! retablir les fils en rapport de melange par rapport a l'air:
-  ! ! On calcule q entre ijb+1 et ije -> on fait pareil pour ratio
-  ! ! puis on boucle en longitude
+  ! On calcule q entre ijb+1 et ije -> on fait pareil pour ratio
+  ! puis on boucle en longitude
   do ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
@@ -495 +495 @@
     if (pole_sud)  ije = ij_end - iip1
 
-    ! ! on a besoin de q entre ij_begin-2*iip1 et ij_end+2*iip1
-    ! ! Si pole sud, entre ij_begin-2*iip1 et ij_end
-    ! ! Si pole Nord, entre ij_begin et ij_end+2*iip1
+    ! on a besoin de q entre ij_begin-2*iip1 et ij_end+2*iip1
+    ! Si pole sud, entre ij_begin-2*iip1 et ij_end
+    ! Si pole Nord, entre ij_begin et ij_end+2*iip1
     DO ij = ijb, ije
       dyqv(ij) = q(ij, l, iq) - q(ij + iip1, l, iq)
@@ -706 +706 @@
     !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l = 1, llm
-      ! ! modif des bornes: CRisi 16 nov 2020
-      ! ! d'abord masse avec bornes corrigées
+      ! modif des bornes: CRisi 16 nov 2020
+      ! d'abord masse avec bornes corrigées
       DO ij = ijbm, ijem
-        ! !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+        !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
 
-      ! ! ensuite Ratio avec anciennes bornes
+      ! ensuite Ratio avec anciennes bornes
       DO ij = ijb, ije
-        ! !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+        !MVals: veiller a ce qu'on n'ait pas de denominateur nul
         if (q(ij, l, iq)>min_qParent) then ! modif 13 nov 2020
           Ratio(ij, l, iq2) = q(ij, l, iq2) / q(ij, l, iq)
@@ -859 +859 @@
   !$OMP THREADPRIVATE(first)
 
-  ! !REAL masseq(ijb_u:ije_u,llm,nqtot),Ratio(ijb_u:ije_u,llm,nqtot) ! CRisi
-  ! ! Ces varibles doivent être déclarées en pointer et en save dans
-  ! ! vlz_loc si on veut qu'elles soient vues par tous les threads.
+  !REAL masseq(ijb_u:ije_u,llm,nqtot),Ratio(ijb_u:ije_u,llm,nqtot) ! CRisi
+  ! Ces varibles doivent être déclarées en pointer et en save dans
+  ! vlz_loc si on veut qu'elles soient vues par tous les threads.
   INTEGER :: ifils, iq2 ! CRisi
 
@@ -870 +870 @@
   !    sens de W
 
-  ! !write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
+  !write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
 
   ijb = ijb_x
@@ -954 +954 @@
     ENDDO
     !$OMP END DO NOWAIT
-    ! !write(*,*) 'vlz 1001'
+    !write(*,*) 'vlz 1001'
 
   ELSE ! countcfl>=1
@@ -1013 +1013 @@
               lorig(ij, l) = lorig(ij, l) - 1
             ENDIF
-            ! ! CRisi 24nov2020: ajout d'un message d'erreur clair au lieu d'un plantage
-            ! ! pour seg fault
+            ! CRisi 24nov2020: ajout d'un message d'erreur clair au lieu d'un plantage
+            ! pour seg fault
             if (lorig(ij, l)==0) then
               CALL abort_gcm("vlz in vlsplt_loc", &
@@ -1064 +1064 @@
     DO l = 1, llm
       DO ij = ijb, ije
-        ! !MVals: veiller a ce qu'on n'ait pas de denominateur nul
+        !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)>min_qParent) then
@@ -1071 +1071 @@
           Ratio(ij, l, iq2) = min_ratio
         endif
-        ! !wq(ij,l,iq2)=wq(ij,l,iq) ! correction bug le 15mai2015
+        !wq(ij,l,iq2)=wq(ij,l,iq) ! correction bug le 15mai2015
         w(ij, l, iq2) = wq(ij, l, iq)
       enddo