00001 subroutine printout_quad_pole(iseq)
00002 use phys_constant, only : c, g, msol, pi
00003 use def_quantities, only : gravmass_sph
00004 use def_matter_parameter, only : ome, radi
00005 use def_quantities, only : Iij, Itf, dt1Itf, dt2Itf, dt3Itf, &
00006 & LGW, dJdt, hplus, hcross
00007 implicit none
00008 real(8) :: ff, rr, hplus30k, hcross30k, hplus30M, hcross30M
00009 integer :: i, j
00010 integer, intent(in) :: iseq
00011
00012 if (iseq.eq.1) then
00013 open(300,file='rnsquadpole.dat',status='unknown')
00014 else
00015 open(300,file='rnsquadpole.dat',status='old',position='append')
00016 end if
00017
00018 do i=1,3
00019 do j=1,3
00020 write (300, '(a36, i1, a1, i1, a4, 1es23.15)') &
00021 & 'Mass quadrupole moment I(', i, ',', j, ') = ', Iij(i, j)
00022 end do
00023 end do
00024 write (300, *)
00025 do i=1,3
00026 do j=1,3
00027 write (300, '(a36, i1, a1, i1, a4, 1es23.15)') &
00028 & 'Reduced quadrupole moment Itf(',i,',',j, ') = ', Itf(i,j)
00029 end do
00030 end do
00031 write (300, *)
00032 do i=1,3
00033 do j=1,3
00034 write (300, '(a36, i1, a1, i1, a4, 1es23.15)') &
00035 & 'First time derivative of QP dt1Itf(',i,',',j, ') = ', dt1Itf(i,j)
00036 end do
00037 end do
00038 write (300, *)
00039 do i=1,3
00040 do j=1,3
00041 write (300, '(a36, i1, a1, i1, a4, 1es23.15)') &
00042 & 'Second time derivative of QP dt2Itf(',i,',',j, ') = ', dt2Itf(i,j)
00043 end do
00044 end do
00045 write (300, *)
00046 do i=1,3
00047 do j=1,3
00048 write (300, '(a36, i1, a1, i1, a4, 1es23.15)') &
00049 & 'Third time derivative of QP dt3Itf(',i,',',j, ') = ', dt3Itf(i,j)
00050 end do
00051 end do
00052
00053 write (300, *)
00054 write(300, '(a10, 1es23.15)') 'LGW = ', LGW
00055 do i = 1, 3
00056 write(300, '(a5, i1, a4, 1es23.15)') 'dJdt(', i, ') = ', dJdt(i)
00057 end do
00058
00059 write(300, *)
00060
00061
00062
00063 rr = 30.0d+8/(3.24d-14*147670.0d0)
00064 hplus30k = hplus/rr
00065 hcross30k = hcross/rr
00066 write(300, '(a16, 1es23.15)') 'h+ (r=30kpc) = ', hplus30k
00067 write(300, '(a16, 1es23.15)') 'hx (r=30kpc) = ', hcross30k
00068
00069 rr = 30.0d+11/(3.24d-14*147670.0d0)
00070 hplus30M = hplus/rr
00071 hcross30M = hcross/rr
00072 write(300, '(a16, 1es23.15)') 'h+ (r=30Mpc) = ', hplus30M
00073 write(300, '(a16, 1es23.15)') 'hx (r=30Mpc) = ', hcross30M
00074
00075 ff = (ome/radi) *c**3 / (pi*g*msol)
00076 write(300, '(a16, 1es23.15)') 'frequency[Hz] = ', ff
00077
00078 close(300)
00079
00080 end subroutine printout_quad_pole
