64 integer :: itime,i,ix,jy,kz,ks,kp,sp_count_i,sp_count_r
65 real :: sp_fact,fact_recept
66 real,
parameter :: smallnum = tiny(0.0)
68 character(len=3) :: anspec
78 write(anspec,
'(i3.3)') ks
79 open(97,file=path(2)(1:length(2))//
'grid_initial'// &
80 '_'//anspec,form=
'unformatted')
83 do kp=1,maxpointspec_act
85 if (ind_rel.eq.1)
then
86 fact_recept=rho_rel(kp)
104 write(97) (sparse_dump_i(i),i=1,sp_count_i)
106 write(97) (sparse_dump_r(i),i=1,sp_count_r)
108 write(97) (sparse_dump_i(i),i=1,sp_count_i)
110 write(97) (sparse_dump_r(i),i=1,sp_count_r)
121 if (init_cond(ix,jy,kz,ks,kp).gt.smallnum)
then
122 if (sp_zer.eqv..true.)
then
123 sp_count_i=sp_count_i+1
124 sparse_dump_i(sp_count_i)= &
125 ix+jy*numxgrid+kz*numxgrid*numygrid
127 sp_fact=sp_fact*(-1.)
129 sp_count_r=sp_count_r+1
130 sparse_dump_r(sp_count_r)=sp_fact* &
131 init_cond(ix,jy,kz,ks,kp)/xmass(kp,ks)*fact_recept
139 write(97) (sparse_dump_i(i),i=1,sp_count_i)
141 write(97) (sparse_dump_r(i),i=1,sp_count_r)
subroutine initial_cond_output(itime)