Changeset 1992 for LMDZ5/trunk/libf/phylmd/plevel.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (10 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/plevel.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/plevel.F) (1 diff)

LMDZ5/trunk/libf/phylmd/plevel.F90

@@ -1 @@
-!
+
 ! $Header$
-!
-c================================================================
-c================================================================
-      SUBROUTINE plevel(ilon,ilev,lnew,pgcm,pres,Qgcm,Qpres)
-c================================================================
-c================================================================
-      USE netcdf
-      USE dimphy
-      IMPLICIT none
 
-cym#include "dimensions.h"
-cy#include "dimphy.h"
+! ================================================================
+! ================================================================
+SUBROUTINE plevel(ilon, ilev, lnew, pgcm, pres, qgcm, qpres)
+  ! ================================================================
+  ! ================================================================
+  USE netcdf
+  USE dimphy
+  IMPLICIT NONE
 
-c================================================================
-c
-c Interpoler des champs 3-D u, v et g du modele a un niveau de
-c pression donnee (pres)
-c
-c INPUT:  ilon ----- nombre de points
-c         ilev ----- nombre de couches
-c         lnew ----- true si on doit reinitialiser les poids
-c         pgcm ----- pressions modeles
-c         pres ----- pression vers laquelle on interpolle
-c         Qgcm ----- champ GCM
-c         Qpres ---- champ interpolle au niveau pres
-c
-c================================================================
-c
-c   arguments :
-c   -----------
+  ! ym#include "dimensions.h"
+  ! y#include "dimphy.h"
 
-      INTEGER ilon, ilev
-      logical lnew
+  ! ================================================================
 
-      REAL pgcm(ilon,ilev)
-      REAL Qgcm(ilon,ilev)
-      real pres
-      REAL Qpres(ilon)
+  ! Interpoler des champs 3-D u, v et g du modele a un niveau de
+  ! pression donnee (pres)
 
-c   local :
-c   -------
+  ! INPUT:  ilon ----- nombre de points
+  ! ilev ----- nombre de couches
+  ! lnew ----- true si on doit reinitialiser les poids
+  ! pgcm ----- pressions modeles
+  ! pres ----- pression vers laquelle on interpolle
+  ! Qgcm ----- champ GCM
+  ! Qpres ---- champ interpolle au niveau pres
 
-cym      INTEGER lt(klon), lb(klon)
-cym      REAL ptop, pbot, aist(klon), aisb(klon)
+  ! ================================================================
 
