Ignore:
Timestamp:
Aug 2, 2024, 9:58:25 PM (6 months ago)
Author:
abarral
Message:

Put dimensions.h and paramet.h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90

    r5158 r5159  
    22! Copyright (c) 2015, Regents of the University of Colorado
    33! All rights reserved.
    4 !
     4
    55! Redistribution and use in source and binary forms, with or without modification, are
    66! permitted provided that the following conditions are met:
    7 !
     7
    88! 1. Redistributions of source code must retain the above copyright notice, this list of
    99!    conditions and the following disclaimer.
    10 !
     10
    1111! 2. Redistributions in binary form must reproduce the above copyright notice, this list
    1212!    of conditions and the following disclaimer in the documentation and/or other
    1313!    materials provided with the distribution.
    14 !
     14
    1515! 3. Neither the name of the copyright holder nor the names of its contributors may be
    1616!    used to endorse or promote products derived from this software without specific prior
    1717!    written permission.
    18 !
     18
    1919! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
    2020! EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
     
    2626! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
    2727! OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
    28 !
     28
    2929! History:
    3030! July 2006: John Haynes      - Initial version
     
    4949    !   evalues the integral (f ds) between f(index=i1) and f(index=i2)
    5050    !   using the AVINT procedure
    51     !
     51
    5252    ! Inputs:
    5353    !   [f]    functional values
     
    5555    !   [i1]   index of lower limit
    5656    !   [i2]   index of upper limit
    57     !
     57
    5858    ! Returns:
    5959    !   result of path integral
    60     !
     60
    6161    ! Notes:
    6262    !   [s] may be in forward or reverse numerical order
    63     !
     63
    6464    ! Requires:
    6565    !   mrgrnk package
    66     !
     66
    6767    ! Created:
    6868    !   02/02/06  John Haynes (haynes@atmos.colostate.edu)
     
    114114    ! Purpose:
    115115    !   estimate the integral of unevenly spaced data
    116     !
     116
    117117    ! Inputs:
    118118    !   [ftab]     functional values
     
    121121    !   [a]        lower limit of integration
    122122    !   [b]        upper limit of integration
    123     !
     123
    124124    ! Outputs:
    125125    !   [result]   approximate value of integral
    126     !
     126
    127127    ! Reference:
    128128    !   From SLATEC libraries, in public domain
    129     !
     129
    130130    !***********************************************************************
    131     !
     131
    132132    !  AVINT estimates the integral of unevenly spaced data.
    133     !
     133
    134134    !  Discussion:
    135     !
     135
    136136    !    The method uses overlapping parabolas and smoothing.
    137     !
     137
    138138    !  Modified:
    139     !
     139
    140140    !    30 October 2000
    141141    !    4 January 2008, A. Bodas-Salcedo. Error control for XTAB taken out of
    142142    !                    loop to allow vectorization.
    143     !
     143
    144144    !  Reference:
    145     !
     145
    146146    !    Philip Davis and Philip Rabinowitz,
    147147    !    Methods of Numerical Integration,
    148148    !    Blaisdell Publishing, 1967.
    149     !
     149
    150150    !    P E Hennion,
    151151    !    Algorithm 77,
     
    153153    !    Communications of the Association for Computing Machinery,
    154154    !    Volume 5, page 96, 1962.
    155     !
     155
    156156    !  Parameters:
    157     !
     157
    158158    !    Input, real ( kind = 8 ) FTAB(NTAB), the function values,
    159159    !    FTAB(I) = F(XTAB(I)).
    160     !
     160
    161161    !    Input, real ( kind = 8 ) XTAB(NTAB), the abscissas at which the
    162162    !    function values are given.  The XTAB's must be distinct
    163163    !    and in ascending order.
    164     !
     164
    165165    !    Input, integer NTAB, the number of entries in FTAB and
    166166    !    XTAB.  NTAB must be at least 3.
    167     !
     167
    168168    !    Input, real ( kind = 8 ) A, the lower limit of integration.  A should
    169169    !    be, but need not be, near one endpoint of the interval
    170170    !    (X(1), X(NTAB)).
    171     !
     171
    172172    !    Input, real ( kind = 8 ) B, the upper limit of integration.  B should
    173173    !    be, but need not be, near one endpoint of the interval
    174174    !    (X(1), X(NTAB)).
    175     !
     175
    176176    !    Output, real ( kind = 8 ) RESULT, the approximate value of the integral.
    177177    ! ########################################################################## 
     
    328328  ! Purpose:
    329329  !   Returns the gamma function
    330   !
     330
    331331  ! Input:
    332332  !   [x]   value to compute gamma function of
    333   !
     333
    334334  ! Returns:
    335335  !   gamma(x)
    336   !
     336
    337337  ! Coded:
    338338  !   02/02/06  John Haynes (haynes@atmos.colostate.edu)
Note: See TracChangeset for help on using the changeset viewer.