00001 subroutine read_parameter_mpt(impt)
00002 use phys_constant, only : long
00003 use def_matter_parameter, only : emdc, pinx
00004 use def_quantities, only : restmass_sph, gravmass_sph, &
00005 & MoverR_sph, schwarz_radi_sph
00006 use grid_parameter
00007 implicit none
00008 integer,intent(in) :: impt
00009 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00010 real(long) :: emdc_ini
00011 open(1,file='rnspar_mpt'//np(impt)//'.dat',status='old')
00012 read(1,'(4i5)') nrg, ntg, npg, nlg
00013 read(1,'(4i5)') nrf, ntf, npf, nlf
00014 read(1,'(2i5,2(3x,a2))') nrf_deform, nrgin, NS_shape, EQ_point
00015 read(1,'(1p,3e10.3)') rgin, rgmid, rgout
00016 read(1,'(/,1i5,2(4x,a1))') iter_max, sw_mass_iter, sw_art_deform
00017 read(1,'(1p,2e10.3)') conv0_gra, conv_ini
00018 read(1,'(1p,2e10.3)') conv0_den, conv0_vep
00019 read(1,'(2(3x,a2),3x,3a1)') indata_type, outdata_type, chrot, chgra, chope
00020 read(1,'(1p,2e10.3)') eps, mass_eps
00021 read(1,'(/,2i5)') num_sol_seq, deform_par
00022 read(1,'(1p,2e14.6)') emdc_ini, pinx
00023 read(1,'(1p,2e14.6)') restmass_sph, gravmass_sph
00024 read(1,'(1p,2e14.6)') MoverR_sph, schwarz_radi_sph
00025 close(1)
00026 ratio = dble(nrf_deform)/dble(nrf)
00027 emdc = emdc_ini
00028 ntgpolp = 0; ntgpolm = ntg; ntgeq = ntg/2; ntgxy = ntg/2
00029 npgxzp = 0; npgxzm = npg/2; npgyzp = npg/4; npgyzm = 3*(npg/4)
00030 ntfpolp = 0; ntfpolm = ntf; ntfeq = ntf/2; ntfxy = ntf/2
00031 npfxzp = 0; npfxzm = npf/2; npfyzp = npf/4; npfyzm = 3*(npf/4)
00032 end subroutine read_parameter_mpt
00033
00034
00035
00036 subroutine read_parameter_mpt_cactus(impt, dir_path)
00037 use phys_constant, only : long
00038 use def_matter_parameter, only : emdc, pinx
00039 use def_quantities, only : restmass_sph, gravmass_sph, &
00040 & MoverR_sph, schwarz_radi_sph
00041 use grid_parameter
00042 implicit none
00043 integer,intent(in) :: impt
00044 character*400, intent(in) :: dir_path
00045 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00046 real(long) :: emdc_ini
00047
00048 open(1,file=trim(dir_path)//'/'//'rnspar_mpt'//np(impt)//'.dat',status='old')
00049 read(1,'(4i5)') nrg, ntg, npg, nlg
00050 read(1,'(4i5)') nrf, ntf, npf, nlf
00051 read(1,'(2i5,2(3x,a2))') nrf_deform, nrgin, NS_shape, EQ_point
00052 read(1,'(1p,3e10.3)') rgin, rgmid, rgout
00053 read(1,'(/,1i5,2(4x,a1))') iter_max, sw_mass_iter, sw_art_deform
00054 read(1,'(1p,2e10.3)') conv0_gra, conv_ini
00055 read(1,'(1p,2e10.3)') conv0_den, conv0_vep
00056 read(1,'(2(3x,a2),3x,3a1)') indata_type, outdata_type, chrot, chgra, chope
00057 read(1,'(1p,2e10.3)') eps, mass_eps
00058 read(1,'(/,2i5)') num_sol_seq, deform_par
00059 read(1,'(1p,2e14.6)') emdc_ini, pinx
00060 read(1,'(1p,2e14.6)') restmass_sph, gravmass_sph
00061 read(1,'(1p,2e14.6)') MoverR_sph, schwarz_radi_sph
00062 close(1)
00063 ratio = dble(nrf_deform)/dble(nrf)
00064 emdc = emdc_ini
00065 ntgpolp = 0; ntgpolm = ntg; ntgeq = ntg/2; ntgxy = ntg/2
00066 npgxzp = 0; npgxzm = npg/2; npgyzp = npg/4; npgyzm = 3*(npg/4)
00067 ntfpolp = 0; ntfpolm = ntf; ntfeq = ntf/2; ntfxy = ntf/2
00068 npfxzp = 0; npfxzm = npf/2; npfyzp = npf/4; npfyzm = 3*(npf/4)
00069 end subroutine read_parameter_mpt_cactus
00070
