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
kanchukaitis · Mar 20, 2020 · 36adf9e · 36adf9e
1 parent f7699ef
commit 36adf9e
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git 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
@@ -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;