Changeset 2132

Timestamp:

2010-09-29T17:31:36+02:00 (14 years ago)

Author:

cbricaud

Message:

add change from DEV_r1784_GLS

Location:

branches/devmercator2010_1

Files:

• CONFIG/ORCA2_LIM/EXP00/namelist (modified) (1 diff)
• DOC/TexFiles/Biblio/Biblio.bib (modified) (9 diffs)
• DOC/TexFiles/Chapters/Chap_ZDF.tex (modified) (1 diff)
• DOC/TexFiles/Figures/tabgls.png (copied) (copied from branches/DEV_r1784_GLS/DOC/TexFiles/Figures/tabgls.png)
• DOC/TexFiles/Namelist/namgls (copied) (copied from branches/DEV_r1784_GLS/DOC/TexFiles/Namelist/namgls)
• NEMO/OPA_SRC/DYN/dynzdf.F90 (modified) (2 diffs)
• NEMO/OPA_SRC/DYN/dynzdf_imp.F90 (modified) (4 diffs)
• NEMO/OPA_SRC/IOM/restart.F90 (modified) (2 diffs)
• NEMO/OPA_SRC/TRA/trazdf.F90 (modified) (2 diffs)
• NEMO/OPA_SRC/ZDF/zdfbfr.F90 (modified) (1 diff)
• NEMO/OPA_SRC/ZDF/zdfgls.F90 (copied) (copied from branches/DEV_r1784_GLS/NEMO/OPA_SRC/ZDF/zdfgls.F90)
• NEMO/OPA_SRC/ZDF/zdfini.F90 (modified) (4 diffs)
• NEMO/OPA_SRC/step.F90 (modified) (3 diffs)

branches/devmercator2010_1/CONFIG/ORCA2_LIM/EXP00/namelist

@@ -607 +607 @@
    nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
 /
+-----------------------------------------------------------------------
+&namzdf_gls                !   GLS vertical diffusion  ("key_zdfgls")
+!-----------------------------------------------------------------------
+   rn_emin     =   1.e-6   !  minimum value of e   [m2/s2]
+   rn_epsmin   =   1.e-12  !  minimum value of eps [m2/s3]
+   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
+   clim_galp   =   0.53    !  galperin limit
+   ln_crban = .TRUE.       !  Use Craig & Banner (1994) surface wave mixing parametrisation
+   ln_sigpsi = .TRUE.      !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
+   rn_crban = 100.         !  Craig and Banner 1994 constant for wb tke flux
+   rn_charn =  70000.      !  Charnock constant for wb induced roughness length
+   nn_tkebc_surf  =   1    !  surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
+   nn_tkebc_bot   =   1    !  bottom tke condition (0/1=Dir/Neum)
+   nn_psibc_surf  =   1    !  surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
+   nn_psibc_bot   =   1    !  bottom psi condition (0/1=Dir/Neum)
+   nn_stab_func   =   2    !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
+   nn_clos        =   1    !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
+/
 !-----------------------------------------------------------------------
 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")

branches/devmercator2010_1/DOC/TexFiles/Biblio/Biblio.bib

@@ -85 +85 @@
 
 @STRING{Tellus = {Tellus}}
+
+@STRING{RGSP = {Rev. Geophys. Space Phys.}
 
 @ARTICLE{Adcroft_Campin_OM04,
@@ -523 +525 @@
 }
 
+@ARTICLE{Canuto_2001,
+  author = {V. M. Canuto and A. Howard and Y. Cheng and M. S. Dubovikov},
+  title = {Ocean turbulence. PartI: One-point closure model-momentum and heat vertical diffusivities},
+  journal = JPO,
+  year = {2001},
+  volume = {24, 12},
+  pages = {2546-2559},
+  owner = {gr},
+  timestamp = {2010.09.09}
+}
+
 @ARTICLE{Cox1987,
   author = {M. Cox},
@@ -532 +545 @@
   owner = {gm},
   timestamp = {2007.08.03}
+}
+
+@ARTICLE{Craig_Banner_1994,
+  author = {P. D. Banner and M. L. Banner},
+  title = {Modeling wave-enhanced turbulence in the ocean surface layer},
+  journal = JPO,
+  year = {1994},
+  volume = {24, 12},
+  pages = {2546-2559},
+  owner = {g5},
+  timestamp = {2010.09.09}
 }
 
@@ -705 +729 @@
   owner = {gm},
   timestamp = {2007.08.04}
+}
+
+@ARTICLE{Galperin_1988,
+  author = {B. Galperin and L. H. Kantha and S. Hassid and A. Rosati},
+  title = {A quasi-equilibrium turbulent energy model for geophysical flows},
+  journal = JAS,
+  year = {1988},
+  volume = {45},
+  pages = {55-62},
+  owner = {gr},
+  timestamp = {2010.09.09}
 }
 
