From 6310fcf6193e44524fa16ae077e51ad95b0d854b Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Wed, 30 Jun 2021 10:40:38 -0400 Subject: [PATCH] Declare array c0 in spec64() with its active length rather than its maximum length. --- lib/qra/q65/q65_loops.f90 | 2 +- lib/spec64.f90 | 5 ++--- 2 files changed, 3 insertions(+), 4 deletions(-) diff --git a/lib/qra/q65/q65_loops.f90 b/lib/qra/q65/q65_loops.f90 index f6a80a2a9..ec93adb09 100644 --- a/lib/qra/q65/q65_loops.f90 +++ b/lib/qra/q65/q65_loops.f90 @@ -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 diff --git a/lib/spec64.f90 b/lib/spec64.f90 index b3b7d54d6..73be44d18 100644 --- a/lib/spec64.f90 +++ b/lib/spec64.f90 @@ -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)