diff --git a/lib/fst280_decode.f90 b/lib/fst280_decode.f90 index 9a9958a9e..bfbbded60 100644 --- a/lib/fst280_decode.f90 +++ b/lib/fst280_decode.f90 @@ -49,7 +49,7 @@ contains integer*1 message101(101),message74(74) logical badsync,unpk77_success integer*2 iwave(300*12000) - + hmod=1 !### pass as arg ### Keff=91 ndeep=3 diff --git a/lib/jt9.f90 b/lib/jt9.f90 index 3c5a3f9c2..49bfbea2a 100644 --- a/lib/jt9.f90 +++ b/lib/jt9.f90 @@ -202,20 +202,18 @@ program jt9 go to 2 1 nutc=0 2 nsps=0 + npts=TRperiod*12000.d0 if(TRperiod.eq.60.d0) then nsps=6912 shared_data%params%nzhsym=181 + npts=(TRperiod-6.d0)*12000.d0 endif - if(mode.eq.280) then - nsps=6912 - npts=TRperiod*12000.d0 - endif + if(mode.eq.280) nsps=6912 if(nsps.eq.0) stop 'Error: bad TRperiod' - + kstep=nsps/2 k=0 nhsym0=-999 - npts=(TRperiod-6.d0)*12000.d0 if(iarg .eq. offset + 1) then call init_timer (trim(data_dir)//'/timer.out') call timer('jt9 ',0) @@ -250,7 +248,7 @@ program jt9 close(unit=wav%lun) shared_data%params%nutc=nutc shared_data%params%ndiskdat=.true. - shared_data%params%ntr=60 + shared_data%params%ntr=TRperiod shared_data%params%nfqso=nrxfreq shared_data%params%newdat=.true. shared_data%params%npts8=74736