Commit 871e3677 authored by sbraccia's avatar sbraccia

Unit stdout (set in Modules/io_global.f90) is used to write on standard output

instead of 6 or *.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@374 c92efa57-630b-4861-b058-cf58834340f0
parent 22d01804
......@@ -74,7 +74,7 @@ subroutine addusddens (drhoscf, dbecsum, irr, mode0, npe, iflag)
allocate (ylmk0( ngm , lqx * lqx))
allocate (qmod ( ngm))
if (.not.lgamma) allocate (qpg( 3 , ngm))
! write(6,*) aux, ylmk0, qmod
! WRITE( stdout,*) aux, ylmk0, qmod
!
! And then we compute the additional charge in reciprocal space
!
......
......@@ -137,7 +137,7 @@ subroutine cgsolve_all (h_psi, cg_psi, e, d0psi, dpsi, h_diag, &
allocate (rhoold ( nbnd))
allocate (auxr ( ndmx , nbnd))
allocate (eu ( nbnd))
! write(6,*) g,gp,t,h,hold
! WRITE( stdout,*) g,gp,t,h,hold
kter = 0
kter_eff = 0.d0
......@@ -189,7 +189,7 @@ subroutine cgsolve_all (h_psi, cg_psi, e, d0psi, dpsi, h_diag, &
call reduce (1, rho (ibnd) )
#endif
anorm = sqrt (rho (ibnd) )
! write(6,'(2i5,e20.5)') iter,ibnd,anorm
! WRITE( stdout,'(2i5,e20.5)') iter,ibnd,anorm
#ifdef FLUSH
! call flush(6)
#endif
......
......@@ -100,7 +100,7 @@ enddo
! nt=ityp(na)
! do ijh=1,nh(nt)*(nh(nt)+1)/2
! do ipol=1,3
! write(6,'(3i5,f20.10)') na, ijh, ipol,
! WRITE( stdout,'(3i5,f20.10)') na, ijh, ipol,
! + alphasum(ijh,ipol,na,1)
! enddo
! enddo
......
......@@ -82,7 +82,7 @@ enddo
! do na=1,nat
! nt=ityp(na)
! do ijh=1,nh(nt)*(nh(nt)+1)/2
! write(6,'(2i5,f20.10)') na, ijh, becsum(ijh,na,1)
! WRITE( stdout,'(2i5,f20.10)') na, ijh, becsum(ijh,na,1)
! enddo
! enddo
! call stop_ph(.true.)
......
......@@ -90,7 +90,7 @@ do ipol = 1,3
if (mb.gt.nb) call errore('compute_qdipol',' mb.gt.nb', 1)
dpqq(ih,jh,ipol,nt)=fact*ap(lp,ivl,jvl)*qrad2(mb,nb,nt)
dpqq(jh,ih,ipol,nt)=dpqq(ih,jh,ipol,nt)
! write(6,'(3i5,2f15.9)') ih,jh,ipol,dpqq(ih,jh,ipol,nt)
! WRITE( stdout,'(3i5,2f15.9)') ih,jh,ipol,dpqq(ih,jh,ipol,nt)
enddo
enddo
endif
......
......@@ -68,7 +68,7 @@ subroutine compute_weight (wgg)
enddo
! do ibnd=1,nbnd
! do jbnd=1,nbnd
! write(6,'(3i5,f20.10)') ibnd, jbnd, ik,wgg(ibnd,jbnd,ik)
! WRITE( stdout,'(3i5,f20.10)') ibnd, jbnd, ik,wgg(ibnd,jbnd,ik)
! enddo
! enddo
......
......@@ -21,6 +21,7 @@ subroutine d2ionq (nat, ntyp, ityp, zv, tau, alat, omega, q, at, &
! April 1997: parallel stuff added (SdG)
!
#include "machine.h"
USE io_global, ONLY : stdout
use parameters, only : DP
implicit none
!
......@@ -140,7 +141,7 @@ subroutine d2ionq (nat, ntyp, ityp, zv, tau, alat, omega, q, at, &
if (upperbound.gt.1.d-9) goto 11
write (6, '(/5x,"Alpha used in Ewald sum = ",f8.4)') alpha
WRITE( stdout, '(/5x,"Alpha used in Ewald sum = ",f8.4)') alpha
call setv (2 * 3 * nat * nmodes, 0.d0, dy1, 1)
call setv (2 * 3 * nat * nmodes, 0.d0, dy2, 1)
call setv (2 * 3 * nat * nmodes, 0.d0, dy3, 1)
......
......@@ -14,7 +14,7 @@ subroutine dielec
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -59,8 +59,8 @@ call poolreduce (9, epsilon)
!
! symmetrize
!
! write(6,'(/,10x,"Unsymmetrized in crystal axis ",/)')
! write(6,'(10x,"(",3f15.5," )")') ((epsilon(ipol,jpol),
! WRITE( stdout,'(/,10x,"Unsymmetrized in crystal axis ",/)')
! WRITE( stdout,'(10x,"(",3f15.5," )")') ((epsilon(ipol,jpol),
! + ipol=1,3),jpol=1,3)
......@@ -68,8 +68,8 @@ call symtns (epsilon, nsym, s)
!
! pass to cartesian axis
!
! write(6,'(/,10x,"Symmetrized in crystal axis ",/)')
! write(6,'(10x,"(",3f15.5," )")') ((epsilon(ipol,jpol),
! WRITE( stdout,'(/,10x,"Symmetrized in crystal axis ",/)')
! WRITE( stdout,'(10x,"(",3f15.5," )")') ((epsilon(ipol,jpol),
! + ipol=1,3),jpol=1,3)
call trntns (epsilon, at, bg, 1)
!
......@@ -81,9 +81,9 @@ enddo
!
! and print the result
!
write (6, '(/,10x,"Dielectric constant in cartesian axis ",/)')
WRITE( stdout, '(/,10x,"Dielectric constant in cartesian axis ",/)')
write (6, '(10x,"(",3f18.9," )")') ((epsilon(ipol,jpol), ipol=1,3), jpol=1,3)
WRITE( stdout, '(10x,"(",3f18.9," )")') ((epsilon(ipol,jpol), ipol=1,3), jpol=1,3)
call stop_clock ('dielec')
return
......
......@@ -126,7 +126,7 @@ subroutine drhodv (nu_i0, nper, drhoscf)
!
! add to the rest of the dynamical matrix
!
! write(6,*) 'drhodv dyn, wdyn'
! WRITE( stdout,*) 'drhodv dyn, wdyn'
! call tra_write_matrix('drhodv dyn',dyn,u,nat)
! call tra_write_matrix('drhodv wdyn',wdyn,u,nat)
......
......@@ -80,7 +80,7 @@ subroutine drhodvus (irr, imode0, dvscfin, npe)
call reduce (18 * nat * nat, dyn1)
#endif
! write(6,*) 'drhodvus dyn1, dyn'
! WRITE( stdout,*) 'drhodvus dyn1, dyn'
! call tra_write_matrix('drhodvus dyn1',dyn1,u,nat)
! call tra_write_matrix('drhodvus dyn',dyn,u,nat)
! call stop_ph(.true.)
......
......@@ -235,7 +235,7 @@ subroutine dvanqq
! do ih=1,nh(1)
! do jh=1,nh(1)
! do ipol=1,3
! write(6,'(3i5,2f20.10)') ipol,ih,jh,int2(ih,jh,ipol,1,1)
! WRITE( stdout,'(3i5,2f20.10)') ipol,ih,jh,int2(ih,jh,ipol,1,1)
! enddo
! enddo
! enddo
......
......@@ -18,7 +18,7 @@ subroutine dvpsi_e (kpoint, ipol)
!
#include "machine.h"
!
USE io_global, ONLY : stdout
use pwcom
USE wavefunctions, ONLY: evc
use parameters, only : DP
......@@ -199,7 +199,7 @@ subroutine dvpsi_e (kpoint, ipol)
h_diag, npwx, npw, thresh, kpoint, lter, conv_root, anorm, &
nbnd_occ (kpoint) )
if (.not.conv_root) write (6, '(5x,"kpoint",i4," ibnd",i4, &
if (.not.conv_root) WRITE( stdout, '(5x,"kpoint",i4," ibnd",i4, &
& " linter: root not converged ",e10.3)') &
kpoint, ibnd, anorm
#ifdef FLUSH
......
......@@ -16,6 +16,7 @@ subroutine dyndia (xq, nmodes, nat, ntyp, ityp, amass, iudyn, dyn, &
!
!
#include "machine.h"
USE io_global, ONLY : stdout
use parameters, only : DP
implicit none
!
......@@ -94,7 +95,7 @@ subroutine dyndia (xq, nmodes, nat, ntyp, ityp, amass, iudyn, dyn, &
!
! Writes on output the displacements and the normalized frequencies.
!
write (6, 9000) (xq (i), i = 1, 3)
WRITE( stdout, 9000) (xq (i), i = 1, 3)
if (iudyn.ne.0) write (iudyn, 9000) (xq (i), i = 1, 3)
9000 format(/,5x,'Diagonalizing the dynamical matrix', &
......@@ -103,7 +104,7 @@ subroutine dyndia (xq, nmodes, nat, ntyp, ityp, amass, iudyn, dyn, &
do nu_i = 1, nmodes
w1 = sqrt (abs (w2 (nu_i) ) )
if (w2 (nu_i) .lt.0.d0) w1 = - w1
write (6, 9010) nu_i, w1 * rydthz, w1 * rydcm1
WRITE( stdout, 9010) nu_i, w1 * rydthz, w1 * rydcm1
if (iudyn.ne.0) write (iudyn, 9010) nu_i, w1 * rydthz, w1 * &
rydcm1
9010 format (5x,'omega(',i2,') =',f15.6,' [THz] =',f15.6,' [cm-1]')
......@@ -119,7 +120,7 @@ subroutine dyndia (xq, nmodes, nat, ntyp, ityp, amass, iudyn, dyn, &
if (iudyn.ne.0) write (iudyn, '(" (",6f10.6," ) ")') (dyn (mu, &
nu_i) / sqrt (unorm) , mu = 1, 3 * nat)
enddo
write (6, '(1x,74("*"))')
WRITE( stdout, '(1x,74("*"))')
if (iudyn.ne.0) write (iudyn, '(1x,74("*"))')
return
......
......@@ -206,10 +206,10 @@ subroutine dynmat_us
#endif
! do na = 1,nat
! do nb = 1,nat
! write(6, '(2i3)') na,nb
! WRITE( stdout, '(2i3)') na,nb
! do icart = 1,3
! na_icart = 3*(na-1)+icart
! write(6,'(6f13.8)')
! WRITE( stdout,'(6f13.8)')
! + (dynwrk(na_icart,3*(nb-1)+jcart), jcart=1,3)
! end do
! end do
......
......@@ -15,7 +15,7 @@ subroutine dynmatrix
! output.
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -77,7 +77,7 @@ subroutine dynmatrix
!
if (iswitch.eq. - 4) then
do nu = 1, 3 * nat
write (6, '(2i5,2f10.6)') modenum, nu, dyn (modenum, nu)
WRITE( stdout, '(2i5,2f10.6)') modenum, nu, dyn (modenum, nu)
enddo
call stop_ph (.false.)
......
......@@ -13,7 +13,7 @@ subroutine ef_shift (drhoscf, ldos, ldoss, dos_ef, irr, npe, flag)
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
USE wavefunctions, ONLY: evc
use parameters, only : DP
......@@ -66,7 +66,7 @@ external w0gauss
!
call start_clock ('ef_shift')
if (.not.flag) then
write (6, * )
WRITE( stdout, * )
do ipert = 1, npert (irr)
delta_n = (0.d0, 0.d0)
do is = 1, nspin
......@@ -81,7 +81,7 @@ if (.not.flag) then
! symmetrizes the Fermi energy shift
!
call sym_def (def, irr)
write (6, '(5x,"Pert. #",i3,": Fermi energy shift (Ryd) =", 2f10.4)') &
WRITE( stdout, '(5x,"Pert. #",i3,": Fermi energy shift (Ryd) =", 2f10.4)') &
(ipert, def (ipert) , ipert = 1, npert (irr) )
!
! corrects the density response accordingly...
......
......@@ -77,7 +77,7 @@ subroutine incdrhoscf (drhoscf, weight, ik, dbecsum, mode)
call cft3s (dpsic, nr1s, nr2s, nr3s, nrx1s, nrx2s, nrx3s, + 2)
do ir = 1, nrxxs
drhoscf (ir) = drhoscf (ir) + wgt * conjg (psi (ir) ) * dpsic (ir)
! if (ir.lt.20) write (6,*) drhoscf(ir)
! if (ir.lt.20) WRITE( stdout,*) drhoscf(ir)
enddo
enddo
......
......@@ -8,6 +8,7 @@
!---------------------------------------------------------------------
subroutine io_pattern (fildrho,nirr,npert,u,iflag)
!---------------------------------------------------------------------
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
#ifdef __PARA
......@@ -38,12 +39,12 @@ subroutine io_pattern (fildrho,nirr,npert,u,iflag)
call seqopn(iunit,filname,'formatted',exst)
if (iflag.gt.0) then
write (6,*) 'WRITING PATTERNS ON FILE ', filname
WRITE( stdout,*) 'WRITING PATTERNS ON FILE ', filname
write(iunit,*) nirr
write(iunit,*) (npert(i),i=1,nirr)
write(iunit,*) u
else
write (6,*) 'READING PATTERNS FROM FILE ', filname
WRITE( stdout,*) 'READING PATTERNS FROM FILE ', filname
read(iunit,*) nirr
read(iunit,*) (npert(i),i=1,nirr)
read(iunit,*) u
......
......@@ -175,7 +175,7 @@ subroutine localdos (ldos, ldoss, dos_ef)
! check = check + omega*DREAL(ldos(nl(1),is))
! call cft3(ldos(1,is),nr1,nr2,nr3,nrx1,nrx2,nrx3,+1)
! end do
! write (*,*) ' check ', check, dos_ef
! WRITE( stdout,*) ' check ', check, dos_ef
!check
!
deallocate(becsum1)
......
......@@ -16,7 +16,7 @@ program phonon
! computed. In the case q=0 the dielectric constant and the effective
! charges are computed.
!
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -36,7 +36,7 @@ program phonon
call start_clock ('PHONON')
gamma_only = .false.
call startup (nd_nmbr, code, version_number)
write (6, '(/5x,"Ultrasoft (Vanderbilt) Pseudopotentials")')
WRITE( stdout, '(/5x,"Ultrasoft (Vanderbilt) Pseudopotentials")')
!
! and begin with the initialization part
!
......@@ -55,7 +55,7 @@ program phonon
call print_clock ('PHONON')
if (epsil.and.irr0.le.0) then
write (6, '(/,5x," Computing electric fields")')
WRITE( stdout, '(/,5x," Computing electric fields")')
call solve_e
if (convt) then
......
......@@ -31,7 +31,7 @@ subroutine phq_init
! second part of the dynamical matrix.
!
#include"machine.h"
USE io_global, ONLY : stdout
use pwcom
USE wavefunctions, ONLY: evc
use parameters, only : DP
......@@ -128,10 +128,10 @@ subroutine phq_init
if (abs (xq (1) - (xk (1, ikq) - xk (1, ikk) ) ) .gt.1.d-8 .or. &
abs (xq (2) - (xk (2, ikq) - xk (2, ikk) ) ) .gt.1.d-8 .or. &
abs (xq (3) - (xk (3, ikq) - xk (3, ikk) ) ) .gt.1.d-8) then
write (6, * ) ikk, ikq, nksq
write (6, * ) (xq (ipol), ipol = 1, 3)
write (6, * ) (xk (ipol, ikq), ipol = 1, 3)
write (6, * ) (xk (ipol, ikk), ipol = 1, 3)
WRITE( stdout, * ) ikk, ikq, nksq
WRITE( stdout, * ) (xq (ipol), ipol = 1, 3)
WRITE( stdout, * ) (xk (ipol, ikq), ipol = 1, 3)
WRITE( stdout, * ) (xk (ipol, ikk), ipol = 1, 3)
call errore ('phq_init', 'wrong order of k points', 1)
endif
endif
......
......@@ -15,7 +15,7 @@ subroutine phq_recover
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -68,26 +68,26 @@ subroutine phq_recover
endif
if (irr0.eq. - 1) then
write (6, '(/,4x," Reading only int1 and int2" )')
WRITE( stdout, '(/,4x," Reading only int1 and int2" )')
close (unit = iunrec, status = 'keep')
elseif (irr0.eq. - 2) then
call errore ('phq_recover', 'recovering from wrong case', 1)
write (6, '(/,4x," Restart from Iteration #",i5, &
WRITE( stdout, '(/,4x," Restart from Iteration #",i5, &
& " of Elect. Field")') iter0 + 1
elseif (irr0.gt.0) then
if (iter0.ne.0) then
write (6, '(/,4x," Restart from Iteration #",i5, &
WRITE( stdout, '(/,4x," Restart from Iteration #",i5, &
& " of Representation #",i5)') iter0 + 1, irr0
else
if (irr0.eq.nirr) then
write (6, '(/,4x," From recover: Dynamical ", &
WRITE( stdout, '(/,4x," From recover: Dynamical ", &
& " Matrix calculation ")')
else
do irr = 1, nirr
if ( (comp_irr (irr) .eq.1) .and. (done_irr (irr) &
.eq.0) ) then
write (6, '(/,4x," Restart from first iteration", &
WRITE( stdout, '(/,4x," Restart from first iteration", &
& " of Representation # ",i5)') irr
goto 1000
endif
......
......@@ -43,7 +43,7 @@ subroutine phq_setup
! Oct-Nov 1998: minor stuff added (SdG)
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
USE constants, ONLY: degspin
use parameters, only : DP
......@@ -160,7 +160,7 @@ subroutine phq_setup
do ibnd = 1, nbnd
if (et (ibnd, ik) .lt.target) nbnd_occ (ik) = ibnd
enddo
if (nbnd_occ (ik) .eq.nbnd) write (6, '(5x,/,&
if (nbnd_occ (ik) .eq.nbnd) WRITE( stdout, '(5x,/,&
&"Possibly too few bands at point ", i4,3f10.5)') &
ik, (xk (ipol, ik) , ipol = 1, 3)
enddo
......
This diff is collapsed.
......@@ -17,7 +17,7 @@ subroutine phqscf
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -61,10 +61,10 @@ subroutine phqscf
imode0 = imode0 + npert (irr1)
enddo
if (npert (irr) .eq.1) then
write (6, '(//,5x,"Representation #", i3," mode # ",i3)') &
WRITE( stdout, '(//,5x,"Representation #", i3," mode # ",i3)') &
irr, imode0 + 1
else
write (6, '(//,5x,"Representation #", i3," modes # ",3i3)') &
WRITE( stdout, '(//,5x,"Representation #", i3," modes # ",3i3)') &
irr, (imode0+irr1, irr1=1,npert(irr))
endif
!
......@@ -85,11 +85,11 @@ subroutine phqscf
endif
if (convt) then
write (6, '(/,5x,"Convergence has been achieved ")')
WRITE( stdout, '(/,5x,"Convergence has been achieved ")')
done_irr (irr) = 1
iter0 = 0
else
write (6, '(/,5x,"No convergence has been achieved ")')
WRITE( stdout, '(/,5x,"No convergence has been achieved ")')
call stop_ph (.false.)
endif
!
......@@ -110,14 +110,14 @@ subroutine phqscf
! if (tcpu.gt.time_max) then
! temporary fix: recover does not work if program stop here
!
write (6, '(/,5x,"Stopping for time limit ",2f10.0)') &
WRITE( stdout, '(/,5x,"Stopping for time limit ",2f10.0)') &
tcpu, time_max
call stop_ph (.false.)
endif
if (irrc.ge.maxirr) then
write (6, '(/,5x,"Stopping at Representation #",i6)') irr
WRITE( stdout, '(/,5x,"Stopping at Representation #",i6)') irr
#ifdef DEBUG
if (me.ne.1.or.mypool.ne.1) close (6)
#endif
......
......@@ -9,17 +9,17 @@
subroutine print_clock_ph
!-----------------------------------------------------------------------
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
implicit none
write (6, * )
WRITE( stdout, * )
call print_clock ('PHONON')
write (6, * ) ' INITIALIZATION: '
WRITE( stdout, * ) ' INITIALIZATION: '
call print_clock ('phq_setup')
call print_clock ('phq_init')
write (6, * )
WRITE( stdout, * )
call print_clock ('phq_init')
if (nlcc_any) call print_clock ('set_drhoc')
call print_clock ('init_vloc')
......@@ -31,13 +31,13 @@ call print_clock ('drho')
if ((epsil.or.zue).and.okvan) call print_clock ('cmpt_qdipol')
if(epsil) then
write (6, * )
write (6, * ) ' DIELECTRIC CONSTANT AND EFFECTIVE CHARGES:'
WRITE( stdout, * )
WRITE( stdout, * ) ' DIELECTRIC CONSTANT AND EFFECTIVE CHARGES:'
call print_clock ('solve_e')
call print_clock ('dielec')
call print_clock ('zstar_eu')
#ifdef TIMING_ZSTAR_US
write (6, * )
WRITE( stdout, * )
call print_clock ('zstar_eu_us')
call print_clock ('zstar_us_1')
call print_clock ('zstar_us_2')
......@@ -46,7 +46,7 @@ if(epsil) then
call print_clock ('zstar_us_5')
#endif
#ifdef TIMING_ADD_DKMDS
write (6, * )
WRITE( stdout, * )
call print_clock ('add_dkmds')
call print_clock ('add_dkmds1')
call print_clock ('add_dkmds2')
......@@ -57,12 +57,12 @@ if(epsil) then
#endif
endif
if(trans) then
write (6, * )
write (6, * ) ' DYNAMICAL MATRIX:'
WRITE( stdout, * )
WRITE( stdout, * ) ' DYNAMICAL MATRIX:'
call print_clock ('dynmat0')
call print_clock ('phqscf')
call print_clock ('dynmatrix')
write (6, * )
WRITE( stdout, * )
call print_clock ('phqscf')
call print_clock ('solve_linter')
call print_clock ('drhodv')
......@@ -70,19 +70,19 @@ if(trans) then
if (zue) call print_clock('add_zstar_1')
if (zue.and.okvan) call print_clock('add_zstar_us')
endif
write (6, * )
WRITE( stdout, * )
call print_clock ('dynmat0')
call print_clock ('dynmat_us')
call print_clock ('addusdynmat1')
call print_clock ('d2ionq')
if (nlcc_any) call print_clock ('dynmatcc')
write (6, * )
WRITE( stdout, * )
call print_clock ('dynmat_us')
call print_clock ('addusdynmat')
write (6, * )
WRITE( stdout, * )
call print_clock ('phqscf')
call print_clock ('solve_linter')
write (6, * )
WRITE( stdout, * )
call print_clock ('solve_linter')
call print_clock ('dvqpsi_us')
call print_clock ('ortho')
......@@ -104,46 +104,46 @@ call print_clock ('adddvscf')
call print_clock ('drhodvus')
write (6, * )
WRITE( stdout, * )
call print_clock ('dvqpsi_us')
call print_clock ('dvqpsi_us_on')
write (6, * )
WRITE( stdout, * )
call print_clock ('cgsolve')
call print_clock ('ch_psi')
write (6, * )
WRITE( stdout, * )
call print_clock ('ch_psi')
call print_clock ('first')
call print_clock ('h_psiq')
call print_clock ('last')
write (6, * )
WRITE( stdout, * )
call print_clock ('h_psiq')
call print_clock ('firstfft')
call print_clock ('product')
call print_clock ('secondfft')
call print_clock ('add_vuspsi')
write (6, * )
WRITE( stdout, * )
call print_clock ('incdrhoscf')
call print_clock ('addusdbec')
write (6, * )
WRITE( stdout, * )
call print_clock ('drhodvus')
call print_clock ('addusddort')
write (6, * )
write (6, * ) ' General routines'
WRITE( stdout, * )
WRITE( stdout, * ) ' General routines'
call print_clock ('ccalbec')
call print_clock ('cft3')
call print_clock ('cft3s')
call print_clock ('cinterpolate')
call print_clock ('davcio')
call print_clock ('write_rec')
write (6, * )
WRITE( stdout, * )
#ifdef __PARA
write (6, * ) ' Parallel routines'
WRITE( stdout, * ) ' Parallel routines'
call print_clock ('reduce')
call print_clock ('poolreduce')
#endif
......
......@@ -21,7 +21,7 @@ subroutine punch_plot_e
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -57,8 +57,8 @@ subroutine punch_plot_e
#endif
if (fildrho.eq.' ') return
write (6, '(/5x,"Calling punch_plot_e" )')
write (6, '(5x,"Writing on file ",a)') fildrho
WRITE( stdout, '(/5x,"Calling punch_plot_e" )')
WRITE( stdout, '(5x,"Writing on file ",a)') fildrho
!
! reads drho from the file
!
......
......@@ -20,7 +20,7 @@ subroutine punch_plot_ph
!
#include "machine.h"
USE io_global, ONLY : stdout
use pwcom
use parameters, only : DP
use phcom
......@@ -59,8 +59,8 @@ subroutine punch_plot_ph
#endif
if (fildrho.eq.' ') return
write (6, '(/5x,"Calling punch_plot_ph" )')
write (6, '(5x,"Writing on file ",a)') fildrho
WRITE( stdout, '(/5x,"Calling punch_plot_ph" )')
WRITE( stdout, '(5x,"Writing on file ",a)') fildrho
!
! reads drho from the file
!
......
......@@ -30,6 +30,8 @@ subroutine set_irr (nat, at, bg, xq, s, invs, nsym, rtau, irt, &
! April 1997: parallel stuff added (SdG)
!
#include "machine.h"
USE io_global, ONLY : stdout
use parameters, only : DP
#ifdef __PARA
use mp, only: mp_bcast
......@@ -185,10 +187,10 @@ subroutine set_irr (nat, at, bg, xq, s, invs, nsym, rtau, irt, &
!
if (iverbosity.eq.1) then
do imode=1,3*nat
write(6, '(2x,"autoval = ", e10.4)') eigen(imode)
write(6, '(2x,"Real(aut_vet)= ( ",6f10.5,")")') &
WRITE( stdout, '(2x,"autoval = ", e10.4)') eigen(imode)
WRITE( stdout, '(2x,"Real(aut_vet)= ( ",6f10.5,")")') &
( DREAL(u(na,imode)), na=1,3*nat )
write(6, '(2x,"Imm(aut_vet)= ( ",6f10.5,")")') &
WRITE( stdout, '(2x,"Imm(aut_vet)= ( ",6f10.5,")")') &
( DIMAG(u(na,imode)), na=1,3*nat )
end do
end if
......@@ -309,9 +311,9 @@ subroutine set_irr (nat, at, bg, xq, s, invs, nsym, rtau, irt, &
! u(na+3*nat/2,1)=(0.d0,0.d0)
! enddo
! u(1,1)=(-1.d0,0.d0)
! write(6,'(" Setting mode for testing ")')
! WRITE( stdout,'(" Setting mode for testing ")')
! do na=1,3*nat
! write(6,*) u(na,1)
! WRITE( stdout,*) u(na,1)
! enddo
! nsymq=1
! minus_q=.false.
......
......@@ -198,10 +198,10 @@ subroutine set_irr_mode (nat, at, bg, xq, s, invs, nsym, rtau, &
enddo
enddo
! write(6,*) 'nsymq',nsymq
! WRITE( stdout,*) 'nsymq',nsymq
! do isymq=1,nsymq
! irot=irgq(isymq)