提交: 451a4306fe74c0b9641118379daaaeeb948dae62 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 451a4306fe74c0b9641118379daaaeeb948dae62 (tree)
• 时间: 2007-11-21 01:42:35
• 作者: iselllo
• Commiter: iselllo

Log Message

I added the code read_test.tcl which specifically reads the output
(velocities and positions) of test.tcl.

更改概述

• add: TCL-codes/read_test.tcl

差异

diff -r 5651476b69b5 -r 451a4306fe74 TCL-codes/read_test.tcl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/TCL-codes/read_test.tcl	Tue Nov 20 16:42:35 2007 +0000
@@ -0,0 +1,153 @@
+
+# first a test on reading a single configuration
+
+
+
+set f [open "config_vel_7" "r"]
+while { [blockfile $f read auto] != "eof" } {}
+close $f
+
+puts "ok reading the block file"
+
+set Mpos2 {}
+set Mvel2 {}
+for {set p 0} {$p <= [setmd max_part]} {incr p} {
+       set pos2 [part $p print pos]
+       set vel2  [part $p print v]
+       lappend Mpos2 $pos2
+       lappend Mvel2 $vel2
+}
+
+puts "Mpos is, $Mpos2"
+puts "M_vel is, $Mvel2"
+
+
+puts "OK reading single configuration for velocity and position"
+
+
+
+
+
+
+#The following 2 parameters are needed to read the saved configurations
+
+set my_increase 1
+
+# number of snapshots of the system
+set tot_config 10
+
+# number of particles
+
+set n_part 10
+
+set new_iter [expr $n_part*$tot_config]
+
+#counter
+set cnt 0
+
+
+#lists containing the velocity and position data
+set Mpos {}
+set Mvel {}
+
+
+set j 0
+while {$j< $tot_config} {
+puts "j is $j"
+set f [open "config_vel_$j" "r"]
+while { [blockfile $f read auto] != "eof" } {}
+close $f
+
+
+
+
+for {set p 0} {$p <= [setmd max_part]} {incr p} {
+       set pos [part $p print pos]
+       set vel  [part $p print v]
+       lappend Mpos $pos
+       lappend Mvel $vel
+}
+
+
+
+incr cnt
+incr j $my_increase
+}
+
+
+#puts "Mpos is, $Mpos"
+
+#puts "Mvel is, $Mvel"
+
+
+
+
+###########################################################################################################
+#############################################################################################################
+#I tested I can read the positions and velocities correctly from the input files, but I also need to save them into 
+# two large files
+
+
+
+puts "the total number of read configurations is, $cnt"
+
+puts "ok reading the configurations"
+
+
+#I define wlist which I'll be using to write the results to a file
+
+set wlist "" 
+
+
+
+
+set j 0
+while {$j< $new_iter} {
+
+foreach value [lindex $Mpos $j] {
+
+lappend wlist [lindex $value ]
+
+}
+incr j $my_increase
+}
+
+
+set obs3 [open "total_configuration.dat" "w"]
+
+foreach r $wlist  { puts $obs3 "$r" }
+close $obs3
+
+
+#Now I save the particle velocities as a 1D array
+#(same as above with the positions)
+
+set wlist2 ""
+
+
+set j 0
+while {$j< $new_iter} {
+
+foreach value [lindex $Mvel $j] {
+
+
+
+
+lappend wlist2 [lindex $value ]
+
+}
+incr j $my_increase
+}
+
+
+
+
+set obs4 [open "total_configuration_vel.dat" "w"]
+
+foreach r $wlist2  { puts $obs4 "$r" }
+close $obs4
+
+
+
+
+puts "So far so good"
\ No newline at end of file