Fixed colnames(state_C) <- colnames(tmpfirstrun) in ReactTrans

R/Rphree_ReactTrans.R

 ##  Functions for dealing with surrogate simulations
 
 ### Marco De Lucia, delucia@gfz-potsdam.de, 2009-2018
-### Time-stamp: "Last modified 2018-08-14 20:59:02 delucia"
+### Time-stamp: "Last modified 2019-05-08 17:18:23 delucia"
 
 
 ##' Todo
@@ -321,12 +321,12 @@ ReactTranspBalanceEq <- function(setup, init, maxtime, step=c("time","iter","fix
         msg("assuming homogeneous medium")
         msg("Names gathered after first run:", paste(colnames(tmpfirstrun), collapse=", "))
         state_C <- matrix(rep(tmpfirstrun,n), byrow=TRUE, ncol=length(prop))
+        colnames(state_C) <- colnames(tmpfirstrun)
     } else {
-        msg("given initial state")
+        msg("given initial state; assuming correct colnames")
         state_C <- init
     }
     
-    colnames(state_C) <- colnames(tmpfirstrun)
     attr(state_C,"immobile") <- immobile
 
     ## ssC <<- state_C