Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_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/vlspltqs_loc.f90 (modified) (11 diffs)

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

@@ -48 +48 @@
   INTEGER :: ijb,ije,ijb_x,ije_x
 
-  ! !write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
+  !write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
   ! &   iq,ijb_x
 
@@ -253 +253 @@
            ENDDO
            niju=iju
-           ! !PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
+           !PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
 
   !  traitement des mailles
@@ -309 +309 @@
   ! 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,nqChildren(iq)=',
+  !write(*,*) 'vlspltqs 336: iq,ijb_x,nqChildren(iq)=',
   ! &     iq,ijb_x,tracers(iq)%nqChildren
 
@@ -317 +317 @@
     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 ! modif 13 nov 2020
@@ -330 +330 @@
   do ifils=1,tracers(iq)%nqChildren
     iq2=tracers(iq)%iqDescen(ifils)
-    ! !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
+    !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
 
-  ! !write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x
+  !write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x
 
   !   calcul des tendances
@@ -341 +341 @@
   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)+ &
@@ -356 +356 @@
 !$OMP END DO NOWAIT
 
-  ! !write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
+  !write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
 
   ! retablir les fils en rapport de melange par rapport a l'air:
@@ -373 +373 @@
   enddo
 
-  ! !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
+  !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
 
   ! CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1)
@@ -453 +453 @@
   l=3
   if ((ij>=ijb).and.(ij<=ije)) then
-    ! !write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
+    !write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
   ! &             ij,l,iq,ijb,q(ij,l,:)
   endif
@@ -713 +713 @@
   if (pole_sud)  ijem=ij_end
 
-  ! !write(lunout,*) 'vlspltqs 737: iq,ijb,ije=',iq,ijb,ije
-  ! !write(lunout,*) 'ij_begin,ij_end=',ij_begin,ij_end
-  ! !write(lunout,*) 'pole_nord,pole_sud=',pole_nord,pole_sud
+  !write(lunout,*) 'vlspltqs 737: iq,ijb,ije=',iq,ijb,ije
+  !write(lunout,*) 'ij_begin,ij_end=',ij_begin,ij_end
+  !write(lunout,*) 'pole_nord,pole_sud=',pole_nord,pole_sud
   do ifils=1,tracers(iq)%nqDescen
     iq2=tracers(iq)%iqDescen(ifils)
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l=1,llm
-      ! ! modif des bornes: CRisi 16 nov 2020
-      ! ! d'abord masse avec bornes corrigees
+      ! modif des bornes: CRisi 16 nov 2020
+      ! d'abord masse avec bornes corrigees
       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 !DO ij=ijbm,ijem
 
-      ! ! 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
-        ! !write(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq)
+        !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)>min_qParent) then ! modif 13 nov 2020
           Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
@@ -742 +742 @@
   do ifils=1,tracers(iq)%nqChildren
     iq2=tracers(iq)%iqDescen(ifils)
-    ! !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2
+    !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2
     CALL vly_loc(Ratio,pente_max,masse,qbyv,iq2)
   enddo