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

Timestamp:

Jul 24, 2024, 2:54:37 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(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

@@ -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
@@ -82 +82 @@
       enddo
 
-!     if (iflag_thermals_ed.ge.1) then
-      if (1==0) then
+!     if (iflag_thermals_ed.ge.1) 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)))  &
@@ -134 +134 @@
           -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
           *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig)))))>1e-10)  &
-          then
-
+          THEN
             zmix(ig)=((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
           *((zlev(ig,lmix(ig)))**2-(zlev(ig,lmix(ig)+1))**2)  &
@@ -152 +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'
@@ -162 +161 @@
          do l=1,nlay
             if (zmix(ig)>=zlev(ig,l).and.  &
-            zmix(ig)<zlev(ig,l+1)) then
+            zmix(ig)<zlev(ig,l+1)) THEN
               lmix(ig)=l
              endif