From 4891eb093f61d37056c50c572e669349dd49a65a Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Tue, 31 May 2016 19:33:33 +0000
Subject: [PATCH] Nucleus translates with the vesicle. Dump and restore includes optional <nucleus> tag within <trisurf> tag if the nucleus is present at all.

---
 src/initial_distribution.c |    6 ++++--
 1 files changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/initial_distribution.c b/src/initial_distribution.c
index 7dfd23d..ae7de0b 100644
--- a/src/initial_distribution.c
+++ b/src/initial_distribution.c
@@ -52,7 +52,9 @@
 	ts_vertex *vtx;
 	ts_tape *tape=vesicle->tape;
 	vesicle->R_nucleus=tape->R_nucleus*tape->R_nucleus;
-
+	vesicle->R_nucleusX=tape->R_nucleusX*tape->R_nucleusX;
+	vesicle->R_nucleusY=tape->R_nucleusY*tape->R_nucleusY;
+	vesicle->R_nucleusZ=tape->R_nucleusZ*tape->R_nucleusZ;
 	vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies;
 
 	//Initialize grafted polymers (brush):
@@ -93,7 +95,7 @@
 	vesicle->dmax=tape->dmax*tape->dmax; /* dmax^2 in the vesicle dmax variable */
 	vesicle->bending_rigidity=tape->xk0;
 	vtx_set_global_values(vesicle); /* make xk0 default value for every vertex */ 
-	ts_fprintf(stdout, "Tape setting: xk0=%e\n",tape->xk0);
+//	ts_fprintf(stdout, "Tape setting: xk0=%e\n",tape->xk0);
 	vesicle->stepsize=tape->stepsize;
 	vesicle->clist->ncmax[0]=tape->ncxmax;
 	vesicle->clist->ncmax[1]=tape->ncymax;

--
Gitblit v1.9.3