1 |
module comconst |
module comconst |
2 |
|
|
3 |
use nr_util, only: pi |
use nr_util, only: twopi |
4 |
|
|
5 |
implicit none |
implicit none |
6 |
|
|
7 |
real, parameter:: daysec = 86400. ! number of seconds per day |
real, parameter:: daysec = 86400. ! number of seconds per day |
8 |
|
|
9 |
REAL dtvr ! time step for dynamics, in s |
REAL, protected:: dtvr ! time step for dynamics, in s |
10 |
REAL dtphys ! time step for physics, in s |
REAL, protected:: dtphys ! time step for physics, in s |
11 |
|
|
12 |
real, parameter:: rad = 6371229. ! radius of the Earth (in m) |
real, parameter:: rad = 6371229. ! radius of the Earth (in m) |
13 |
|
|
21 |
|
|
22 |
real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2) |
real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2) |
23 |
|
|
24 |
real, parameter:: omeg = 2 * pi / daysec |
real, parameter:: omeg = twopi / daysec |
25 |
! angular speed of rotation of the Earth (in rad s-1) |
! angular speed of rotation of the Earth (in rad s-1) |
26 |
|
|
27 |
private pi |
private twopi |
28 |
|
|
29 |
contains |
contains |
30 |
|
|
44 |
|
|
45 |
print *, 'dtvr = ', dtvr |
print *, 'dtvr = ', dtvr |
46 |
print *, 'dtphys = ', dtphys |
print *, 'dtphys = ', dtphys |
|
PRINT *, 'cpp = ', cpp |
|
47 |
PRINT *, 'R = ', r |
PRINT *, 'R = ', r |
|
PRINT *, 'kappa = ', kappa |
|
48 |
|
|
49 |
END SUBROUTINE iniconst |
END SUBROUTINE iniconst |
50 |
|
|