@@ -1083 +1118 @@
   owner = {gm},
   timestamp = {2008.08.31}
+}
+
+@ARTICLE{Kantha_Clayson_1994,
+  author = {L. H. Kantha and C. A. Clayson},
+  title = {An improved mixed layer model for geophysical applications},
+  journal = JGR,
+  year = {1994},
+  volume = {99},
+  pages = {25235-25266},
+  owner = {gr},
+  timestamp = {2010.09.09}
 }
 
@@ -1584 +1630 @@
 }
 
+@ARTICLE{Mellor_Yamada_1982,
+  author = {G. L. Mellor and T. Yamada},
+  title = {Development of a turbulence closure model for geophysical fluid problems},
+  journal = RGSP,
+  year = {1982},
+  volume = {20},
+  pages = {851-875},
+  owner = {gr},
+  timestamp = {2010.09.09}
+}
+
 @ARTICLE{Merryfield1999,
   author = {W. J. Merryfield and G. Holloway and A. E. Gargett},
@@ -1791 +1848 @@
   owner = {gm},
   timestamp = {2007.08.04}
+}
+
+@ARTICLE{Rodi_1987,
+  author = {W. Rodi},
+  title = {Examples of calculation methods for flow and mixing in stratified fluids},
+  journal = JGR,
+  year = {1987},
+  volume = {92, C5},
+  pages = {5305-5328},
+  owner = {gr},
+  timestamp = {2010.09.09}
 }
 
@@ -2214 +2282 @@
 }
 
+@ARTICLE{Umlauf_Burchard_2003,
+  author = {L. Umlauf and H. Burchard},
+  title = {A generic length-scale equation for geophysical turbulence models},
+  journal = {Journal of Marine Systems}, 
+  year = {2003},
+  volume = {61},
+  pages = {235-265},
+  number = {2},
+  owner = {gr},
+  timestamp = {2010.09.09}
+}
+
 @BOOK{UNESCO1983,
   title = {Algorithms for computation of fundamental property of sea water},
@@ -2308 +2388 @@
   owner = {gm},
   timestamp = {2007.08.04}
+}
+
+@ARTICLE{Wilcox_1988,
+  author = {D. C. Wilcox},
+  title = {Reassessment of the scale-determining equation for advanced turbulence models},
+  journal = {AIAA journal},
+  year = {1988},
+  volume = {26, 11},
+  pages = {1299-1310},
+  owner = {gr},
+  timestamp = {2010.09.09}
 }
 

branches/devmercator2010_1/DOC/TexFiles/Chapters/Chap_ZDF.tex

@@ -233 +233 @@
 
 % -------------------------------------------------------------------------------------------------------------
