Commit 7d289e40 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

ts/wf/evelo [FIX]: fix Fortran compilation warnings

parent 723f1c18
......@@ -79,16 +79,16 @@ c commandline
integer maxopt, lastarg, iargc
character*80 argument
parameter(maxopt=11)
character*2 optid(maxopt)
character*3 optid(maxopt)
character*40 optarg(maxopt)
logical optset(maxopt), opthasarg(maxopt)
c debugging
logical debug, verbose
c here are the keys to our commandline options
data optid/2h-d, 2h-v, 2h-H, 2h-R, 2h-I, 2h-A, 2h-P, '-4',
data optid/'-d', '-v', '-H', '-R', '-I', '-A', '-P', '-4',
& '-ty', '-it', '-ot'/
data opthasarg/8*.FALSE.,3*.TRUE./
data optarg/8*1h-,3*'sff'/
data optarg/8*'-',3*'sff'/
c
c------------------------------------------------------------------------------
c basic information
......@@ -193,7 +193,8 @@ c
if (hilbert) comment='calculate Hilbert transform'
if (pi4shift) comment='shift by pi/4'
if (realspec) comment='calculate real part of Fourier spectrum'
if (imagspec) comment='calculate imaginary part of Fourier spectrum'
if (imagspec) comment=
& 'calculate imaginary part of Fourier spectrum'
if (ampspec) comment='calculate Fourier amplitude spectrum'
if (phasespec) comment='calculate Fourier phase spectrum'
print *,comment
......@@ -229,7 +230,7 @@ c
c
do i=1,nspa
ftarray(i)=(0.d0,0.d0)
if (i.le.nsamp) ftarray(i)=cmplx(fdata(i), 0.d0)
if (i.le.nsamp) ftarray(i)=cmplx(fdata(i), 0.)
enddo
c
call tf_dfork(nspa, ftarray, -1.d0)
......@@ -265,31 +266,32 @@ c
scalefact=dt*sqrt(float(nspa))
if (hilbert.or.pi4shift) then
do i=1,nsamp
fdata(i)=real(ftarray(i))
fdata(i)=sngl(real(ftarray(i)))
enddo
elseif (envelope) then
do i=1,nsamp
fdata(i)=2.d0*abs(ftarray(i))
fdata(i)=sngl(2.d0*abs(ftarray(i)))
enddo
elseif (realspec) then
nsamp=nspa
do i=1,nsamp
fdata(i)=scalefact*real(ftarray(i))
fdata(i)=sngl(scalefact*real(ftarray(i)))
enddo
elseif (imagspec) then
nsamp=nspa
do i=1,nsamp
fdata(i)=scalefact*imag(ftarray(i))
fdata(i)=sngl(scalefact*imag(ftarray(i)))
enddo
elseif (ampspec) then
nsamp=nspa
do i=1,nsamp
fdata(i)=scalefact*abs(ftarray(i))
fdata(i)=sngl(scalefact*abs(ftarray(i)))
enddo
elseif (phasespec) then
nsamp=nspa
do i=1,nsamp
fdata(i)=atan2(imag(ftarray(i)),real(ftarray(i)))*180./pi
fdata(i)=sngl(atan2(imag(ftarray(i)),real(ftarray(i))))*
& 180./pi
enddo
else
stop 'ERROR in selection code'
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment