Revision | 93e904e0c512de5fb83a29684987abed670ba955 (tree) |
---|---|
Zeit | 2008-06-16 21:47:03 |
Autor | iselllo |
Commiter | iselllo |
Now the code is able to initialize its own random number generator, so I
do not need any longer any input file with the initial particle
positions.
@@ -2,11 +2,11 @@ | ||
2 | 2 | |
3 | 3 | # number of monomers in each cluster |
4 | 4 | |
5 | -set n_part 5 | |
5 | +set n_part 50 | |
6 | 6 | |
7 | 7 | #number of clusters |
8 | 8 | |
9 | -set type_of_part 10 | |
9 | +set type_of_part 200 | |
10 | 10 | |
11 | 11 | #total number of monomers |
12 | 12 |
@@ -17,7 +17,7 @@ | ||
17 | 17 | #Now the situation is different: the box size is chosen and density is not relevant. |
18 | 18 | |
19 | 19 | |
20 | -set box_l 200. | |
20 | +set box_l 10000. | |
21 | 21 | |
22 | 22 | # set density 0.03 |
23 | 23 |
@@ -47,7 +47,7 @@ | ||
47 | 47 | |
48 | 48 | # Now I am going to read a simple text file, NOT a block file |
49 | 49 | #set filename "clusters_ini_$j" |
50 | -set filename "initial_state_folded_and_shifted" | |
50 | +set filename "copy_of_cluster50_200" | |
51 | 51 | #set filename "initial_state_folded" |
52 | 52 | |
53 | 53 | set fp [open $filename "r"] |
@@ -105,7 +105,7 @@ | ||
105 | 105 | #parameters of the simulation |
106 | 106 | |
107 | 107 | |
108 | -set tot_time 300. | |
108 | +set tot_time 400. | |
109 | 109 | |
110 | 110 | |
111 | 111 | set my_step 0.00125 |
@@ -120,7 +120,7 @@ | ||
120 | 120 | |
121 | 121 | # number of configurations I want to save |
122 | 122 | |
123 | -set N_save 150 | |
123 | +set N_save 4000 | |
124 | 124 | #I save a configuration every time steps |
125 | 125 | |
126 | 126 | set every [expr $N_step/$N_save] |
@@ -176,6 +176,15 @@ | ||
176 | 176 | set obs5 [open "time.dat" "w"] |
177 | 177 | |
178 | 178 | |
179 | +#initialization random numbers | |
180 | + | |
181 | +expr srand(200) | |
182 | +set l ""; for {set i 0} {$i < [setmd n_nodes]} {incr i} { | |
183 | + lappend l [expr int(32768*rand())] | |
184 | +} | |
185 | +eval t_random seed $l | |
186 | + | |
187 | + | |
179 | 188 | |
180 | 189 | |
181 | 190 | #puts "the simulation time now is, [setmd time]" |