Changeset 1537 for trunk/NEMO/OPA_SRC/ZDF/zdfini.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/zdfini.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfini.F90

@@ -22 +22 @@
    USE tranpc          ! convection: non penetrative adjustment
    USE ldfslp          ! iso-neutral slopes
+   USE restart         ! ocean restart
 
    USE in_out_manager  ! I/O manager
+   USE iom             ! IOM library
 
    IMPLICIT NONE
@@ -51 +53 @@
       !!   9.0  !  05-06  (C. Ethe) KPP parameterization
       !!----------------------------------------------------------------------
-      !! * Local declarations
       INTEGER ::   ioptio       ! temporary scalar
-
-      !! * Namelist
-      NAMELIST/namzdf/ ln_zdfexp, n_zdfexp, ln_zdfevd, ln_zdfnpc,   &
-         &             avm0     , avt0    , avevd    , n_evdm
-      !!----------------------------------------------------------------------
-      !!  OPA 9.0, LODYC-IPSL (2003)
-      !! 
+      !!
+      NAMELIST/nam_zdf/ rn_avm0, rn_avt0, nn_avb, nn_havtb, ln_zdfexp, nn_zdfexp,   &
+         &              ln_zdfevd, nn_evdm, rn_avevd, ln_zdfnpc, nn_npc, nn_npcp
       !!----------------------------------------------------------------------
 
-      ! Read namzdf namelist : vertical mixing parameters
-      ! --------------------
-      REWIND( numnam )
-      READ  ( numnam, namzdf )
+      REWIND( numnam )               ! Read nam_zdf namelist : vertical mixing parameters
+      READ  ( numnam, nam_zdf )
 
-      ! Parameter print
-      ! ---------------
-      IF(lwp) THEN
+      IF(lwp) THEN                   ! Parameter print
          WRITE(numout,*)
          WRITE(numout,*) 'zdf_init: vertical physics'
          WRITE(numout,*) '~~~~~~~~'
-         WRITE(numout,*) '   Namelist namzdf : set vertical mixing mixing parameters'
+         WRITE(numout,*) '   Namelist nam_zdf : set vertical mixing mixing parameters'
+         WRITE(numout,*) '      vertical eddy viscosity             rn_avm0   = ', rn_avm0
+         WRITE(numout,*) '      vertical eddy diffusivity           rn_avt0   = ', rn_avt0
+         WRITE(numout,*) '      constant background or profile      nn_avb    = ', nn_avb
+         WRITE(numout,*) '      horizontal variation for avtb       nn_havtb  = ', nn_havtb
          WRITE(numout,*) '      time splitting / backward scheme    ln_zdfexp = ', ln_zdfexp
-         WRITE(numout,*) '      number of time step                  n_zdfexp = ', n_zdfexp
+         WRITE(numout,*) '      number of time step                 nn_zdfexp = ', nn_zdfexp
          WRITE(numout,*) '      enhanced vertical diffusion         ln_zdfevd = ', ln_zdfevd
-         WRITE(numout,*) '      non-penetrative convection          ln_zdfnpc = ', ln_zdfnpc
-         WRITE(numout,*) '      vertical eddy viscosity             avm0      = ', avm0
-         WRITE(numout,*) '      vertical eddy diffusivity           avt0      = ', avt0
-         WRITE(numout,*) '      vertical coefficient for evd        avevd     = ', avevd
-         WRITE(numout,*) '         applied on momentum (=1/0)       n_evdm    = ', n_evdm
+         WRITE(numout,*) '         applied on momentum (=1/0)       nn_evdm   = ', nn_evdm
+         WRITE(numout,*) '      vertical coefficient for evd        rn_avevd  = ', rn_avevd
+         WRITE(numout,*) '      non-penetrative convection (npc)    ln_zdfnpc = ', ln_zdfnpc
+         WRITE(numout,*) '      npc call  frequency                 nn_npc    = ', nn_npc
+         WRITE(numout,*) '      npc print frequency                 nn_npcp   = ', nn_npcp
       ENDIF
 
@@ -144 +141 @@
          &              ' required: ln_zdfevd or ln_zdfnpc logicals' )
 
+
+      !                               !* Background eddy viscosity and diffusivity profil
+      IF( nn_avb == 0 ) THEN                ! Define avmb, avtb from namelist parameter
+         avmb(:) = rn_avm0
+         avtb(:) = rn_avt0                     
+      ELSE                                  ! Background profile of avt (fit a theoretical/observational profile (Krauss 1990)
+         avmb(:) = rn_avm0
+         avtb(:) = rn_avt0 + ( 3.0e-4 - 2 * rn_avt0 ) * 1.0e-4 * gdepw_0(:)   ! m2/s
+         IF(ln_sco .AND. lwp)   CALL ctl_warn( '          avtb profile not valid in sco' )
+      ENDIF
+      !
+      IF( ln_rstart ) THEN                  !  Read avmb, avtb in restart (if exist)
+         ! if ln_traadv_cen, avmb, avtb have been modified in traadv_cen2 module. 
+         ! To ensure the restartability, avmb & avtb are written in the restart 
+         ! file in traadv_cen2 end read here. 
+         IF( iom_varid( numror, 'avmb', ldstop = .FALSE. ) > 0 ) THEN
+            CALL iom_get( numror, jpdom_unknown, 'avmb', avmb )
+            CALL iom_get( numror, jpdom_unknown, 'avtb', avtb )
+         ENDIF
+      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
+      !
    END SUBROUTINE zdf_init