Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (2 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/dyn3d_common/disvert.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90

@@ -3 +3 @@
 SUBROUTINE disvert()
 
-  use ioipsl, only: getin
-  use new_unit_m, only: new_unit
-  use lmdz_assert, only: assert
+  use ioipsl, ONLY: getin
+  use new_unit_m, ONLY: new_unit
+  use lmdz_assert, ONLY: assert
   USE comvert_mod, ONLY: ap, bp, aps, bps, nivsigs, nivsig, dpres, presnivs, &
                          pseudoalt, pa, preff, scaleheight, presinter
@@ -47 +47 @@
   REAL alpha, beta, deltaz
   REAL x
-  character(len=*),parameter :: modname="disvert"
-
-  character(len=24):: vert_sampling
+  CHARACTER(LEN=*),parameter :: modname="disvert"
+
+  CHARACTER(LEN=24):: vert_sampling
   ! (allowed values are "param", "tropo", "strato" and "read")
 
@@ -62 +62 @@
   CALL getin('vert_sampling', vert_sampling)
   WRITE(lunout,*) TRIM(modname)//' vert_sampling = ' // vert_sampling
-  if (llm==39 .and. vert_sampling=="strato") then
+  if (llm==39 .and. vert_sampling=="strato") THEN
      dsigmin=0.3 ! Vieille option par défaut pour CMIP5
   else
@@ -82 +82 @@
      CLOSE(99)
      alpha=deltaz/(llm*scaleheight)
-     write(lunout, *)trim(modname),':scaleheight, alpha, k0, k1, beta', &
+     WRITE(lunout, *)trim(modname),':scaleheight, alpha, k0, k1, beta', &
                                scaleheight, alpha, k0, k1, beta
 
@@ -96 +96 @@
         dzk1=alpha*tanh(l/k0)
         dzk2=alpha*tanh((llm-k1)/k0)*beta**(l-(llm-k1))/log(beta)
-        write(lunout, *)l, sig(l+1), zk, zk-zkm1, dzk1, dzk2
+        WRITE(lunout, *)l, sig(l+1), zk, zk-zkm1, dzk1, dzk2
         zkm1=zk
      enddo
@@ -332 +332 @@
   ENDDO
 
-  write(lunout, *)  trim(modname),': BP '
-  write(lunout, *) bp
-  write(lunout, *)  trim(modname),': AP '
-  write(lunout, *) ap
-
-  write(lunout, *) 'Niveaux de pressions approximatifs aux centres des'
-  write(lunout, *)'couches calcules pour une pression de surface =', preff
-  write(lunout, *) 'et altitudes equivalentes pour une hauteur d echelle de '
-  write(lunout, *) scaleheight,' km'
+  WRITE(lunout, *)  trim(modname),': BP '
+  WRITE(lunout, *) bp
+  WRITE(lunout, *)  trim(modname),': AP '
+  WRITE(lunout, *) ap
+
+  WRITE(lunout, *) 'Niveaux de pressions approximatifs aux centres des'
+  WRITE(lunout, *)'couches calcules pour une pression de surface =', preff
+  WRITE(lunout, *) 'et altitudes equivalentes pour une hauteur d echelle de '
+  WRITE(lunout, *) scaleheight,' km'
   DO l = 1, llm
      dpres(l) = bp(l) - bp(l+1)
@@ -347 +347 @@
      presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff )
      pseudoalt(l) = log(preff/presnivs(l))*scaleheight
-     write(lunout, *)'PRESNIVS(', l, ')=', presnivs(l), ' Z ~ ', &
+     WRITE(lunout, *)'PRESNIVS(', l, ')=', presnivs(l), ' Z ~ ', &
           pseudoalt(l) &
           , ' DZ ~ ', scaleheight*log((ap(l)+bp(l)*preff)/ &
@@ -354 +354 @@
   DO l=1, llmp1
      presinter(l)= ( ap(l)+bp(l)*preff)
-     write(lunout, *)'PRESINTER(', l, ')=', presinter(l)
+     WRITE(lunout, *)'PRESINTER(', l, ')=', presinter(l)
   ENDDO
 
-  write(lunout, *) trim(modname),': PRESNIVS '
-  write(lunout, *) presnivs
+  WRITE(lunout, *) trim(modname),': PRESNIVS '
+  WRITE(lunout, *) presnivs
 
 CONTAINS