00001
00002 include '../Module/phys_constant.f90'
00003 include '../Module/def_quantities.f90'
00004 include '../Module/def_quantities_derived.f90'
00005 include '../Subroutine/IO_input_physq_plot.f90'
00006 include '../Subroutine/IO_output_physq_plot.f90'
00007 include '../Analysis/Subroutine/IO_output_physq_paper.f90'
00008
00009
00010
00011
00012
00013 PROGRAM pick_up_seq_data
00014
00015 implicit none
00016 integer :: flag_restmass
00017 integer :: total_num_sol, pickup_num_sol, sol_id(40)
00018 integer :: ii, jj
00019
00020 open(20,file='rns_pickup_seq.dat',status='old')
00021 read(20,'(2i5)') total_num_sol, pickup_num_sol
00022 do ii = 1, pickup_num_sol
00023 read(20,'(i5)') sol_id(ii)
00024 end do
00025
00026 do jj = 1, total_num_sol
00027 call IO_input_physq_plot(jj,flag_restmass)
00028 do ii = 1, pickup_num_sol
00029 if (jj.eq.sol_id(ii)) then
00030 call IO_output_physq_plot(ii,flag_restmass)
00031 call IO_output_physq_paper(ii)
00032 end if
00033 end do
00034 end do
00035
00036 END PROGRAM pick_up_seq_data