;+
; NAME:
;  vty16_temp_term0_bare
; PURPOSE: (one line)
;  Return temperature in balance with mean absorbed solar flux
; DESCRIPTION:
;  Return temperature in balance with mean absorbed solar flux
; CATEGORY:
;  Volatile Transport
; CALLING SEQUENCE:
;  temp_0 = vty16_temp_term0_bare(sol_0, flux_int, emis)
; INPUTS:
;  sol_0 : mean absorbed solar flux, erg/cm^2/s : float[n_term+1] or float[n_loc,n_term+1]
;  flux_int : internal heat flux, erg/cm^2/s, float or float[n_loc]
;  emis : thermal emissivity, unitless, float or float[n_loc]
; OUTPUTS:
;  temp_0: temperature in balance with mean absorbed solar flux, float or float[n_loc]
; RESTRICTIONS:
; PROCEDURE:
;  Young, L. A. 2016, Volatile transport on inhomogeneous surfaces: II. Numerical calculations (VT3D)
;   Resubmitted to Icarus.
;   Eq 3.2-8
; MODIFICATION HISTORY:
;  Written 2011 Dec 31, by Leslie Young, SwRI
;  2016 Mar 24 LAY. Modified for inclusion in vty16 library.
;-
function vty16_temp_term0_bare, sol_0, flux_int, emis

    physconstants
    
    temp_0 = ( (sol_0 + flux_int) / (emis * !phys.sigma) ) ^ 0.25

    return, temp_0

end