trunk/dyn3d/comconst.f

module comconst

  use nr_util, only: twopi

  implicit none

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

  REAL, protected:: dtvr ! time step for dynamics, in s
  REAL, protected:: 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 = twopi / daysec
  ! angular speed of rotation of the Earth (in rad s-1)

  private twopi

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 *, 'R = ', r

  END SUBROUTINE iniconst

end module comconst
1	module comconst
2
3	use nr_util, only: twopi
4
5	implicit none
6
7	real, parameter:: daysec = 86400. ! number of seconds per day
8
9	REAL, protected:: dtvr ! time step for dynamics, in s
10	REAL, protected:: 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 = twopi / daysec
25	! angular speed of rotation of the Earth (in rad s-1)
26
27	private twopi
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 *, 'R = ', r
48
49	END SUBROUTINE iniconst
50
51	end module comconst