Changeset 12019


Ignore:
Timestamp:
2019-11-30T15:48:32+01:00 (7 weeks ago)
Author:
laurent
Message:

Started to re-write the "sbcblk" par ot the doc…

Location:
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/cfgs/SHARED/namelist_ref

    r11831 r12019  
    250250/ 
    251251!----------------------------------------------------------------------- 
    252 &namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     252&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T) 
    253253!----------------------------------------------------------------------- 
    254254   !                    !  bulk algorithm : 
     
    258258   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    259259      ! 
    260       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    261       rn_zu       = 10.       !  Wind vector reference height (m) 
    262       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    263       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    264       ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    265       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    266       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    267       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    268       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    269       ln_skin_cs = .FALSE.  !  use the cool-skin parameterization (only available in ECMWF and COARE algorithms) !LB 
    270       ln_skin_wl = .FALSE.  !  use the warm-layer        "               "                    " 
    271       ! 
    272       ln_humi_sph = .true.     !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    273       ln_humi_dpt = .false.    !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    274       ln_humi_rlh = .false.    !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
     260      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
     261      rn_zu      = 10.      !  Wind vector reference height (m) 
     262      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
     263      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
     264      ln_taudif  = .false.  !  HF tau contribution: use "mean of stress module  
     265      !                     !  - module of the mean stress" data 
     266      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
     267      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
     268      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity  
     269      !                     !  used to calculate the wind stress 
     270      !                     ! (0. => absolute or 1. => relative winds) 
     271      ln_skin_cs = .false.  !  use the cool-skin parameterization 
     272      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     273      !                     !   ==> only available in ECMWF and COARE algorithms 
     274      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     275      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
     276      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
    275277   ! 
    276278   cn_dir      = './'      !  root directory for the bulk data location 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/doc/latex/NEMO/subfiles/chap_SBC.tex

    r11831 r12019  
    11\documentclass[../main/NEMO_manual]{subfiles} 
     2\usepackage{fontspec} 
     3\usepackage{fontawesome} 
    24 
    35\begin{document} 
     
    504506\label{sec:SBC_flx} 
    505507 
     508% Laurent: DO NOT mix up ``bulk formulae'' (the classic equation) and the ``bulk 
     509% parameterization'' (i.e NCAR, COARE, ECMWF...) 
     510 
    506511\begin{listing} 
    507512  \nlst{namsbc_flx} 
     
    520525See \autoref{subsec:SBC_ssr} for its specification. 
    521526 
    522 %% ================================================================================================= 
     527 
     528 
     529 
     530 
     531 
     532 
     533%% ================================================================================================= 
     534\pagebreak 
     535\newpage 
    523536\section[Bulk formulation (\textit{sbcblk.F90})]{Bulk formulation (\protect\mdl{sbcblk})} 
    524537\label{sec:SBC_blk} 
     
    530543\end{listing} 
    531544 
    532 In the bulk formulation, the surface boundary condition fields are computed with bulk formulae using atmospheric fields 
    533 and ocean (and sea-ice) variables averaged over \np{nn_fsbc}{nn\_fsbc} time-step. 
    534  
    535 The atmospheric fields used depend on the bulk formulae used. 
    536 In forced mode, when a sea-ice model is used, a specific bulk formulation is used. 
    537 Therefore, different bulk formulae are used for the turbulent fluxes computation 
    538 over the ocean and over sea-ice surface. 
    539 For the ocean, four bulk formulations are available thanks to the \href{https://brodeau.github.io/aerobulk/}{Aerobulk} package (\citet{brodeau.barnier.ea_JPO16}): 
    540 the NCAR (formerly named CORE), COARE 3.0, COARE 3.5 and ECMWF bulk formulae. 
    541 The choice is made by setting to true one of the following namelist variable: 
    542  \np{ln_NCAR}{ln\_NCAR}, \np{ln_COARE_3p0}{ln\_COARE\_3p0},  \np{ln_COARE_3p5}{ln\_COARE\_3p5} and  \np{ln_ECMWF}{ln\_ECMWF}. 
    543 For sea-ice, three possibilities can be selected: 
    544 a constant transfer coefficient (1.4e-3; default value), \citet{lupkes.gryanik.ea_JGR12} (\np{ln_Cd_L12}{ln\_Cd\_L12}), and \citet{lupkes.gryanik_JGR15} (\np{ln_Cd_L15}{ln\_Cd\_L15}) parameterizations 
     545In the bulk formulation, the surface boundary condition fields are computed with 
     546bulk formulae using prescribed atmospheric fields and prognostic ocean (and 
     547sea-ice) surface variables averaged over \np{nn_fsbc}{nn\_fsbc} time-step. 
     548 
     549% Turbulent air-sea fluxes are computed using the sea surface properties and 
     550% atmospheric SSVs at height $z$ above the sea surface, with the traditional 
     551% aerodynamic bulk formulae: 
     552 
     553 
     554%%% Bulk formulae are this: 
     555\subsection{Bulk formulae} 
     556% 
     557In NEMO, when the bulk formulation is selected, surface fluxes are computed by means of the traditional bulk formulae: 
     558% 
     559\begin{subequations}\label{eq_bulk} 
     560  \begin{eqnarray} 
     561    \mathbf{\tau} &=& \rho~ C_D ~ \mathbf{U}_z  ~ U_B \label{eq_b_t} \\ 
     562    Q_H           &=& \rho~C_H~C_P~\big[ \theta_z - T_s \big] ~ U_B \label{eq_b_qh} \\ 
     563    E             &=& \rho~C_E    ~\big[    q_s   - q_z \big] ~ U_B \label{eq_b_e}  \\ 
     564    Q_L           &=& -L_v \, E  \label{eq_b_qe} \\ 
     565    % 
     566    Q_{sr}        &=& (1 - a) Q_{sw\downarrow} \\ 
     567    Q_{ir}        &=& \delta (Q_{lw\downarrow} -\sigma T_s^4) 
     568  \end{eqnarray} 
     569\end{subequations} 
     570%lulu 
     571% 
     572From which, the the non-solar heat flux is \[ Q_{ns} = Q_L + Q_H + Q_{ir} \] 
     573% 
     574   \[ \theta_z \simeq T_z+\gamma z \] 
     575   \[  q_s \simeq 0.98\,q_{sat}(T_s,p_a ) \] 
     576 
     577 
     578 
     579where $\mathbf{\tau}$ is the wind stress vector, $Q_H$ the sensible heat flux, 
     580$E$ the evaporation, $Q_L$ the latent heat flux, and $Q_{ir}$ the net longwave 
     581flux. 
     582% 
     583$Q_{sw\downarrow}$ and $Q_{lw\downarrow}$ are the surface downwelling shortwave 
     584and longwave radiative fluxes, respectively. 
     585% 
     586Note: a positive sign of $\mathbf{\tau}$, $Q_H$, and $Q_L$ means a gain of the 
     587relevant quantity for the ocean, while a positive $E$ implies a freshwater loss 
     588for the ocean. 
     589% 
     590$\rho$ is the density of air. $C_D$, $C_H$ and $C_E$ are the BTCs for momentum, 
     591sensible heat, and moisture, respectively.  $C_P$ is the heat capacity of moist 
     592air, and $L_v$ is the latent heat of vaporization of water.  $\theta_z$, $T_z$ 
     593and $q_z$ are the potential temperature, temperature, and specific humidity of 
     594air at height $z$, respectively. $\gamma z$ is a temperature correction term 
     595which accounts for the adiabatic lapse rate and approximates the potential 
     596temperature at height $z$ \citep{Josey_al_2013}.  $\mathbf{U}_z$ is the wind 
     597speed vector at height $z$ (possibly referenced to the surface current 
     598$\mathbf{u_0}$, section \ref{s_res1}.\ref{ss_current}). The bulk scalar wind 
     599speed, $U_B$, is the scalar wind speed, $|\mathbf{U}_z|$, with the potential 
     600inclusion of a gustiness contribution (section 
     601\ref{s_res2}.\ref{ss_calm}). 
     602$P_0$ is the mean sea-level pressure (SLP). 
     603$T_s$ is the sea surface temperature. $q_s$ is the saturation specific humidity 
     604of air at temperature $T_s$ and includes a 2\% reduction to account for the 
     605presence of salt in seawater \citep{Sverdrup_al_1942,Kraus_Businger_1996}. 
     606Depending on the bulk parameterization used, $T_s$ can be the temperature at the 
     607air-sea interface (skin temperature, hereafter SSST) or at a few tens of 
     608centimeters below the surface (bulk sea surface temperature, hereafter SST). 
     609The SSST differs from the SST due to the contributions of two effects of 
     610opposite sign, the \emph{cool skin} and \emph{warm layer} (hereafter CSWL). The 
     611\emph{cool skin} refers to the cooling of the millimeter-scale uppermost layer 
     612of the ocean, in which the net upward flux of heat to the atmosphere is 
     613ineffectively sustained by molecular diffusion. As such, a steep vertical 
     614gradient of temperature must exist to ensure the heat flux continuity with 
     615underlying layers in which the same flux is sustained by turbulence. 
     616The \emph{warm layer} refers to the warming of the upper few meters of the ocean 
     617under sunny conditions. 
     618The CSWL effects are most significant under weak wind conditions due to the 
     619absence of substancial surface vertical mixing (caused by \eg breaking waves). 
     620The impact of the CSWL on the computed TASFs is discussed in section 
     621\ref{s_res1}.\ref{ss_skin}. 
     622 
     623 
     624%%%% Second set of equations (rad): 
     625where $a$ and $\delta$ are the albedo and emissivity of the sea surface, 
     626respectively. 
     627Thus, we use the computed $Q_L$ and $Q_H$ and the 3-hourly surface downwelling 
     628shortwave and longwave radiative fluxes ($Q_{sw\downarrow}$ and 
     629$Q_{lw\downarrow}$, respectively) from ERA-Interim to correct the daily SST 
     630every 3 hours. Due to the implicitness of the problem implied by the dependence 
     631of $Q_{nsol}$ on $T_s$, this correction is done iteratively during the 
     632computation of the TASFs. 
     633 
     634 
     635\subsection{Bulk parameterizations} 
     636 
     637Accuracy of the estimate of surface turbulent fluxes by means of bulk formulae 
     638strongly relies on that of the bulk transfer coefficients: $C_D$, $C_H$ and 
     639$C_E$. They are estimated with what we refer to as a \emph{bulk 
     640parameterization} algorithm. 
     641 
     642... also to adjust humidity and temperature of air to the wind reference measurement 
     643height (generally 10\,m). 
     644 
     645Over the open ocean, four bulk parameterization algorithms are available: 
     646\begin{itemize} 
     647\item NCAR, formerly known as CORE, \citep{large.yeager_rpt04} 
     648\item COARE 3.0 \citep{fairall.bradley.ea_JC03} 
     649\item COARE 3.6 \citep{edson.jampana.ea_JPO13} 
     650\item ECMWF (IFS documentation, cy41) 
     651\end{itemize} 
     652 
     653~ 
     654 
     655% In a typical bulk algorithm, the BTCs under neutral stability conditions are 
     656% defined using \emph{in-situ} flux measurements while their dependence on the 
     657% stability is accounted through the \emph{Monin-Obukhov Similarity Theory} and 
     658% the \emph{flux-profile} relationships \citep[\eg{}][]{Paulson_1970}. BTCs are 
     659% functions of the wind speed and the near-surface stability of the atmospheric 
     660% surface layer (hereafter ASL), and hence, depend on $U_B$, $T_s$, $T_z$, $q_s$ 
     661% and $q_z$. 
     662 
     663 
     664 
     665 
     666\subsection{Cool-skin and warm-layer parameterizations} 
     667 
     668As oposed to the NCAR bulk parameterization, more advanced bulk 
     669parameterizations such as COARE3.x and ECMWF are meant to be used with the skin 
     670temperature $T_s$ rather than the bulk SST (which, in NEMO is the temperature at 
     671the first T-point level). 
     672% 
     673So that, technically, the cool-skin and warm-layer parameterization must be 
     674activated (XXX) to use COARE3.x and ECMWF in a consistant way. 
     675 
     676 
     677\subsection{Air humidity} 
     678 
     679Air humidity can be provided as three different parameters: specific humidity 
     680[kg/kg], relative humidity [\%], or dew-point temperature [K] (LINK to namelist 
     681parameters)... 
     682 
     683 
     684~\\ 
     685 
     686 
     687 
     688 
     689The atmospheric fields used depend on the bulk formulae used.  In forced mode, 
     690when a sea-ice model is used, a specific bulk formulation is used.  Therefore, 
     691different bulk formulae are used for the turbulent fluxes computation over the 
     692ocean and over sea-ice surface. 
     693% 
     694 
     695 
     696thanks to the \href{https://brodeau.github.io/aerobulk/}{Aerobulk} package 
     697(\citet{brodeau.barnier.ea_JPO16}): 
     698 
     699The choice is made by setting to true one of the following namelist 
     700variable: \np{ln_NCAR}{ln\_NCAR}, \np{ln_COARE_3p0}{ln\_COARE\_3p0}, \np{ln_COARE_3p6}{ln\_COARE\_3p6} 
     701and \np{ln_ECMWF}{ln\_ECMWF}.  For sea-ice, three possibilities can be selected: 
     702a constant transfer coefficient (1.4e-3; default 
     703value), \citet{lupkes.gryanik.ea_JGR12} (\np{ln_Cd_L12}{ln\_Cd\_L12}), 
     704and \citet{lupkes.gryanik_JGR15} (\np{ln_Cd_L15}{ln\_Cd\_L15}) parameterizations 
    545705 
    546706Common options are defined through the \nam{sbc_blk}{sbc\_blk} namelist variables. 
     
    557717    j-component of the 10m air velocity  & vtau           & $m.s^{-1}$         & T     \\ 
    558718    \hline 
    559     10m air temperature                  & tair           & \r{}$K$            & T     \\ 
     719    10m air temperature                  & tair           & $K$               & T     \\ 
    560720    \hline 
    561     Specific humidity                    & humi           & \%                 & T     \\ 
     721    Specific humidity                    & humi           & $-$               & T     \\ 
     722    Relative humidity                    & ~              & $\%$              & T     \\ 
     723    Dew-point temperature                & ~              & $K$               & T     \\     
    562724    \hline 
    563     Incoming long wave radiation         & qlw            & $W.m^{-2}$         & T     \\ 
     725    Downwelling longwave radiation       & qlw            & $W.m^{-2}$         & T     \\ 
    564726    \hline 
    565     Incoming short wave radiation        & qsr            & $W.m^{-2}$         & T     \\ 
     727    Downwelling shortwave radiation      & qsr            & $W.m^{-2}$         & T     \\ 
    566728    \hline 
    567729    Total precipitation (liquid + solid) & precip         & $Kg.m^{-2}.s^{-1}$ & T     \\ 
     
    595757Its range must be between zero and one, and it is recommended to set it to 0 at low-resolution (ORCA2 configuration). 
    596758 
    597 As for the flux formulation, information about the input data required by the model is provided in 
     759As for the flux parameterization, information about the input data required by the model is provided in 
    598760the namsbc\_blk namelist (see \autoref{subsec:SBC_fldread}). 
    599761 
    600762%% ================================================================================================= 
    601 \subsection[Ocean-Atmosphere Bulk formulae (\textit{sbcblk\_algo\_coare.F90, sbcblk\_algo\_coare3p5.F90, sbcblk\_algo\_ecmwf.F90, sbcblk\_algo\_ncar.F90})]{Ocean-Atmosphere Bulk formulae (\mdl{sbcblk\_algo\_coare}, \mdl{sbcblk\_algo\_coare3p5}, \mdl{sbcblk\_algo\_ecmwf}, \mdl{sbcblk\_algo\_ncar})} 
     763\subsection[Ocean-Atmosphere Bulk formulae (\textit{sbcblk\_algo\_coare.F90, sbcblk\_algo\_coare3p6.F90, sbcblk\_algo\_ecmwf.F90, sbcblk\_algo\_ncar.F90})]{Ocean-Atmosphere Bulk formulae (\mdl{sbcblk\_algo\_coare}, \mdl{sbcblk\_algo\_coare3p6}, \mdl{sbcblk\_algo\_ecmwf}, \mdl{sbcblk\_algo\_ncar})} 
    602764\label{subsec:SBC_blk_ocean} 
    603765 
    604766Four different bulk algorithms are available to compute surface turbulent momentum and heat fluxes over the ocean. 
    605 COARE 3.0, COARE 3.5 and ECMWF schemes mainly differ by their roughness lenghts computation and consequently 
     767COARE 3.0, COARE 3.6 and ECMWF schemes mainly differ by their roughness lenghts computation and consequently 
    606768their neutral transfer coefficients relationships with neutral wind. 
    607769\begin{itemize} 
     
    615777  This is the so-called DRAKKAR Forcing Set (DFS) \citep{brodeau.barnier.ea_OM10}. 
    616778\item COARE 3.0 (\np[=.true.]{ln_COARE_3p0}{ln\_COARE\_3p0}): See \citet{fairall.bradley.ea_JC03} for more details 
    617 \item COARE 3.5 (\np[=.true.]{ln_COARE_3p5}{ln\_COARE\_3p5}): See \citet{edson.jampana.ea_JPO13} for more details 
     779\item COARE 3.6 (\np[=.true.]{ln_COARE_3p6}{ln\_COARE\_3p6}): See \citet{edson.jampana.ea_JPO13} for more details 
    618780\item ECMWF (\np[=.true.]{ln_ECMWF}{ln\_ECMWF}): Based on \href{https://www.ecmwf.int/node/9221}{IFS (Cy31)} implementation and documentation. 
    619781  Surface roughness lengths needed for the Obukhov length are computed following \citet{beljaars_QJRMS95}. 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/doc/latex/global/packages.tex

    r11702 r12019  
    1818%% Issue with fontawesome pkg: path to FontAwesome.otf has to be hard-coded 
    1919\defaultfontfeatures{ 
    20     Path = /usr/local/texlive/2019/texmf-dist/fonts/opentype/public/fontawesome/ 
     20      Path = /usr/share/texlive/texmf-dist/tex/latex/fontawesome/opentype/ 
    2121} 
    2222\usepackage{academicons, fontawesome, newtxtext} 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/doc/namelists/namsbc_blk

    r11831 r12019  
    11!----------------------------------------------------------------------- 
    2 &namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     2&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T) 
    33!----------------------------------------------------------------------- 
    44   !                    !  bulk algorithm : 
    5    ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     5   ln_NCAR      = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
    66   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    7    ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    8    ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     7   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
     8   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    99      ! 
    10       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    11       rn_zu       = 10.       !  Wind vector reference height (m) 
    12       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    13       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    14       ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    15       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    16       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    17       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    18       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    19  
     10      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
     11      rn_zu      = 10.      !  Wind vector reference height (m) 
     12      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
     13      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
     14      ln_taudif  = .false.  !  HF tau contribution: use "mean of stress module  
     15      !                     !  - module of the mean stress" data 
     16      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
     17      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
     18      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity  
     19      !                     !  used to calculate the wind stress 
     20      !                     ! (0. => absolute or 1. => relative winds) 
     21      ln_skin_cs = .false.  !  use the cool-skin parameterization 
     22      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     23      !                     !   ==> only available in ECMWF and COARE algorithms 
     24      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     25      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
     26      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     27   ! 
    2028   cn_dir      = './'      !  root directory for the bulk data location 
    2129   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/README.md

    r11930 r12019  
    66## Objectives 
    77 
    8 ```STATION_ASF``` is a demonstration case that mimics an in-situ station (buoy, platform) dedicated to the estimation of surface air-sea fluxes by means of the measurement of traditional meteorological surface parameters. 
     8```STATION_ASF``` is a demonstration case that mimics an in-situ station (buoy, platform) dedicated to the estimation of surface air-sea fluxes by means of the more widely measured traditional meteorological surface parameters (sea and atmosphere). 
    99 
    1010```STATION_ASF``` is based on the merging of the "single column" and the "standalone surface module" configurations of NEMO. In short, it coulb defined as "SAS meets C1D". As such, the spatial domain of ```STATION_ASF``` is punctual (1D, well actually 3 x 3 as in C1D). 
Note: See TracChangeset for help on using the changeset viewer.