Ignore:
Timestamp:
Jul 23, 2024, 7:14:34 PM (7 months ago)
Author:
abarral
Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

Location:
LMDZ6/branches/Amaury_dev/libf/grid
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h

    r5099 r5105  
    22! $Header$
    33
    4 c--------------------------------------------------------------
    5          REAL ripx
    6          REAL fx,fxprim,fy,fyprim,ri,rj,bigy
    7 c
    8 c....stretching in x...
    9 c
    10         ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim)
    11         fx  (  ri )= ripx(ri) + transx  +
    12      *         alphax * SIN( ripx(ri)+transx-pxo ) - pi
    13         fxprim(ri) = 2.*pi/REAL(iim)  *
    14      *        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
     4!--------------------------------------------------------------
     5   REAL :: ripx
     6   REAL :: fx,fxprim,fy,fyprim,ri,rj,bigy
     7!
     8!....stretching in x...
     9!
     10  ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim)
     11  fx  (  ri )= ripx(ri) + transx  + &
     12        alphax * SIN( ripx(ri)+transx-pxo ) - pi
     13  fxprim(ri) = 2.*pi/REAL(iim)  * &
     14        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
    1515
    16 c....stretching in y...
    17 c
    18         bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
    19         fy(rj)     =  ( bigy(rj) + transy  +
    20      *        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
    21         fyprim(rj) = ( pi/jjm ) * ( 1.+
    22      *           alphay * COS( bigy(rj)+transy-pyo ) )
     16!....stretching in y...
     17!
     18  bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
     19  fy(rj)     =  ( bigy(rj) + transy  + &
     20        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
     21  fyprim(rj) = ( pi/jjm ) * ( 1.+ &
     22        alphay * COS( bigy(rj)+transy-pyo ) )
    2323
    24 c      fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
    25 c    *  transy ))
    26 c       fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy))
    27 c--------------------------------------------------------------
     24  ! fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
     25! *  transy ))
     26!   fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy))
     27!--------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/grid/fxy_reg.h

    r5099 r5105  
    22! $Header$
    33
    4 c-----------------------------------------------------------------------
    5 c INCLUDE 'fxyprim.h'
    6 c
    7 c    ................................................................
    8 c    ................  Fonctions in line  ...........................
    9 c    ................................................................
    10 c
    11       REAL fy, fx, fxprim, fyprim
    12       REAL ri, rj
    13 c
    14 c
    15       fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
    16       fyprim( rj ) =    pi/REAL(jjm)
     4!-----------------------------------------------------------------------
     5! INCLUDE 'fxyprim.h'
     6!
     7!    ................................................................
     8!    ................  Fonctions in line  ...........................
     9!    ................................................................
     10!
     11REAL :: fy, fx, fxprim, fyprim
     12REAL :: ri, rj
     13!
     14!
     15fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
     16fyprim( rj ) =    pi/REAL(jjm)
    1717
    18 c    fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    19 c    fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
     18! fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
     19! fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    2020
    21       fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
    22 c    fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
    23       fxprim( ri ) = 2.*pi/REAL(iim)
    24 c
    25 c
    26 c    La valeur de pi est passee par le common/const/ou /const2/ .
    27 c    Sinon, il faut la calculer avant d'appeler ces fonctions .
    28 c
    29 c   ----------------------------------------------------------------
    30 c    Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
    31 c   -----------------------------------------------------------------
    32 c
    33 c    .....  ici, on a l'application particuliere suivante   ........
    34 c
    35 c                **************************************
    36 c                **     x = 2. * pi/iim *  X         **
    37 c                **     y =      pi/jjm *  Y         **
    38 c                **************************************
    39 c
    40 c   ..................................................................
    41 c   ..................................................................
    42 c
    43 c
    44 c
    45 c-----------------------------------------------------------------------
     21fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     22! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     23fxprim( ri ) = 2.*pi/REAL(iim)
     24!
     25!
     26!    La valeur de pi est passee par le common/const/ou /const2/ .
     27!    Sinon, il faut la calculer avant d'appeler ces fonctions .
     28!
     29!   ----------------------------------------------------------------
     30! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
     31!   -----------------------------------------------------------------
     32!
     33!    .....  ici, on a l'application particuliere suivante   ........
     34!
     35!            **************************************
     36!            **     x = 2. * pi/iim *  X         **
     37!            **     y =      pi/jjm *  Y         **
     38!            **************************************
     39!
     40!   ..................................................................
     41!   ..................................................................
     42!
     43!
     44!
     45!-----------------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/grid/fxy_sin.h

    r5099 r5105  
    22! $Header$
    33
    4 c-----------------------------------------------------------------------
    5 c INCLUDE 'fxyprim.h'
    6 c
    7 c    ................................................................
    8 c    ................  Fonctions in line  ...........................
    9 c    ................................................................
    10 c
    11       REAL fy, fx, fxprim, fyprim
    12       REAL ri, rj
    13 c
    14 c
    15       fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    16       fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
     4!-----------------------------------------------------------------------
     5! INCLUDE 'fxyprim.h'
     6!
     7!    ................................................................
     8!    ................  Fonctions in line  ...........................
     9!    ................................................................
     10!
     11REAL :: fy, fx, fxprim, fyprim
     12REAL :: ri, rj
     13!
     14!
     15fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
     16fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    1717
    18       fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
    19 c    fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
    20       fxprim( ri ) = 2.*pi/REAL(iim)
    21 c
    22 c
    23 c    La valeur de pi est passee par le common/const/ou /const2/ .
    24 c    Sinon, il faut la calculer avant d'appeler ces fonctions .
    25 c
    26 c   ----------------------------------------------------------------
    27 c    Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
    28 c   -----------------------------------------------------------------
    29 c
    30 c    .....  ici, on a l'application particuliere suivante   ........
    31 c
    32 c                **************************************
    33 c                **     x = 2. * pi/iim *  X         **
    34 c                **     y =      pi/jjm *  Y         **
    35 c                **************************************
    36 c
    37 c   ..................................................................
    38 c   ..................................................................
    39 c
    40 c
    41 c
    42 c-----------------------------------------------------------------------
     18fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     19! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     20fxprim( ri ) = 2.*pi/REAL(iim)
     21!
     22!
     23!    La valeur de pi est passee par le common/const/ou /const2/ .
     24!    Sinon, il faut la calculer avant d'appeler ces fonctions .
     25!
     26!   ----------------------------------------------------------------
     27! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
     28!   -----------------------------------------------------------------
     29!
     30!    .....  ici, on a l'application particuliere suivante   ........
     31!
     32!            **************************************
     33!            **     x = 2. * pi/iim *  X         **
     34!            **     y =      pi/jjm *  Y         **
     35!            **************************************
     36!
     37!   ..................................................................
     38!   ..................................................................
     39!
     40!
     41!
     42!-----------------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/grid/fxyprim.h

    r5099 r5105  
    22! $Header$
    33
    4 c-----------------------------------------------------------------------
    5 c INCLUDE 'fxyprim.h'
    6 c
    7 c    ................................................................
    8 c    ................  Fonctions in line  ...........................
    9 c    ................................................................
    10 c
    11       REAL fy, fx, fxprim, fyprim
    12       REAL ri, rj
    13 c
    14 c
    15       fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
    16       fyprim( rj ) =    pi/REAL(jjm)
     4!-----------------------------------------------------------------------
     5! INCLUDE 'fxyprim.h'
     6!
     7!    ................................................................
     8!    ................  Fonctions in line  ...........................
     9!    ................................................................
     10!
     11REAL :: fy, fx, fxprim, fyprim
     12REAL :: ri, rj
     13!
     14!
     15fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
     16fyprim( rj ) =    pi/REAL(jjm)
    1717
    18 c    fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    19 c    fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
     18! fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
     19! fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    2020
    21       fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
    22 c    fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
    23       fxprim( ri ) = 2.*pi/REAL(iim)
    24 c
    25 c
    26 c    La valeur de pi est passee par le common/const/ou /const2/ .
    27 c    Sinon, il faut la calculer avant d'appeler ces fonctions .
    28 c
    29 c   ----------------------------------------------------------------
    30 c    Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
    31 c   -----------------------------------------------------------------
    32 c
    33 c    .....  ici, on a l'application particuliere suivante   ........
    34 c
    35 c                **************************************
    36 c                **     x = 2. * pi/iim *  X         **
    37 c                **     y =      pi/jjm *  Y         **
    38 c                **************************************
    39 c
    40 c   ..................................................................
    41 c   ..................................................................
    42 c
    43 c
    44 c
    45 c-----------------------------------------------------------------------
     21fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     22! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     23fxprim( ri ) = 2.*pi/REAL(iim)
     24!
     25!
     26!    La valeur de pi est passee par le common/const/ou /const2/ .
     27!    Sinon, il faut la calculer avant d'appeler ces fonctions .
     28!
     29!   ----------------------------------------------------------------
     30! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
     31!   -----------------------------------------------------------------
     32!
     33!    .....  ici, on a l'application particuliere suivante   ........
     34!
     35!            **************************************
     36!            **     x = 2. * pi/iim *  X         **
     37!            **     y =      pi/jjm *  Y         **
     38!            **************************************
     39!
     40!   ..................................................................
     41!   ..................................................................
     42!
     43!
     44!
     45!-----------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.