Changeset 12126
- Timestamp:
- 2019-12-09T18:11:20+01:00 (4 years ago)
- 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 65 65 ln_clobber = .true. ! clobber (overwrite) an existing file 66 66 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) 68 68 nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 69 69 / … … 257 257 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 258 258 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) 260 260 ! 261 261 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 557 557 558 558 Note: all the NEMO Fortran routines involved in the present section have been 559 initially developed (and are still develop ped in parallel) in559 initially developed (and are still developed in parallel) in 560 560 the \href{https://brodeau.github.io/aerobulk/}{\texttt{AeroBulk}} open-source project 561 561 \citep{brodeau.barnier.ea_JPO17}. … … 568 568 % 569 569 \begin{subequations}\label{eq_bulk} 570 \label{eq:SBC_bulk_form} 570 571 \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 \\ 575 576 % 576 577 Q_{sr} &=& (1 - a) Q_{sw\downarrow} \\ … … 592 593 and longwave radiative fluxes, respectively. 593 594 % 594 Note: a positive sign of $\mathbf{\tau}$, the various fluxes of heat implies a595 gain of the relevant quantity for the ocean, while a positive $E$ implies a 596 freshwater loss for the ocean.595 Note: a positive sign for $\mathbf{\tau}$, $Q_H$, $Q_L$, $Q_{sr}$ or $Q_{ir}$ 596 implies a gain of the relevant quantity for the ocean, while a positive $E$ 597 implies a freshwater loss for the ocean. 597 598 % 598 599 $\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). 600 coefficients for momentum, sensible heat, and moisture, respectively. 601 601 % 602 602 $C_P$ is the heat capacity of moist air, and $L_v$ is the latent heat of … … 614 614 % 615 615 The 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. 617 617 % 618 618 $a$ and $\delta$ are the albedo and emissivity of the sea surface, respectively.\\ … … 621 621 % 622 622 $T_s$ is the sea surface temperature. $q_s$ is the saturation specific humidity 623 of air at temperature $T_s$ andincludes a 2\% reduction to account for the623 of air at temperature $T_s$; it includes a 2\% reduction to account for the 624 624 presence of salt in seawater \citep{sverdrup.johnson.ea_1942,kraus.businger_QJRMS96}. 625 625 Depending on the bulk parametrization used, $T_s$ can either be the temperature … … 637 637 (\np[=.true.]{ln_NCAR}{ln\_NCAR}) for which $T_s$ is the bulk SST (\ie~temperature 638 638 at first T-point level). 639 640 639 641 640 For more details on all these aspects the reader is invited to refer … … 656 655 657 656 658 For the open ocean, four bulk parametrization algorithms are available :657 For the open ocean, four bulk parametrization algorithms are available in NEMO: 659 658 \begin{itemize} 660 659 \item NCAR, formerly known as CORE, \citep{large.yeager_rpt04,large.yeager_CD09} … … 665 664 666 665 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. 666 With respect to version 3, the principal advances in version 3.6 of the COARE 667 bulk parametrization are built around improvements in the representation of the 668 effects of waves on 669 fluxes \citep{edson.jampana.ea_JPO13,brodeau.barnier.ea_JPO17}. This includes 670 improved relationships of surface roughness, and whitecap fraction on wave 671 parameters. It is therefore recommended to chose version 3.6 over 3. 671 672 672 673 … … 708 709 \subsubsection{NCAR} 709 710 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: 711 NCAR bulk parametrizations (formerly known as CORE) is meant to be used with the 712 CORE II atmospheric forcing \citep{large.yeager_CD09}. The expected sea surface 713 temperature is the bulk SST. Hence the following namelist parameters must be 714 set: 713 715 % 714 716 \begin{verbatim} … … 728 730 \subsubsection{ECMWF} 729 731 % 730 With a DFS* or any ECMWF-based type of atmospheric forcing, we strongly731 recommend to use the ECMWF bulk parametrizations with the cool-skin and 732 warm-layer parametrizations turned on. In ECMWF reanalyzes, since air733 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 as735 follows:732 With an atmospheric forcing based on a reanalysis of the ECMWF, such as the 733 Drakkar Forcing Set \citep{brodeau.barnier.ea_OM10}, we strongly recommend to 734 use the ECMWF bulk parametrizations with the cool-skin and warm-layer 735 parametrizations activated. In ECMWF reanalyzes, since air temperature and 736 humidity are provided at the 2\,m height, and given that the humidity is 737 distributed as the dew-point temperature, the namelist must be tuned as follows: 736 738 % 737 739 \begin{verbatim} … … 758 760 % 759 761 Since the ECMWF parametrization is largely based on the COARE* parametrization, 760 the two algorithms are very similar in terms of structure and closure approach761 (see \citet{brodeau.barnier.ea_JPO17} for the differences). As such, the 762 namelist tuning for COARE 3.x is identical to that ofECMWF:762 the two algorithms are very similar in terms of structure and closure 763 approach. As such, the namelist tuning for COARE 3.x is identical to that of 764 ECMWF: 763 765 % 764 766 \begin{verbatim} … … 777 779 778 780 779 ~ 781 %lulu 780 782 781 783 … … 845 847 846 848 \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 HightFrequency (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.) 848 850 \np{rn_zqt}{rn\_zqt}: is the height of humidity and temperature measurements (m) 849 851 \np{rn_zu}{rn\_zu}: is the height of wind measurements (m) … … 856 858 Its range must be between zero and one, and it is recommended to set it to 0 at low-resolution (ORCA2 configuration). 857 859 858 As for the flux paramet erization, information about the input data required by the model is provided in860 As for the flux parametrization, information about the input data required by the model is provided in 859 861 the namsbc\_blk namelist (see \autoref{subsec:SBC_fldread}). 860 862 -
NEMO/branches/2019/dev_ASINTER-01-05_merged/doc/namelists/namsbc_blk
r12081 r12126 6 6 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 7 7 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) 9 9 ! 10 10 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 123 123 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 124 124 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) 126 126 rn_zqt = 2.5 ! Air temperature & humidity reference height (m) [PAPA station] 127 127 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 123 123 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 124 124 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) 126 126 rn_zqt = 2.5 ! Air temperature & humidity reference height (m) [PAPA station] 127 127 rn_zu = 4. ! Wind vector reference height (m) [PAPA station] … … 135 135 ln_humi_rlh = .true. ! humidity specified below in "sn_humi" is relative humidity [%] if .true. 136 136 ! 137 cn_dir = './'! root directory for the bulk data location137 cn_dir = './' ! root directory for the bulk data location 138 138 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 139 139 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 164 164 ln_read_frq = .false. ! specify whether we must read frq or not 165 165 166 cn_dir = './' ! root directory for the ocean data location166 cn_dir = './' ! root directory for the ocean data location 167 167 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 168 168 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 216 216 &nameos ! ocean Equation Of Seawater (default: NO selection) 217 217 !----------------------------------------------------------------------- 218 ln_eos80 = .true. 218 ln_eos80 = .true. ! = Use EOS80 219 219 / 220 220 !!====================================================================== -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_ecmwf-noskin_cfg
r12084 r12126 123 123 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 124 124 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) 126 126 rn_zqt = 2.5 ! Air temperature & humidity reference height (m) [PAPA station] 127 127 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 123 123 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 124 124 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) 126 126 rn_zqt = 2.5 ! Air temperature & humidity reference height (m) [PAPA station] 127 127 rn_zu = 4. ! Wind vector reference height (m) [PAPA station] … … 135 135 ln_humi_rlh = .true. ! humidity specified below in "sn_humi" is relative humidity [%] if .true. 136 136 ! 137 cn_dir = './'! root directory for the bulk data location137 cn_dir = './' ! root directory for the bulk data location 138 138 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 139 139 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 164 164 ln_read_frq = .false. ! specify whether we must read frq or not 165 165 166 cn_dir = './' ! root directory for the ocean data location166 cn_dir = './' ! root directory for the ocean data location 167 167 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 168 168 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 216 216 &nameos ! ocean Equation Of Seawater (default: NO selection) 217 217 !----------------------------------------------------------------------- 218 ln_eos80 = .true. 218 ln_eos80 = .true. ! = Use EOS80 219 219 / 220 220 !!====================================================================== -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/STATION_ASF/EXPREF/namelist_ncar_cfg
r12084 r12126 135 135 ln_humi_rlh = .true. ! humidity specified below in "sn_humi" is relative humidity [%] if .true. 136 136 ! 137 cn_dir = './'! root directory for the bulk data location137 cn_dir = './' ! root directory for the bulk data location 138 138 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 139 139 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 164 164 ln_read_frq = .false. ! specify whether we must read frq or not 165 165 166 cn_dir = './' ! root directory for the ocean data location166 cn_dir = './' ! root directory for the ocean data location 167 167 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 168 168 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 216 216 &nameos ! ocean Equation Of Seawater (default: NO selection) 217 217 !----------------------------------------------------------------------- 218 ln_eos80 = .true. 218 ln_eos80 = .true. ! = Use EOS80 219 219 / 220 220 !!======================================================================
Note: See TracChangeset
for help on using the changeset viewer.