\hypertarget{namespacesflx}{}\doxysubsection{sflx Module Reference}
\label{namespacesflx}\index{sflx@{sflx}}
\doxysubsubsection*{Functions/\+Subroutines}
\begin{DoxyCompactItemize}
\item 
subroutine \mbox{\hyperlink{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_ga36e98f3ac91f839b816fd167eab30901}{gfssflx}} (nsoil, couple, icein, ffrozp, dt, zlvl, sldpth, swdn, swnet, lwdn, sfcems, sfcprs, sfctmp, sfcspd, prcp, q2, q2sat, dqsdt2, th2, ivegsrc, vegtyp, soiltyp, slopetyp, shdmin, alb, snoalb, rhonewsn, exticeden, bexpp, xlaip, lheatstrg, tbot, cmc, t1, stc, smc, sh2o, sneqv, ch, cm, z0, nroot, shdfac, snowh, albedo, eta, sheat, ec, edir, et, ett, esnow, drip, dew, beta, etp, ssoil, flx1, flx2, flx3, runoff1, runoff2, runoff3, snomlt, sncovr, rc, pc, rsmin, xlai, rcs, rct, rcq, rcsoil, soilw, soilm, smcwlt, smcdry, smcref, smcmax, errmsg, errflg)
\begin{DoxyCompactList}\small\item\em This is the entity of GFS Noah LSM model of physics subroutines. It is a soil/veg/snowpack land-\/surface model to update soil moisture, soil ice, soil temperature, skin temperature, snowpack water content, snowdepth, and all terms of the surface energy balance and surface water balance (excluding input atmospheric forcings of downward radiation and precipitation). \end{DoxyCompactList}\end{DoxyCompactItemize}


