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

Timestamp:

2005-09-02T17:06:49+02:00 (19 years ago)

Author:

opalod

Message:

nemo_v1_update_002 : CT : Integration of the KPP turbulent closure scheme

File:

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

trunk/NEMO/OPA_SRC/ZDF/zdfini.F90

@@ -10 +10 @@
    !! * Modules used
    USE par_oce         ! mesh and scale factors
-   USE ldftra_oce      ! ???
-   USE ldfdyn_oce      ! ???
+   USE ldftra_oce      ! ocean active tracers: lateral physics
+   USE ldfdyn_oce      ! ocean dynamics lateral physics
    USE zdf_oce         ! TKE vertical mixing          
-   USE zdftke          ! TKE vertical mixing          
+   USE lib_mpp         ! distribued memory computing
+   USE zdftke          ! TKE vertical mixing  
+   USE zdfkpp          ! KPP vertical mixing          
    USE zdfddm          ! double diffusion mixing      
    USE zdfevd          ! enhanced vertical diffusion  
    USE zdfric          ! Richardson vertical mixing   
    USE tranpc          ! convection: non penetrative adjustment
-   USE ldfslp          ! ???
+   USE ldfslp          ! iso-neutral slopes
 
    USE in_out_manager  ! I/O manager
@@ -46 +48 @@
       !!        !  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
       !!----------------------------------------------------------------------
       !! * Local declarations
@@ -84 +87 @@
       ! ------------------------
       ! ... vertical mixing
-      ! time stepping scheme (N.B. TKE scheme => force the use of implicit scheme)
-      IF( ln_zdfexp .AND. .NOT.lk_zdftke ) THEN  
+      ! time stepping scheme (N.B. TKE or KPP schemes => force the use of implicit scheme)
+      IF( ( ln_zdfexp .AND. .NOT.lk_zdftke ) .OR. ( ln_zdfexp .AND. .NOT.lk_zdfkpp ) ) THEN  
          l_trazdf_exp = .TRUE.           ! use explicit scheme
          l_trazdf_imp = .FALSE.
@@ -135 +138 @@
          ioptio = ioptio+1
       ENDIF
+      IF( lk_zdfkpp ) THEN
+         IF(lwp) WRITE(numout,*) '             KPP dependent eddy coef.'
+         ioptio = ioptio+1
+      ENDIF
       IF( ioptio == 0 .OR. ioptio > 1 .AND. .NOT. lk_esopa ) THEN
           IF(lwp) WRITE(numout,cform_err)
@@ -156 +163 @@
          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(lk_mpp) THEN
+            IF(lwp) WRITE(numout,cform_err)
+            IF(lwp) WRITE(numout,*) '             The KPP scheme is not ready to run in MPI'
+         ENDIF
+      ENDIF
       IF ( ioptio > 1 .AND. .NOT. lk_esopa ) THEN
           IF(lwp) WRITE(numout,cform_err)