-cym      save lt,lb,ptop,pbot,aist,aisb
-      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: lt,lb
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: aist,aisb
-c$OMP THREADPRIVATE(lt,lb,aist,aisb)      
-      REAL,SAVE :: ptop, pbot
-c$OMP THREADPRIVATE(ptop, pbot)      
-      LOGICAL,SAVE :: first = .true.
-c$OMP THREADPRIVATE(first)
-      INTEGER i, k
-c
-      REAL missing_val
-c
-      missing_val=nf90_fill_real
-c
-      if (first) then
-        allocate(lt(klon),lb(klon),aist(klon),aisb(klon))
-        first=.false.
-      endif
-      
-c=====================================================================
-      if (lnew) then
-c   on r�nitialise les r�ndicages et les poids
-c=====================================================================
+  ! arguments :
+  ! -----------
+
+  INTEGER ilon, ilev
+  LOGICAL lnew
+
+  REAL pgcm(ilon, ilev)
+  REAL qgcm(ilon, ilev)
+  REAL pres
+  REAL qpres(ilon)
+
+  ! local :
+  ! -------
+
+  ! ym      INTEGER lt(klon), lb(klon)
+  ! ym      REAL ptop, pbot, aist(klon), aisb(klon)
+
+  ! ym      save lt,lb,ptop,pbot,aist,aisb
+  INTEGER, ALLOCATABLE, SAVE, DIMENSION (:) :: lt, lb
+  REAL, ALLOCATABLE, SAVE, DIMENSION (:) :: aist, aisb
+  !$OMP THREADPRIVATE(lt,lb,aist,aisb)
+  REAL, SAVE :: ptop, pbot
+  !$OMP THREADPRIVATE(ptop, pbot)
+  LOGICAL, SAVE :: first = .TRUE.
+  !$OMP THREADPRIVATE(first)
+  INTEGER i, k
+
+  REAL missing_val
+
+  missing_val = nf90_fill_real
+
+  IF (first) THEN
+    ALLOCATE (lt(klon), lb(klon), aist(klon), aisb(klon))
+    first = .FALSE.
+  END IF
+
+  ! =====================================================================
+  IF (lnew) THEN
+    ! on r�nitialise les r�ndicages et les poids
+    ! =====================================================================
 
 
-c Chercher les 2 couches les plus proches du niveau a obtenir
-c
-c Eventuellement, faire l'extrapolation a partir des deux couches
-c les plus basses ou les deux couches les plus hautes:
-      DO 130 i = 1, klon
-         IF ( ABS(pres-pgcm(i,ilev) ) .LT.
-     .        ABS(pres-pgcm(i,1)) ) THEN
-            lt(i) = ilev     ! 2
-            lb(i) = ilev-1   ! 1
-         ELSE
-            lt(i) = 2
-            lb(i) = 1
-         ENDIF
-  130 CONTINUE
-      DO 150 k = 1, ilev-1
-         DO 140 i = 1, klon
-            pbot = pgcm(i,k)
-            ptop = pgcm(i,k+1)
-            IF (ptop.LE.pres .AND. pbot.GE.pres) THEN
-               lt(i) = k+1
-               lb(i) = k
-            ENDIF
-  140    CONTINUE
-  150 CONTINUE
-c
-c Interpolation lineaire:
-c
+    ! Chercher les 2 couches les plus proches du niveau a obtenir
+
+    ! Eventuellement, faire l'extrapolation a partir des deux couches
+    ! les plus basses ou les deux couches les plus hautes:
+    DO i = 1, klon
+      IF (abs(pres-pgcm(i,ilev))<abs(pres-pgcm(i,1))) THEN
+        lt(i) = ilev ! 2
+        lb(i) = ilev - 1 ! 1
+      ELSE
+        lt(i) = 2
+        lb(i) = 1
+      END IF
+    END DO
+    DO k = 1, ilev - 1
       DO i = 1, klon
-c interpolation en logarithme de pression:
-c
-c ...   Modif . P. Le Van    ( 20/01/98) ....
-c       Modif Fr��ic Hourdin (3/01/02)
+        pbot = pgcm(i, k)
+        ptop = pgcm(i, k+1)
+        IF (ptop<=pres .AND. pbot>=pres) THEN
+          lt(i) = k + 1
+          lb(i) = k
+        END IF
+      END DO
+    END DO
 
-        aist(i) = LOG( pgcm(i,lb(i))/ pres )
-     .       / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i)) )
-        aisb(i) = LOG( pres / pgcm(i,lt(i)) )
-     .       / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i)))
-      enddo
+    ! Interpolation lineaire:
+
+    DO i = 1, klon
+      ! interpolation en logarithme de pression:
+
+      ! ...   Modif . P. Le Van    ( 20/01/98) ....
+      ! Modif Fr��ic Hourdin (3/01/02)
+
+      aist(i) = log(pgcm(i,lb(i))/pres)/log(pgcm(i,lb(i))/pgcm(i,lt(i)))
+      aisb(i) = log(pres/pgcm(i,lt(i)))/log(pgcm(i,lb(i))/pgcm(i,lt(i)))
+    END DO
 
 
-      endif ! lnew
+  END IF ! lnew
 
-c======================================================================
-c    inteprollation
-c======================================================================
+  ! ======================================================================
+  ! inteprollation
+  ! ======================================================================
 
-      do i=1,klon
-         Qpres(i)= Qgcm(i,lb(i))*aisb(i)+Qgcm(i,lt(i))*aist(i)
-      enddo
-c
-c Je mets les vents a zero quand je rencontre une montagne
-      do i = 1, klon
-         if (pgcm(i,1).LT.pres) THEN
-            Qpres(i)=missing_val
-         endif
-      enddo
+  DO i = 1, klon
+    qpres(i) = qgcm(i, lb(i))*aisb(i) + qgcm(i, lt(i))*aist(i)
+  END DO
 
-c
-      RETURN
-      END
+  ! Je mets les vents a zero quand je rencontre une montagne
+  DO i = 1, klon
+    IF (pgcm(i,1)<pres) THEN
+      qpres(i) = missing_val
+    END IF
+  END DO
+
+
+  RETURN
+END SUBROUTINE plevel