\doxysubsubsection{Function/\+Subroutine Documentation}
\mbox{\Hypertarget{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_ga36e98f3ac91f839b816fd167eab30901}\label{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_ga36e98f3ac91f839b816fd167eab30901}} 
\index{sflx@{sflx}!gfssflx@{gfssflx}}
\index{gfssflx@{gfssflx}!sflx@{sflx}}
\doxyparagraph{\texorpdfstring{gfssflx()}{gfssflx()}}
{\footnotesize\ttfamily subroutine sflx\+::gfssflx (\begin{DoxyParamCaption}\item[{integer, intent(in)}]{nsoil,  }\item[{integer, intent(in)}]{couple,  }\item[{integer, intent(in)}]{icein,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{ffrozp,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{dt,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{zlvl,  }\item[{real (kind=kind\+\_\+phys), dimension(nsoil), intent(in)}]{sldpth,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{swdn,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{swnet,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{lwdn,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{sfcems,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{sfcprs,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{sfctmp,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{sfcspd,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{prcp,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{q2,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{q2sat,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{dqsdt2,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{th2,  }\item[{integer, intent(in)}]{ivegsrc,  }\item[{integer, intent(in)}]{vegtyp,  }\item[{integer, intent(in)}]{soiltyp,  }\item[{integer, intent(in)}]{slopetyp,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{shdmin,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{alb,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{snoalb,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{rhonewsn,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{exticeden,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{bexpp,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{xlaip,  }\item[{real (kind=kind\+\_\+phys), intent(in)}]{lheatstrg,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{tbot,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{cmc,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{t1,  }\item[{real (kind=kind\+\_\+phys), dimension(nsoil), intent(inout)}]{stc,  }\item[{real (kind=kind\+\_\+phys), dimension(nsoil), intent(inout)}]{smc,  }\item[{real (kind=kind\+\_\+phys), dimension(nsoil), intent(inout)}]{sh2o,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{sneqv,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{ch,  }\item[{real (kind=kind\+\_\+phys), intent(inout)}]{cm,  }\item[{real (kind=kind\+\_\+phys)}]{z0,  }\item[{integer, intent(out)}]{nroot,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{shdfac,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{snowh,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{albedo,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{eta,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{sheat,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{ec,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{edir,  }\item[{real (kind=kind\+\_\+phys), dimension(nsoil), intent(out)}]{et,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{ett,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{esnow,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{drip,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{dew,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{beta,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{etp,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{ssoil,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{flx1,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{flx2,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{flx3,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{runoff1,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{runoff2,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{runoff3,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{snomlt,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{sncovr,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rc,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{pc,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rsmin,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{xlai,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rcs,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rct,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rcq,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{rcsoil,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{soilw,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{soilm,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{smcwlt,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{smcdry,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{smcref,  }\item[{real (kind=kind\+\_\+phys), intent(out)}]{smcmax,  }\item[{character(len=$\ast$), intent(out)}]{errmsg,  }\item[{integer, intent(out)}]{errflg }\end{DoxyParamCaption})}

The land-\/surface model component was substantially upgraded from the Oregon State University (OSU) land surface model to EMC\textquotesingle{}s new Noah Land Surface Model (Noah LSM) during the major implementation in the NCEP Global Forecast System (GFS) on May 31, 2005. Forecast System (GFS). The Noah LSM embodies about 10 years of upgrades (see \cite{chen_et_al_1996}, \cite{koren_et_al_1999}, \cite{ek_et_al_2003}) to its ancestor, the OSU LSM. The Noah LSM upgrade includes\+:
\begin{DoxyItemize}
\item An increase from two (10, 190 cm thick) to four soil layers (10, 30, 60, 100 cm thick)
\item Addition of frozen soil physics
\item Add glacial ice treatment
\item Two snowpack states (SWE, density)
\item New formulations for infiltration and runoff account for sub-\/grid variability in precipitation and soil moisture
\item Revised physics of the snowpack and its influence on surface heat fluxes and albedo
\item Higher canopy resistance
\item Spatially varying root depth
\item Surface fluxes weighted by snow cover fraction
\item Improved thermal conduction in soil/snow
\item Improved seasonality of green vegetation cover.
\item Improved evaporation treatment over bare soil and snowpack
\end{DoxyItemize}


\begin{DoxyParams}[1]{Parameters}
\mbox{\texttt{ in}}  & {\em nsoil} & integer, number of soil layers (\texorpdfstring{$>$}{>}=2 but \texorpdfstring{$<$}{<}=nsold) \\
\hline
\mbox{\texttt{ in}}  & {\em couple} & integer, =0\+:uncoupled (land model only), =1\+:coupled with parent atmos model \\
\hline
\mbox{\texttt{ in}}  & {\em icein} & integer, sea-\/ice flag (=1\+: sea-\/ice, =0\+: land) \\
\hline
\mbox{\texttt{ in}}  & {\em ffrozp} & real, flag for snow-\/rain detection (1.=snow, 0.=rain) \\
\hline
\mbox{\texttt{ in}}  & {\em dt} & real, time step (\texorpdfstring{$<$}{<}3600 sec) \\
\hline
\mbox{\texttt{ in}}  & {\em zlvl} & real, height abv atmos ground forcing vars ( $m$) \\
\hline
\mbox{\texttt{ in}}  & {\em sldpth} & real, thickness of each soil layer ( $m$), nsoil \\
\hline
\mbox{\texttt{ in}}  & {\em swdn} & real, downward SW radiation flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ in}}  & {\em swnet} & real, downward SW net (dn-\/up) flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ in}}  & {\em lwdn} & real, downward LW radiation flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ in}}  & {\em sfcems} & real, sfc LW emissivity (fractional) \\
\hline
\mbox{\texttt{ in}}  & {\em sfcprs} & real, pressure at height zlvl above ground( $Pa$) \\
\hline
\mbox{\texttt{ in}}  & {\em sfctmp} & real, air temp at height zlvl above ground ( $K$) \\
\hline
\mbox{\texttt{ in}}  & {\em sfcspd} & real, wind speed at height zlvl above ground ( $m s^{-1}$) \\
\hline
\mbox{\texttt{ in}}  & {\em prcp} & real, precipitation rate ( $kgm^{-2}s^{-1}$) \\
\hline
\mbox{\texttt{ in}}  & {\em q2} & real, mixing ratio at hght zlvl above ground ( $kgkg^{-1}$) \\
\hline
\mbox{\texttt{ in}}  & {\em q2sat} & real, sat mixing ratio at zlvl above ground ( $kgkg^{-1}$) \\
\hline
\mbox{\texttt{ in}}  & {\em dqsdt2} & real, slope of sat specific humidity curve at t=sfctmp ( $kgkg^{-1}k^{-1}$) \\
\hline
\mbox{\texttt{ in}}  & {\em th2} & real, air potential temperature at zlvl above ground ( $K$) \\
\hline
\mbox{\texttt{ in}}  & {\em ivegsrc} & integer, sfc veg type data source UMD or IGBP \\
\hline
\mbox{\texttt{ in}}  & {\em vegtyp} & integer, vegetation type (integer index) \\
\hline
\mbox{\texttt{ in}}  & {\em soiltyp} & integer, soil type (integer index) \\
\hline
\mbox{\texttt{ in}}  & {\em slopetyp} & integer, class of sfc slope (integer index) \\
\hline
\mbox{\texttt{ in}}  & {\em shdmin} & real, min areal coverage of green veg (fraction) \\
\hline
\mbox{\texttt{ in}}  & {\em alb} & real, background snow-\/free sfc albedo (fraction) \\
\hline
\mbox{\texttt{ in}}  & {\em snoalb} & real, max albedo over deep snow (fraction) \\
\hline
\mbox{\texttt{ in}}  & {\em bexpp} & real, perturbation of soil type \char`\"{}b\char`\"{} parameter (perturbation) \\
\hline
\mbox{\texttt{ in}}  & {\em xlaip} & real, perturbation of leave area index (perturbation) \\
\hline
\mbox{\texttt{ in}}  & {\em lheatstrg} & logical, flag for canopy heat storage parameterization \\
\hline
\mbox{\texttt{ in,out}}  & {\em tbot} & real, bottom soil temp ( $K$) (local yearly-\/mean sfc air temp) \\
\hline
\mbox{\texttt{ in,out}}  & {\em cmc} & real, canopy moisture content ( $m$) \\
\hline
\mbox{\texttt{ in,out}}  & {\em t1} & real, ground/canopy/snowpack eff skin temp ( $K$) \\
\hline
\mbox{\texttt{ in,out}}  & {\em stc} & real, soil temp ( $K$) \\
\hline
\mbox{\texttt{ in,out}}  & {\em smc} & real, total soil moisture (vol fraction) \\
\hline
\mbox{\texttt{ in,out}}  & {\em sh2o} & real, unfrozen soil moisture (vol fraction), note\+: frozen part = smc-\/sh2o \\
\hline
\mbox{\texttt{ in,out}}  & {\em sneqv} & real, water-\/equivalent snow depth ( $m$), note\+: snow density = snwqv/snowh \\
\hline
\mbox{\texttt{ in,out}}  & {\em ch} & real, sfc exchange coeff for heat \& moisture ( $ms^{-1}$), note\+: conductance since it\textquotesingle{}s been mult by wind \\
\hline
\mbox{\texttt{ in,out}}  & {\em cm} & real, sfc exchange coeff for momentum ( $ms^{-1}$), note\+: conductance since it\textquotesingle{}s been mult by wind \\
\hline
\mbox{\texttt{ in,out}}  & {\em z0} & real, roughness length ( $m$) \\
\hline
\mbox{\texttt{ out}}  & {\em nroot} & integer, number of root layers \\
\hline
\mbox{\texttt{ out}}  & {\em shdfac} & real, aeral coverage of green veg (fraction) \\
\hline
\mbox{\texttt{ out}}  & {\em snowh} & real, snow depth ( $m$) \\
\hline
\mbox{\texttt{ out}}  & {\em albedo} & real, sfc albedo incl snow effect (fraction) \\
\hline
\mbox{\texttt{ out}}  & {\em eta} & real, downward latent heat flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em sheat} & real, downward sensible heat flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em ec} & real, canopy water evaporation ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em edir} & real, direct soil evaporation ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em et} & real, plant transpiration ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em ett} & real, total plant transpiration ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em esnow} & real, sublimation from snowpack ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em drip} & real, through-\/fall of precip and/or dew in excess of canopy water-\/holding capacity ( $m$) \\
\hline
\mbox{\texttt{ out}}  & {\em dew} & real, dewfall (or frostfall for t\texorpdfstring{$<$}{<}273.\+15) ( $m$) \\
\hline
\mbox{\texttt{ out}}  & {\em beta} & real, ratio of actual/potential evap \\
\hline
\mbox{\texttt{ out}}  & {\em etp} & real, potential evaporation ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em ssoil} & real, upward soil heat flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em flx1} & real, precip-\/snow sfc flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em flx2} & real, freezing rain latent heat flux ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em flx3} & real, phase-\/change heat flux from snowmelt ( $W/m^2$) \\
\hline
\mbox{\texttt{ out}}  & {\em runoff1} & real, surface runoff ( $ms^{-1}$) not infiltrating sfc \\
\hline
\mbox{\texttt{ out}}  & {\em runoff2} & real, sub sfc runoff ( $ms^{-1}$) (baseflow) \\
\hline
\mbox{\texttt{ out}}  & {\em runoff3} & real, excess of porosity for a given soil layer \\
\hline
\mbox{\texttt{ out}}  & {\em snomlt} & real, snow melt ( $m$) (water equivalent) \\
\hline
\mbox{\texttt{ out}}  & {\em sncovr} & real, fractional snow cover \\
\hline
\mbox{\texttt{ out}}  & {\em rc} & real, canopy resistance (s/m) \\
\hline
\mbox{\texttt{ out}}  & {\em pc} & real, plant coeff (fraction) where pc$\ast$etp=transpi \\
\hline
\mbox{\texttt{ out}}  & {\em rsmin} & real, minimum canopy resistance (s/m) \\
\hline
\mbox{\texttt{ out}}  & {\em xlai} & real, leaf area index (dimensionless) \\
\hline
\mbox{\texttt{ out}}  & {\em rcs} & real, incoming solar rc factor (dimensionless) \\
\hline
\mbox{\texttt{ out}}  & {\em rct} & real, air temperature rc factor (dimensionless) \\
\hline
\mbox{\texttt{ out}}  & {\em rcq} & real, atoms vapor press deficit rc factor \\
\hline
\mbox{\texttt{ out}}  & {\em rcsoil} & real, soil moisture rc factor (dimensionless) \\
\hline
\mbox{\texttt{ out}}  & {\em soilw} & real, available soil moisture in root zone \\
\hline
\mbox{\texttt{ out}}  & {\em soilm} & real, total soil column moisture (frozen+unfrozen) ( $m$) \\
\hline
\mbox{\texttt{ out}}  & {\em smcwlt} & real, wilting point (volumetric) \\
\hline
\mbox{\texttt{ out}}  & {\em smcdry} & real, dry soil moisture threshold (volumetric) \\
\hline
\mbox{\texttt{ out}}  & {\em smcref} & real, soil moisture threshold (volumetric) \\
\hline
\mbox{\texttt{ out}}  & {\em smcmax} & real, porosity (sat val of soil mois) \\
\hline
\end{DoxyParams}
\hypertarget{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_general_sflx}{}\doxysubsubsection{GFS Noah LSM General Algorithm}\label{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_general_sflx}

\begin{DoxyItemize}
\item Set ice = -\/1 and green vegetation fraction (shdfac) = 0 for glacial-\/ice land.
\item Calculate soil layer depth below ground.
\item For ice, set green vegetation fraction (shdfac) = 0. and set sea-\/ice layers of equal thickness and sum to 3 meters
\item Otherwise, calculate depth (negative) below ground from top skin sfc to bottom of each soil layer.
\item Call \mbox{\hyperlink{group___noah___l_s_m_ga1c19b4bae5de2a391c019b6778a58dad_ga1c19b4bae5de2a391c019b6778a58dad}{redprm()}} to set the land-\/surface paramters, including soil-\/type and veg-\/type dependent parameters.
\item Calculate perturbated soil type \char`\"{}b\char`\"{} parameter. Following Gehne et al. (2019) \cite{Gehne_2019} , a perturbation of LAI \char`\"{}leaf area index\char`\"{} (xlaip) and a perturbation of the empirical exponent parameter b in the soil hydraulic conductivity calculation (bexpp) are added to account for the uncertainties of LAI and b associated with different vegetation types and soil types using a linear scaling. The spatial pattern of xlaip is drawn from a normal distribution with a standard deviation of 0.\+25 while makes the LAI between 0 and 8. The spatial pattern of bexpp is drawn from a normal distribution with a standard deviation of 0.\+4, and is bounded between -\/1 and 1.
\item Calculate perturbated leaf area index.
\item Initialize precipitation logicals.
\item Over sea-\/ice or glacial-\/ice, if water-\/equivalent snow depth ({\itshape sneqv}) below threshold lower bound (0.\+01 m for sea-\/ice, 0.\+10 m for glacial-\/ice), then set at lower bound and store the source increment in subsurface runoff/baseflow (runoff2).
\item For sea-\/ice and glacial-\/ice cases, set smc and sh2o values = 1.\+0 as a flag for non-\/soil medium.
\item If input snowpack ({\itshape sneqv}) is nonzero, then call \mbox{\hyperlink{group___noah___l_s_m_ga028fc448a3afe77726222a337008ffa4_ga028fc448a3afe77726222a337008ffa4}{csnow()}} to compute snow density ({\itshape sndens}) and snow thermal conductivity ({\itshape sncond}).
\item Determine if it\textquotesingle{}s precipitating and what kind of precipitation it is. if it\textquotesingle{}s precipitating and the air temperature is colder than $0^oC$, it\textquotesingle{}s snowing! if it\textquotesingle{}s precipitating and the air temperature is warmer than $0^oC$, but the ground temperature is colder than $0^oC$, freezing rain is presumed to be falling.
\item If either precipitation flag ({\itshape snowng}, {\itshape frzgra}) is set as true\+:
\item Since all precip is added to snowpack, no precip infiltrates into the soil so that {\itshape prcp1} is set to zero.
\item Call \mbox{\hyperlink{group___noah___l_s_m_gab816dbe84c4db91178251b08553c6d02_gab816dbe84c4db91178251b08553c6d02}{snow\+\_\+new()}} to update snow density based on new snowfall, using old and new snow.
\item Call \mbox{\hyperlink{group___noah___l_s_m_ga028fc448a3afe77726222a337008ffa4_ga028fc448a3afe77726222a337008ffa4}{csnow()}} to update snow thermal conductivity.
\item If precipitation is liquid (rain), hence save in the precip variable that later can wholely or partially infiltrate the soil (along with any canopy \char`\"{}drip\char`\"{} added to this later).
\item Determine snowcover fraction and albedo fraction over sea-\/ice, glacial-\/ice, and land.
\item Call \mbox{\hyperlink{group___noah___l_s_m_gafca324f1a0f2bcca12116291bfc8ad68_gafca324f1a0f2bcca12116291bfc8ad68}{snfrac()}} to calculate snow fraction cover.
\item Call \mbox{\hyperlink{group___noah___l_s_m_ga18a0f02e4a1a32b6cde35cbe67b09908_ga18a0f02e4a1a32b6cde35cbe67b09908}{alcalc()}} to calculate surface albedo modification due to snowdepth state.
\item Calculate thermal diffusivity ({\itshape df1})\+:
\begin{DoxyItemize}
\item For sea-\/ice case and glacial-\/ice case, this is constant( $df1=2.2$).
\end{DoxyItemize}
\item For non-\/glacial land case, call \mbox{\hyperlink{group___noah___l_s_m_ga8bc3154096e6364c5fcddf884c003bef_ga8bc3154096e6364c5fcddf884c003bef}{tdfcnd()}} to calculate the thermal diffusivity of top soil layer (\cite{peters-lidard_et_al_1998}).
\item Add subsurface heat flux reduction effect from the overlying green canopy, adapted from section 2.\+1.\+2 of \cite{peters-lidard_et_al_1997}.
\item Calculate subsurface heat flux, {\itshape ssoil}, from final thermal diffusivity of surface mediums,{\itshape df1} above, and skin temperature and top mid-\/layer soil temperature.
\item For uncoupled mode, call \mbox{\hyperlink{group___noah___l_s_m_ga95d889f436d70014fe640fa1b2420545_ga95d889f436d70014fe640fa1b2420545}{snowz0()}} to calculate surface roughness ({\itshape z0}) over snowpack using snow condition from the previous timestep.
\item Calculate virtual temps and virtual potential temps needed by subroutines sfcdif and penman.
\item Calculate the total downward radiation ({\itshape fdown}) = net solar ({\itshape swnet}) + downward longwave ({\itshape lwdn}) as input of \mbox{\hyperlink{group___noah___l_s_m_ga25b90f8a8049eb17cc6cb6a3fecf3fdc_ga25b90f8a8049eb17cc6cb6a3fecf3fdc}{penman()}} and other surface energy budget calculations.
\item Call \mbox{\hyperlink{group___noah___l_s_m_ga25b90f8a8049eb17cc6cb6a3fecf3fdc_ga25b90f8a8049eb17cc6cb6a3fecf3fdc}{penman()}} to calculate potential evaporation ({\itshape etp}), and other partial products and sums for later calculations.
\item Call \mbox{\hyperlink{group___noah___l_s_m_gaa292bfa0bd38e58ab9389cb38e4d09c0_gaa292bfa0bd38e58ab9389cb38e4d09c0}{canres()}} to calculate the canopy resistance and convert it into pc if nonzero greenness fraction.
\item Now decide major pathway branch to take depending on whether snowpack exists or not\+:
\item For no snowpack is present, call \mbox{\hyperlink{group___noah___l_s_m_ga92994e4737ecfd692a6704778558c9c2_ga92994e4737ecfd692a6704778558c9c2}{nopac()}} to calculate soil moisture and heat flux values and update soil moisture contant and soil heat content values.
\item For a snowpack is present, call \mbox{\hyperlink{group___noah___l_s_m_gad65a889bf3f431ade2a312674671ab41_gad65a889bf3f431ade2a312674671ab41}{snopac()}}.
\item Noah LSM post-\/processing\+:
\begin{DoxyItemize}
\item Calculate sensible heat (h) for return to parent model.
\end{DoxyItemize}
\item Convert units and/or sign of total evap (eta), potential evap (etp), subsurface heat flux (s), and runoffs for what parent model expects.
\item Convert the sign of soil heat flux so that\+:
\begin{DoxyItemize}
\item ssoil\texorpdfstring{$>$}{>}0\+: warm the surface (night time)
\item ssoil\texorpdfstring{$<$}{<}0\+: cool the surface (day time)
\end{DoxyItemize}
\item For the case of land (but not glacial-\/ice)\+: convert runoff3 (internal layer runoff from supersat) from $m$ to $ms^-1$ and add to subsurface runoff/baseflow (runoff2). runoff2 is already a rate at this point.
\item For the case of sea-\/ice (ice=1) or glacial-\/ice (ice=-\/1), add any snowmelt directly to surface runoff (runoff1) since there is no soil medium, and thus no call to subroutine smflx (for soil moisture tendency).
\item Calculate total column soil moisture in meters (soilm) and root-\/zone soil moisture availability (fraction) relative to porosity/saturation.
\end{DoxyItemize}

References \mbox{\hyperlink{group___noah___l_s_m_ga18a0f02e4a1a32b6cde35cbe67b09908_ga18a0f02e4a1a32b6cde35cbe67b09908}{alcalc()}}, \mbox{\hyperlink{group___noah___l_s_m_gaa292bfa0bd38e58ab9389cb38e4d09c0_gaa292bfa0bd38e58ab9389cb38e4d09c0}{canres()}}, \mbox{\hyperlink{group___noah___l_s_m_ga028fc448a3afe77726222a337008ffa4_ga028fc448a3afe77726222a337008ffa4}{csnow()}}, \mbox{\hyperlink{group___noah___l_s_m_ga92994e4737ecfd692a6704778558c9c2_ga92994e4737ecfd692a6704778558c9c2}{nopac()}}, \mbox{\hyperlink{group___noah___l_s_m_ga25b90f8a8049eb17cc6cb6a3fecf3fdc_ga25b90f8a8049eb17cc6cb6a3fecf3fdc}{penman()}}, \mbox{\hyperlink{group___noah___l_s_m_ga1c19b4bae5de2a391c019b6778a58dad_ga1c19b4bae5de2a391c019b6778a58dad}{redprm()}}, \mbox{\hyperlink{group___noah___l_s_m_gafdc2dc9815673befca7228be47616e41_gafdc2dc9815673befca7228be47616e41}{sfcdif()}}, \mbox{\hyperlink{group___noah___l_s_m_gafca324f1a0f2bcca12116291bfc8ad68_gafca324f1a0f2bcca12116291bfc8ad68}{snfrac()}}, \mbox{\hyperlink{group___noah___l_s_m_gad65a889bf3f431ade2a312674671ab41_gad65a889bf3f431ade2a312674671ab41}{snopac()}}, \mbox{\hyperlink{group___noah___l_s_m_gab816dbe84c4db91178251b08553c6d02_gab816dbe84c4db91178251b08553c6d02}{snow\+\_\+new()}}, \mbox{\hyperlink{group___noah___l_s_m_ga95d889f436d70014fe640fa1b2420545_ga95d889f436d70014fe640fa1b2420545}{snowz0()}}, and \mbox{\hyperlink{group___noah___l_s_m_ga8bc3154096e6364c5fcddf884c003bef_ga8bc3154096e6364c5fcddf884c003bef}{tdfcnd()}}.



Referenced by \mbox{\hyperlink{group___noah___l_s_m_ga5d624f72a0da3ffcdf7a01fb9862e0e8_ga5d624f72a0da3ffcdf7a01fb9862e0e8}{lsm\+\_\+noah\+::lsm\+\_\+noah\+\_\+run()}}.

Here is the call graph for this function\+:
\nopagebreak
\begin{figure}[H]
\begin{center}
\leavevmode
\includegraphics[width=350pt]{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_ga36e98f3ac91f839b816fd167eab30901_cgraph}
\end{center}
\end{figure}
Here is the caller graph for this function\+:
\nopagebreak
\begin{figure}[H]
\begin{center}
\leavevmode
\includegraphics[width=308pt]{group___noah___l_s_m_ga36e98f3ac91f839b816fd167eab30901_ga36e98f3ac91f839b816fd167eab30901_icgraph}
\end{center}
\end{figure}
