;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;The code is to transfer the normalized translational and rotational energy of the total particles, sma-;
;ll, medium and large particles in PFC to the files in '.log' type, respectively                         ;
;This is a case of granular flows with fractal dimension 3.5 at the normalized time of 6.0.             ;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

new

restore 3p5_normalized_6.sav

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Write the recorded data into the specified 'table'.;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

table 1 erase
table 2 erase
table 3 erase
table 4 erase

hist write 1 sk 10 table 1
hist write 2 sk 10 table 2
hist write 3 sk 10 table 3
hist write 4 sk 10 table 4

def datadeal

 nn = table_size(2)

 loop n (1,nn)
  aa2 = ytable(2,n)
  aa3 = ytable(3,n)
  aa4 = ytable(4,n)

  xtable(12,n) = ytable(1,n)/2.61338 ;2.61338 is the value of the characteristic time scale.
  xtable(13,n) = ytable(1,n)/2.61338
  xtable(14,n) = ytable(1,n)/2.61338

  ytable(12,n) = aa2
  ytable(13,n) = aa3
  ytable(14,n) = aa4
 endloop

end
datadeal

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;The rotation velocity of the small, medium and large particles in 'table' is transferred to the files in '.log' type, respectively.;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

set log on
set logfile 3p5_rotation_x  ;The rotation velocity of the small particles.
pr tab 12
set log off

set log on
set logfile 3p5_rotation_z  ;The rotation velocity of the medium particles
pr tab 13
set log off

set log on
set logfile 3p5_rotation_d  ;The rotation velocity energy of the large particles
pr tab 14
set log off