--- trunk/libf/dyn3d/comconst.f90	2008/02/27 13:16:39	3
+++ trunk/dyn3d/comconst.f	2014/03/05 14:57:53	82
@@ -1,32 +1,53 @@
 module comconst
 
-  use dimens_m, only: jjm
+  use nr_util, only: pi
 
   implicit none
 
-  INTEGER im, jm, lllm, imp1
-  integer, parameter:: jmp1 = jjm + 1
-  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 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 r
-  real, parameter:: cpp = 1004.70885, kappa = 0.2857143
-  REAL cotot, unsim
+
+  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 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)
+
+  private pi
 
 contains
 
-  subroutine initialize
+  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 *, "Call sequence information: initialize"
-    pi     =  Acos(-1.)
-    omeg   = 2 * pi / daysec
+    print *, 'dtvr = ', dtvr
+    print *, 'dtphys = ', dtphys
+    PRINT *, 'cpp = ', cpp
+    PRINT *, 'R = ', r
+    PRINT *, 'kappa = ', kappa
 
-  end subroutine initialize
+  END SUBROUTINE iniconst
 
 end module comconst