00001 subroutine read_surf_parameter_mpt(impt)
00002 use phys_constant, only : long
00003 use grid_parameter
00004 implicit none
00005 integer,intent(in) :: impt
00006 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00007
00008 if(impt==1 .or. impt==2) then
00009 open(1,file='rnspar_surf_mpt'//np(impt)//'.dat',status='old')
00010 read(1,'(i5,4x,a1,3i5)') nrg_1, sw_L1_iter, sw_sepa, sw_quant, sw_spin
00011
00012 read(1,'(1p,1e23.15)') r_surf
00013 read(1,'(1p,1e23.15)') sepa_km
00014 read(1,'(1p,1e23.15)') target_qt
00015 read(1,'(1p,1e23.15)') target_sx
00016 read(1,'(1p,1e23.15)') target_sy
00017 read(1,'(1p,1e23.15)') target_sz
00018 close(1)
00019 end if
00020 end subroutine read_surf_parameter_mpt
00021
00022
00023 subroutine read_surf_parameter_mpt_cactus(impt, dir_path)
00024 use phys_constant, only : long
00025 use grid_parameter
00026 implicit none
00027 integer,intent(in) :: impt
00028 character*400, intent(in) :: dir_path
00029 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00030
00031 if(impt==1 .or. impt==2) then
00032 open(1,file=trim(dir_path)//'/'//'rnspar_surf_mpt'//np(impt)//'.dat',status='old')
00033 read(1,'(i5,4x,a1,3i5)') nrg_1, sw_L1_iter, sw_sepa, sw_quant, sw_spin
00034
00035 read(1,'(1p,1e23.15)') r_surf
00036 read(1,'(1p,1e23.15)') sepa_km
00037 read(1,'(1p,1e23.15)') target_qt
00038 read(1,'(1p,1e23.15)') target_sx
00039 read(1,'(1p,1e23.15)') target_sy
00040 read(1,'(1p,1e23.15)') target_sz
00041 close(1)
00042 end if
00043 end subroutine read_surf_parameter_mpt_cactus
00044
