Index: LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h	(revision 5099)
+++ LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h	(revision 5105)
@@ -2,26 +2,26 @@
 ! $Header$
 
-c--------------------------------------------------------------
-         REAL ripx
-         REAL fx,fxprim,fy,fyprim,ri,rj,bigy
-c
-c....stretching in x...
-c
-        ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim) 
-        fx  (  ri )= ripx(ri) + transx  +
-     *         alphax * SIN( ripx(ri)+transx-pxo ) - pi
-        fxprim(ri) = 2.*pi/REAL(iim)  *
-     *        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
+!--------------------------------------------------------------
+   REAL :: ripx
+   REAL :: fx,fxprim,fy,fyprim,ri,rj,bigy
+!
+!....stretching in x...
+!
+  ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim)
+  fx  (  ri )= ripx(ri) + transx  + &
+        alphax * SIN( ripx(ri)+transx-pxo ) - pi
+  fxprim(ri) = 2.*pi/REAL(iim)  * &
+        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
 
-c....stretching in y...
-c
-        bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
-        fy(rj)     =  ( bigy(rj) + transy  +
-     *        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
-        fyprim(rj) = ( pi/jjm ) * ( 1.+
-     *           alphay * COS( bigy(rj)+transy-pyo ) )
+!....stretching in y...
+!
+  bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
+  fy(rj)     =  ( bigy(rj) + transy  + &
+        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
+  fyprim(rj) = ( pi/jjm ) * ( 1.+ &
+        alphay * COS( bigy(rj)+transy-pyo ) )
 
-c       fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
-c     *  transy ))
-c       fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy)) 
-c--------------------------------------------------------------
+  ! fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
+! *  transy ))
+!   fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy))
+!--------------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/grid/fxy_reg.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/grid/fxy_reg.h	(revision 5099)
+++ LMDZ6/branches/Amaury_dev/libf/grid/fxy_reg.h	(revision 5105)
@@ -2,44 +2,44 @@
 ! $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         **
+!            **************************************
+!
+!   ..................................................................
+!   ..................................................................
+!
+!
+!
+!-----------------------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/grid/fxy_sin.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/grid/fxy_sin.h	(revision 5099)
+++ LMDZ6/branches/Amaury_dev/libf/grid/fxy_sin.h	(revision 5105)
@@ -2,41 +2,41 @@
 ! $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)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
-      fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
+!-----------------------------------------------------------------------
+! INCLUDE 'fxyprim.h'
+!
+!    ................................................................
+!    ................  Fonctions in line  ...........................
+!    ................................................................
+!
+REAL :: fy, fx, fxprim, fyprim
+REAL :: ri, 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         **
+!            **************************************
+!
+!   ..................................................................
+!   ..................................................................
+!
+!
+!
+!-----------------------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/grid/fxyprim.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/grid/fxyprim.h	(revision 5099)
+++ LMDZ6/branches/Amaury_dev/libf/grid/fxyprim.h	(revision 5105)
@@ -2,44 +2,44 @@
 ! $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         **
+!            **************************************
+!
+!   ..................................................................
+!   ..................................................................
+!
+!
+!
+!-----------------------------------------------------------------------
