Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.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_noterre.f90 (modified) (17 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90

@@ -30 +30 @@
   REAL :: snorm
   REAL :: alpha,beta,gama,delta,deltaz
-  real :: quoi,quand
+  REAL :: quoi,quand
   REAL :: zsig(llm),sig(llm+1)
   INTEGER :: np,ierr
-  integer :: ierr1,ierr2,ierr3,ierr4
+  INTEGER :: ierr1,ierr2,ierr3,ierr4
   REAL :: x
 
   REAL :: SSUM
   EXTERNAL SSUM
-  real :: newsig
+  REAL :: newsig
   REAL :: dz0,dz1,nhaut,sig1,esig,csig,zz
-  real :: tt,rr,gg, prevz
-  real :: s(llm),dsig(llm)
-
-  integer :: iz
-  real :: z, ps,p
-  character(len=*),parameter :: modname="disvert_noterre"
+  REAL :: tt,rr,gg, prevz
+  REAL :: s(llm),dsig(llm)
+
+  INTEGER :: iz
+  REAL :: z, ps,p
+  CHARACTER(LEN=*),parameter :: modname="disvert_noterre"
 
   !
@@ -55 +55 @@
   hybrid=.TRUE. ! default value for hybrid (ie: use hybrid coordinates)
   CALL getin('hybrid',hybrid)
-  write(lunout,*) trim(modname),': hybrid=',hybrid
+  WRITE(lunout,*) trim(modname),': hybrid=',hybrid
 
   ! Ouverture possible de fichiers typiquement E.T.
@@ -61 +61 @@
      open(99,file="esasig.def",status='old',form='formatted', &
            iostat=ierr2)
-     if(ierr2/=0) then
+     IF(ierr2/=0) THEN
           close(99)
           open(99,file="z2sig.def",status='old',form='formatted', &
@@ -71 +71 @@
   !   ----------------------------------------
 
-  IF(ierr2==0) then
-
+  IF(ierr2==0) THEN
      ! Lecture de esasig.def :
      ! Systeme peu souple, mais qui respecte en theorie
@@ -78 +77 @@
      ! <-> energie cinetique, d'apres la note de Frederic Hourdin...
 
-     write(lunout,*)'*****************************'
-     write(lunout,*)'WARNING reading esasig.def'
-     write(lunout,*)'*****************************'
+     WRITE(lunout,*)'*****************************'
+     WRITE(lunout,*)'WARNING reading esasig.def'
+     WRITE(lunout,*)'*****************************'
      READ(99,*) scaleheight
      READ(99,*) dz0
@@ -126 +125 @@
   !   ----------------------------------------
 
-  ELSE IF(ierr4==0) then
-     write(lunout,*)'****************************'
-     write(lunout,*)'Reading z2sig.def'
-     write(lunout,*)'****************************'
+  ELSE IF(ierr4==0) THEN
+     WRITE(lunout,*)'****************************'
+     WRITE(lunout,*)'Reading z2sig.def'
+     WRITE(lunout,*)'****************************'
 
      READ(99,*) scaleheight
@@ -146 +145 @@
   !-----------------------------------------------------------------------
   ELSE
-     write(lunout,*) 'didn t you forget something ??? '
-     write(lunout,*) 'We need file  z2sig.def ! (OR esasig.def)'
+     WRITE(lunout,*) 'didn t you forget something ??? '
+     WRITE(lunout,*) 'We need file  z2sig.def ! (OR esasig.def)'
      stop
   ENDIF
@@ -170 +169 @@
 
   if (hybrid) then  ! use hybrid coordinates
-     write(lunout,*) "*********************************"
-     write(lunout,*) "Using hybrid vertical coordinates"
-     write(lunout,*)
+     WRITE(lunout,*) "*********************************"
+     WRITE(lunout,*) "Using hybrid vertical coordinates"
+     WRITE(lunout,*)
      ! Coordonnees hybrides avec mod
      DO l = 1, llm
@@ -183 +182 @@
      ap(llmp1) = 0.
   else ! use sigma coordinates
-     write(lunout,*) "********************************"
-     write(lunout,*) "Using sigma vertical coordinates"
-     write(lunout,*)
+     WRITE(lunout,*) "********************************"
+     WRITE(lunout,*) "Using sigma vertical coordinates"
+     WRITE(lunout,*)
      ! Pour ne pas passer en coordonnees hybrides
      DO l = 1, llm
@@ -196 +195 @@
   bp(llmp1) =   0.
 
-  write(lunout,*) trim(modname),': BP '
-  write(lunout,*)  bp
-  write(lunout,*) trim(modname),': AP '
-  write(lunout,*)  ap
+  WRITE(lunout,*) trim(modname),': BP '
+  WRITE(lunout,*)  bp
+  WRITE(lunout,*) trim(modname),': AP '
+  WRITE(lunout,*)  ap
 
   ! Calcul au milieu des couches :
@@ -214 +213 @@
   ENDDO
 
-  if (hybrid) then
+  if (hybrid) THEN
      aps(llm) = aps(llm-1)**2 / aps(llm-2)
      bps(llm) = 0.5*(bp(llm) + bp(llm+1))
@@ -222 +221 @@
   end if
 
-  write(lunout,*) trim(modname),': BPs '
-  write(lunout,*)  bps
-  write(lunout,*) trim(modname),': APs'
-  write(lunout,*)  aps
+  WRITE(lunout,*) trim(modname),': BPs '
+  WRITE(lunout,*)  bps
+  WRITE(lunout,*) trim(modname),': APs'
+  WRITE(lunout,*)  aps
 
   DO l = 1, llm
@@ -232 +231 @@
   ENDDO
 
-  write(lunout,*)trim(modname),' : PRESNIVS'
-  write(lunout,*)presnivs
-  write(lunout,*)'Pseudo altitude of Presnivs : (for a scale ', &
+  WRITE(lunout,*)trim(modname),' : PRESNIVS'
+  WRITE(lunout,*)presnivs
+  WRITE(lunout,*)'Pseudo altitude of Presnivs : (for a scale ', &
         'height of ',scaleheight,' km)'
-  write(lunout,*)pseudoalt
+  WRITE(lunout,*)pseudoalt
 
   ! --------------------------------------------------
@@ -252 +251 @@
   !   do l=2,llm
   ! approximation of scale height for Venus
-  !      if (zsig(l-1).le.55.) then
+  !      if (zsig(l-1).le.55.) THEN
   !         scaleheight = 15.5 - zsig(l-1)/55.*10.
   !      else
@@ -260 +259 @@
   !    .    log((aps(l) + bps(l)*ps)/(aps(l-1) + bps(l-1)*ps))
   !   END DO
-  !   write(53,'(I3,50F10.5)') iz, zsig
+  !   WRITE(53,'(I3,50F10.5)') iz, zsig
   !  END DO
   !  close(53)
@@ -296 +295 @@
 
   IMPLICIT NONE
-  real :: x1, x2, sig,pa,preff, newsig, F
-  integer :: j
+  REAL :: x1, x2, sig,pa,preff, newsig, F
+  INTEGER :: j
 
   newsig = sig
   x1=0
   x2=1
-  if (sig>=1) then
+  if (sig>=1) THEN
         newsig= sig
-  else if (sig*preff/pa>=0.25) then
+  else if (sig*preff/pa>=0.25) THEN
     DO J=1,9999  ! nombre d''iteration max
       F=((1 -pa/preff)*exp(1-1./newsig**2)+(pa/preff)*newsig)/sig
-      ! write(0,*) J, ' newsig =', newsig, ' F= ', F
-      if (F>1) then
+      ! WRITE(0,*) J, ' newsig =', newsig, ' F= ', F
+      if (F>1) THEN
           X2 = newsig
           newsig=(X1+newsig)*0.5
@@ -319 +318 @@
       IF(abs(10.*log(F))<1.E-5) goto 999
     END DO
-   else   !    if (sig*preff/pa.le.0.25) then
+   else   !    if (sig*preff/pa.le.0.25) THEN
          newsig= sig*preff/pa
    end if