+%        GLS Generic Length Scale Scheme 
+% -------------------------------------------------------------------------------------------------------------
+\subsection{GLS Generic Length Scale (\key{zdfgls})}
+\label{ZDF_gls}
+
+%--------------------------------------------namgls---------------------------------------------------------
+\namdisplay{namgls}
+%--------------------------------------------------------------------------------------------------------------
+
+The model allows to resolve two prognostic equations for turbulent 
+kinetic energy $\bar {e}$ and a generic length scale \citep{Umlauf_Burchard_2003}. Thanks to the latter, commonly 
+used closures can be retrieved: $k-kl$ \citep{Mellor_Yamada_1982}, $k-{\epsilon }$ \citep{Rodi_1987} and $k-{\omega }$ 
+\citep{Wilcox_1988}. These equations could be written in a generic form with the incorporation 
+of a new variable : ${\psi} = (C^{0}_{\mu})^{p} \ {\bar{e}}^{m} \ l^{n}$.
+
+\begin{equation} \label{Eq_zdfgls_e}
+\frac{\partial \bar{e}}{\partial t} = 
+\frac{A^{vm}}{{\sigma_e} {e_3} }\;\left[ {\left( {\frac{\partial u}{\partial k}} \right)^2
+                                                        +\left( {\frac{\partial v}{\partial k}} \right)^2} \right]
+-A^{vT}\,N^2
++\frac{1}{e_3}  \;\frac{\partial }{\partial k}\left[ {\frac{A^{vm}}{e_3 }
+                                \;\frac{\partial \bar{e}}{\partial k}} \right]
+- \epsilon \;
+\end{equation}
+
+\begin{equation} \label{Eq_zdfgls_psi}
+\frac{\partial \psi}{\partial t} = \frac{\psi}{\bar{e}} 
+(\frac{{C_1}A^{vm}}{{\sigma_{\psi}} {e_3} }\;\left[ {\left( {\frac{\partial u}{\partial k}} \right)^2
+                                                        +\left( {\frac{\partial v}{\partial k}} \right)^2} \right]
+-{C_3}A^{vT}\,N^2- C_2{\epsilon}Fw) 
++\frac{1}{e_3}  \;\frac{\partial }{\partial k}\left[ {\frac{A^{vm}}{e_3 }
+                                \;\frac{\partial \psi}{\partial k}} \right]\;
+\end{equation}
+
+\begin{equation} \label{Eq_zdfgls_kz}
+   \begin{split}
+         A^{vm} &= C_{\mu} \ \sqrt {\bar{e}} \ l         \\
+         A^{vT} &= C_{\mu'}\ \sqrt {\bar{e}} \ l
+   \end{split}
+\end{equation}
+
+\begin{equation} \label{Eq_zdfgls_eps}
+{\epsilon} = (C^{0}_{\mu}) \frac{\bar {e}^{3/2}}{l} \;
+\end{equation}
+where $N$ is the local Brunt-Vais\"{a}l\"{a} frequency (see \S\ref{TRA_bn2}) and $\epsilon$ the dissipation rate.
+In function of the parameters k, m and n, common turbulent closure could be retrieved.
+The constants C1, C2, C3, ${\sigma_e}$, ${\sigma_{\psi}}$ and the wall function (Fw) depends of the choice of the turbulence model.
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+\begin{figure}[!h]
+\centering
+\includegraphics[scale=0.7]{./TexFiles/Figures/tabgls.png}
+\caption {Values of the parameters in function of the model of turbulence.}
+\label{tabgls}
+\end{figure}
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+About the Mellor-Yamada model, the negativity of n allows to use a wall function to force
+the convergence of the mixing length towards Kzb (K: Kappa and zb: rugosity length) value
+near physical boundaries (logarithmic boundary layer law). $C_{\mu}$ and $C_{\mu'}$ are calculated from stability functions 
+of \citet{Galperin_1988}, \citet{Kantha_Clayson_1994} or \citet{Canuto_2001}.
+$C^{0}_{\mu}$ depends of the choice of the stability function.
+
+The boundary condition at the surface and the bottom could be calculated thanks to Diriclet or Neumann condition.
+The wave effect on the mixing could be also being considered \citep{Craig_Banner_1994}.
+
+-------------------------------------------------------------------------------------------------------------
 %        K Profile Parametrisation (KPP) 
 % -------------------------------------------------------------------------------------------------------------

branches/devmercator2010_1/NEMO/OPA_SRC/DYN/dynzdf.F90

@@ -107 +107 @@
       USE zdftke_old
       USE zdftke
+      USE zdfgls
       USE zdfkpp
       !!----------------------------------------------------------------------
@@ -116 +117 @@
 
       ! Force implicit schemes
-      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp )   nzdf = 1   ! TKE or KPP physics
-      IF( ln_dynldf_iso                               )   nzdf = 1   ! iso-neutral lateral physics
-      IF( ln_dynldf_hor .AND. ln_sco                  )   nzdf = 1   ! horizontal lateral physics in s-coordinate
+      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfgls .OR. lk_zdfkpp )   nzdf = 1   ! TKE, GLS or KPP physics
+      IF( ln_dynldf_iso                                              )   nzdf = 1   ! iso-neutral lateral physics
+      IF( ln_dynldf_hor .AND. ln_sco                                 )   nzdf = 1   ! horizontal lateral physics in s-coordinate
 
       IF( lk_esopa )    nzdf = -1                   ! Esopa key: All schemes used

