00001 subroutine IO_input_initial_1D_qeos
00002 use phys_constant, only : long
00003 use grid_parameter
00004 use coordinate_grav_r, only : rg
00005 use def_metric
00006 use def_matter
00007 use def_matter_parameter, only : ome, ber, radi
00008 implicit none
00009 real(long) :: tmp, tmp1, dummy
00010 integer :: nrftmp, ir, nrgtmp
00011
00012
00013
00014 rs = 1.0d0
00015 rhog = 0.0d0
00016 ome = 1.0d-2
00017 omef(0:nrf, 0:ntf, 0:npf) = ome
00018 omeg(0:nrf, 0:ntf, 0:npf) = ome
00019
00020 open(2,file='rnsflu_1D.ini',status='old')
00021
00022 read(2,'(5i5)') nrftmp
00023 do ir = 0, nrftmp
00024 read(2,'(1p,6e20.12)') dummy, tmp
00025 rhog(ir, 0:ntg, 0:npg) = tmp
00026 end do
00027 rhof(0:nrf, 0:ntf, 0:npf) = rhog(0:nrf, 0:ntf, 0:npf)
00028 read(2,'(1p,6e20.12)') ber, radi
00029
00030 close(2)
00031
00032
00033
00034 open(3,file='rnsgra_1D.ini',status='old')
00035
00036 read(3,'(5i5)') nrgtmp
00037 do ir = 0, nrgtmp
00038 read(3,'(1p,6e20.12)') dummy, tmp, tmp1
00039 psi(ir, 0:ntg, 0:npg) = tmp
00040 alph(ir, 0:ntg, 0:npg) = tmp1
00041 alps = tmp * tmp1
00042 end do
00043
00044 close(3)
00045 bvxd = 0.0d0
00046 bvyd = 0.0d0
00047 bvzd = 0.0d0
00048 tfkij = 0.0d0
00049 tfkijkij = 0.0d0
00050 end subroutine IO_input_initial_1D_qeos
