Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.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/phylmd/lmdz_thermcell_plume.f90 (modified) (10 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.f90

@@ -169 +169 @@
     ! dans une couche l>1
     !------------------------------------------------------------------------------
-    do ig = 1, ngrid
+    DO ig = 1, ngrid
       ! Le panache va prendre au debut les caracteristiques de l'air contenu
       ! dans cette couche.
@@ -188 +188 @@
     !boucle de calcul de la vitesse verticale dans le thermique
     !==============================================================================
-    do l = 2, nlay - 1
+    DO l = 2, nlay - 1
       !==============================================================================
 
@@ -194 +194 @@
       ! On decide si le thermique est encore actif ou non
       ! AFaire : Il faut sans doute ajouter entr_star a alim_star dans ce test
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         active(ig) = active(ig) &
                 &                 .AND. zw2(ig, l)>1.e-10 &
@@ -214 +214 @@
       ztemp(:) = zpspsk(:, l) * ztla(:, l - 1)
       CALL thermcell_qsat(ngrid, active, pplev(:, l), ztemp, zqta(:, l - 1), zqsat(:))
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         !       PRINT*,'active',active(ig),ig,l
         IF(active(ig)) THEN
@@ -250 +250 @@
           zdz2 = zlev(ig, lt) - zlev(ig, l)
 
-          do while (lmel>zdz2)
+          DO while (lmel>zdz2)
             lt = lt + 1
             zlt = zlev(ig, lt)
@@ -283 +283 @@
       !-------------------------------------------------
 
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         IF (active(ig)) THEN
           !          zw2m=max(0.5*(w_est(ig,l)+w_est(ig,l+1)),0.1)
@@ -331 +331 @@
 
       activetmp(:) = active(:) .AND. f_star(:, l + 1)>1.e-10
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         IF (activetmp(ig)) THEN
           Zsat = .FALSE.
@@ -346 +346 @@
       ztemp(:) = zpspsk(:, l) * ztla(:, l)
       CALL thermcell_qsat(ngrid, activetmp, pplev(:, l), ztemp, zqta(:, l), zqsatth(:, l))
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         IF (activetmp(ig)) THEN
           ! on ecrit de maniere conservative (sat ou non)
@@ -376 +376 @@
 
       nbpb = 0
-      do ig = 1, ngrid
+      DO ig = 1, ngrid
         IF (zw2(ig, l + 1)>0. .AND. zw2(ig, l + 1)<1.e-10) THEN
           !               stop 'On tombe sur le cas particulier de thermcell_dry'
@@ -420 +420 @@
 
     !on recalcule alim_star_tot
-    do ig = 1, ngrid
+    DO ig = 1, ngrid
       alim_star_tot(ig) = 0.
     enddo
-    do ig = 1, ngrid
-      do l = 1, lalim(ig) - 1
+    DO ig = 1, ngrid
+      DO l = 1, lalim(ig) - 1
         alim_star_tot(ig) = alim_star_tot(ig) + alim_star(ig, l)
       enddo