Index: trunk/NEMO/OPA_SRC/phycst.F90
===================================================================
--- trunk/NEMO/OPA_SRC/phycst.F90	(revision 699)
+++ trunk/NEMO/OPA_SRC/phycst.F90	(revision 703)
@@ -4,5 +4,14 @@
    !!     Definition of of both ocean and ice parameters used in the code
    !!=====================================================================
-   !! * Modules used
+   !! History :        !  90-10  (C. Levy - G. Madec)  Original code
+   !!                  !  91-11  (G. Madec)
+   !!                  !  91-12  (M. Imbard)
+   !!             8.5  !  02-08  (G. Madec, C. Ethe)  F90, add ice constants
+   !!             9.0  !  06-08  (G. Madec)  style 
+   !!----------------------------------------------------------------------
+
+   !!----------------------------------------------------------------------
+   !!   phy_cst  : define and print physical constant and domain parameters
+   !!----------------------------------------------------------------------
    USE par_oce          ! ocean parameters
    USE in_out_manager   ! I/O manager
@@ -11,18 +20,9 @@
    PRIVATE
 
-   !! * Routine accessibility
-   PUBLIC phy_cst          ! routine called by inipar.F90
+   PUBLIC   phy_cst     ! routine called by inipar.F90
 
-   !! * Shared module variables
-   INTEGER, PUBLIC, DIMENSION(12) ::   &  !:
-      nbiss = (/ 31, 29, 31, 30, 31, 30,      &  !: number of days per month
-         &       31, 31, 30, 31, 30, 31 /) ,  &  !  (leap-year)
-      nobis = (/ 31, 28, 31, 30, 31, 30,      &  !: number of days per month
-         &       31, 31, 30, 31, 30, 31 /)       !  (365 days a year)
-   
-   REAL(wp), PUBLIC ::                        &  !:
-      rpi = 3.141592653589793_wp           ,  &  !: pi
-      rad = 3.141592653589793_wp / 180._wp ,  &  !: conversion from degre into radian
-      rsmall = 0.5 * EPSILON( 1. )               !: smallest real computer value
+   REAL(wp), PUBLIC ::   rpi = 3.141592653589793_wp             !: pi
+   REAL(wp), PUBLIC ::   rad = 3.141592653589793_wp / 180._wp   !: conversion from degre into radian
+   REAL(wp), PUBLIC ::   rsmall = 0.5 * EPSILON( 1. )           !: smallest real computer value
    
    REAL(wp), PUBLIC ::          & !:
@@ -61,6 +61,6 @@
       xlic    = 300.33e+6_wp  ,   &  !: volumetric latent heat fusion of ice
       xsn     =   2.8e+6      ,   &  !: latent heat of sublimation of snow
-      rhoic   = 900._wp       ,   &  !: density of sea ice (kg/m3)
-      rhosn   = 330._wp       ,   &  !: density of snow (kg/m3)
+      rhoic   = 900._wp       ,   &  !: volumic mass of sea ice (kg/m3)
+      rhosn   = 330._wp       ,   &  !: volumic mass of snow (kg/m3)
       emic    =   0.97_wp     ,   &  !: emissivity of snow or ice
       sice    =   6.0_wp      ,   &  !: salinity of ice (psu)
@@ -70,9 +70,10 @@
       vkarmn  =   0.4_wp      ,   &  !: von Karman constant
       stefan  =   5.67e-8_wp         !: Stefan-Boltzmann constant 
-      !!----------------------------------------------------------------------
-      !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-      !! $Id$
-      !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
-      !!----------------------------------------------------------------------
+
+   !!----------------------------------------------------------------------
+   !!  OPA 9.0 , LOCEAN-IPSL (2006) 
+   !! $Id$
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
    
 CONTAINS
@@ -83,15 +84,6 @@
       !!
       !! ** Purpose :   Print model parameters and set and print the constants
-      !!
-      !! ** Method  :   no
-      !!
-      !! History :
-      !!        !  90-10  (C. Levy - G. Madec)  Original code
-      !!        !  91-11  (G. Madec)
-      !!        !  91-12  (M. Imbard)
-      !!   8.5  !  02-08  (G. Madec, C. Ethe)  F90, add ice constants 
       !!----------------------------------------------------------------------
-      !! * Local variables
-      CHARACTER (len=64) ::   cform = "(A9, 3(A13, I7) )" 
+      CHARACTER (len=64) ::   cform = "(A12, 3(A13, I7) )" 
       !!----------------------------------------------------------------------
 
@@ -103,22 +95,17 @@
       ! ----------------
       IF(lwp) THEN
-         WRITE(numout,*) '       parameter file'
-         WRITE(numout,*)
+         WRITE(numout,*) '       Domain info'
          WRITE(numout,*) '          dimension of model'
-         WRITE(numout,*) '              Local domain      Global domain       Data domain '
-         WRITE(numout,cform) '         ','   jpi     : ', jpi, '   jpiglo  : ', jpiglo, '   jpidta  : ', jpidta
-         WRITE(numout,cform) '         ','   jpj     : ', jpj, '   jpjglo  : ', jpjglo, '   jpjdta  : ', jpjdta
-         WRITE(numout,cform) '         ','   jpk     : ', jpk, '   jpk     : ', jpk   , '   jpkdta  : ', jpkdta
-         WRITE(numout,*)      '        ','   jpij    : ', jpij
-         WRITE(numout,*)
+         WRITE(numout,*) '                 Local domain      Global domain       Data domain '
+         WRITE(numout,cform) '            ','   jpi     : ', jpi, '   jpiglo  : ', jpiglo, '   jpidta  : ', jpidta
+         WRITE(numout,cform) '            ','   jpj     : ', jpj, '   jpjglo  : ', jpjglo, '   jpjdta  : ', jpjdta
+         WRITE(numout,cform) '            ','   jpk     : ', jpk, '   jpk     : ', jpk   , '   jpkdta  : ', jpkdta
+         WRITE(numout,*)      '           ','   jpij    : ', jpij
          WRITE(numout,*) '          mpp local domain info (mpp)'
          WRITE(numout,*) '             jpni    : ', jpni, '   jpreci  : ', jpreci
          WRITE(numout,*) '             jpnj    : ', jpnj, '   jprecj  : ', jprecj
          WRITE(numout,*) '             jpnij   : ', jpnij
-
-         WRITE(numout,*)
          WRITE(numout,*) '          lateral domain boundary condition type : jperio  = ', jperio
-         WRITE(numout,*) '          domain island (use in rigid-lid case)  : jpisl   = ', jpisl 
-         WRITE(numout,*) '                                                   jpnisl  = ', jpnisl
+         WRITE(numout,*) '          domain island (use in rigid-lid case)  : jpisl   = ', jpisl, '   jpnisl  = ', jpnisl
       ENDIF
 
@@ -126,5 +113,5 @@
       ! ----------------
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '       constants'
+      IF(lwp) WRITE(numout,*) '       Constants'
 
       IF(lwp) WRITE(numout,*)
