00001 subroutine printout_NS_shape_seq(iseq)
00002 use trigonometry_grav_theta, only : sinthg, costhg
00003 use trigonometry_grav_phi, only : sinphig, cosphig
00004 use grid_parameter, only : ntf, npf, ntfeq, ntfxy, &
00005 & npfxzp, npfxzm, npfyzp, npfyzm
00006 use def_matter, only : rs
00007 implicit none
00008 integer :: iseq
00009 integer :: it, ip
00010
00011 if (iseq.eq.1) open(20,file='rnsshape_seq_xy.dat',status='unknown')
00012 if (iseq.eq.1) open(21,file='rnsshape_seq_xz.dat',status='unknown')
00013 if (iseq.eq.1) open(22,file='rnsshape_seq_zy.dat',status='unknown')
00014 if (iseq.eq.1) open(23,file='rnsshape_seq_yz.dat',status='unknown')
00015
00016 it = ntfeq
00017 do ip = 0, npf
00018 write(20,'(1p,2e14.6)') &
00019 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*sinthg(it)*sinphig(ip)
00020 end do
00021 write(20,'(/)')
00022
00023 ip = npfxzp
00024 do it = 0, ntf
00025 write(21,'(1p,2e14.6)') &
00026 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*costhg(it)
00027 end do
00028 ip = npfxzm
00029 do it = ntf-1, 0, -1
00030 write(21,'(1p,2e14.6)') &
00031 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*costhg(it)
00032 end do
00033 write(21,'(/)')
00034
00035 ip = npfyzp
00036 do it = 0, ntf
00037 write(22,'(1p,2e14.6)') &
00038 & rs(it,ip)*costhg(it), rs(it,ip)*sinthg(it)*sinphig(ip)
00039 end do
00040 ip = npfyzm
00041 do it = ntf-1, 0, -1
00042 write(22,'(1p,2e14.6)') &
00043 & rs(it,ip)*costhg(it), rs(it,ip)*sinthg(it)*sinphig(ip)
00044 end do
00045 write(22,'(/)')
00046
00047 ip = npfyzp
00048 do it = 0, ntf
00049 write(23,'(1p,2e14.6)') &
00050 & rs(it,ip)*sinthg(it)*sinphig(ip), rs(it,ip)*costhg(it)
00051 end do
00052 ip = npfyzm
00053 do it = ntf-1, 0, -1
00054 write(23,'(1p,2e14.6)') &
00055 & rs(it,ip)*sinthg(it)*sinphig(ip), rs(it,ip)*costhg(it)
00056 end do
00057 write(23,'(/)')
00058
00059 end subroutine printout_NS_shape_seq
