Changeset 1559 for trunk/NEMO/OPA_SRC/ZDF/zdfini.F90

Timestamp:

2009-07-29T16:03:14+02:00 (15 years ago)

Author:

ctlod

Message:

only cosmetic changes

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfini.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfini.F90

@@ -1 +1 @@
 MODULE zdfini
    !!======================================================================
-   !!              ***  MODULE  zdfini  ***
-   !! Ocean physics : define vertical mixing variables
-   !!=====================================================================
+   !!                      ***  MODULE  zdfini  ***
+   !! Ocean physics :   read vertical mixing namelist and check consistancy
+   !!======================================================================
+   !! History :  8.0  ! 1997-06  (G. Madec)  Original code from inimix
+   !!            1.0  ! 2002-08  (G. Madec)  F90 : free form
+   !!             -   ! 2005-06  (C. Ethe) KPP parameterization
+   !!             -   ! 2009-07  (G. Madec) add avmb, avtb in restart for cen2 advection
+   !!----------------------------------------------------------------------
 
    !!----------------------------------------------------------------------
    !!   zdf_init    : initialization, namelist read, and parameters control
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE par_oce         ! mesh and scale factors
    USE ldftra_oce      ! ocean active tracers: lateral physics
@@ -30 +34 @@
    PRIVATE
 
-   !! *  Routine accessibility
-   PUBLIC zdf_init          ! routine called by opa.F90
+   PUBLIC   zdf_init   ! routine called by opa.F90
+   
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
+   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 
    !! $Id$
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
@@ -47 +51 @@
       !!
       !! ** Method  :   Read namelist namzdf, control logicals 
-      !!
-      !! History :
-      !!        !  97-06  (G. Madec)  Original code from inimix
-      !!   8.5  !  02-08  (G. Madec)  F90 : free form
-      !!   9.0  !  05-06  (C. Ethe) KPP parameterization
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio       ! temporary scalar
@@ -59 +58 @@
       !!----------------------------------------------------------------------
 
-      REWIND( numnam )               ! Read nam_zdf namelist : vertical mixing parameters
+      REWIND( numnam )           !* Read nam_zdf namelist : vertical mixing parameters
       READ  ( numnam, nam_zdf )
 
-      IF(lwp) THEN                   ! Parameter print
+      IF(lwp) THEN               !* Parameter print
          WRITE(numout,*)
          WRITE(numout,*) 'zdf_init: vertical physics'
@@ -81 +80 @@
       ENDIF
 
-      ! Parameter & logicals controls
-      ! -----------------------------
-      ! ... check of vertical mixing scheme on tracers
-      !           ==> will be done in trazdf module
-
-      ! ... check of mixing coefficient
+      !                          !* Parameter & logical controls
+      !                          !  ----------------------------
+      !
+      !                               ! ... check of vertical mixing scheme on tracers
+      !                                              ==> will be done in trazdf module
+      !
+      !                               ! ... check of mixing coefficient
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          vertical mixing option :'
+      IF(lwp) WRITE(numout,*) '   vertical mixing option :'
       ioptio = 0
       IF( lk_zdfcst ) THEN
-         IF(lwp) WRITE(numout,*) '             constant eddy diffusion coefficients'
+         IF(lwp) WRITE(numout,*) '      constant eddy diffusion coefficients'
          ioptio = ioptio+1
       ENDIF
       IF( lk_zdfric ) THEN
-         IF(lwp) WRITE(numout,*) '             Richardson dependent eddy coefficients'
+         IF(lwp) WRITE(numout,*) '      Richardson dependent eddy coefficients'
          ioptio = ioptio+1
       ENDIF
       IF( lk_zdftke_old ) THEN
-         IF(lwp) WRITE(numout,*) '             TKE dependent eddy coefficients'
+         IF(lwp) WRITE(numout,*) '      TKE dependent eddy coefficients'
          ioptio = ioptio+1
       ENDIF
       IF( lk_zdftke ) THEN
-         IF(lwp) WRITE(numout,*) '             TKE dependent eddy coefficients'
+         IF(lwp) WRITE(numout,*) '      TKE dependent eddy coefficients'
          ioptio = ioptio+1
       ENDIF
       IF( lk_zdfkpp ) THEN
-         IF(lwp) WRITE(numout,*) '             KPP dependent eddy coefficients'
+         IF(lwp) WRITE(numout,*) '      KPP dependent eddy coefficients'
          ioptio = ioptio+1
       ENDIF
       IF( ioptio == 0 .OR. ioptio > 1 .AND. .NOT. lk_esopa ) &
-         CALL ctl_stop( ' one and only one vertical diffusion option has to be defined ' )
-
-      ! ... Convection
+         &   CALL ctl_stop( ' one and only one vertical diffusion option has to be defined ' )
+      !
+      !                               ! ... Convection
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          convection :'
+      IF(lwp) WRITE(numout,*) '   convection :'
       ioptio = 0
       IF( ln_zdfnpc ) THEN
-         IF(lwp) WRITE(numout,*) '             use non penetrative convective scheme'
+         IF(lwp) WRITE(numout,*) '      use non penetrative convective scheme'
          ioptio = ioptio+1
       ENDIF
       IF( ln_zdfevd ) THEN
-         IF(lwp) WRITE(numout,*) '             use enhanced vertical dif. scheme'
+         IF(lwp) WRITE(numout,*) '      use enhanced vertical dif. scheme'
          ioptio = ioptio+1
       ENDIF
       IF( lk_zdftke_old .OR. lk_zdftke ) THEN
-         IF(lwp) WRITE(numout,*) '             use the 1.5 turbulent closure'
+         IF(lwp) WRITE(numout,*) '      use the 1.5 turbulent closure'
       ENDIF
       IF( lk_zdfkpp ) THEN
-         IF(lwp) WRITE(numout,*) '             use the KPP closure scheme'
+         IF(lwp) WRITE(numout,*) '      use the KPP closure scheme'
          IF(lk_mpp) THEN
             IF(lwp) WRITE(numout,cform_err)
-            IF(lwp) WRITE(numout,*) '             The KPP scheme is not ready to run in MPI'
+            IF(lwp) WRITE(numout,*) 'The KPP scheme is not ready to run in MPI'
          ENDIF
       ENDIF
-      IF ( ioptio > 1 .AND. .NOT. lk_esopa ) &
-           CALL ctl_stop( ' chose between ln_zdfnpc', '           and ln_zdfevd' )
+      IF ( ioptio > 1 .AND. .NOT. lk_esopa )   CALL ctl_stop( ' chose between ln_zdfnpc and ln_zdfevd' )
       IF( ioptio == 0 .AND. .NOT.( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp ) ) &
          CALL ctl_stop( ' except for TKE sor KPP physics, a convection scheme is', &