;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;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 file 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
table 5 erase
table 6 erase
table 7 erase
table 8 erase
table 9 erase
table 10 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
hist write 5 sk 10 table 5
hist write 6 sk 10 table 6
hist write 7 sk 10 table 7
hist write 8 sk 10 table 8
hist write 9 sk 10 table 9
hist write 10 sk 10 table 10

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Obtain the numbers of small, medium and large particles;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

def findballnum
  bshineng_sum = 1.8e9 ;the value of the initial potential energy
  bp = ball_head
  loop while bp # null

    if  b_color(bp) = 0 then
      nsmall = nsmall + 1
    endif

    if  b_color(bp) = 1 then
      nmedium = nmedium + 1
    endif

    if  b_color(bp) = 2 then
      nlarge = nlarge + 1
    endif

  bp = b_next(bp)
  end_loop

end
findballnum

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Obtain the numbers of small, medium and large particles;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

def datadeal

 nn = table_size(2)

 loop n (1,nn)
  aa2 = ytable(2,n)
  aa3 = ytable(3,n)
  aa4 = ytable(4,n)
  aa5 = ytable(5,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
  xtable(15,n) = ytable(1,n)/2.61338
  ytable(12,n) = aa2/bshineng_sum
  ytable(13,n) = aa3/bshineng_sum
  ytable(14,n) = aa4/bshineng_sum
  ytable(15,n) = aa5/bshineng_sum
  
  aa6 = ytable(6,n)
  aa7 = ytable(7,n)
  aa8 = ytable(8,n)
  aa9 = ytable(9,n)
  aa10 = ytable(10,n)
  xtable(16,n) = ytable(1,n)/2.61338
  xtable(17,n) = ytable(1,n)/2.61338
  xtable(18,n) = ytable(1,n)/2.61338
  xtable(19,n) = ytable(1,n)/2.61338
  xtable(20,n) = ytable(1,n)/2.61338

  ytable(16,n) = aa6/bshineng_sum
  ytable(17,n) = aa7/bshineng_sum
  ytable(18,n) = aa8/bshineng_sum
  ytable(19,n) = aa9/bshineng_sum
  ytable(20,n) = aa10/bshineng_sum

 endloop

end
datadeal

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;The data of the normalized rotational energy of the total particles, small, medium and large par-;
;ticles in 'table' are transferred to the files in '.log' type, respectively.                     ;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

set log on
set logfile 3p5_xuanzhuan_x  ;The normalized rotational energy of the small particles.
pr tab 12
set log off

set log on
set logfile 3p5_xuanzhuan_z  ;The normalized rotational energy of the medium particles.
pr tab 13
set log off

set log on
set logfile 3p5_xuanzhuan_d  ;The normalized rotational energy of the large particles.
pr tab 14
set log off

set log on
set logfile 3p5_xuanzhuan_sum  ;The normalized rotational energy of the total particles.
pr tab 15
set log off

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;The data of the normalized translational energy of the total particles, small, medium and large  ;
;particles in 'table' are transferred to the files in '.log' type, respectively.                  ;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

set log on
set logfile 3p5_pingdong_x  ;The normalized translational energy of the small particles.
pr tab 16
set log off

set log on
set logfile 3p5_pingdong_z  ;The normalized translational energy of the medium particles.
pr tab 17
set log off

set log on
set logfile 3p5_pingdong_d  ;The normalized translational energy of the large particles.
pr tab 18
set log off

set log on
set logfile 3p5_xpingdong_sum  ;The normalized translational energy of the total particles.
pr tab 19
set log off