1 |
guez |
3 |
module comconst |
2 |
|
|
|
3 |
guez |
39 |
use nr_util, only: pi |
4 |
guez |
3 |
|
5 |
|
|
implicit none |
6 |
|
|
|
7 |
|
|
real, parameter:: daysec = 86400. ! number of seconds per day |
8 |
guez |
79 |
|
9 |
|
|
REAL dtvr ! time step for dynamics, in s |
10 |
|
|
REAL dtphys ! time step for physics, in s |
11 |
|
|
|
12 |
guez |
3 |
real, parameter:: rad = 6371229. ! radius of the Earth (in m) |
13 |
guez |
79 |
|
14 |
|
|
real, parameter:: cpp = 1004.70885 |
15 |
|
|
! specific heat capacity at constant pressure of dry air, in J K-1 kg-1 |
16 |
|
|
|
17 |
|
|
real, parameter:: kappa = 0.2857143 ! r / cpp |
18 |
|
|
|
19 |
|
|
real, parameter:: r = cpp * kappa |
20 |
|
|
! specific ideal gas constant for dry air, in J K-1 kg-1 |
21 |
|
|
|
22 |
guez |
3 |
real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2) |
23 |
|
|
|
24 |
guez |
39 |
real, parameter:: omeg = 2 * pi / daysec |
25 |
|
|
! angular speed of rotation of the Earth (in rad s-1) |
26 |
guez |
3 |
|
27 |
guez |
79 |
private pi |
28 |
guez |
3 |
|
29 |
guez |
79 |
contains |
30 |
|
|
|
31 |
|
|
SUBROUTINE iniconst |
32 |
|
|
|
33 |
|
|
! From dyn3d/iniconst.F,v 1.1.1.1 2004/05/19 12:53:05 |
34 |
|
|
! P. Le Van |
35 |
|
|
|
36 |
|
|
USE conf_gcm_m, ONLY: day_step, iphysiq |
37 |
|
|
|
38 |
|
|
IMPLICIT NONE |
39 |
|
|
|
40 |
|
|
!----------------------------------------------------------------------- |
41 |
|
|
|
42 |
|
|
dtvr = daysec / real(day_step) |
43 |
|
|
dtphys = iphysiq * dtvr |
44 |
|
|
|
45 |
|
|
print *, 'dtvr = ', dtvr |
46 |
|
|
print *, 'dtphys = ', dtphys |
47 |
|
|
PRINT *, 'cpp = ', cpp |
48 |
|
|
PRINT *, 'R = ', r |
49 |
|
|
PRINT *, 'kappa = ', kappa |
50 |
|
|
|
51 |
|
|
END SUBROUTINE iniconst |
52 |
|
|
|
53 |
guez |
3 |
end module comconst |