#Init
	variable f index Cu-100
	log     $f.log
	variable freq equal 1
	
	units metal
	dimension 3
	newton on
	boundary p p p
	atom_style atomic
	
#Model
	read_data $f.xrd
	
#Set
	#Force
			pair_style eam
			pair_coeff * * Cu_u3.eam
				
	#Param
			neighbor	0.5 bin
			neigh_modify	delay 10

	
#Proc
	#Equi
		#Fix
		  #velocity all create tmin 2349852 mom yes rot yes sum yes dist gaussian 

		#Out
			thermo_style custom step temp pe pzz lz
			dump		1 all xyz ${freq} $f.xyz		
			timestep	0.001
		#Run
			min_style cg
			
			
			thermo ${freq}
			fix 1 all deform 1 z delta 0.0 -21.688E-2 units box
			run 1 post no
			unfix 1
			
			label begin
			variable t loop 21
			
			minimize	1.0e-6 0.001 1000 10000
			#fix 1 all nvt ${tmin} ${tmax} 0.2
			#run 10000 post no
			#unfix 1
			
			# crd_chg 
			fix 1 all deform 1 z delta 0.0 21.688E-3 units box
			run 1 post no
			unfix 1
			
			next t
			jump Young-100.in begin