Contained here is the two-season water storage model with isotopic fluxes. It can be run without any packages. It reproduces the Figure 12 in the manusript "The seasonality of deuterium excess in non-polar precipitation". If you have any question about the code, please contact the lead author (Zhengyu Xia, [email protected]).