libsim Versione 7.2.6
phys_const.f90
1! Copyright (C) 2010 ARPA-SIM <urpsim@smr.arpa.emr.it>
2! authors:
3! Davide Cesari <dcesari@arpa.emr.it>
4! Paolo Patruno <ppatruno@arpa.emr.it>
5
6! This program is free software; you can redistribute it and/or
7! modify it under the terms of the GNU General Public License as
8! published by the Free Software Foundation; either version 2 of
9! the License, or (at your option) any later version.
10
11! This program is distributed in the hope that it will be useful,
12! but WITHOUT ANY WARRANTY; without even the implied warranty of
13! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14! GNU General Public License for more details.
15
16! You should have received a copy of the GNU General Public License
17! along with this program. If not, see <http://www.gnu.org/licenses/>.
18!> \brief Costanti fisiche (SINGOLA PRECISIONE).
19!!
20!! Questo modulo definisce delle costanti fisiche di uso comune
21!! in meteorologia e discipline affini. Attualmente esse sono definite
22!! solo in singola precisione, sarebbe opportuno trovare un metodo
23!! elegante per definirle semiautomaticamente anche in doppia precisione
24!! affinché possano essere usate, ad esempio, dalle routine di
25!! conversione geo-UTM della classe geo_coord_class.
26!! \ingroup base
28IMPLICIT NONE
29
30! pi = 4.*ATAN(1.)
31REAL,PARAMETER :: pi =3.141593 !< pi greco (\f$\pi\f$)
32REAL,PARAMETER :: rearth = 6370997. !< raggio medio della Terra (\f$m\f$)
33REAL,PARAMETER :: gearth = 9.80665 !< accelerazione gravitazionale media alla superficie terrestre (\f$m/s^2\f$)
34REAL,PARAMETER :: omearth = 2.*pi/86164. !< velocità angolare terrestre (\f$s^{-1}\f$)
35REAL,PARAMETER :: degrad = pi/180. !< fattore di conversione da gradi a radianti
36REAL,PARAMETER :: raddeg = 180./pi !< fattore di conversione da radianti a gradi
37REAL,PARAMETER :: t0c = 273.15 !< 0 gradi celsius in kelvin
38REAL,PARAMETER :: mwater = 18.0153 !< molar mass of water (mol)
39REAL,PARAMETER :: rgas = 8314.472 !< universal gas constant R
40REAL,PARAMETER :: mdry = 28.9644 !< molar mass of dry air
41! rd= 287.05831986852826 ma si trova rd=287.05
42REAL,PARAMETER :: rd = rgas/mdry !< costante dei gas per l'aria secca
43! rv= 461.522816717 ma in vecchi testi 8314.3/18.016=461.495337478
44REAL,PARAMETER :: rv = rgas/mwater !< costante dei gas per il vapore acqueo (\f$JK^{-1}Kg^{-1}\f$)
45REAL,PARAMETER :: eps0 = rd/rv !< \f$rd/rv\f$
46REAL,PARAMETER :: epsy = rv/rd-1. !< \f$rv/rd -1\f$
47REAL,PARAMETER :: rcp = 2./7. !< rapporto tra calore specifico...?
48REAL,PARAMETER :: cpd = rd/rcp !< calore specifico dell'aria secca a pressione costante
49REAL,PARAMETER :: cvd = cpd-rd !< calore specifico dell'aria secca a volume costante
50REAL,PARAMETER :: lvw = 2.5e+6 !< calore latente...?
51REAL,PARAMETER :: convff = 1.94 !< fattore di conversione da m/s a KT
52! c1 era parameter(c1=0.378)
53REAL,PARAMETER :: c1 = 1.- mwater/mdry !< 1. - ( mwater/mdry )
54! coefficents for Tetens formula, from Murray, 1967, J.Appl.Meteor.
55REAL,PARAMETER :: tetens_aw = 17.2694 !< a coefficient over water
56REAL,PARAMETER :: tetens_bw = 35.86 !< b coefficient over water (K)
57REAL,PARAMETER :: tetens_ai = 21.8746 !< a coefficient over ice
58REAL,PARAMETER :: tetens_bi = 7.66 !< b coefficient over ice (K)
59REAL,PARAMETER :: tetens_e0_hpa = 6.1078 !< saturated water vapour pressure at 0C in hPa
60
61END MODULE phys_const
62
63!> \brief Costanti fisiche (DOUBLEPRECISION).
64!!
65!! Questo modulo definisce delle costanti fisiche di uso comune
66!! in meteorologia e discipline affini. Attualmente esse sono definite
67!! solo in singola precisione, sarebbe opportuno trovare un metodo
68!! elegante per definirle semiautomaticamente anche in doppia precisione
69!! affinché possano essere usate, ad esempio, dalle routine di
70!! conversione geo-UTM della classe geo_coord_class.
71!! \ingroup base
73IMPLICIT NONE
74
75! pi = 4.*ATAN(1.)
76DOUBLEPRECISION,PARAMETER :: pi = 3.1415926535897932d0 !< pi greco (\f$\pi\f$)
77
78DOUBLEPRECISION,PARAMETER :: rearth = 6370997.0d0 !< raggio medio della Terra (\f$m\f$)
79DOUBLEPRECISION,PARAMETER :: gearth = 9.80665d0 !< accelerazione gravitazionale media alla superficie terrestre (\f$m/s^2\f$)
80DOUBLEPRECISION,PARAMETER :: omearth = 2.0d0*pi/86164.0d0 !< velocità angolare terrestre (\f$s^{-1}\f$)
81DOUBLEPRECISION,PARAMETER :: degrad = pi/180.0d0 !< fattore di conversione da gradi a radianti
82DOUBLEPRECISION,PARAMETER :: raddeg = 180.0d0/pi !< fattore di conversione da radianti a gradi
83DOUBLEPRECISION,PARAMETER :: t0c = 273.15d0 !< 0 gradi celsius in kelvin
84DOUBLEPRECISION,PARAMETER :: mwater = 18.0153d0 !< molar mass of water (mol)
85DOUBLEPRECISION,PARAMETER :: rgas = 8314.472d0 !< universal gas constant R
86DOUBLEPRECISION,PARAMETER :: mdry = 28.9644d0 !< molar mass of dry air
87! rd= 287.05831986852826 ma si trova rd=287.05
88DOUBLEPRECISION,PARAMETER :: rd = rgas/mdry !< costante dei gas per l'aria secca
89! rv= 461.522816717 ma in vecchi testi 8314.3/18.016=461.495337478
90DOUBLEPRECISION,PARAMETER :: rv = rgas/mwater !< costante dei gas per il vapore acqueo (\f$JK^{-1}Kg^{-1}\f$)
91DOUBLEPRECISION,PARAMETER :: eps0 = rd/rv !< \f$rd/rv\f$
92DOUBLEPRECISION,PARAMETER :: epsy = rv/rd-1.0d0 !< \f$rv/rd -1\f$
93DOUBLEPRECISION,PARAMETER :: rcp = 2.0d0/7.0d0 !< rapporto tra calore specifico...?
94DOUBLEPRECISION,PARAMETER :: cpd = rd/rcp !< calore specifico dell'aria secca a pressione costante
95DOUBLEPRECISION,PARAMETER :: cvd = cpd-rd !< calore specifico dell'aria secca a volume costante
96DOUBLEPRECISION,PARAMETER :: lvw = 2.5d+6 !< calore latente...?
97DOUBLEPRECISION,PARAMETER :: convff = 1.94d0 !< fattore di conversione da KT a m/s
98! c1 era parameter(c1=0.378)
99DOUBLEPRECISION,PARAMETER :: c1 = 1.0d0- mwater/mdry !< 1. - ( mwater/mdry )
100! coefficents for Tetens formula, from Murray, 1967, J.Appl.Meteor.
101DOUBLEPRECISION,PARAMETER :: tetens_aw = 17.2694d0 !< a coefficient over water
102DOUBLEPRECISION,PARAMETER :: tetens_bw = 35.86d0 !< b coefficient over water (K)
103DOUBLEPRECISION,PARAMETER :: tetens_ai = 21.8746d0 !< a coefficient over ice
104DOUBLEPRECISION,PARAMETER :: tetens_bi = 7.66d0 !< b coefficient over ice (K)
105DOUBLEPRECISION,parameter :: tetens_e0_hpa = 6.1078d0 !< saturated water vapour pressure at 0C in hPa
106
108
Costanti fisiche (DOUBLEPRECISION).
Costanti fisiche (SINGOLA PRECISIONE).

Generated with Doxygen.