Declare array c0 in spec64() with its active length rather than its maximum length.

This commit is contained in:
Joe Taylor 2021-06-30 10:40:38 -04:00
parent f543cc21e7
commit 6310fcf619
2 changed files with 3 additions and 4 deletions

View File

@ -61,7 +61,7 @@ subroutine q65_loops(c00,npts2,nsps2,nsubmode,ndepth,jpk0, &
jpk=jpk0 + nsps2*ndt/16 !tsym/16
jpk=max(0,jpk)
jpk=min(29000,jpk)
call spec64(c0,nsps2,mode_q65,jpk,s3,LL,NN)
call spec64(c0,npts2,nsps2,mode_q65,jpk,s3,LL,NN)
call pctile(s3,LL*NN,40,base)
s3=s3/base
where(s3(1:LL*NN)>s3lim) s3(1:LL*NN)=s3lim

View File

@ -1,8 +1,7 @@
subroutine spec64(c0,nsps,mode_q65,jpk,s3,LL,NN)
subroutine spec64(c0,npts,nsps,mode_q65,jpk,s3,LL,NN)
parameter (MAXFFT=20736)
!### Fix this:
complex c0(0:1800000-1) !Complex spectrum of dd()
complex c0(0:npts-1) !Complex spectrum of dd()
complex cs(0:MAXFFT-1) !Complex symbol spectrum
real s3(LL,NN) !Synchronized symbol spectra
real xbase0(LL),xbase(LL)