trunk/dyn3d/comconst.f

module comconst

  use nr_util, only: pi

  implicit none

  real, parameter:: daysec = 86400. ! number of seconds per day

  REAL dtvr ! time step for dynamics, in s
  REAL dtphys ! time step for physics, in s

  real, parameter:: rad = 6371229. ! radius of the Earth (in m)

  real, parameter:: cpp = 1004.70885
  ! specific heat capacity at constant pressure of dry air, in J K-1 kg-1

  real, parameter:: kappa = 0.2857143 ! r / cpp

  real, parameter:: r = cpp * kappa
  ! specific ideal gas constant for dry air, in J K-1 kg-1

  real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2)

  real, parameter:: omeg = 2 * pi / daysec
  ! angular speed of rotation of the Earth (in rad s-1)

  private pi

contains

  SUBROUTINE iniconst

    ! From dyn3d/iniconst.F,v 1.1.1.1 2004/05/19 12:53:05
    ! P. Le Van

    USE conf_gcm_m, ONLY: day_step, iphysiq

    IMPLICIT NONE

    !-----------------------------------------------------------------------

    dtvr = daysec / real(day_step)
    dtphys  = iphysiq * dtvr

    print *, 'dtvr = ', dtvr
    print *, 'dtphys = ', dtphys
    PRINT *, 'cpp = ', cpp
    PRINT *, 'R = ', r
    PRINT *, 'kappa = ', kappa

  END SUBROUTINE iniconst

end module comconst
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