Changeset 1537 for trunk/NEMO/OPA_SRC/ZDF/zdftke_old.F90

Timestamp:

2009-07-24T16:22:23+02:00 (15 years ago)

Author:

ctlod

Message:

ensure the restartability of the 2nd order advection scheme,see ticket: 489

File:

• trunk/NEMO/OPA_SRC/ZDF/zdftke_old.F90 (modified) (9 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdftke_old.F90

@@ -65 +65 @@
 #endif
 
-   !                                       !!! ** Namelist  namtke  **
+   !                                       !!! ** Namelist  nam_tke  **
    LOGICAL  ::   ln_rstke = .FALSE.         ! =T restart with tke from a run without tke
    LOGICAL  ::   ln_mxl0  = .FALSE.         ! mixing length scale surface value as function of wind stress or not
@@ -73 +73 @@
    INTEGER  ::   nn_pdl   =  1              ! Prandtl number or not (ratio avt/avm) (=0/1)
    INTEGER  ::   nn_ave   =  1              ! horizontal average or not on avt, avmu, avmv (=0/1)
-   INTEGER  ::   nn_avb   =  0              ! constant or profile background on avt (=0/1)
    REAL(wp) ::   rn_ediff = 0.1_wp          ! coefficient for avt: avt=rn_ediff*mxl*sqrt(e)
    REAL(wp) ::   rn_ediss = 0.7_wp          ! coefficient of the Kolmogoroff dissipation 
@@ -81 +80 @@
    REAL(wp) ::   rn_emin0 = 1.e-4_wp        ! surface minimum value of tke (m2/s2)
    REAL(wp) ::   rn_cri   = 2._wp / 9._wp   ! critic Richardson number
-   INTEGER  ::   nn_havtb = 1               ! horizontal shape or not for avtb (=0/1)
    INTEGER  ::   nn_etau  = 0               ! type of depth penetration of surface tke (=0/1/2)
    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1/2)
@@ -88 +86 @@
    REAL(wp) ::   rn_efr   = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean
    REAL(wp) ::   rn_lc    = 0.15_wp         ! coef to compute vertical velocity of Langmuir cells
-
-   REAL(wp), DIMENSION (jpi,jpj) ::   avtb_2d   ! set in tke_init, for other modif than ice
 
    !! * Substitutions
@@ -714 +710 @@
       !!      viscosity when using a tke turbulent closure scheme
       !!
-      !! ** Method  :   Read the namtke namelist and check the parameters
+      !! ** Method  :   Read the nam_tke namelist and check the parameters
       !!      called at the first timestep (nit000)
       !!
-      !! ** input   :   Namlist namtke
+      !! ** input   :   Namlist nam_tke
       !!
       !! ** Action  :   Increase by 1 the nstop flag is setting problem encounter
@@ -731 +727 @@
 # endif
       !!
-      NAMELIST/namtke/ ln_rstke, rn_ediff, rn_ediss, rn_ebb  , rn_efave, rn_emin,   &
-         &             rn_emin0, rn_cri  , nn_itke , nn_mxl  , nn_pdl  , nn_ave ,   &
-         &             nn_avb  , ln_mxl0 , rn_lmin , rn_lmin0, nn_havtb, nn_etau,   &
-         &             nn_htau , rn_efr  , ln_lc   , rn_lc 
+      NAMELIST/nam_tke/ ln_rstke, rn_ediff, rn_ediss, rn_ebb  , rn_efave, rn_emin,   &
+         &              rn_emin0, rn_cri  , nn_itke , nn_mxl  , nn_pdl  , nn_ave ,   &
+         &              ln_mxl0 , rn_lmin , rn_lmin0, nn_etau,   &
+         &              nn_htau , rn_efr  , ln_lc   , rn_lc 
       !!----------------------------------------------------------------------
 
-      ! Read Namelist namtke : Turbulente Kinetic Energy
+      ! Read Namelist nam_tke : Turbulente Kinetic Energy
       ! --------------------
       REWIND ( numnam )
-      READ   ( numnam, namtke )
+      READ   ( numnam, nam_tke )
 
       ! Compute boost associated with the Richardson critic
@@ -756 +752 @@
          WRITE(numout,*) 'zdf_tke_init : tke turbulent closure scheme (old scheme)'
          WRITE(numout,*) '~~~~~~~~~~~~'
-         WRITE(numout,*) '          Namelist namtke : set tke mixing parameters'
+         WRITE(numout,*) '          Namelist nam_tke : set tke mixing parameters'
          WRITE(numout,*) '             restart with tke from no tke              ln_rstke = ', ln_rstke
          WRITE(numout,*) '             coef. to compute avt                      rn_ediff = ', rn_ediff
@@ -770 +766 @@
          WRITE(numout,*) '             critic Richardson nb                      rn_cri   = ', rn_cri
          WRITE(numout,*) '                and its associated coeff.              eboost   = ', eboost
-         WRITE(numout,*) '             constant background or profile            nn_avb   = ', nn_avb
          WRITE(numout,*) '             surface mixing length = F(stress) or not  ln_mxl0  = ', ln_mxl0
          WRITE(numout,*) '             surface  mixing length minimum value      rn_lmin0 = ', rn_lmin0
          WRITE(numout,*) '             interior mixing length minimum value      rn_lmin0 = ', rn_lmin0
-         WRITE(numout,*) '             horizontal variation for avtb             nn_havtb = ', nn_havtb
          WRITE(numout,*) '             test param. to add tke induced by wind    nn_etau  = ', nn_etau
          WRITE(numout,*) '             flag for computation of exp. tke profile  nn_htau  = ', nn_htau
@@ -852 +846 @@
          !
       END SELECT
-
-
-      ! Background eddy viscosity and diffusivity profil
-      ! ------------------------------------------------
-      IF( nn_avb == 0 ) THEN      ! Define avmb, avtb from namelist parameter
-         avmb(:) = avm0
-         avtb(:) = avt0
-      ELSE                      ! Background profile of avt (fit a theoretical/observational profile (Krauss 1990) 
-         avmb(:) = avm0
-         avtb(:) = avt0 + ( 3.0e-4 - 2 * avt0 ) * 1.0e-4 * gdepw_0(:)   ! m2/s
-         IF(ln_sco .AND. lwp)   CALL ctl_warn( '          avtb profile not valid in sco' )
-      ENDIF
-      !
-      !                         ! 2D shape of the avtb
-      avtb_2d(:,:) = 1.e0             ! uniform 
-      !
-      IF( nn_havtb == 1 ) THEN      ! decrease avtb in the equatorial band
-           !  -15S -5S : linear decrease from avt0 to avt0/10.
-           !  -5S  +5N : cst value avt0/10.
-           !   5N  15N : linear increase from avt0/10, to avt0
-           WHERE(-15. <= gphit .AND. gphit < -5 )   avtb_2d = (1.  - 0.09 * (gphit + 15.))
-           WHERE( -5. <= gphit .AND. gphit <  5 )   avtb_2d =  0.1
-           WHERE(  5. <= gphit .AND. gphit < 15 )   avtb_2d = (0.1 + 0.09 * (gphit -  5.))
-      ENDIF
 
       ! Initialization of vertical eddy coef. to the background value