From 36adf9ed6fb28909fd7e9bd32baa2a2861a57491 Mon Sep 17 00:00:00 2001
From: Kevin Anchukaitis <kevin.anchukaitis@gmail.com>
Date: Fri, 20 Mar 2020 13:51:38 -0700
Subject: [PATCH] updated README with citation, additional model output in
 structure

updared readme with full citation, created output space for the number of total cells (cambium+xylem) per day to compare to wood anatomical and microcoring measurements
---
 README.md | 2 +-
 vsm.m     | 6 +++++-
 2 files changed, 6 insertions(+), 2 deletions(-)

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;