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 2019-05-08 17:18:23 delucia"
+### Time-stamp: "Last modified 2019-05-08 17:19:57 delucia"
 
 
 ##' Todo
@@ -624,12 +624,12 @@ ReactTranspBalanceKin <- function(setup, init, maxtime, step=c("time","iter","fi
         tmpfirstrun <- RunPQC(initsim)
         msg("assuming homogeneous medium")
         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
 
     ## we save with pe/pH