Changeset 12126


Ignore:
Timestamp:
2019-12-09T18:11:20+01:00 (6 months ago)
Author:
laurent
Message:

Writing SBCBLK doc…

Location:
NEMO/branches/2019/dev_ASINTER-01-05_merged
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/cfgs/SHARED/namelist_ref

    r12081 r12126  
    6565   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    6666   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    67    ln_xios_read = .FALSE.  !  use XIOS to read restart file (only for a single file restart) 
     67   ln_xios_read = .false.  !  use XIOS to read restart file (only for a single file restart) 
    6868   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 
    6969/ 
     
    257257   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    258258   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    259    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     259   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    260260      ! 
    261261      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/doc/latex/NEMO/subfiles/chap_SBC.tex

    r12063 r12126  
    557557 
    558558Note: all the NEMO Fortran routines involved in the present section have been 
    559  initially developed (and are still developped in parallel) in 
     559 initially developed (and are still developed in parallel) in 
    560560 the \href{https://brodeau.github.io/aerobulk/}{\texttt{AeroBulk}} open-source project 
    561561\citep{brodeau.barnier.ea_JPO17}. 
     
    568568% 
    569569\begin{subequations}\label{eq_bulk} 
     570  \label{eq:SBC_bulk_form} 
    570571  \begin{eqnarray} 
    571     \mathbf{\tau} &=& \rho~ C_D ~ \mathbf{U}_z  ~ U_B \label{eq_b_t} \\ 
    572     Q_H           &=& \rho~C_H~C_P~\big[ \theta_z - T_s \big] ~ U_B \label{eq_b_qh} \\ 
    573     E             &=& \rho~C_E    ~\big[    q_s   - q_z \big] ~ U_B \label{eq_b_e}  \\ 
    574     Q_L           &=& -L_v \, E  \label{eq_b_qe} \\ 
     572    \mathbf{\tau} &=& \rho~ C_D ~ \mathbf{U}_z  ~ U_B \\ 
     573    Q_H           &=& \rho~C_H~C_P~\big[ \theta_z - T_s \big] ~ U_B \\ 
     574    E             &=& \rho~C_E    ~\big[    q_s   - q_z \big] ~ U_B \\ 
     575    Q_L           &=& -L_v \, E \\ 
    575576    % 
    576577    Q_{sr}        &=& (1 - a) Q_{sw\downarrow} \\ 
     
    592593and longwave radiative fluxes, respectively. 
    593594% 
    594 Note: a positive sign of $\mathbf{\tau}$, the various fluxes of heat implies a 
    595 gain of the relevant quantity for the ocean, while a positive $E$ implies a 
    596 freshwater loss for the ocean. 
     595Note: a positive sign for $\mathbf{\tau}$, $Q_H$, $Q_L$, $Q_{sr}$ or $Q_{ir}$ 
     596implies a gain of the relevant quantity for the ocean, while a positive $E$ 
     597implies a freshwater loss for the ocean. 
    597598% 
    598599$\rho$ is the density of air. $C_D$, $C_H$ and $C_E$ are the bulk transfer 
    599 coefficients for momentum, sensible heat, and moisture, respectively (hereafter 
    600 referred to as BTCs). 
     600coefficients for momentum, sensible heat, and moisture, respectively. 
    601601% 
    602602$C_P$ is the heat capacity of moist air, and $L_v$ is the latent heat of 
     
    614614% 
    615615The bulk scalar wind speed, namely $U_B$, is the scalar wind speed, 
    616 $|\mathbf{U}_z|$, with the potential inclusion of a gustiness contribution . 
     616$|\mathbf{U}_z|$, with the potential inclusion of a gustiness contribution. 
    617617% 
    618618$a$ and $\delta$ are the albedo and emissivity of the sea surface, respectively.\\ 
     
    621621% 
    622622$T_s$ is the sea surface temperature. $q_s$ is the saturation specific humidity 
    623 of air at temperature $T_s$ and includes a 2\% reduction to account for the 
     623of air at temperature $T_s$; it includes a 2\% reduction to account for the 
    624624presence of salt in seawater \citep{sverdrup.johnson.ea_1942,kraus.businger_QJRMS96}. 
    625625Depending on the bulk parametrization used, $T_s$ can either be the temperature 
     
    637637(\np[=.true.]{ln_NCAR}{ln\_NCAR}) for which $T_s$ is the bulk SST (\ie~temperature 
    638638at first T-point level). 
    639  
    640639 
    641640For more details on all these aspects the reader is invited to refer 
     
    656655 
    657656 
    658 For the open ocean, four bulk parametrization algorithms are available: 
     657For the open ocean, four bulk parametrization algorithms are available in NEMO: 
    659658\begin{itemize} 
    660659\item NCAR, formerly known as CORE, \citep{large.yeager_rpt04,large.yeager_CD09} 
     
    665664 
    666665 
    667 Differences between versions 3.0 and 3.6 of the COARE algorithm mainly ... wind 
    668 stress BLABLA \citep{edson.jampana.ea_JPO13,brodeau.barnier.ea_JPO17}. 
    669 Therefore it is recommanded to use version 3.6 of the COARE algorithms rather 
    670 than version 3. 
     666With respect to version 3, the principal advances in version 3.6 of the COARE 
     667bulk parametrization are built around improvements in the representation of the 
     668effects of waves on 
     669fluxes \citep{edson.jampana.ea_JPO13,brodeau.barnier.ea_JPO17}. This includes 
     670improved relationships of surface roughness, and whitecap fraction on wave 
     671parameters. It is therefore recommended to chose version 3.6 over 3. 
    671672 
    672673 
     
    708709\subsubsection{NCAR} 
    709710 
    710 NCAR bulk parametrizations (formerly know as CORE) is meant to be used with the 
    711 CORE II atmospheric forcing \citep{large.yeager_CD09}. Hence the following 
    712 namelist parameters must be set: 
     711NCAR bulk parametrizations (formerly known as CORE) is meant to be used with the 
     712CORE II atmospheric forcing \citep{large.yeager_CD09}. The expected sea surface 
     713temperature is the bulk SST. Hence the following namelist parameters must be 
     714set: 
    713715% 
    714716\begin{verbatim} 
     
    728730\subsubsection{ECMWF} 
    729731% 
    730 With a DFS* or any ECMWF-based type of atmospheric forcing, we strongly 
    731 recommend to use the ECMWF bulk parametrizations with the cool-skin and 
    732 warm-layer parametrizations turned on. In ECMWF reanalyzes, since air 
    733 temperature and humidity are provided at the 2\,m height, and given that the 
    734 humidity is provided as the dew-point temperature, the namelist must be tuned as 
    735 follows: 
     732With an atmospheric forcing based on a reanalysis of the ECMWF, such as the 
     733Drakkar Forcing Set \citep{brodeau.barnier.ea_OM10}, we strongly recommend to 
     734use the ECMWF bulk parametrizations with the cool-skin and warm-layer 
     735parametrizations activated. In ECMWF reanalyzes, since air temperature and 
     736humidity are provided at the 2\,m height, and given that the humidity is 
     737distributed as the dew-point temperature, the namelist must be tuned as follows: 
    736738% 
    737739\begin{verbatim} 
     
    758760% 
    759761Since the ECMWF parametrization is largely based on the COARE* parametrization, 
    760 the two algorithms are very similar in terms of structure and closure approach 
    761 (see \citet{brodeau.barnier.ea_JPO17} for the differences). As such, the 
    762 namelist tuning for COARE 3.x is identical to that of ECMWF: 
     762the two algorithms are very similar in terms of structure and closure 
     763approach. As such, the namelist tuning for COARE 3.x is identical to that of 
     764ECMWF: 
    763765% 
    764766\begin{verbatim} 
     
    777779 
    778780 
    779 ~ 
     781%lulu 
    780782 
    781783 
     
    845847 
    846848\np{cn_dir}{cn\_dir} is the directory of location of bulk files 
    847 \np{ln_taudif}{ln\_taudif} is the flag to specify if we use Hight Frequency (HF) tau information (.true.) or not (.false.) 
     849%\np{ln_taudif}{ln\_taudif} is the flag to specify if we use High Frequency (HF) tau information (.true.) or not (.false.) 
    848850\np{rn_zqt}{rn\_zqt}: is the height of humidity and temperature measurements (m) 
    849851\np{rn_zu}{rn\_zu}: is the height of wind measurements (m) 
     
    856858Its range must be between zero and one, and it is recommended to set it to 0 at low-resolution (ORCA2 configuration). 
    857859 
    858 As for the flux parameterization, information about the input data required by the model is provided in 
     860As for the flux parametrization, information about the input data required by the model is provided in 
    859861the namsbc\_blk namelist (see \autoref{subsec:SBC_fldread}). 
    860862 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/doc/namelists/namsbc_blk

    r12081 r12126  
    66   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    77   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    8    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     8   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    99      ! 
    1010      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_coare3p6-noskin_cfg

    r12084 r12126  
    123123   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    124124   ln_COARE_3p6 = .true.    ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     125   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    126126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
    127127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

    r12084 r12126  
    123123   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    124124   ln_COARE_3p6 = .true.    ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     125   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    126126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
    127127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
     
    135135      ln_humi_rlh = .true.  !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    136136   ! 
    137    cn_dir = './'  !  root directory for the bulk data location 
     137   cn_dir      = './'      !  root directory for the bulk data location 
    138138   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    139139   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
     
    164164      ln_read_frq = .false.   !  specify whether we must read frq or not 
    165165 
    166    cn_dir      = './'     !  root directory for the ocean data location 
     166   cn_dir      = './'      !  root directory for the ocean data location 
    167167   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    168168   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     
    216216&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    217217!----------------------------------------------------------------------- 
    218    ln_eos80    = .true.          !  = Use EOS80 
     218   ln_eos80    = .true.         !  = Use EOS80 
    219219/ 
    220220!!====================================================================== 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_ecmwf-noskin_cfg

    r12084 r12126  
    123123   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    124124   ln_COARE_3p6 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125    ln_ECMWF     = .true.    ! "ECMWF"     algorithm   (IFS cycle 31) 
     125   ln_ECMWF     = .true.    ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    126126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
    127127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

    r12084 r12126  
    123123   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    124124   ln_COARE_3p6 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125    ln_ECMWF     = .true.    ! "ECMWF"     algorithm   (IFS cycle 31) 
     125   ln_ECMWF     = .true.    ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    126126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
    127127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
     
    135135      ln_humi_rlh = .true.  !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    136136   ! 
    137    cn_dir = './'  !  root directory for the bulk data location 
     137   cn_dir      = './'      !  root directory for the bulk data location 
    138138   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    139139   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
     
    164164      ln_read_frq = .false.   !  specify whether we must read frq or not 
    165165 
    166    cn_dir      = './'     !  root directory for the ocean data location 
     166   cn_dir      = './'      !  root directory for the ocean data location 
    167167   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    168168   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     
    216216&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    217217!----------------------------------------------------------------------- 
    218    ln_eos80    = .true.          !  = Use EOS80 
     218   ln_eos80    = .true.         !  = Use EOS80 
    219219/ 
    220220!!====================================================================== 
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

    r12084 r12126  
    135135      ln_humi_rlh = .true.  !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    136136   ! 
    137    cn_dir = './'  !  root directory for the bulk data location 
     137   cn_dir      = './'      !  root directory for the bulk data location 
    138138   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    139139   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
     
    164164      ln_read_frq = .false.   !  specify whether we must read frq or not 
    165165 
    166    cn_dir      = './'     !  root directory for the ocean data location 
     166   cn_dir      = './'      !  root directory for the ocean data location 
    167167   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    168168   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     
    216216&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    217217!----------------------------------------------------------------------- 
    218    ln_eos80    = .true.          !  = Use EOS80 
     218   ln_eos80    = .true.         !  = Use EOS80 
    219219/ 
    220220!!====================================================================== 
Note: See TracChangeset for help on using the changeset viewer.