Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dv2.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_dv2.F90 (modified) (6 diffs)

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

@@ -73 +73 @@
          do ig=1,ngrid
             ltherm(ig,k)=(fm(ig,k+1)+detr(ig,k))*ptimestep > 1.e-5*masse(ig,k)
-            if(ltherm(ig,k).and.larga(ig)>0.) then
+            IF(ltherm(ig,k).and.larga(ig)>0.) THEN
                gamma0(ig,k)=masse(ig,k)  &
            *sqrt( 0.5*(fraca(ig,k+1)+fraca(ig,k)) )  &
@@ -90 +90 @@
 
          do ig=1,ngrid
-            if (ltherm(ig,k)) then
+            if (ltherm(ig,k)) THEN
                dua(ig,k)=ua(ig,k-1)-u(ig,k-1)
                dva(ig,k)=va(ig,k-1)-v(ig,k-1)
@@ -116 +116 @@
 !  par le module du vent dans la couche en dessous.
 !  Mais pourquoi donc ???
-               if (ltherm(ig,k)) then
+               if (ltherm(ig,k)) THEN
 !   On choisit une relaxation lineaire.
 !                 gamma(ig,k)=gamma0(ig,k)
@@ -177 +177 @@
 !----------------------
 
-   if(prt_level>=10) then
+   IF(prt_level>=10) THEN
       do k=1,nlay
          do ig=1,ngrid
@@ -187 +187 @@
    endif
 
-     if (nlarga0>0) then
+     if (nlarga0>0) THEN
           PRINT*,'WARNING !!!!!! DANS THERMCELL_DV2 '
           PRINT*,nlarga0,' points pour lesquels laraga=0. dans un thermique'
@@ -193 +193 @@
      endif
 
-      return
+      RETURN
       end
 END MODULE lmdz_thermcell_dv2