11 |
real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1 |
real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1 |
12 |
|
|
13 |
! A1.1 Astronomical constants |
! A1.1 Astronomical constants |
14 |
REAL RSIYEA, RSIDAY, ROMEGA |
|
15 |
|
REAL ROMEGA |
16 |
real, parameter:: RDAY = 86400. |
real, parameter:: RDAY = 86400. |
17 |
real, parameter:: REA = 149597870000. |
|
18 |
real, parameter:: REPSM = 0.409093 |
REAL, parameter:: RSIYEA = 365.25 * RDAY * 2. * PI / 6.283076 |
19 |
|
! sideral year, in s |
20 |
|
|
21 |
|
REAL, parameter:: RSIDAY = RDAY / (1. + RDAY / RSIYEA) ! sideral day, in s |
22 |
|
|
23 |
! A1.2 Geoide |
! A1.2 Geoide |
24 |
real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2 |
real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2 |
81 |
|
|
82 |
! 2. DEFINE ASTRONOMICAL CONSTANTS |
! 2. DEFINE ASTRONOMICAL CONSTANTS |
83 |
|
|
|
RSIYEA = 365.25*RDAY*2.*PI/6.283076 |
|
|
RSIDAY = RDAY/(1.+RDAY/RSIYEA) |
|
84 |
ROMEGA = 2.*PI/RSIDAY |
ROMEGA = 2.*PI/RSIDAY |
85 |
|
|
86 |
print *, 'Astronomical constants ' |
print *, 'Astronomical constants ' |
|
print '('' day = '', E13.7, '' s'')', RDAY |
|
|
print '('' half g. axis = '', E13.7, '' m'')', REA |
|
|
print '('' mean anomaly = '', E13.7, '' -'')', REPSM |
|
|
print '('' sideral year = '', E13.7, '' s'')', RSIYEA |
|
|
print '('' sideral day = '', E13.7, '' s'')', RSIDAY |
|
87 |
print '('' omega = '', E13.7, '' s-1'')', ROMEGA |
print '('' omega = '', E13.7, '' s-1'')', ROMEGA |
88 |
|
|
89 |
! 3. DEFINE GEOIDE. |
! 3. DEFINE GEOIDE. |