Commit: 99fae860e52b3bfc1532a174bdabde783f954f23 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisão: 99fae860e52b3bfc1532a174bdabde783f954f23 (tree)
• Hora: 2007-10-05 22:09:46
• Autor: iselllo
• Commiter: iselllo

Mensagem de Log

I corrected a few typos and went on developing the code. The aim is to
calculate the time-dependent radial density function.

Mudança Sumário

• modified: TCL-codes/mytest2.tcl

Diff

diff -r 333392d74b58 -r 99fae860e52b TCL-codes/mytest2.tcl
--- a/TCL-codes/mytest2.tcl	Fri Oct 05 10:44:41 2007 +0000
+++ b/TCL-codes/mytest2.tcl	Fri Oct 05 13:09:46 2007 +0000
@@ -65,10 +65,10 @@
 
 #prepare the saving of the results
 
-set obs [open "rg_lorenzo.dat" "w"] 
+#set obs [open "rg_lorenzo.dat" "w"] 
 
 
-for {set i 0} { $i < 20 } { incr i} {
+for {set i 0} { $i < 3 } { incr i} {
 set temp [expr [analyze energy kinetic]/(1.5*$n_part)]
 puts "t=[setmd time] E=[analyze energy total], T=$temp"
 integrate $integ_steps
@@ -90,11 +90,25 @@
 
 }
 
-close $obs
+#close $obs
 puts "end of integration" 
 
-plotObs "rg.dat" { 1:2 } labels { "time" "rg" } out "rg"
-exec gv rg.ps
+#plotObs "rg.dat" { 1:2 } labels { "time" "rg" } out "rg"
+#exec gv rg.ps
+
+set f [open "config_1" "r"]
+while { [blockfile $f read auto] != "eof" } {}
+close $f
+
+puts "ok reading the block file"
+
+set rdf [analyze rdf 0 0 0.9 [expr $box_l/2] 100]
+set rlist ""
+set rdflist ""
+foreach value [lindex $rdf 1] {
+lappend rlist [lindex $value 0]
+lappend rdflist [lindex $value 1]
+}
 
 
 puts "So far so good"
\ No newline at end of file