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

LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_height.F90

@@ -28 +28 @@
 
 !calcul de la hauteur max du thermique
-      do ig=1,ngrid
+      DO ig=1,ngrid
          lmax(ig)=lalim(ig)
       enddo
-      do ig=1,ngrid
-         do l=nlay,lalim(ig)+1,-1
+      DO ig=1,ngrid
+         DO l=nlay,lalim(ig)+1,-1
             IF (zw2(ig,l)<=1.e-10) THEN
                lmax(ig)=l-1
@@ -41 +41 @@
 ! On traite le cas particulier qu'il faudrait éviter ou le thermique
 ! atteind le haut du modele ...
-      do ig=1,ngrid
+      DO ig=1,ngrid
       IF ( zw2(ig,nlay) > 1.e-10 ) THEN
           PRINT*,'WARNING !!!!! W2 thermiques non nul derniere couche '
@@ -49 +49 @@
 
 ! pas de thermique si couche 1 stable
-      do ig=1,ngrid
+      DO ig=1,ngrid
          IF (lmin(ig)>1) THEN
              lmax(ig)=1
@@ -58 +58 @@
 
 ! Determination de zw2 max
-      do ig=1,ngrid
+      DO ig=1,ngrid
          wmax(ig)=0.
       enddo
 
-      do l=1,nlay
-         do ig=1,ngrid
+      DO l=1,nlay
+         DO ig=1,ngrid
             IF (l<=lmax(ig)) THEN
                 IF (zw2(ig,l)<0.)THEN
@@ -77 +77 @@
 
 !   Longueur caracteristique correspondant a la hauteur des thermiques.
-      do  ig=1,ngrid
+      DO  ig=1,ngrid
          zmax(ig)=0.
          zlevinter(ig)=zlev(ig,1)
@@ -87 +87 @@
          num(:)=0.
          denom(:)=0.
-         do ig=1,ngrid
-          do l=1,nlay
+         DO ig=1,ngrid
+          DO l=1,nlay
              num(ig)=num(ig)+zw2(ig,l)*zlev(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
              denom(ig)=denom(ig)+zw2(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
           enddo
        enddo
-       do ig=1,ngrid
+       DO ig=1,ngrid
        IF (denom(ig)>1.e-10) THEN
           zmax(ig)=2.*num(ig)/denom(ig)
@@ -102 +102 @@
       else
 !CR:Calcul de zmax continu via le linter      
-      do  ig=1,ngrid
+      DO  ig=1,ngrid
 ! calcul de zlevinter
           zlevinter(ig)=(zlev(ig,lmax(ig)+1)-zlev(ig,lmax(ig)))*  &
@@ -114 +114 @@
 
 !CR:calcul du zcong
-      do  ig=1,ngrid
+      DO  ig=1,ngrid
 ! calcul de zlevintercong
          zlevintercong(ig)=(zlev(ig,lcong(ig)+1)-zlev(ig,lcong(ig)))*  &
@@ -127 +127 @@
 
 ! def de  zmix continu (profil parabolique des vitesses)
-      do ig=1,ngrid
+      DO ig=1,ngrid
            IF (lmix(ig)>1) THEN
 ! test 
@@ -158 +158 @@
 
 ! calcul du nouveau lmix correspondant
-      do ig=1,ngrid
-         do l=1,nlay
+      DO ig=1,ngrid
+         DO l=1,nlay
             IF (zmix(ig)>=zlev(ig,l).AND.  &
             zmix(ig)<zlev(ig,l+1)) THEN