| 1 |
module comconst |
module comconst |
| 2 |
|
|
|
use dimens_m, only: jjm |
|
| 3 |
use nr_util, only: pi |
use nr_util, only: pi |
| 4 |
|
|
| 5 |
implicit none |
implicit none |
| 6 |
|
|
|
INTEGER im, jm, lllm, imp1 |
|
|
integer, parameter:: jmp1 = jjm + 1 |
|
|
integer lllmm1, lllmp1, lcl |
|
|
REAL dtvr ! time step for dynamics (in s) |
|
| 7 |
real, parameter:: daysec = 86400. ! number of seconds per day |
real, parameter:: daysec = 86400. ! number of seconds per day |
| 8 |
REAL dtphys ! time step of integration |
|
| 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) |
real, parameter:: rad = 6371229. ! radius of the Earth (in m) |
| 13 |
real r |
|
| 14 |
real, parameter:: cpp = 1004.70885 ! J K-1 kg-1 |
real, parameter:: cpp = 1004.70885 |
| 15 |
real, parameter:: kappa = 0.2857143 |
! specific heat capacity at constant pressure of dry air, in J K-1 kg-1 |
| 16 |
REAL cotot, unsim |
|
| 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) |
real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2) |
| 23 |
|
|
| 24 |
real, parameter:: omeg = 2 * pi / daysec |
real, parameter:: omeg = 2 * pi / 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 jjm, pi |
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 |
end module comconst |
Parent Directory
| 
Revision Log
| 
Patch