diff --git a/README.md b/README.md
index 398e074..5357b23 100644
--- a/README.md
+++ b/README.md
@@ -6,7 +6,7 @@ This repository provides the Vaganov-Shashkin tree-ring growth model (VSM) in MA
 
 Please cite this work when using this model:
 
-Anchukaitis, K.J., M.N. Evans, M. K. Hughes, and E. Vaganov, An interpreted language implementation of the Vaganov-Shashkin tree-ring proxy system model, submitted to *Dendrochronologia*, 2019
+Anchukaitis, K.J., M.N. Evans, M. K. Hughes, and E. Vaganov, An interpreted language implementation of the Vaganov-Shashkin tree-ring proxy system model, *Dendrochronologia*, 60, 125677, doi:[10.1016/j.dendro.2020.125677](https://doi.org/10.1016/j.dendro.2020.125677), 2020
 
 A preprint is available at [EartharXiv](https://doi.org/10.31223/osf.io/ruyad). 
 
diff --git a/vsm.m b/vsm.m
index f4bf141..6400dfe 100755
--- a/vsm.m
+++ b/vsm.m
@@ -376,6 +376,8 @@
             if sm(1,cyear+1) >= parameters.Wmax; sm(1,cyear+1) = parameters.Wmax; end % error catching
             if isnan(sm(1,cyear+1)); sm(1,cyear+1) = parameters.Wmin; end  % error catching
         end
+    
+        cellCount(t,cyear) = nring(cyear);
         
     end % end day (t) cycle
     
@@ -427,4 +429,6 @@
 output.fday             = fday;
 output.startDay         = sday;
 output.endDay           = eday;
-output.parameters       = parameters; 
+output.parameters       = parameters;
+output.difftime         = RT;
+output.cellCount        = cellCount;