;+
; NAME:
;  vt3d_temp_term0_local
; 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 = vt3d_temp_term0_local(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]
;  latheat : latent heat of sublimation, erg/g, loat or float[n_loc]
;  mflux_esc : mass flux due to escape, g/cm^2/s, 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 4.2-3
; MODIFICATION HISTORY:
;  Written 2011 Dec 31, by Leslie Young, SwRI
;  2016 Mar 24 LAY. Modified for inclusion in vt3d library.
;-
function vt3d_temp_term0_local, sol_0, flux_int, emis, latheat, mflux_esc

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

    return, temp_0

end