Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F90

Timestamp:

Jul 24, 2024, 1:17:08 PM (8 weeks 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/dyn3d/vlsplt.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F90

@@ -89 +89 @@
      ENDDO
   ENDDO
-  ! ! CRisi: aussi pour les fils
+  ! CRisi: aussi pour les fils
   do ifils=1,tracers(iq)%nqDescen
     iq2=tracers(iq)%iqDescen(ifils)
@@ -143 +143 @@
   REAL :: u_mq(ip1jmp1,llm)
 
-  ! ! CRisi
+  ! CRisi
   REAL :: masseq(ip1jmp1,llm,nqtot),Ratio(ip1jmp1,llm,nqtot)
   INTEGER :: ifils,iq2 ! CRisi
@@ -361 +361 @@
   ! 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,nqDesc(iq)=',iq,tracers(iq)%nqDescen
+  !write(*,*) 'vlsplt 326: iq,nqDesc(iq)=',iq,tracers(iq)%nqDescen
 
   do ifils=1,tracers(iq)%nqDescen
@@ -367 +367 @@
     DO l=1,llm
       DO ij=iip2,ip1jm
-        ! ! On a besoin de q et masse seulement entre iip2 et ip1jm
-        ! !masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq)
+        ! On a besoin de q et masse seulement entre iip2 et ip1jm
+        !masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq)
   !           !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
-        ! !Mvals: veiller a ce qu'on n'ait pas de denominateur nul
+        !Mvals: veiller a ce qu'on n'ait pas de denominateur nul
         masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass)
         if (q(ij,l,iq)>min_qParent) then
@@ -391 +391 @@
   DO l=1,llm
      DO ij=iip2+1,ip1jm
-        ! !MVals: veiller a ce qu'on ait pas de denominateur nul
+        !MVals: veiller a ce qu'on 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)+ &
@@ -404 +404 @@
   ENDDO
 
-  ! ! retablir les fils en rapport de melange par rapport a l'air:
-  ! ! On calcule q entre iip2+1,ip1jm -> on fait pareil pour ratio
-  ! ! puis on boucle en longitude
+  ! retablir les fils en rapport de melange par rapport a l'air:
+  ! On calcule q entre iip2+1,ip1jm -> on fait pareil pour ratio
+  ! puis on boucle en longitude
   do ifils=1,tracers(iq)%nqDescen
     iq2=tracers(iq)%iqDescen(ifils)