Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_inicp.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (7 weeks 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/cv3_inicp.F90 (modified) (5 diffs)

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

@@ -10 +10 @@
   ! **************************************************************
   USE lmdz_abort_physic, ONLY: abort_physic
+  USE lmdz_yomcst2
+
   IMPLICIT NONE
-  include "YOMCST2.h"
 
   INTEGER iflag_clos
-  CHARACTER (LEN=20) :: modname = 'cv3_inicp'
-  CHARACTER (LEN=80) :: abort_message
+  CHARACTER (LEN = 20) :: modname = 'cv3_inicp'
+  CHARACTER (LEN = 80) :: abort_message
 
   ! --   Mixing probability distribution functions
 
   REAL qcoef1, qcoef2, qff, qfff, qmix, rmix, qmix1, rmix1, qmix2, rmix2, f
-  REAL :: sumcoef,sigma,aire,pdf,mu,df,ff
+  REAL :: sumcoef, sigma, aire, pdf, mu, df, ff
 
-  qcoef1(f) = tanh(f/gammas)
-  qcoef2(f) = (tanh(f/gammas)+gammas*log(cosh((1.-f)/gammas)/cosh(f/gammas)))
-  qff(f) = max(min(f,1.), 0.)
+  qcoef1(f) = tanh(f / gammas)
+  qcoef2(f) = (tanh(f / gammas) + gammas * log(cosh((1. - f) / gammas) / cosh(f / gammas)))
+  qff(f) = max(min(f, 1.), 0.)
   qfff(f) = min(qff(f), scut)
-  qmix1(f) = (tanh((qff(f)-fmax)/gammas)+qcoef1max)/qcoef2max
-  rmix1(f) = (gammas*log(cosh((qff(f)-fmax)/gammas))+qff(f)*qcoef1max)/ &
-    qcoef2max
-  qmix2(f) = -log(1.-qfff(f))/scut
-  rmix2(f) = (qfff(f)+(1.-qff(f))*log(1.-qfff(f)))/scut
-  qmix(f) = qqa1*qmix1(f) + qqa2*qmix2(f)
-  rmix(f) = qqa1*rmix1(f) + qqa2*rmix2(f)
+  qmix1(f) = (tanh((qff(f) - fmax) / gammas) + qcoef1max) / qcoef2max
+  rmix1(f) = (gammas * log(cosh((qff(f) - fmax) / gammas)) + qff(f) * qcoef1max) / &
+          qcoef2max
+  qmix2(f) = -log(1. - qfff(f)) / scut
+  rmix2(f) = (qfff(f) + (1. - qff(f)) * log(1. - qfff(f))) / scut
+  qmix(f) = qqa1 * qmix1(f) + qqa2 * qmix2(f)
+  rmix(f) = qqa1 * rmix1(f) + qqa2 * rmix2(f)
 
   ! cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -76 +77 @@
 
     sumcoef = qqa1 + qqa2
-    qqa1 = qqa1/sumcoef
-    qqa2 = qqa2/sumcoef
+    qqa1 = qqa1 / sumcoef
+    qqa2 = qqa2 / sumcoef
 
     qcoef1max = qcoef1(fmax)
@@ -90 +91 @@
     ! do ff = 0.0 + df, 1.0 - 2.*df, df
     ff = df
-    DO WHILE (ff<=1.0-2.*df)
-      pdf = (qmix(ff+df)-qmix(ff))*(1.-ff)/df
-      aire = aire + (qmix(ff+df)-qmix(ff))*(1.-ff)
-      mu = mu + pdf*ff*df
+    DO WHILE (ff<=1.0 - 2. * df)
+      pdf = (qmix(ff + df) - qmix(ff)) * (1. - ff) / df
+      aire = aire + (qmix(ff + df) - qmix(ff)) * (1. - ff)
+      mu = mu + pdf * ff * df
       ! c              WRITE(*,*) pdf,  Qmix(ff), aire, ff
       ff = ff + df
@@ -100 +101 @@
     ! do ff=0.0+df,1.0 - 2.*df,df
     ff = df
-    DO WHILE (ff<=1.0-2.*df)
-      pdf = (qmix(ff+df)-qmix(ff))*(1.-ff)/df
-      sigma = sigma + pdf*(ff-mu)*(ff-mu)*df
+    DO WHILE (ff<=1.0 - 2. * df)
+      pdf = (qmix(ff + df) - qmix(ff)) * (1. - ff) / df
+      sigma = sigma + pdf * (ff - mu) * (ff - mu) * df
       ff = ff + df
     END DO
     sigma = sqrt(sigma)
 
-    IF (abs(aire-1.0)>0.02) THEN
+    IF (abs(aire - 1.0)>0.02) THEN
       PRINT *, 'WARNING:: AREA OF MIXING PDF IS::', aire
       abort_message = ''
@@ -116 +117 @@
   END IF !  (iflag_mix .gt. 0)
 
-
 END SUBROUTINE cv3_inicp