program main
	use calcvector
	use constants
	implicit none
	
	call fileinit
	call fileread
	call calclength(6,21) !*********Bの大きさで最後に値を切り落とす*****
	call unitvector(6,15) !b=B/|B|
	call rot(6,9) !J=rotB
	call hasihosei(9)
	call hasihosei(10)
	call hasihosei(11)
	call vpro(6,9,18) !BxJ
	call unitvector(18,12) !a=BxJ/|BxJ|
	call vpro(12,15,18) !c=axb
	call crosseta(9,9)
	call crosseta(10,10)
	call crosseta(11,11) !J->etaJ
	call spro(9,15,1) !J b
	call spro(9,18,5) !J c
	call vprominus(2,6,9) !-VxB
	call spro(9,12,6) !(-VxB)a
	call spro(9,18,7) !(-VxB)c
	
!**********ここから渦度****************
	call rot(2,9) !omega=rotV
	call hasihosei(9)
	call hasihosei(10)
	call hasihosei(11)
	call spro(9,15,2) !omega para
	call scalarcrossvector(2,15,12) !omegapara(vec)
	call vecsubtraction(9,12,18)
	call calclength(18,3)
	
	!*********************ここから切り落とし用
	call kiriotosi(21,1,2*1.3946E-3)
	call kiriotosi(21,2,2*1.3946E-3)
	call kiriotosi(21,3,2*1.3946E-3)
	call kiriotosi(21,5,2*1.3946E-3)
	call kiriotosi(21,6,2*1.3946E-3)
	call kiriotosi(21,7,2*1.3946E-3)
	!*******************************************
	call filewrite(6)
	call filewrite(1)
	call filewrite(7)
	call filewrite(5)
	call filewrite(2)
	call filewrite(3)
end program main