branches/devmercator2010_1/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

@@ -24 +24 @@
    USE phycst          ! physical constants
    USE in_out_manager  ! I/O manager
+#if defined key_zdfgls
+   USE zdfbfr, ONLY : bfrua, bfrva, wbotu, wbotv ! bottom stresses
+   USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
+#endif
 
    IMPLICIT NONE
@@ -75 +79 @@
       REAL(wp) ::   zzwi                               ! temporary scalars
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwi        ! temporary workspace arrays
+#if defined key_zdfgls
+      INTEGER :: ikbu, ikbv, ikbum1, ikbvm1
+      REAL(wp) :: zcbcu, zcbcv
+#endif
       !!----------------------------------------------------------------------
 
@@ -173 +181 @@
          END DO
       END DO
+
+#if defined key_zdfgls
+      ! Save bottom stress for next time step
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ikbu   = MIN( mbathy(ji+1,jj  ), mbathy(ji,jj) )
+            ikbum1 = MAX( ikbu-1, 1 )
+            wbotu(ji,jj) = bfrua(ji,jj) * ua(ji,jj,ikbum1) * umask(ji,jj,ikbum1)
+         END DO
+      END DO
+      CALL lbc_lnk( wbotu(:,:), 'U', -1. )
+#endif
 
       ! Normalization to obtain the general momentum trend ua
@@ -272 +292 @@
       END DO
 
+#if defined key_zdfgls
+      ! Save bottom stress for next time step
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ikbv   = MIN( mbathy(ji,jj+1), mbathy(ji,jj) )
+            ikbvm1 = MAX( ikbv-1, 1 )
+            wbotv(ji,jj) = bfrva(ji,jj) * va(ji,jj,ikbvm1) * vmask(ji,jj,ikbvm1)
+         END DO
+      END DO
+      CALL lbc_lnk( wbotv(:,:), 'V', -1. )
+#endif
+
       ! Normalization to obtain the general momentum trend va
       DO jk = 1, jpkm1

branches/devmercator2010_1/NEMO/OPA_SRC/IOM/restart.F90

@@ -26 +26 @@
    USE zdfmxl          ! mixed layer depth
    USE trdmld_oce      ! ocean active mixed layer tracers trends variables
+#if defined key_zdfgls
+   USE zdfbfr, ONLY : wbotu, wbotv ! bottom stresses
+   USE zdf_oce
+#endif
 
    IMPLICIT NONE
@@ -136 +140 @@
 #endif
 
+#if defined key_zdfgls
+      ! Save bottom stresses
+      CALL iom_rstput( kt, nitrst, numrow, 'wbotu' , wbotu )
+      CALL iom_rstput( kt, nitrst, numrow, 'wbotv' , wbotv )
+#endif
+
       IF( kt == nitrst ) THEN
          CALL iom_close( numrow )     ! close the restart file (only at last time step)

branches/devmercator2010_1/NEMO/OPA_SRC/TRA/trazdf.F90

@@ -125 +125 @@
       USE zdftke_old
       USE zdftke
+      USE zdfgls
       USE zdfkpp
       !!----------------------------------------------------------------------
@@ -139 +140 @@
 
       ! Force implicit schemes
