1 |
module comconst |
2 |
|
3 |
use nr_util, only: pi |
4 |
|
5 |
implicit none |
6 |
|
7 |
real, parameter:: daysec = 86400. ! number of seconds per day |
8 |
|
9 |
REAL dtvr ! time step for dynamics, in s |
10 |
REAL dtphys ! time step for physics, in s |
11 |
|
12 |
real, parameter:: rad = 6371229. ! radius of the Earth (in m) |
13 |
|
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 |
real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2) |
23 |
|
24 |
real, parameter:: omeg = 2 * pi / daysec |
25 |
! angular speed of rotation of the Earth (in rad s-1) |
26 |
|
27 |
private pi |
28 |
|
29 |
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 |
end module comconst |