Protect against jpk being too large (xdt > 4.8 s ?) in q65_loops.

This commit is contained in:
Joe Taylor 2021-06-29 17:35:07 -04:00
parent 915571827c
commit f543cc21e7
1 changed files with 2 additions and 1 deletions

View File

@ -59,7 +59,8 @@ subroutine q65_loops(c00,npts2,nsps2,nsubmode,ndepth,jpk0, &
ndt=idt/2
if(mod(idt,2).eq.0) ndt=-ndt
jpk=jpk0 + nsps2*ndt/16 !tsym/16
if(jpk.lt.0) jpk=0
jpk=max(0,jpk)
jpk=min(29000,jpk)
call spec64(c0,nsps2,mode_q65,jpk,s3,LL,NN)
call pctile(s3,LL*NN,40,base)
s3=s3/base