-      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp )   nzdf = 1      ! TKE or KPP physics
-      IF( ln_traldf_iso                               )   nzdf = 1      ! iso-neutral lateral physics
-      IF( ln_traldf_hor .AND. ln_sco                  )   nzdf = 1      ! horizontal lateral physics in s-coordinate
+      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfgls .OR. lk_zdfkpp )   nzdf = 1      ! TKE, GLS or KPP physics
+      IF( ln_traldf_iso                                              )   nzdf = 1      ! iso-neutral lateral physics
+      IF( ln_traldf_hor .AND. ln_sco                                 )   nzdf = 1      ! horizontal lateral physics in s-coordinate
 
       IF( ln_zdfexp .AND. nzdf == 1 )   THEN

branches/devmercator2010_1/NEMO/OPA_SRC/ZDF/zdfbfr.F90

@@ -29 +29 @@
    
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   bfrua , bfrva   !: Bottom friction coefficients set in zdfbfr
+#if defined key_zdfgls
+   REAL(wp), PUBLIC                     ::   rn_hbro =  0.003_wp  ! Bottom roughness (m)
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wbotu, wbotv    !  Bottom stresses
+#endif
 
    !                                    !!* Namelist nambfr: bottom friction namelist *

branches/devmercator2010_1/NEMO/OPA_SRC/ZDF/zdfini.F90

@@ -20 +20 @@
    USE zdftke_old      ! TKE vertical mixing  (old scheme)
    USE zdftke          ! TKE vertical mixing
+   USE zdfgls          ! GLS vertical mixing
    USE zdfkpp          ! KPP vertical mixing          
    USE zdfddm          ! double diffusion mixing      
@@ -106 +107 @@
          ioptio = ioptio+1
       ENDIF
+      IF( lk_zdfgls ) THEN
+         IF(lwp) WRITE(numout,*) '      GLS dependent eddy coefficients'
+         ioptio = ioptio+1
+      ENDIF
       IF( lk_zdfkpp ) THEN
          IF(lwp) WRITE(numout,*) '      KPP dependent eddy coefficients'
@@ -128 +133 @@
          IF(lwp) WRITE(numout,*) '      use the 1.5 turbulent closure'
       ENDIF
+      IF( lk_zdfgls ) THEN
+         IF(lwp) WRITE(numout,*) '      use the GLS closure scheme'
+      ENDIF
       IF( lk_zdfkpp ) THEN
          IF(lwp) WRITE(numout,*) '      use the KPP closure scheme'
@@ -136 +144 @@
       ENDIF
       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', &
+      IF( ioptio == 0 .AND. .NOT.( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfgls .OR. lk_zdfkpp ) ) &
+         CALL ctl_stop( ' except for TKE, GLS or KPP physics, a convection scheme is', &
          &              ' required: ln_zdfevd or ln_zdfnpc logicals' )
 

branches/devmercator2010_1/NEMO/OPA_SRC/step.F90

@@ -90 +90 @@
    USE zdftke_old      ! old TKE vertical mixing      (zdf_tke_old routine)
    USE zdftke          ! TKE vertical mixing              (zdf_tke routine)
+   USE zdfgls          ! GLS vertical mixing              (zdf_gls routine)  
    USE zdfkpp          ! KPP vertical mixing              (zdf_kpp routine)
    USE zdfddm          ! double diffusion mixing          (zdf_ddm routine)
@@ -204 +205 @@
       IF( lk_zdftke_old )   CALL zdf_tke_old( kstp )  ! TKE closure scheme for Kz (old scheme)
       IF( lk_zdftke     )   CALL zdf_tke    ( kstp )  ! TKE closure scheme for Kz
+      IF( lk_zdfgls     )   CALL zdf_gls    ( kstp )  ! GLS closure scheme for Kz
       IF( lk_zdfkpp     )   CALL zdf_kpp    ( kstp )  ! KPP closure scheme for Kz
       IF( lk_zdfcst     )   THEN                      ! Constant Kz (reset avt, avm[uv] to the background value)
@@ -223 +225 @@
                                                       ! write tke information in the restart file
       IF( lrst_oce .AND. lk_zdftke )   CALL tke_rst( kstp, 'WRITE' )
+                                                      ! write gls information in the restart file
+      IF( lrst_oce .AND. lk_zdfgls )   CALL gls_rst( kstp, 'WRITE' )
       !
       !  LATERAL  PHYSICS