Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.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/yamada_c.F90 (modified) (7 diffs)

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

@@ -54 +54 @@
 
       REAL timestep
-      real plev(klon,klev+1)
-      real play(klon,klev)
-      real ustar(klon)
-      real kmin,qmin,pblhmin(klon),coriol(klon)
+      REAL plev(klon,klev+1)
+      REAL play(klon,klev)
+      REAL ustar(klon)
+      REAL kmin,qmin,pblhmin(klon),coriol(klon)
       REAL zlev(klon,klev+1)
       REAL zlay(klon,klev)
@@ -74 +74 @@
       REAL kn(klon,klev)
       REAL kq(klon,klev)
-      real ff(klon,klev+1),delta(klon,klev+1)
-      real aa(klon,klev+1),aa0,aa1
-      integer iflag_pbl,ngrid
-      integer nlay,nlev
-
-      logical first
-      integer ipas
+      REAL ff(klon,klev+1),delta(klon,klev+1)
+      REAL aa(klon,klev+1),aa0,aa1
+      INTEGER iflag_pbl,ngrid
+      INTEGER nlay,nlev
+
+      LOGICAL first
+      INTEGER ipas
       save first,ipas
 !FH/IM     data first,ipas/.TRUE.,0/
@@ -89 +89 @@
 
 
-      integer ig,k
-
-
-      real ri,zrif,zalpha,zsm,zsn
-      real rif(klon,klev+1),sm(klon,klev+1),alpha(klon,klev)
-
-      real m2(klon,klev+1),dz(klon,klev+1),zq,n2(klon,klev+1)
+      INTEGER ig,k
+
+
+      REAL ri,zrif,zalpha,zsm,zsn
+      REAL rif(klon,klev+1),sm(klon,klev+1),alpha(klon,klev)
+
+      REAL m2(klon,klev+1),dz(klon,klev+1),zq,n2(klon,klev+1)
       REAL, DIMENSION(klon,klev+1) :: km2,kn2,sqrtq
-      real dtetadz(klon,klev+1)
-      real m2cstat,mcstat,kmcstat
-      real l(klon,klev+1)
-      real leff(klon,klev+1)
-      real,allocatable,save :: l0(:)
+      REAL dtetadz(klon,klev+1)
+      REAL m2cstat,mcstat,kmcstat
+      REAL l(klon,klev+1)
+      REAL leff(klon,klev+1)
+      REAL,ALLOCATABLE,save :: l0(:)
 !$OMP THREADPRIVATE(l0)      
-      real sq(klon),sqz(klon),zz(klon,klev+1)
-      integer iter
-
-      real ric,rifc,b1,kap
+      REAL sq(klon),sqz(klon),zz(klon,klev+1)
+      INTEGER iter
+
+      REAL ric,rifc,b1,kap
       save ric,rifc,b1,kap
       data ric,rifc,b1,kap/0.195,0.191,16.6,0.4/
 !$OMP THREADPRIVATE(ric,rifc,b1,kap)
-      real frif,falpha,fsm
-      real fl,zzz,zl0,zq2,zn2
-
-      real rino(klon,klev+1),smyam(klon,klev),styam(klon,klev)
-      real lyam(klon,klev),knyam(klon,klev)
-      real w2yam(klon,klev),t2yam(klon,klev)
+      REAL frif,falpha,fsm
+      REAL fl,zzz,zl0,zq2,zn2
+
+      REAL rino(klon,klev+1),smyam(klon,klev),styam(klon,klev)
+      REAL lyam(klon,klev),knyam(klon,klev)
+      REAL w2yam(klon,klev),t2yam(klon,klev)
       logical,save :: firstcall=.TRUE.
 !$OMP THREADPRIVATE(firstcall)       
@@ -134 +134 @@
 
       okiophys=klon==1
-      if (firstcall) THEN
+      IF (firstcall) THEN
         CALL getin_p('iflag_tke_diff',iflag_tke_diff)
         allocate(l0(klon))
@@ -277 +277 @@
 !        n2(ig,k)=0.
          ri=n2(ig,k)/max(m2(ig,k),1.e-10)
-         if (ri<ric) THEN
+         IF (ri<ric) THEN
             rif(ig,k)=frif(ri)
          else
@@ -300 +300 @@
 
 !   Mise a jour de l0
-      if (iflag_pbl==8.or.iflag_pbl==10) THEN
+      IF (iflag_pbl==8.OR.iflag_pbl==10) THEN
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Iterative computation of l0
@@ -438 +438 @@
 
 ! Computation of turbulent diffusivities
-!  IF (25<=iflag_pbl.and.iflag_pbl<=28) THEN
+!  IF (25<=iflag_pbl.AND.iflag_pbl<=28) THEN
 !    DO k=2,klev
 !       sqrtq(:,k)=sqrt(0.5*(q2(:,k)+q2(:,k-1)))