The authors provide a very comprehensive description of TARTES model, which serves as a very good technical documentation for the model. They also compared the TARTES model results with a few other widely-used snow radiative transfer models (AART, DISORT-Mie, SNICAR-ADv3), which generally show reasonable agreement (within 0.02 for snow albedo). Overall, the manuscript is well written. I have a few minor comments for the authors to address.

Minor comments:

1. It is not clear what the specific improvements in TARTES version 2 are compared to version 1, which needs to be clarified.
2. It would be good to also summarize/discuss some underlying assumptions that may limit the TARTES model applications or the cautions users need to take when applying the model.
3. Lines 53-78: These widely-used snow albedo/radiative transfer models have been reviewed and described in details in He and Flanner (2020; https://doi.org/10.1007/978-3-030-38696-2_3), which would be a resource to mention and refer the interested readers to.
4. Equation (4): please introduce the delta function and F0 variable.
5. It would be good to have a table listing all the input and output variables (with units) for TARTES.
6. Equations (36-37): please give the mathematical expressions of “A” and “B”.
7. Line 244: “… we incorporate the exponential terms in the vector X”. This sentence is not very clear to me.
8. Equations (53-55): It seems that B is not included in these equations, which however should be, right?
9. Line 300: For the negative albedo case, did the authors also reset other relevant quantities (e.g., internal light absorption energy, actinic flux, etc.)? If so, reset to what values?
10. Line 320: Some studies (e.g., Peltoniemi, 2007: https://doi.org/10.1016/j.jqsrt.2007.05.009; He et al., 2017: https://doi.org/10.1002/2017GL072916) have quantified the impact of treating snow grains as densely packed medium on snow albedo, which is worth discussing briefly here.
11. Section 3.2: It is not very clear what the key differences in model physics and/or parameters between the Python and Fortran versions, which needs to be clarified.
12. Line 506: It should note that this conclusion here is applied to semi-infinite snowpack tested in this study.
13. It would be good to add a short paragraph to discuss future plans for TARTES model improvements/developments.

This manuscript provides a comprehensive technical description of the TARTES v2.0 radiative transfer model for snowpack.  TARTES is widely used by the community and is embedded in the CROCUS snow thermodynamic model.  A unique and valuable feature of TARTES is the representation of ice particle asphericity via two parameters that can vary continuously and are not tied to any particular shapes, thus enabling the representation of the "optical shape" of snow via a continuum. The manuscript is well-written and appropriately includes both technical descriptions and comparisons against other snowpack radiative transfer models.  The introduction is well-referenced and provides useful background to the topic.  I have only minor comments and am happy to recommend publication of the manuscript in GMD.

Minor comments:

line 30: "understand" -> "understanding"

lines 34-36: "SSA ... advantageously replaces the grain size as it can be rigorously defined and calculated for any porous medium" - While I agree about the advantages of SSA, "calculating" the surface area of complex shapes and porous media is often non-trivial.  To me, the main advantage of SSA is simply that it is a well-defined physical quantity, whereas the meaning of effective radius can be unclear for complex shapes and porous media.

line 107: "an" -> "a"

line 123: The latter part of the sentence also refers to diffuse illumination, so perhaps change the first party to "... illuminated by a beam source and diffuse light..."

line 163: "into" -> "in"

line 244: The appearance and application of these exponential terms is not immediately clear to me.  Could you please clarify or elaborate on this point?

line 373: "When impurities are added in realistic low quantities, we assume the extinction coefficient of snow is unchanged..." - Although this is certainly a valid approximation for most realistic impurity mixing ratios, there may be situations where high impurity loads appreciably affect the total extinction coefficient, particularly in the near-IR spectrum.  For example, very high dust loads can flatten out the 1.03um ice absorption feature (e.g., Fair et al, 2022, https://tc.copernicus.org/articles/16/3801/2022/) . Hence, it might be useful to provide a rough upper limit to the "low quantity" that applies for this assumption to hold, or in general to define the limits of applicability for model users.

line 466: Briefly, what are the user-controlled inputs to this "atmospheric_incident_spectrum" function?

line 492: The accuracy of the delta-Eddington technique in handling diffuse incident light was also assessed more recently by Dang et al (2019, https://doi.org/10.5194/tc-13-2325-2019)

line 576: "benefice" -> "benefit"

p.29 and Fig. 11: The comparison provided is useful because it shows how the default representations of BC in each model affect the simulated albedo.  Because the BC optical properties are slightly different in each model, one could further explore the sources of differences between TARTES and SNICAR by imposing identical BC properties.  This could be accomplished, e.g., by directly importing the BC MAE values from the SNICAR optics library into TARTES, similar to how dust MAE is used in TARTES.  I am not requesting this for manuscript revisions, but merely highlighting it as an informative sensitivity study for the future.

line 605: "67um": Typo, should be "67nm"

line 605: 67nm appears to be close to the monodisperse radius of maximum MAE.  Because the SNICAR BC properties are mass-weighted from a lognormal distribution, the MAE value at the effective radius corresponding to the maximum MAE for a monodisperse distribution (i.e., the trough in Fig. 12) will always be less than the monodisperse MAE at that radius, since the weighted average incorporates lower MAE values from both sides of that radius.  This could explain the phenomenon described in this part of the text.  But other model factors could also contribute to differences in the absolute simulated albedos shown in Fig 12.  For example, are identical ice refractive indices used in each model simulation?

line 624: "This rates" -> "This rate"

line 638: "model as TARTES" -> "model such as TARTES"

line 654: For which wavelength or wavelength range does g=0.82 apply?

line 655: "These values are the defaults in TARTES v2.0" - Actually, it might be helpful to include a table of all of the default parameter settings in TARTES v2.0, but I leave this to the authors to decide.

Finally, I agree with a comment from the other referee that it would be useful to highlight differences and improvements between versions 2.0 and 1.0 of TARTES.