\over 2}(\cos(\theta_j) + \cos(\theta_{j+1}))\right],
\label{thetacenter}\end{aligned}$$ and the zone widths are $(\Delta r)_{i'} = r_{i+1} - r_i$ and $(\Delta \theta)_{j'} = \theta_{j+1} - \theta_j$. The definitions of $(\cos\vartheta)_{\beta'}$, $(\sin\vartheta)_{\beta'}$, and $(\cos\varphi)_{\gamma'}$ in (\[spaceDivergence\]) will be given below. Finally, the values of ${\cal N}$ on zone surfaces in (\[spaceDivergence\]) are given by a particular linear interpolation of zone center values. This linear interpolation—which depends on the neutrino mean free paths—has the effect of shifting from “diamond” differencing in diffusive regimes to “upwind” (or “donor-cell”) differencing in free-streaming regions; see [@liebendoerfer02] for details.

A conservative differencing of the momentum space divergence in (\[example\]) is $$\begin{aligned}
-{1\over {\cal V}_{\alpha',\beta',\gamma'}}\left(1\over r\right)_{i'} \left[({\cal A}_\vartheta\,{\cal N})_{\beta+1,\gamma'} -
({\cal A}_\vartheta\,{\cal N})_{\beta,\gamma'} \right] &-&
\nonumber \\
{1\over  {\cal V}_{\alpha',\beta',\gamma'}}{(\cot\theta)_{j'} \over r_{i'}}
\left[({\cal A}_\varphi\,{\cal N})_{\beta',\gamma+1} -
({\cal A}_\varphi\,{\cal N})_{\beta',\gamma} \right]&,& \label{momentumDivergence}\end{aligned}$$ where the momentum space “geometric” factors are $$\begin{aligned}
{\cal V}_{\alpha',\beta',\gamma'}&=& (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta'}) (\Delta \epsilon)_{\alpha'}  (\Delta\vartheta)_{\beta'} (\Delta\varphi)_{\gamma'}, \\
({\cal A}_\vartheta)_{\beta,\gamma'} &=& (\epsilon_{\alpha'})^2 [\sin(\vartheta_{\beta})]^2 (\Delta \epsilon)_{\alpha'} (\Delta\varphi)_{\gamma'}, \\
({\cal A}_\varphi)_{\beta',\gamma} &=&
\over 2}(\cos(\theta_j) + \cos(\theta_{j+1}))\right ], 
 \label{thetacenter}\end{aligned}$$ and the zone widths are $ (\Delta r)_{i' } = r_{i+1 } - r_i$ and $ (\Delta \theta)_{j' } = \theta_{j+1 } - \theta_j$. The definitions of $ (\cos\vartheta)_{\beta'}$, $ (\sin\vartheta)_{\beta'}$, and $ (\cos\varphi)_{\gamma'}$ in (\[spaceDivergence\ ]) will be given below. ultimately, the value of $ { \cal N}$ on zone surfaces in (\[spaceDivergence\ ]) are given by a particular analogue interpolation of zone center values. This analogue interjection — which depends on the neutrino mean free paths — has the impression of shifting from “ diamond ” differencing in diffusing regimes to “ upwind ” (or “ donor - cell ”) differencing in free - streaming regions; see [ @liebendoerfer02 ] for details. 

 A cautious differencing of the momentum space divergence in (\[example\ ]) is $ $ \begin{aligned } 
 -{1\over { \cal V}_{\alpha',\beta',\gamma'}}\left(1\over r\right)_{i' } \left[({\cal A}_\vartheta\,{\cal N})_{\beta+1,\gamma' } - 
 ({ \cal A}_\vartheta\,{\cal N})_{\beta,\gamma' } \right ] & - & 
 \nonumber \\ 
 { 1\over   { \cal V}_{\alpha',\beta',\gamma'}}{(\cot\theta)_{j' } \over r_{i' } } 
 \left[({\cal A}_\varphi\,{\cal N})_{\beta',\gamma+1 } - 
 ({ \cal A}_\varphi\,{\cal N})_{\beta',\gamma } \right ] &, & \label{momentumDivergence}\end{aligned}$$ where the momentum space “ geometric ” agent are $ $ \begin{aligned } 
 { \cal V}_{\alpha',\beta',\gamma'}&= & (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta' }) (\Delta \epsilon)_{\alpha' }   (\Delta\vartheta)_{\beta' } (\Delta\varphi)_{\gamma' }, \\ 
 ({ \cal A}_\vartheta)_{\beta,\gamma' } & = & (\epsilon_{\alpha'})^2 [ \sin(\vartheta_{\beta})]^2 (\Delta \epsilon)_{\alpha' } (\Delta\varphi)_{\gamma' }, \\ 
 ({ \cal A}_\varphi)_{\beta',\gamma } & = &
\oveg 2}(\cos(\theta_j) + \cos(\theta_{j+1}))\rinht],
\label{thetacenjee}\end{almgned}$$ ahd the zune widths are $(\Delta r)_{i'} = r_{i+1} - r_u$ and $(\Delta \theta)_{j'} = \theta_{j+1} - \theta_j$. Tje definutiois of $(\cos\varthete)_{\geta'}$, $(\sik\rarthsba)_{\betc'}$, end $(\cos\varphi)_{\gakma'}$ in (\[spaweDivergence\]) whlu ye given below. Finally, the values of ${\cal N}$ pn zone surfaces in (\[x[aceSpvtrgence\]) are given by a particular linear interpolatiom of zone center values. Thls llnear interpolatioj—which depebds jb the neutrivo mean frte paths—has jhe effect of shifting from “diamund” dnfferencing ib dlxfusive regmmes tj “upwind” (or “donor-celn”) diffetencing in frec-streemint regions; see [@liebendmerfer02] for details.

W conservdtnve differencing of tye momettum spazw dkvedgxncs in (\[edam'le\]) is $$\begih{aligned}
-{1\ovee {\cal V}_{\alpha',\beta',\gamka'}}\jvgt(1\over r\righf)_{i'} \lefe[({\cwl A}_\vartheta\,{\cal N})_{\beta+1,\gamma'} -
({\cal A}_\varthetd\,{\cam N})_{\beta,\gamma'} \right] &-&
\nonunber \\
{1\over  {\cal V}_{\alpha',\heta',\gamma'}}{(\sot\theta)_{j'} \over r_{i'}}
\left[({\cal A}_\varphi\,{\cal N})_{\beta',\gamma+1} -
({\can A}_\vacpfi\,{\cco N})_{\begq',\gwmma} \right]&,& \label{momentumDivergence}\end{aligned}$$ rgete the momentum fpace “geomeyrlc” sactors are $$\bggin{aliyhes}
{\cal V}_{\alpha',\beta',\gamla'}&=& (\epsijon_{\alpha'})^2 \sin(\vawtheya_{\beta'}) (\Delta \epsilon)_{\alpha'}  (\Eelta\varthetc)_{\bera'} (\Delta\varphi)_{\gammc'}, \\
({\cal A}_\varthzta)_{\bets,\gamms'} &=& (\epsilon_{\alpha'})^2 [\sin(\vartketa_{\befa})]^2 (\Delta \epdilon)_{\alphz'} (\Delta\varphi)_{\gammx'}, \\
({\csl A}_\varphi)_{\beta',\gamma} &=&
\over 2}(\cos(\theta_j) + \cos(\theta_{j+1}))\right], \label{thetacenter}\end{aligned}$$ and the are r)_{i'} = - r_i$ and \theta_j$. definitions of $(\cos\vartheta)_{\beta'}$, and $(\cos\varphi)_{\gamma'}$ in will be given below. Finally, the of ${\cal N}$ on zone surfaces in (\[spaceDivergence\]) are given by a particular interpolation of zone center values. This linear interpolation—which depends on the neutrino mean paths—has effect shifting “diamond” differencing in diffusive regimes to “upwind” (or “donor-cell”) differencing in free-streaming regions; see [@liebendoerfer02] for A conservative differencing of the momentum space divergence (\[example\]) is $$\begin{aligned} -{1\over V}_{\alpha',\beta',\gamma'}}\left(1\over r\right)_{i'} \left[({\cal A}_\vartheta\,{\cal N})_{\beta+1,\gamma'} ({\cal N})_{\beta,\gamma'} \right] \nonumber {1\over V}_{\alpha',\beta',\gamma'}}{(\cot\theta)_{j'} \over r_{i'}} A}_\varphi\,{\cal N})_{\beta',\gamma+1} - ({\cal A}_\varphi\,{\cal N})_{\beta',\gamma} \right]&,& \label{momentumDivergence}\end{aligned}$$ where the momentum space “geometric” factors are $$\begin{aligned} {\cal (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta'}) (\Delta\vartheta)_{\beta'} (\Delta\varphi)_{\gamma'}, ({\cal &=& [\sin(\vartheta_{\beta})]^2 (\Delta \epsilon)_{\alpha'} ({\cal A}_\varphi)_{\beta',\gamma} &=&
\over 2}(\cos(\theta_j) + \cos(\theta_{j+1}))\rigHt],
\label{theTacenTer}\End{AlIgneD}$$ and The zone widths aRE $(\DelTa r)_{i'} = r_{i+1} - r_i$ and $(\Delta \theta)_{j'} = \Theta_{J+1} - \tHEta_j$. tHe DefinItions oF $(\CoS\VArtHeTa)_{\BetA'}$, $(\sIN\vArtheTa)_{\bEta'}$, and $(\cOs\varphi)_{\gaMma'}$ In (\[SpaceDivergeNCe\]) Will be giveN beLow. Finally, thE vaLues of ${\CaL N}$ oN Zone sUrfAces iN (\[spacedIvergeNce\]) are givEn BY a partICular liNEAr InteRpolation of zone ceNTeR Values. This lineAr inteRpOLaTIOn—wHicH depends on ThE neutRIno mean FReE PAThs—HAs the effect of Shifting froM “DiaMond” diFfEreNCing in DiffuSiVE reGimes to “upwiNd” (or “Donor-cell”) DifferENcing in FRee-streAming rEgiOns; See [@lIEbEnDoeRfER02] foR DeTaiLS.

A cOnservatIvE dIfferEnciNG OF The mOmeNtum Space Divergence in (\[eXamPle\]) iS $$\BegIn{aliGned}
-{1\oVer {\cAl v}_{\alphA',\beta',\gAmma'}}\lEfT(1\over r\right)_{i'} \lefT[({\cal a}_\vartheta\,{\Cal n})_{\bEta+1,\GaMma'} -
({\caL a}_\varthEta\,{\Cal n})_{\beta,\gaMma'} \righT] &-&
\NonUmBER \\
{1\OvEr  {\cal V}_{\alpha',\beta',\gamMa'}}{(\COT\tHeta)_{j'} \oveR r_{i'}}
\lefT[({\CaL A}_\VArphi\,{\cal n})_{\bEta',\GammA+1} -
({\CAl A}_\vaRphi\,{\CAl n})_{\beta',\gamMa} \righT]&,& \LaBeL{momentUmdivergEnCe}\eNd{aLigneD}$$ WherE the moMentum spAce “geOMetric” factors aRE $$\begin{aligned}
{\CAl v}_{\ALpHA',\betA',\gaMma'}&=& (\epsilon_{\aLpha'})^2 \SIn(\vaRtheTA_{\bEta'}) (\dElta \ePsiloN)_{\aLPhA'}  (\delta\vartheta)_{\beta'} (\DeLtA\varphI)_{\gammA'}, \\
({\cal A}_\vartheta)_{\Beta,\gamma'} &=& (\ePSILon_{\alpha'})^2 [\Sin(\vARtHEta_{\beta})]^2 (\Delta \epSilon)_{\Alpha'} (\Delta\VArphi)_{\gamMa'}, \\
({\cal a}_\varphi)_{\bEta',\gamma} &=&
\over 2}(\cos(\theta_j) +\cos(\thet a_{j+ 1}) )\r ig ht], 
\la bel{thetacente r }\en d{aligned}$$ and the z one w id t hs a r e$(\De lta r)_ { i' } = r _{ i+ 1}-r _i $ and  $( \Delta\theta)_{j '}=\theta_{j+1} -\theta_j$.  Th e definition s o f $(\c os \va r theta )_{ \beta '}$, $ ( \sin\v artheta)_ {\ b eta'}$ ,  and $( \ c os \var phi)_{\gamma'}$ i n  ( \ [spaceDivergen ce\])wi l lb e  gi ven  below. Fi na lly,t he valu e so f ${\ c al N}$ on zon e surfacesi n ( \[spac eD ive r gence\ ]) ar eg ive n by a part icul ar linear  inter p olation of zone  cente r v alu es.T hi slin ea r  in t er pol a tio n—whichde pe nds o n th e n e utri nomean  free  paths—has th e e ffec t  of  shif tingfrom  “ diamo nd” di ffere nc ing in diffusiv e re gimes to“up wi nd”  ( or “d o nor-ce ll” ) d ifferen cing in fre e- s t r ea ming regions; see[@ l i eb endoerfe r02] f o rde t ails.

A  c ons erva t i ve di ffer e nc ing of t he mom e nt um  spacedi vergen ce  in  (\ [exam p le\] ) is $ $\begin{ align e d}
-{1\over {\ c al V}_{\alpha ' ,\ b e ta ' ,\ga mma '}}\left(1\ over r\ri ght) _ {i '}\ left[ ({\ca lA }_ \ vartheta\,{\cal N}) _{ \beta+ 1,\ga mma'} -
({\ca l A}_\vart h e t a\,{\cal  N}) _ {\ b eta,\gamma'} \ right ] &-&
\non u mber \\{1\ov er  {\ca l V}_{\al p h a',\beta ',\ gam ma' }}{ ( \ co t\theta)_{j'} \ over  r _{i'}}\le ft[({\c alA}_ \va rph i\ ,{\cal N} )_{\beta ', \g am ma +1}  -
({ \ cal A}_\ va rph i\ ,{\ cal N } )_{\be ta',\ gamm a}  \ r igh t]&,& \ l ab e l {mom en tu mDiv erg en ce}\e nd{a l ign ed}$$ w here themom e ntum  s pa ce “geo metric” facto rs  are $$\be gi n{a ligned } {\cal V} _{\alpha',\beta',\gamma ' }&=& (\ eps ilon_ {\al pha'})^2\si n(\var the t a_{\be ta'})(\Del ta  \e p s ilon) _ { \a lph a' }  (\Delta \ v art heta) _{ \bet a'} (\D elta\varphi)_{\gam m a'} , \\
({\cal A }_\ vart h e ta )_{ \ be t a,\ ga m ma' } &=& (\epsilon_{ \alpha'})^ 2[ \s in(\varthe t a_{ \b eta})]^ 2 (\Del ta \e p silon)_ {\alpha'}  (\Delta\ va rphi ) _ {\g amma'}, \\ 
({\calA}_\varph i )_{\b e ta ',\ga mma } &=&
\over 2}(\cos(\theta_j)_+ \cos(\theta_{j+1}))\right],
\label{thetacenter}\end{aligned}$$_and the zone widths_are $(\Delta_r)_{i'}_= r_{i+1}_-_r_i$ and $(\Delta_\theta)_{j'} = \theta_{j+1}_- \theta_j$. The definitions_of $(\cos\vartheta)_{\beta'}$, $(\sin\vartheta)_{\beta'}$,_and_$(\cos\varphi)_{\gamma'}$ in (\[spaceDivergence\]) will be given below. Finally, the values of ${\cal N}$ on_zone_surfaces in_(\[spaceDivergence\])_are_given by a particular linear_interpolation of zone center values._This linear_interpolation—which depends on the neutrino mean free paths—has_the_effect of shifting_from “diamond” differencing in diffusive regimes to “upwind” (or_“donor-cell”) differencing in free-streaming regions; see_[@liebendoerfer02] for details.

A_conservative_differencing_of the momentum space_divergence in (\[example\]) is $$\begin{aligned}
-{1\over {\cal_V}_{\alpha',\beta',\gamma'}}\left(1\over r\right)_{i'} \left[({\cal A}_\vartheta\,{\cal N})_{\beta+1,\gamma'} -
({\cal_A}_\vartheta\,{\cal N})_{\beta,\gamma'} \right] &-&
\nonumber \\
{1\over  {\cal_V}_{\alpha',\beta',\gamma'}}{(\cot\theta)_{j'} \over r_{i'}}
\left[({\cal A}_\varphi\,{\cal N})_{\beta',\gamma+1} -
({\cal_A}_\varphi\,{\cal N})_{\beta',\gamma} \right]&,& \label{momentumDivergence}\end{aligned}$$ where_the momentum_space “geometric” factors are $$\begin{aligned}
{\cal_V}_{\alpha',\beta',\gamma'}&=& (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta'})_(\Delta \epsilon)_{\alpha'}_ (\Delta\vartheta)_{\beta'} (\Delta\varphi)_{\gamma'},_\\
({\cal A}_\vartheta)_{\beta,\gamma'} &=& (\epsilon_{\alpha'})^2 [\sin(\vartheta_{\beta})]^2 (\Delta_\epsilon)_{\alpha'} (\Delta\varphi)_{\gamma'}, \\
({\cal_A}_\varphi)_{\beta',\gamma} &=&