00001 subroutine printout_NS_shape
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 :: it, ip
00009
00010 open(20,file='rnsshape.dat',status='unknown')
00011
00012 it = ntfeq
00013 do ip = 0, npf
00014 write(20,'(1p,2e14.6)') &
00015 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*sinthg(it)*sinphig(ip)
00016 end do
00017 write(20,'(/)')
00018
00019 ip = npfxzp
00020 do it = 0, ntf
00021 write(20,'(1p,2e14.6)') &
00022 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*costhg(it)
00023 end do
00024 ip = npfxzm
00025 do it = ntf-1, 0, -1
00026 write(20,'(1p,2e14.6)') &
00027 & rs(it,ip)*sinthg(it)*cosphig(ip), rs(it,ip)*costhg(it)
00028 end do
00029 write(20,'(/)')
00030
00031 ip = npfyzp
00032 do it = 0, ntf
00033 write(20,'(1p,2e14.6)') &
00034 & rs(it,ip)*costhg(it), rs(it,ip)*sinthg(it)*sinphig(ip)
00035 end do
00036 ip = npfyzm
00037 do it = ntf-1, 0, -1
00038 write(20,'(1p,2e14.6)') &
00039 & rs(it,ip)*costhg(it), rs(it,ip)*sinthg(it)*sinphig(ip)
00040 end do
00041 close(20)
00042
00043 end subroutine printout_NS_shape
