Changeset 5246 for LMDZ6/trunk/libf/grid/fxyprim.h

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/grid/fxyprim.h (modified) (1 diff)

LMDZ6/trunk/libf/grid/fxyprim.h

@@ -2 +2 @@
 ! $Header$
 !
-c-----------------------------------------------------------------------
-c INCLUDE 'fxyprim.h'
-c
-c    ................................................................
-c    ................  Fonctions in line  ...........................
-c    ................................................................
-c
-      REAL  fy, fx, fxprim, fyprim
-      REAL  ri, rj
-c
-c
-      fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
-      fyprim( rj ) =    pi/REAL(jjm)
+!-----------------------------------------------------------------------
+! INCLUDE 'fxyprim.h'
+!
+!    ................................................................
+!    ................  Fonctions in line  ...........................
+!    ................................................................
+!
+REAL :: fy, fx, fxprim, fyprim
+REAL :: ri, rj
+!
+!
+fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
+fyprim( rj ) =    pi/REAL(jjm)
 
-c     fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
-c     fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
+! fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
+! fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
 
-      fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
-c     fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
-      fxprim( ri ) = 2.*pi/REAL(iim)
-c
-c
-c    La valeur de pi est passee par le common/const/ou /const2/ .
-c    Sinon, il faut la calculer avant d'appeler ces fonctions .
-c
-c   ----------------------------------------------------------------
-c     Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
-c   -----------------------------------------------------------------
-c
-c    .....  ici, on a l'application particuliere suivante   ........
-c
-c                **************************************
-c                **     x = 2. * pi/iim *  X         **
-c                **     y =      pi/jjm *  Y         **
-c                **************************************
-c
-c   ..................................................................
-c   ..................................................................
-c
-c
-c
-c-----------------------------------------------------------------------
+fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
+! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
+fxprim( ri ) = 2.*pi/REAL(iim)
+!
+!
+!    La valeur de pi est passee par le common/const/ou /const2/ .
+!    Sinon, il faut la calculer avant d'appeler ces fonctions .
+!
+!   ----------------------------------------------------------------
+! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
+!   -----------------------------------------------------------------
+!
+!    .....  ici, on a l'application particuliere suivante   ........
+!
+!            **************************************
+!            **     x = 2. * pi/iim *  X         **
+!            **     y =      pi/jjm *  Y         **
+!            **************************************
+!
+!   ..................................................................
+!   ..................................................................
+!
+!
+!
+!-----------------------------------------------------------------------