File: mean.f90

package info (click to toggle)
flexpart 9.02-23
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 4,944 kB
  • sloc: f90: 14,310; makefile: 28; sh: 18
file content (66 lines) | stat: -rw-r--r-- 3,274 bytes parent folder | download | duplicates (6)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
!**********************************************************************
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann   *
!                                                                     *
! This file is part of FLEXPART.                                      *
!                                                                     *
! FLEXPART is free software: you can redistribute it and/or modify    *
! it under the terms of the GNU General Public License as published by*
! the Free Software Foundation, either version 3 of the License, or   *
! (at your option) any later version.                                 *
!                                                                     *
! FLEXPART is distributed in the hope that it will be useful,         *
! but WITHOUT ANY WARRANTY; without even the implied warranty of      *
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       *
! GNU General Public License for more details.                        *
!                                                                     *
! You should have received a copy of the GNU General Public License   *
! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
!**********************************************************************

subroutine mean(x,xm,xs,number)

  !*****************************************************************************
  !                                                                            *
  !  This subroutine calculates mean and standard deviation of a given element.*
  !                                                                            *
  !      AUTHOR: Andreas Stohl, 25 January 1994                                *
  !                                                                            *
  !*****************************************************************************
  !                                                                            *
  ! Variables:                                                                 *
  ! x(number)           field of input data                                    *
  ! xm                  mean                                                   *
  ! xs                  standard deviation                                     *
  ! number              number of elements of field x                          *
  !                                                                            *
  ! Constants:                                                                 *
  ! eps                 tiny number                                            *
  !                                                                            *
  !*****************************************************************************

  implicit none

  integer :: number,i
  real :: x(number),xm,xs,xl,xq,xaux
  real,parameter :: eps=1.0e-30

  xl=0.
  xq=0.
  do i=1,number
    xl=xl+x(i)
    xq=xq+x(i)*x(i)
  end do

  xm=xl/real(number)

  xaux=xq-xl*xl/real(number)

  if (xaux.lt.eps) then
    xs=0.
  else
    xs=sqrt(xaux/real(number-1))
  endif

end subroutine mean