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

Timestamp:

Aug 2, 2024, 2:12:03 PM (3 months ago)

Author:

abarral

Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F90 (modified) (19 diffs)

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

@@ -71 +71 @@
   CALL SCOPY(ijp1llm, masse, 1, zm(1, 1, iq), 1)
 
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     CALL SCOPY(ijp1llm, q(1, 1, iq2), 1, zq(1, 1, iq2), 1)
@@ -91 +91 @@
   ENDDO
   ! CRisi: aussi pour les fils
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -317 +317 @@
             i = ijq - (j - 1) * iip1
             !   accumulation pour les mailles completements advectees
-            do while(zu_m>masse(ijq, l, iq))
+            DO while(zu_m>masse(ijq, l, iq))
               u_mq(ij, l) = u_mq(ij, l) + q(ijq, l, iq) &
                       * masse(ijq, l, iq)
@@ -332 +332 @@
             i = ijq - (j - 1) * iip1
             !   accumulation pour les mailles completements advectees
-            do while(-zu_m>masse(ijq, l, iq))
+            DO while(-zu_m>masse(ijq, l, iq))
               u_mq(ij, l) = u_mq(ij, l) - q(ijq, l, iq) &
                       * masse(ijq, l, iq)
@@ -361 +361 @@
   !WRITE(*,*) 'vlsplt 326: iq,nqDesc(iq)=',iq,tracers(iq)%nqDescen
 
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -378 +378 @@
     enddo
   enddo
-  do ifils = 1, tracers(iq)%nqChildren
+  DO ifils = 1, tracers(iq)%nqChildren
     iq2 = tracers(iq)%iqDescen(ifils)
     CALL vlx(Ratio, pente_max, masseq, u_mq, iq2)
@@ -405 +405 @@
   ! On calcule q entre iip2+1,ip1jm -> on fait pareil pour ratio
   ! puis on boucle en longitude
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -431 +431 @@
   !    ********************************************************************
   ! q,masse_adv_v,w sont des arguments d'entree  pour le s-pg ....
-  ! dq         sont des arguments de sortie pour le s-pg ....
+  ! dq            sont des arguments de sortie pour le s-pg ....
   !
   !
@@ -479 +479 @@
     PRINT*, 'Shema  Amont nouveau  appele dans  Vanleer   '
     first = .FALSE.
-    do i = 2, iip1
+    DO i = 2, iip1
       coslon(i) = cos(rlonv(i))
       sinlon(i) = sin(rlonv(i))
@@ -673 +673 @@
   ! WRITE(*,*) 'vly 689: iq,nqDesc(iq)=',iq,tracers(iq)%nqDescen
 
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -692 +692 @@
   enddo
 
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     CALL vly(Ratio, pente_max, masseq, qbyv, iq2)
@@ -709 +709 @@
     massepn = ssum(iim, masse(1, l, iq), 1)
     qpn = 0.
-    do ij = 1, iim
+    DO ij = 1, iim
       qpn = qpn + masse(ij, l, iq) * q(ij, l, iq)
     enddo
     qpn = (qpn + convpn) / (massepn + convmpn)
-    do ij = 1, iip1
+    DO ij = 1, iip1
       q(ij, l, iq) = qpn
     enddo
@@ -720 +720 @@
     masseps = ssum(iim, masse(ip1jm + 1, l, iq), 1)
     qps = 0.
-    do ij = ip1jm + 1, ip1jmp1 - 1
+    DO ij = ip1jm + 1, ip1jmp1 - 1
       qps = qps + masse(ij, l, iq) * q(ij, l, iq)
     enddo
     qps = (qps + convps) / (masseps + convmps)
-    do ij = ip1jm + 1, ip1jmp1
+    DO ij = ip1jm + 1, ip1jmp1
       q(ij, l, iq) = qps
     enddo
@@ -730 +730 @@
 
   ! retablir les fils en rapport de melange par rapport a l'air:
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -752 +752 @@
   !    ********************************************************************
   !    q,pbaru,pbarv,w sont des arguments d'entree  pour le s-pg ....
-  ! dq         sont des arguments de sortie pour le s-pg ....
+  ! dq            sont des arguments de sortie pour le s-pg ....
   !   --------------------------------------------------------------------
   IMPLICIT NONE
@@ -821 +821 @@
 
   DO l = 1, llm - 1
-    do  ij = 1, ip1jmp1
+    DO  ij = 1, ip1jmp1
       IF(w(ij, l + 1)>0.) THEN
         sigw = w(ij, l + 1) / masse(ij, l + 1, iq)
@@ -841 +841 @@
   ! Il faut faire ça avant d'avoir mis à jour q et masse
   ! !WRITE(*,*) 'vlsplt 942: iq,nqChildren(iq)=',iq,nqChildren(iq)
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm
@@ -858 +858 @@
   enddo
 
-  do ifils = 1, tracers(iq)%nqChildren
+  DO ifils = 1, tracers(iq)%nqChildren
     iq2 = tracers(iq)%iqDescen(ifils)
     CALL vlz(Ratio, pente_max, masseq, wq, iq2)
@@ -874 +874 @@
 
   ! retablir les fils en rapport de melange par rapport a l'air:
-  do ifils = 1, tracers(iq)%nqDescen
+  DO ifils = 1, tracers(iq)%nqDescen
     iq2 = tracers(iq)%iqDescen(ifils)
     DO l = 1, llm