--- trunk/libf/dyn3d/comconst.f90	2008/02/27 13:16:39	3
+++ trunk/libf/dyn3d/comconst.f90	2013/06/24 15:39:52	70
@@ -1,6 +1,7 @@
 module comconst
 
   use dimens_m, only: jjm
+  use nr_util, only: pi
 
   implicit none
 
@@ -9,24 +10,17 @@
   integer lllmm1, lllmp1, lcl
   REAL dtvr ! time step for dynamics (in s)
   real, parameter:: daysec = 86400. ! number of seconds per day
-  REAL pi, dtphys, dtdiss
+  REAL dtphys ! time step of integration
   real, parameter:: rad = 6371229. ! radius of the Earth (in m)
   real r
-  real, parameter:: cpp = 1004.70885, kappa = 0.2857143
+  real, parameter:: cpp = 1004.70885 ! J K-1 kg-1
+  real, parameter:: kappa = 0.2857143
   REAL cotot, unsim
   real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2)
-  real omeg ! angular speed of rotation of the Earth (in rad s-1)
 
-  private jjm
+  real, parameter:: omeg = 2 * pi / daysec
+  ! angular speed of rotation of the Earth (in rad s-1)
 
-contains
-
-  subroutine initialize
-
-    print *, "Call sequence information: initialize"
-    pi     =  Acos(-1.)
-    omeg   = 2 * pi / daysec
-
-  end subroutine initialize
+  private jjm, pi
 
 end module comconst