00001 subroutine read_omega
00002   use phys_constant, only : long, pi
00003   use def_bh_parameter, only : ome_bh
00004   implicit none
00005   open(1,file='omega_last.txt',status='old')
00006   read(1,'(1p,e20.12)') ome_bh
00007   close(1)
00008 end subroutine read_omega