source: trunk/src/raerod.f90 @ 28

Last change on this file since 28 was 4, checked in by mlanger, 9 years ago
File size: 4.0 KB
Line 
1!**********************************************************************
2! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
3! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
4! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann   *
5!                                                                     *
6! This file is part of FLEXPART.                                      *
7!                                                                     *
8! FLEXPART is free software: you can redistribute it and/or modify    *
9! it under the terms of the GNU General Public License as published by*
10! the Free Software Foundation, either version 3 of the License, or   *
11! (at your option) any later version.                                 *
12!                                                                     *
13! FLEXPART is distributed in the hope that it will be useful,         *
14! but WITHOUT ANY WARRANTY; without even the implied warranty of      *
15! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       *
16! GNU General Public License for more details.                        *
17!                                                                     *
18! You should have received a copy of the GNU General Public License   *
19! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
20!**********************************************************************
21
22function raerod (l,ust,z0)
23
24  !*****************************************************************************
25  !                                                                            *
26  !     Calculation of the aerodynamical resistance ra from ground up to href  *
27  !                                                                            *
28  !     AUTHOR: Matthias Langer, modified by Andreas Stohl (6 August 1993)     *
29  !                                                                            *
30  !     Literature:                                                            *
31  !     [1]  Hicks/Baldocchi/Meyers/Hosker/Matt (1987), A Preliminary          *
32  !             Multiple Resistance Routine for Deriving Dry Deposition        *
33  !             Velocities from Measured Quantities.                           *
34  !             Water, Air and Soil Pollution 36 (1987), pp.311-330.           *
35  !     [2]  Scire/Yamartino/Carmichael/Chang (1989),                          *
36  !             CALGRID: A Mesoscale Photochemical Grid Model.                 *
37  !             Vol II: User's Guide. (Report No.A049-1, June, 1989)           *
38  !                                                                            *
39  !     Variable list:                                                         *
40  !     L     = Monin-Obukhov-length [m]                                       *
41  !     ust   = friction velocity [m/sec]                                      *
42  !     z0    = surface roughness length [m]                                   *
43  !     href  = reference height [m], for which deposition velocity is         *
44  !             calculated                                                     *
45  !                                                                            *
46  !     Constants:                                                             *
47  !     karman    = von Karman-constant (~0.4)                                 *
48  !     ramin = minimum resistence of ra (1 s/m)                               *
49  !                                                                            *
50  !     Subprograms and functions:                                             *
51  !     function psih (z/L)                                                    *
52  !                                                                            *
53  !*****************************************************************************
54
55  use par_mod
56
57  implicit none
58
59  real :: l,psih,raerod,ust,z0
60
61  raerod=(alog(href/z0)-psih(href,l)+psih(z0,l))/(karman*ust)
62
63end function raerod
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG