Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/geopot.f90

Timestamp:

Oct 21, 2024, 2:58:45 PM (23 hours ago)

Author:

abarral

Message:

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)

File:

• LMDZ6/trunk/libf/dyn3d_common/geopot.f90 (moved) (moved from LMDZ6/trunk/libf/dyn3d_common/geopot.F) (1 diff)

LMDZ6/trunk/libf/dyn3d_common/geopot.f90

@@ -2 +2 @@
 ! $Header$
 !
-      SUBROUTINE geopot (ngrid, teta, pk, pks, phis, phi )
-      IMPLICIT NONE
+SUBROUTINE geopot (ngrid, teta, pk, pks, phis, phi )
+  IMPLICIT NONE
 
-c=======================================================================
-c
-c   Auteur:  P. Le Van
-c   -------
-c
-c   Objet:
-c   ------
-c
-c    *******************************************************************
-c    ....   calcul du geopotentiel aux milieux des couches    .....
-c    *******************************************************************
-c
-c     ....   l'integration se fait de bas en haut  ....
-c
-c     .. ngrid,teta,pk,pks,phis sont des argum. d'entree pour le s-pg ..
-c              phi               est un  argum. de sortie pour le s-pg .
-c
-c=======================================================================
-c-----------------------------------------------------------------------
-c   Declarations:
-c   -------------
+  !=======================================================================
+  !
+  !   Auteur:  P. Le Van
+  !   -------
+  !
+  !   Objet:
+  !   ------
+  !
+  !    *******************************************************************
+  !    ....   calcul du geopotentiel aux milieux des couches    .....
+  !    *******************************************************************
+  !
+  ! ....   l'integration se fait de bas en haut  ....
+  !
+  ! .. ngrid,teta,pk,pks,phis sont des argum. d'entree pour le s-pg ..
+  !          phi               est un  argum. de sortie pour le s-pg .
+  !
+  !=======================================================================
+  !-----------------------------------------------------------------------
+  !   Declarations:
+  !   -------------
 
-      include "dimensions.h"
-      include "paramet.h"
+  include "dimensions.h"
+  include "paramet.h"
 
-c   Arguments:
-c   ----------
+  !   Arguments:
+  !   ----------
 
-      INTEGER ngrid
-      REAL teta(ngrid,llm),pks(ngrid),phis(ngrid),pk(ngrid,llm) ,
-     *       phi(ngrid,llm)
+  INTEGER :: ngrid
+  REAL :: teta(ngrid,llm),pks(ngrid),phis(ngrid),pk(ngrid,llm) , &
+        phi(ngrid,llm)
 
 
-c   Local:
-c   ------
+  !   Local:
+  !   ------
 
-      INTEGER  l, ij
+  INTEGER :: l, ij
 
 
-c-----------------------------------------------------------------------
-c     calcul de phi au niveau 1 pres du sol  .....
+  !-----------------------------------------------------------------------
+  ! calcul de phi au niveau 1 pres du sol  .....
 
-      DO   1  ij  = 1, ngrid
-      phi( ij,1 ) = phis( ij ) + teta(ij,1) * ( pks(ij) - pk(ij,1) )
-   1  CONTINUE
+  DO  ij  = 1, ngrid
+  phi( ij,1 ) = phis( ij ) + teta(ij,1) * ( pks(ij) - pk(ij,1) )
+  END DO
 
-c     calcul de phi aux niveaux superieurs  .......
+  ! calcul de phi aux niveaux superieurs  .......
 
-      DO  l = 2,llm
-        DO  ij    = 1,ngrid
-        phi(ij,l) = phi(ij,l-1) + 0.5 * ( teta(ij,l)  + teta(ij,l-1) ) 
-     *                              *   (  pk(ij,l-1) -  pk(ij,l)    )
-        ENDDO
-      ENDDO
+  DO  l = 2,llm
+    DO  ij    = 1,ngrid
+    phi(ij,l) = phi(ij,l-1) + 0.5 * ( teta(ij,l)  + teta(ij,l-1) ) &
+          *   (  pk(ij,l-1) -  pk(ij,l)    )
+    ENDDO
+  ENDDO
 
-      RETURN
-      END
+  RETURN
+END SUBROUTINE geopot