Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_height.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_height.F90 (modified) (10 diffs)

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

@@ -13 +13 @@
 
 ! Entree
-      integer, intent(in) :: ngrid,nlay
-      real, intent(in), dimension(ngrid) :: linter,lintercong
-      real, intent(in), dimension(ngrid,nlay+1) :: zlev
+      INTEGER, INTENT(IN) :: ngrid,nlay
+      REAL, INTENT(IN), DIMENSION(ngrid) :: linter,lintercong
+      REAL, INTENT(IN), DIMENSION(ngrid,nlay+1) :: zlev
 ! Sortie
-      real, intent(out), dimension(ngrid) :: wmax,zmax,zmax0,zmix,zcong
-      integer, intent(out), dimension(ngrid) :: lmax,lcong
+      REAL, INTENT(OUT), DIMENSION(ngrid) :: wmax,zmax,zmax0,zmix,zcong
+      INTEGER, INTENT(OUT), DIMENSION(ngrid) :: lmax,lcong
 ! Les deux
-     integer, intent(inout), dimension(ngrid) :: lmix,lalim,lmin
-     real, intent(inout), dimension(ngrid,nlay+1) :: zw2
+     INTEGER, INTENT(INOUT), DIMENSION(ngrid) :: lmix,lalim,lmin
+     REAL, INTENT(INOUT), DIMENSION(ngrid,nlay+1) :: zw2
 
 ! local
-     real, dimension(ngrid) :: num,denom,zlevinter,zlevintercong
-     integer ig,l
+     REAL, DIMENSION(ngrid) :: num,denom,zlevinter,zlevintercong
+     INTEGER ig,l
 
 !calcul de la hauteur max du thermique
@@ -33 +33 @@
       do ig=1,ngrid
          do l=nlay,lalim(ig)+1,-1
-            if (zw2(ig,l)<=1.e-10) THEN
+            IF (zw2(ig,l)<=1.e-10) THEN
                lmax(ig)=l-1
             endif
@@ -42 +42 @@
 ! atteind le haut du modele ...
       do ig=1,ngrid
-      if ( zw2(ig,nlay) > 1.e-10 ) THEN
+      IF ( zw2(ig,nlay) > 1.e-10 ) THEN
           PRINT*,'WARNING !!!!! W2 thermiques non nul derniere couche '
           lmax(ig)=nlay
@@ -50 +50 @@
 ! pas de thermique si couche 1 stable
       do ig=1,ngrid
-         if (lmin(ig)>1) THEN
+         IF (lmin(ig)>1) THEN
              lmax(ig)=1
              lmin(ig)=1
@@ -64 +64 @@
       do l=1,nlay
          do ig=1,ngrid
-            if (l<=lmax(ig)) THEN
-                if (zw2(ig,l)<0.)THEN
+            IF (l<=lmax(ig)) THEN
+                IF (zw2(ig,l)<0.)THEN
                   PRINT*,'pb2 zw2<0'
                 endif
@@ -83 +83 @@
 
 !     if (iflag_thermals_ed.ge.1) THEN
-      if (1==0) THEN
+      IF (1==0) THEN
 !CR:date de quand le calcul du zmax continu etait buggue 
          num(:)=0.
@@ -94 +94 @@
        enddo
        do ig=1,ngrid
-       if (denom(ig)>1.e-10) THEN
+       IF (denom(ig)>1.e-10) THEN
           zmax(ig)=2.*num(ig)/denom(ig)
           zmax0(ig)=zmax(ig)
@@ -128 +128 @@
 ! def de  zmix continu (profil parabolique des vitesses)
       do ig=1,ngrid
-           if (lmix(ig)>1) THEN
+           IF (lmix(ig)>1) THEN
 ! test 
-              if (((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
+              IF (((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
           *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
           -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
@@ -151 +151 @@
           endif
 !test
-         if ((zmax(ig)-zmix(ig))<=0.) THEN
+         IF ((zmax(ig)-zmix(ig))<=0.) THEN
             zmix(ig)=0.9*zmax(ig)
 !            PRINT*,'pb zmix>zmax'
@@ -160 +160 @@
       do ig=1,ngrid
          do l=1,nlay
-            if (zmix(ig)>=zlev(ig,l).and.  &
+            IF (zmix(ig)>=zlev(ig,l).AND.  &
             zmix(ig)<zlev(ig,l+1)) THEN
               lmix(ig)=l